plot {1/(1+erf(sqrt pi/2 y))-1/(1+y)(1+ pi/12 y^3/(1+y)(1+pi/120 (y^3-9y^2)/(1+y)))},{y,0,1.0}