integrate [(x-a*cos(theta))/((x-a*cos(theta))^2+(y-a*sin(theta))^2)^(3/2),{theta, 0, 2*Pi}]