integrate 1/sqrt(r^2+s^2-2*r*s*cos(phi)+y^2) d(phi)