07-Jan-2022 23:25:51 ozone_ode_test(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Solve the ozone ODE. Parameter values: k2 = 100000 k3 = 1e-16 t0 = 0 y0 = (0,0,5e+11,8e+11) tstop = 172800 ozone_ode15s(): Use ode15s() to solve the ozone ODE. Number of time steps = 366 Average timestep = 473.425 seconds Graphics saved as "ozone_ode15s_y1.png" Graphics saved as "ozone_ode15s_y2.png" Graphics saved as "ozone_ode15s_y3.png" Graphics saved as "ozone_ode15s_y4.png" Graphics saved as "ozone_ode15s_conservation.png" ozone_ode15s: Normal end of execution. ozone_midpoint(): Use midpoint() to solve the ozone ODE. Number of time steps = 3000 Timestep size = 57.6 seconds Graphics saved as "ozone_midpoint_y1.png" Graphics saved as "ozone_midpoint_y2.png" Graphics saved as "ozone_midpoint_y3.png" Graphics saved as "ozone_midpoint_y4.png" Graphics saved as "ozone_midpoint_conservation.png" ozone_midpoint(): Normal end of execution. ozone_ode_test: Normal end of execution. 07-Jan-2022 23:26:08