function y = pdfn (x)
y = 1./sqrt(2*pi) .* exp(-(x.^2)./2);
