rand('seed',23);
n=500;
Sigma=20;
z=rand(n);
z=sqrt(-Sigma*2*log(1-z));
y=z(:,1);
hist(z,60); pause;
end