Six-body problem solution using symplectic integrators

In this paper, an application of symplectic integrators for solving an important problem of non-linear dynamics is considered. Numerical algorithms for solving six-body problem are described, built on explicit second order Runge-Kutta method, implicit Euler method, symplectic Verlet method and D-method. The examined six-body problem time-domain solution is given, as well as the system's full energy dynamics for the implemented numerical algorithm. Conclusions of symplectic integrators efficiency are made based on simulation results.