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.
Download apps, toolboxes, and other File Exchange content using Add-On Explorer in MATLAB.