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