Introduction - If you have any usage issues, please Google them yourself
Programming multi-body dynamics with symplectic methods matlb Fully vectorized n-body equation nbodyVect (x0, dx0, mus, t) integrates the n-body equations of motion using either a symplectic second order integrator, or any of the standard built-in or user supplied first order integrators.