integrate (p + q cos(theta))^3