function x = stdnorminv( p )

x = -sqrt(2)*erfcinv(2*p);