Thu Feb 3 20:47:26 2022 reaction_ode_test(): Python version: 3.6.9 Solve reaction_ode(). Reaction differential equation: da/dt = - k a b db/dt = - k a b dc/dt = + k a b parameters: k = 0.02 , (reaction rate) t0 = 0.0 , (initial time, in seconds s) a0 = 10.0 , (initial amount of species a in Mol/L) b0 = 12.0 , (initial amount of species b in Mol/L) c0 = 3.0 , (initial amount of species c in Mol/L) tstop = 20.0 , (final time, in seconds s) reaction_solve_ivp(): Use solve_ivp() to solve reaction_ode(). Graphics saved as "reaction_solve_ivp_solution.png" Graphics saved as "reaction_solve_ivp_conservation.png" reaction_ode_test(): Normal end of execution. Thu Feb 3 20:47:26 2022