integrate 1/(a+b\cos[theta])d\theta