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

07-Oct-2025 21:22:54

glycolysis_ode_test():
  MATLAB/Octave version 9.11.0.2358333 (R2021b) Update 7
  Solve glycolysis_ode().

  parameters:
    alpha = 0.08
    gamma = 0.6
    t0    = 0
    y0    = (0.9,0.7)
    tstop = 50

glycolysis_equilibrium_test()
  Verify that dy/dt=[0;0] for any parameter values.

  (a,b) y_equi, t_equi, dydt

  (0.814724,0.905792): (0.905792,0.553939)  0.126987  (5.55112e-17,-5.55112e-17)
  (0.913376,0.632359): (0.632359,0.481521)  0.0975404  (5.55112e-17,-5.55112e-17)
  (0.278498,0.546882): (0.546882,0.946854)  0.957507  (-5.55112e-17,5.55112e-17)
  (0.964889,0.157613): (0.157613,0.159248)  0.970593  (-2.08167e-17,2.08167e-17)
  (0.957167,0.485376): (0.485376,0.406936)  0.80028  (-2.77556e-17,2.77556e-17)

glycolysis_ode45()
  Use ode45() to solve glycolysis_ode().

  Number of variable steps = 237
  Graphics saved as "glycolysis_ode45_plot.png"
  Graphics saved as "glycolysis_ode45_phase.png"

glycolysis_ode_test():
  Normal end of execution.

07-Oct-2025 21:23:01