Let mi denote the mass of body i, and (xi, yi, zi) 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.