numerically integrate exp(-x^2)/(1+x^2) from -1000 to 1000