plot y=e^{-x^2/2}(1+erf(-3x))