We compare the performance of the Stormer methods of orders nine and thirteen, two optimized symplectic explicit Runge-Kutta-Nystrom methods of orders four and seven, and the integrators DIVA, STEP, DXRK8, RKSUITE, ODEX, ODEX2, RKNINT on four challenging N-body simulations of the Solar System.
