08-Jan-2022 10:45:29

two_body_ode_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Test two_body_ode().

  Parameter values:
    t0    = 0
    y0    = (1,0,0,0.8)
    tstop = 7.79

initial_orbit():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  This simulation follows a small body for two orbits
  around a relatively massive body - such as Mercury around
  the sun.
  Kepler's equations for a two body system are used.
  Note that the orbit is NOT an ellipse.  But that's OK,
  because the planet is far from its equilibrium orbit.
  Graphics saved as "initial_orbit.png"

initial_orbit:
  Normal end of execution.

orbital_decay()
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  This simulation follows a small body for 20 orbits
  around a relatively massive body - such as Mercury around
  the sun.
  Kepler's equations for a two body system are used.
  Initially, the orbit is NOT an ellipse, but as time passes,
  the orbit decays into an elliptical shape.
  Graphics saved as "orbital_decay"

orbital_decay:
  Normal end of execution.

two_body_ode_test:
  Normal end of execution.

08-Jan-2022 10:45:34