29-May-2021 13:41:40
kepler_ode_test)_:
MATLAB/Octave version 4.2.2
Solve the kepler ODE.
Parameter values:
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 = 259
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.
29-May-2021 13:41:48