07-Jan-2022 20:33:12

glycolysis_ode_test():
  MATLAB/Octave version 9.8.0.1380330 (R2020a) Update 2
  Solve the 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.6,1.36364)  0.126987  (-5.55112e-17,5.55112e-17)
  (0.913376,0.632359): (0.6,1.36364)  0.0975404  (-5.55112e-17,5.55112e-17)
  (0.278498,0.546882): (0.6,1.36364)  0.957507  (-5.55112e-17,5.55112e-17)
  (0.964889,0.157613): (0.6,1.36364)  0.970593  (-5.55112e-17,5.55112e-17)
  (0.957167,0.485376): (0.6,1.36364)  0.80028  (-5.55112e-17,5.55112e-17)

glycolysis_ode45()
  Use ode45() to solve the 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-Jan-2022 20:33:17