int_0^x e^(-t^2)erf(ts)/(2 sqrt(pi)) dt