diff --git a/matlab/recursive_moments.m b/matlab/recursive_moments.m
index 921cee357..88a8ca8c5 100644
--- a/matlab/recursive_moments.m
+++ b/matlab/recursive_moments.m
@@ -1,6 +1,4 @@
 function [mu,sigma,offset] = recursive_moments(m0,s0,data,offset)
-
-% function [mu,sigma,offset] = recursive_moments(m0,s0,data,offset)
 % Recursive estimation of order one and two moments (expectation and
 % covariance matrix). 
 % 
@@ -27,13 +25,14 @@ function [mu,sigma,offset] = recursive_moments(m0,s0,data,offset)
 [T,n] = size(data);
 
 for t = 1:T
-  tt = t+offset;
-  m1 = m0 + (1/tt)*(data(t,:)'-m0);
-  qq = m1*m1';
-  s1 = s0 + (1/tt)*(data(t,:)'*data(t,:)-qq-s0) + ((tt-1)/tt)*(m0*m0'-qq');
-  m0 = m1;
-  s0 = s1;
+    tt = t+offset;
+    m1 = m0 + (data(t,:)'-m0)/tt;
+    qq = m1*m1';
+    s1 = s0 + ( (data(t,:)'*data(t,:)-qq-s0) + (tt-1)*(m0*m0'-qq') )/tt;
+    m0 = m1;
+    s0 = s1;
 end
 
-mu = m1; sigma = s1;
+mu = m1; 
+sigma = s1;
 offset = offset+T;
\ No newline at end of file