henon_heiles_ode_test

henon_heiles_ode_test, an Octave code which calls henon_heiles_ode(), which solves the Henon-Heiles ordinary differential equations (ODE) which model the motion of a star around the galactic center.

Licensing:

The computer code and data files described and made available on this web page are distributed under the MIT license

Related Data and Programs:

henon_heiles_ode, an Octave code which solves the Henon-Heiles ordinary differential equations (ODE) which model the motion of a star around the galactic center.

Source Code:

• henon_heiles_ode_test.m, calls all the tests.
• henon_heiles_ode_test.sh, runs all the tests.
• henon_heiles_ode_test.txt, the output file.

• euler.m, solves an ODE using the Euler method.
• henon_heiles_euler.m, used euler() to solve the Henon-Heiles problem.
• henon_heiles_euler_plot.png, a plot of t versus x, x', y, y'..
• henon_heiles_euler_orbit.png, a orbit plot of x versus y.
• henon_heiles_euler_conserved.png, a plot of the conserved quantity.

• midpoint_fixed.m, solves an ODE using the midpoint method, with fixed point iteration.
• henon_heiles_midpoint.m, used midpoint_fixed() to solve the Henon-Heiles problem.
• henon_heiles_midpoint_plot.png, a plot of t versus x, x', y, y'.
• henon_heiles_midpoint_orbit.png, a orbit plot of x versus y.
• henon_heiles_midpoint_conserved.png, a plot of the conserved quantity.

• henon_heiles_ode45.m, used ode45() to solve the Henon-Heiles problem.
• henon_heiles_ode45_plot.png, a plot of t versus x, x', y, y'.
• henon_heiles_ode45_orbit.png, a orbit plot of x versus y.
• henon_heiles_ode45_conserved.png, a plot of the conserved quantity.