Home License -- for personal use only.  Not for government,
	academic, research, commercial, or other organizational use.

13-May-2025 17:09:33

reaction_twoway_ode_test():
  MATLAB/Octave version 9.11.0.2358333 (R2021b) Update 7
  Solve reaction_twoway_ode().

  Reaction differential equation:
    dW1/dt = - k1 W1 + k2 W2
    dW2/dt = + k1 W1 - k2 W2

  parameters:
    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 reaction_twoway_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_midpoint:
  Solve reaction_twoway_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_ode45():
  Solve reaction_twoway_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_trapezoidal:
  Solve reaction_twoway_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_ode_test():
  Normal end of execution.

13-May-2025 17:09:46