integrate exp(x^2)*(erf(x) - erf(a)) dx