GlycolyticOscillation

Systems · Continuous · Chem bio systems

Dimension: 3

Equations

\[ \begin{aligned} \dot{y_{0}} &= \nu - \frac{s_{1} y_{0} \left(y_{0} + 1\right) \left(y_{1} + 1\right)^{2}}{l_{1} + \left(y_{0} + 1\right)^{2} \left(y_{1} + 1\right)^{2}} \\ \dot{y_{1}} &= \frac{q_{1} s_{1} y_{0} \left(y_{0} + 1\right) \left(y_{1} + 1\right)^{2}}{l_{1} + \left(y_{0} + 1\right)^{2} \left(y_{1} + 1\right)^{2}} - \frac{s_{2} y_{1} \left(y_{2} + 1\right)^{2} \left(d y_{1} + 1\right)}{l_{2} + \left(y_{2} + 1\right)^{2} \left(d y_{1} + 1\right)^{2}} \\ \dot{y_{2}} &= - k y_{2} + \frac{q_{2} s_{2} y_{1} \left(y_{2} + 1\right)^{2} \left(d y_{1} + 1\right)}{l_{2} + \left(y_{2} + 1\right)^{2} \left(d y_{1} + 1\right)^{2}} \end{aligned} \]

Parameters

parameter	default
d	0.0
k	4.422
l1	500000000.0
l2	100
nu	1.0
q1	50
q2	0.02
s1	22.2222
s2	22.2222

Usage

import tsdynamics as ts

sys = ts.GlycolyticOscillation()
traj = sys.integrate(final_time=100.0, dt=0.01)

Back to Chem bio systems