integral 1/2 (1 + EllipticTheta(3, 0, e^(-(π t)/2)))/(sqrt t(t^2+1)) from .001 to infinity