Let m_i denote the mass of body i, and (x_i, y_i, z_i) denote the position of body i. Then, the following differential equation system describes the movement of the N bodies:

If we let:

then we can re-write the system as:

This system can be solved using the Modified Picard Method of Parker and Sochaki.