A MATLAB Script for Predicting the Evolution of Lunar Orbits

Script for propagating spacecraft lunar orbits subject to non-spherical lunar gravity and third-body gravity.

PDF document and MATLAB script named lprop_matlab which can be used to predict the long-term evolution of lunar orbits. The software integrates the selenocentric or moon-centered Cartesian equations of motion using a variable step size Runge-Kutta-Fehlberg 7(8) numerical method. The selenocentric equations of motion include the non-spherical gravity of the moon and (optionally) the point-mass gravity of the earth and sun. The ephemeris source for this script is JPL DE421 and the default lunar gravity model is LP150q.