NIntegrate[(1+Erf[(3*x/Sqrt[1000]-2)/Sqrt[2]]