08-Jan-2022 09:15:47 rigid_body_ode_test(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Solve the rigid body ODE. rigid_body_backward_euler(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Solve the rigid body ODE using backward_euler(). Number of fixed steps = 1000 Graphics saved as "rigid_body_backward_backward euler_plot.png" Graphics saved as "rigid_body_backward_backward_euler_phase.png" Graphics saved as "rigid_body_backward_euler_conservation_h1.png" Graphics saved as "rigid_body_backward_euler_conservation_h2.png" rigid_body_euler(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Solve the rigid body ODE using euler(). Number of fixed steps = 1000 Graphics saved as "rigid_body_euler_plot.png" Graphics saved as "rigid_body_euler_phase.png" Graphics saved as "rigid_body_euler_conservation_h1.png" Graphics saved as "rigid_body_euler_conservation_h2.png" rigid_body_midpoint(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Use midpoint() to solve the rigid body ODE. Number of equal steps n = 20000 Graphics saved as "rigid_body_midpoint_plot.png" Graphics saved as "rigid_body_midpoint_phase.png" Graphics saved as "rigid_body_midpoint_conservation_h1.png" Graphics saved as "rigid_body_midpoint_conservation_h2.png" rigid_body_rk4(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Solve the rigid body ODE using rk4(). Number of fixed time steps = 20000 Graphics saved as "rigid_body_rk4_plot.png" Graphics saved as "rigid_body_rk4_phase.png" Graphics saved as "rigid_body_rk4_conservation_h1.png" Graphics saved as "rigid_body_rk4_conservation_h2.png" rigid_body_ode15s(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Solve the rigid body ODE using ode15s(). Number of (variable) time steps = 25165 Graphics saved as "rigid_body_ode15s_plot.png" Graphics saved as "rigid_body_ode15s_phase.png" Graphics saved as "rigid_body_ode15s_conservation_h1.png" Graphics saved as "rigid_body_ode15s_conservation_h2.png" rigid_body_ode23(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Solve the rigid body ODE using ode23(). Number of (variable) time steps = 21266 Graphics saved as "rigid_body_ode23_plot.png" Graphics saved as "rigid_body_ode23_phase.png" Graphics saved as "rigid_body_ode23_conservation_h1.png" Graphics saved as "rigid_body_ode23_conservation_h2.png" rigid_body_ode23s(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Solve the rigid body ODE using ode23s(). Number of (variable) time steps = 44097 Graphics saved as "rigid_body_ode23s_plot.png" Graphics saved as "rigid_body_ode23s_phase.png" Graphics saved as "rigid_body_ode23s_conservation_h1.png" Graphics saved as "rigid_body_ode23s_conservation_h2.png" rigid_body_ode45(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Solve the rigid body ODE using ode45(). Number of (variable) time steps = 24909 Graphics saved as "rigid_body_ode45_plot.png" Graphics saved as "rigid_body_ode45_phase.png" Graphics saved as "rigid_body_ode45_conservation_h1.png" Graphics saved as "rigid_body_ode45_conservation_h2.png" rigid_body_trapezoidal(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Solve the rigid body ODE using the trapezoidal method. Number of fixed steps = 20001 Graphics saved as "rigid_body_trapezoidal_plot.png" Graphics saved as "rigid_body_trapezoidal_phase.png" Graphics saved as "rigid_body_trapezoidal_conservation_h1.png" Graphics saved as "rigid_body_trapezoidal_conservation_h2.png" rigid_body_ode_test(): Normal end of execution. 08-Jan-2022 09:16:44