Integrate 1/sqrt((d-cos(theta))^2+(sin(theta))^2) from 0 to 2*pi