% Cumulative probability distribution of standard normal distribution
function cumpd=cumpd(u)
arg=u/sqrt(2);
cumpd=0.5+0.5*erf(arg);
