Integrate[(b x^2+c x+d-Erf[x/Sqrt[2]])^2,x]