Hi all, I'm currently working on some n-body simulations and I'm trying to implement higher order integrators. Currently, I have a Leapfrog and a Hermite integrator (From the student series on The Art of Computational Science). The results I get from Leapfrog are excellent, with respect to the accuracy and energy preservation. But, I would like to try implementing higher order symplectic integrators. Does any know of these higher order integrators or of any well written resources on it? Thanks!