two_body_ode_test, an Octave code which calls two_body_ode(), which computes the behavior of two bodies, constrained to lie in a plane, moving under the influence of gravity, with one body much more massive than the other.
The computer code and data files described and made available on this web page are distributed under the MIT license
two_body_ode, an Octave code which simulates the behavior of two bodies, constrained to lie in a plane, moving under the influence of gravity, with one body much more massive than the other.
INITIAL_ORBIT simulates the problem over approximately two orbits.
ORBITAL_DECAY computes about twenty successive orbits, showing how the orbit gradually decays to a more elliptical form.