07-Jan-2022 21:41:01

henon_heiles_ode_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Solve henon heiles ODE using various ODE solvers.

  Parameter values:
    e     = 0.05
    lam   = 1
    t0    = 0
    y0    = (0,0.316228,0,0)
    tstop = 10000

henon_heiles_midpoint_test():
  Use midpoint to solve the Henon Heiles problem.

  Number of fixed steps was 20000
  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_ode23s():
  Use ode23s to solve the Henon Heiles problem.

  Number of variable steps was 72330
  Graphics saved as "henon_heiles_ode23s_plot.png"
  Graphics saved as "henon_heiles_ode23s_orbit.png"
  Graphics saved as "henon_heiles_ode23s_conserved.png"

henon_heiles_ode45_test():
  Use ode45 to solve the Henon Heiles problem.

  Number of variable steps was 59189
  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.

07-Jan-2022 21:41:31