Sat Jan 29 18:55:24 2022

glycolysis_ode_test():
  Python version: 3.6.9
  Solve glycolysis_ode().

  parameters:
    alpha =  0.08
    gamma =  0.6
    t0    =  0.0
    y0    =  [0.9 0.7]
    tstop =  50.0

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

  (a,b) y_equi, t_equi, dydt

  ( 0.7582905939120552 , 0.06693142380987682 ):  [0.06693142 0.0877478 ] , 0.2516366981132845 , [ 3.25260652e-18 -3.25260652e-18]
  ( 0.2379434947351019 , 0.6738943774566316 ):  [0.67389438 0.97372728] , 0.7703803386788162 , [0. 0.]
  ( 0.3673394988588655 , 0.4103606472699515 ):  [0.41036065 0.76597641] , 0.23310218421045892 , [-2.77555756e-17  2.77555756e-17]
  ( 0.0393706807042129 , 0.5565746976564439 ):  [0.5565747  1.59410269] , 0.20713977936296202 , [0. 0.]
  ( 0.9837289591548201 , 0.4864003115256632 ):  [0.48640031 0.39858612] , 0.17351671521742973 , [ 2.77555756e-17 -2.77555756e-17]

glycolysis_solve_ivp()
  Use solve_ivp() to solve glycolysis_ode().
  Graphics saved as "glycolysis_solve_ivp_plot.png"
  Graphics saved as "glycolysis_solve_ivp_phase.png"

glycolysis_ode_test():
  Normal end of execution.
Sat Jan 29 18:55:25 2022