07-Jan-2022 22:37:07

kepler_ode_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Solve the Kepler ODE.

  Parameters:
    e     = 0.6
    t0    = 0
    y0    = (0.4,0,0,2)
    tstop = 120

kepler_euler_test():
  Use euler() to solve the Kepler problem.

  Number of fixed steps = 10000
  Graphics saved as "kepler_euler_plot.png"
  Graphics saved as "kepler_euler_phase.png"
  Graphics saved as "kepler_euler_h1.png"
  Graphics saved as "kepler_euler_h2.png"

kepler_midpoint_test():
  Use midpoint() to solve the Kepler problem.

  Number of fixed steps = 1000
  Graphics saved as "kepler_midpoint_plot.png"
  Graphics saved as "kepler_midpoint_phase.png"
  Graphics saved as "kepler_midpoint_h1.png"
  Graphics saved as "kepler_midpoint_h2.png"

kepler_ode45_test:
  Use ode45() to solve the Kepler ODE.

  Number of fixed steps = 1069
  Graphics saved as "kepler_ode45_plot.png"
  Graphics saved as "kepler_ode45_phase.png"
  Graphics saved as "kepler_ode45_h1.png"
  Graphics saved as "kepler_ode45_h2.png"

kepler_ode_test():
  Normal end of execution.

07-Jan-2022 22:37:17