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