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