Integrate[Exp[x^2 - x^2 (Erf[x] + 1/2) ] , {x,0,Infinity}]