08-Jan-2022 09:14:42 reaction_twoway_ode_test(): MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2 Solve the two-way reaction ODE. Reaction differential equation: dW1/dt = - k1 W1 + k2 W2 dW2/dt = + k1 W1 - k2 W2 Parameter values: k1 = 1, (reaction rate) k2 = 10, (reaction rate) t0 = 0, (initial time, in seconds s) w10 = 0.1, (initial amount of species W1) w20 = 0.9, (initial amount of species W2) tstop = 1, (final time, in seconds s) reaction_twoway_euler(): Solve the two way reaction ODE using euler() Number of steps will be 51 Stepsize will be 0.02 Graphics saved as "reaction_twoway_euler_solution.png" Graphics saved as "reaction_twoway_euler_conservation.png" reaction_twoway_euler: Normal end of execution. reaction_twoway_midpoint: Solve the two way reaction ODE using midpoint() Number of steps will be 51 Stepsize will be 0.02 Graphics saved as "reaction_twoway_midpoint_solution.png" Graphics saved as "reaction_twoway_midpoint_conservation.png" reaction_twoway_midpoint(): Normal end of execution. reaction_twoway_ode45(): Solve the two way reaction ODE using ode45() ode45() used 57 steps, of variable size. Graphics saved as "reaction_twoway_ode45_solution.png" Graphics saved as "reaction_twoway_ode45_conservation.png" reaction_twoway_ode45: Normal end of execution. reaction_twoway_trapezoidal: Solve the two way reaction ODE using trapezoidal() Number of steps will be 51 Stepsize will be 0.02 Graphics saved as "reaction_twoway_trapezoidal_solution.png" Graphics saved as "reaction_twoway_trapezoidal_conservation.png" reaction_twoway_trapezoidal: Normal end of execution. reaction_twoway_ode_test(): Normal end of execution. 08-Jan-2022 09:14:52