function phi=phi(x);
    phi=exp(-x.^2/2)/(2*pi)^(1/2);
end;