function z=varianssi(x)
n=length(x);
a=sum(x)/n;
z=sum((x-a).^2)/(n-1);