integrate (x^2+r^2-2*x*r*cos t)^(-2) dt