function r = autocor(eta,M)

r=NaN(M+1,1);

%varEta = (eta'-mean(eta))*(eta-mean(eta));
N = length(eta);

for(n=0:M)
% r(i) = (eta(1:end-i)-mean(eta))'*(eta(i+1:end) - mean(eta))/varEta;
  r(n+1) = eta(1:end-n)'*eta(n+1:end)/(N-n+1);
end
