Home License -- for personal use only. Not for government, academic, research, commercial, or other organizational use. 13-May-2025 19:09:11 two_body_ode_test(): MATLAB/Octave version 9.11.0.2358333 (R2021b) Update 7 Solve two_body_ode(). parameters: t0 = 0 y0 = (1,0,0,0.8) tstop = 7.79 initial_orbit(): MATLAB/Octave version 9.11.0.2358333 (R2021b) Update 7 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.11.0.2358333 (R2021b) Update 7 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. 13-May-2025 19:09:18