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.