integrate log(1+n/cos(theta))dtheta