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