Mutual gravitational potential and torque of solid bodies via inertia integrals

The mutual gravitational potential and the mutual gravitational torque of two bodies of arbitrary shape are expanded to the fourth order. The derivations are based on Cartesian coordinates, inertia integrals with relation to the principal reference frames of each body, and the relative rotation matrix. The current formulation is convenient to utilize in high precision problems in rotational dynamics.