-1/(8sqrt pi)(7-sqrt2(1+2/pi arctan sqrt 2)+int_0^infty x e^(-x^2) (erfc x)^3/(1-erfc x /2) dx)