This LiveScript is designed to propagate any conic section orbit trajectory about a user-selected spherical celestial object within the solar system using Keplerian orbital mechanics. Using only initial postion and velocity state vectors, one can propagate circular, elliptical, parabolic, or hyperbolic orbits about an orbit body over time using a Cowell Method numerical integrator. The example compares the unperturbed two-body propagation with a purturbed propagated orbit due to orbit body spherical oblateness of the same initial state vectors. Once the orbit state vectors are obtained, it plots the orbit trajectories in three dimensions specific to the selected orbit body and see the differences in trajectories.
Richard DeMark (2019). Orbit Trajectory Propagation (https://www.mathworks.com/matlabcentral/fileexchange/69439-orbit-trajectory-propagation), MATLAB Central File Exchange. Retrieved .