**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 GNU LGPL 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.

- two_body_ode_test.m, calls all the tests.
- two_body_ode_test.sh, runs all the tests.
- two_body_ode_test.txt, the output file.

**INITIAL_ORBIT** simulates the problem over approximately two orbits.

- initial_orbit.m, the source code.
- initial_orbit.png, an image of the trajectory.

**ORBITAL_DECAY** computes about twenty successive orbits, showing
how the orbit gradually decays to a more elliptical form.

- orbital_decay.m, the source code.
- orbital_decay.png, an image of the trajectory.