The normal cumulative distribution function with mean m and standard deviation s, at the value of x :
y=integral form -inf to x of exp(-0.5*((t-m)./s)^2)./(sqrt(2*pi)*s)dt