plot (1/sqrt(2*pi))*e^(-(1/2)*(x-1)^2)*(1/2)*(1+erf((1/sqrt(2))*4*(x-1))) from x = 0 to x = 5 from y = 0 to y = 0.5