integrate (a^2tan^2theta)/(a-b cos theta) dtheta