Vectorized N-Body Equation

Fully vectorized n-body equations of motion in second and first order form, along with C MEX implementations of an RKN 12-10 integrator.

nbodyVect provides fully vectorized n-body gravitational equations for use with any built-in (or compatible) MATLAB ode solver. Version 2 of the code also provides event functions for detecting collisions and conjunctions, as well as 3 C implementations of the RKN1210 integrator from Dormand et al. (1987), based on the original MATLAB implementation by Rody Oldenhuis.

The MEX functions have been tested on Windows and macOS/Linux and should compile with no issue if a compatible compiler is installed (see https://www.mathworks.com/help/matlab/call-mex-files-1.html for details).

nbodyRKN_c and nbodyRKN_c_Events can be compiled on their own, e.g.:

mex nbodyRKN_c

nbodyRKN_c_Conj requires linking to levmar (https://users.ics.forth.gr/~lourakis/levmar/) - see file header for compilation notes.

Script nbodyVect_test provides a basic test of integrator functionality.