function v=moments(x)
% gives mean and standard deviation
v=mean(x);
n=length(x);
s=var(x)*(n-1)/n;
v=[v sqrt(s)];