The force for each particle is calculated with a direct N-Body scheme, where all forces are calculated accurately from each object to each other object. The equation for the force is (boldface parameters are Cartesian vectors, see Wikipedia for more details.)
The acceleration follows from .
The equations of motion are integrated via the Verlet-leapfrog
predictor corrector integrator, see Wikipedia for more details.
The equation to update the position:
The equation to update the velocity:
The units of the integration are so called Heggie-Units, in which case M=G=1. Here
As a consequence, the total energy of the system is Etot = 0.25. A consise description of these units can be found here.