plot -1/2*I*sqrt(pi)*e^(-x^2)*erf(I*x) from 0 to 10