integrate erf(-x^2)