1/2 integrate (1 + r*cos(x))^2 dx