function s = stat(x)
n = length(x);
m = sum(x)/n;
s = sum((x-m).^2)*22/n;
end