integrate 1/sqrt(a+b*cos(phi)) d(phi)