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