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.

Downloads:
Image (Overview, Lights Off)
Image (Pluto, Lights Off)
Image (Sun, Lights Off)
Image (Sun, Lights On)
Image (Earth, Lights On)
Image (Earth, Lights On)
Paper (Modified Picard Method)
Paper (Solution of the N Body Problem)
Video (640x480, WMV, 12Mb)
Video (640x480, AVI, 30Mb)
Video (320x240, MPG, 12Mb)