PDF document and three MATLAB functions that can be used to propagate two body or “unperturbed” satellite orbits. This type of propagation is also called the orbital initial value problem (IVP). All three algorithms are valid for elliptical and hyperbolic orbits, and each numerical method can propagate forward or backward in time. For programming flexibility, the input and output arguments for each routine are identical.
David Eagle (2020). MATLAB Functions for Two Body Orbit Propagation (https://www.mathworks.com/matlabcentral/fileexchange/48723-matlab-functions-for-two-body-orbit-propagation), MATLAB Central File Exchange. Retrieved .