function p = stdnormcdf( x )

p = 0.5*erfc(-x/sqrt(2));