14-Nov-2020 09:27:11

henon_heiles_ode_test:
  MATLAB/Octave version 4.2.2
  Test henon_heiles_ode using various ODE solvers.

  Parameter lambda = 1
  Initial energy E = 0.05
  Initial condition (X,X',Y,Y') = (0,0.316228,0,0)

henon_heiles_euler_test:
  Use euler to solve the Henon Heiles problem.
  Graphics saved as "henon_heiles_euler_plot.png"
  Graphics saved as "henon_heiles_euler_orbit.png"
  Graphics saved as "henon_heiles_euler_conserved.png"

henon_heiles_midpoint_test:
  Use midpoint to solve the Henon Heiles problem.
  Graphics saved as "henon_heiles_midpoint_plot.png"
  Graphics saved as "henon_heiles_midpoint_orbit.png"
  Graphics saved as "henon_heiles_midpoint_conserved.png"

henon_heiles_ode45_test:
  Use ode45 to solve the Henon Heiles problem.
  Graphics saved as "henon_heiles_ode45_plot.png"
  Graphics saved as "henon_heiles_ode45_orbit.png"
  Graphics saved as "henon_heiles_ode45_conserved.png"

henon_heiles_ode_test:
  Normal end of execution.

14-Nov-2020 09:27:43