function  prob = normalCDF(x)
% CDF for N(0,1)

prob = 0.5 + erf(x/sqrt(2))/2;
