function fz = stdnormdens(z)
fz = 1/sqrt(2*pi)*exp(-1/2*z.^2);
