Home License -- for personal use only.  Not for government,
	academic, research, commercial, or other organizational use.

13-May-2025 15:13:14

kepler_ode_test():
  MATLAB/Octave version 9.11.0.2358333 (R2021b) Update 7
  Solve kepler_ode().

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

kepler_euler_test():
  Use euler() to solve kepler_ode().

  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 kepler_ode().

  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():
  kepler_ode() is solved by ode45().

  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.

13-May-2025 15:13:26