plot {z, sqrt(z^2 - log1p((0.25 - erfc(z/sqrt(2))/2) / (exp(-z^2/2)/sqrt(2*pi)) * z))}, z = 0 .. 3