kepler_ode

kepler_ode, a Python code which defines the ordinary differential equations (ODE) for a Kepler two body gravitational problem.

Licensing:

The information on this web page is distributed under the MIT license.

Languages:

kepler_ode is available in a MATLAB version and an Octave version and a Python version.

Related Data and codes:

python_ode, Python codes which sets up various systems of ordinary differential equations (ODE).

References:

1. Ernst Hairer, Christian Lubich, Gerhard Wanner,
   Geometric Numerical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations,
   Springer, 2006,
   ISSN: 0179-3632.

Source Code:

• kepler_ode.py, the source code.
• kepler_ode.sh, runs all the tests.
• kepler_ode.txt, the output file.

• kepler_solve_ivp_plot.png, a plot of t versus p1, p2, q1, q2.
• kepler_solve_ivp_phase.png, a phase plane plot of p1 versus p2.
• kepler_solve_ivp_h1.png, conservation of energy.
• kepler_solve_ivp_h2.png, conservation of angular momentum.