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 information on this web page is 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.