Integrate exp(-a x^2)/(1+exp(x)) from -infinity to +infinity