08-Jan-2022 09:14:31

reaction_ode_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Solve the reaction ODE.

  Reaction differential equation:
    da/dt = - k a b
    db/dt = - k a b
    dc/dt = + k a b

  Parameter values:
    k     = 0.02, (reaction rate)
    t0    = 0, (initial time, in seconds s)
    a0    = 10, (initial amount of species a in Mol/L)
    b0    = 12, (initial amount of species b in Mol/L)
    c0    = 3, (initial amount of species c in Mol/L)
    tstop = 20, (final time, in seconds s)

reaction_ode45():
  Use ode45() to solve the reaction ODE.

  Number of variable size steps = 49
  Graphics saved as "reaction_ode45_solution.png"
  Graphics saved as "reaction_ode45_conservation.png"

reaction_ode_test:
  Normal end of execution.

08-Jan-2022 09:14:37