06-Jan-2023 20:23:41

glycolysis_ode_test():
  MATLAB/Octave version 4.2.2
  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.199287,0.165786): (0.165786,0.731067)  0.518014  (-2.08167e-17,2.08167e-17)
  (0.327928,0.778244): (0.778244,0.833602)  0.290005  (0,0)
  (0.539782,0.320628): (0.320628,0.498967)  0.181475  (4.16334e-17,-4.16334e-17)
  (0.59003,0.107386): (0.107386,0.178512)  0.00265105  (6.93889e-18,-6.93889e-18)
  (0.498811,0.1437): (0.1437,0.276632)  0.587842  (-7.80626e-18,7.80626e-18)

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

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

glycolysis_ode_test():
  Normal end of execution.

06-Jan-2023 20:23:41