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.
You are now following this Submission
- You will see updates in your followed content feed
- You may receive emails, depending on your communication preferences
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.
Cite As
Dmitry Savransky (2026). Vectorized N-Body Equation (https://se.mathworks.com/matlabcentral/fileexchange/27820-vectorized-n-body-equation), MATLAB Central File Exchange. Retrieved .
Acknowledgements
Inspired by: rkn86, RKN1210 - A 12th/10th order Runge-Kutta-Nyström integrator
General Information
- Version 2.0.0.0 (22.9 KB)
MATLAB Release Compatibility
- Compatible with R2016b and later releases
Platform Compatibility
- Windows
- macOS
- Linux
