ExcitableCell

Systems · Continuous · Chem bio systems

Dimension: 3

Equations

\[ \begin{aligned} \dot{y_{0}} &= \frac{5.74594990367292 \cdot 10^{-6} gi \left(vi - y_{0}\right) \left(y_{0} + 25\right)^{3} e^{- 0.05 y_{0}}}{\left(1 - 0.0820849986238988 e^{- 0.1 y_{0}}\right)^{3} \left(0.00574594990367292 e^{- 0.05 y_{0}} + \frac{1}{e^{- 0.1 y_{0} - 2} + 1}\right) \left(4 e^{- \frac{y_{0}}{18} - \frac{25}{9}} + \frac{0.1 \left(y_{0} + 25\right)}{1 - 0.0820849986238988 e^{- 0.1 y_{0}}}\right)^{3}} + \frac{gkc y_{2} \left(vk - y_{0}\right)}{y_{2} + 1} + gkv y_{1}^{4} \left(vk - y_{0}\right) + gl \left(vl - y_{0}\right) \\ \dot{y_{1}} &= 230 \left(- y_{1} + \frac{0.01 \left(y_{0} + 20\right)}{\left(1 - e^{- 0.1 y_{0} - 2}\right) \left(0.125 e^{- \frac{y_{0}}{80} - \frac{3}{8}} + \frac{0.01 \left(y_{0} + 20\right)}{1 - e^{- 0.1 y_{0} - 2}}\right)}\right) \left(0.125 e^{- \frac{y_{0}}{80} - \frac{3}{8}} + \frac{0.01 \left(y_{0} + 20\right)}{1 - e^{- 0.1 y_{0} - 2}}\right) \\ \dot{y_{2}} &= \rho \left(- kc y_{2} + \frac{5.74594990367292 \cdot 10^{-6} \left(vc - y_{0}\right) \left(y_{0} + 25\right)^{3} e^{- 0.05 y_{0}}}{\left(1 - 0.0820849986238988 e^{- 0.1 y_{0}}\right)^{3} \left(0.00574594990367292 e^{- 0.05 y_{0}} + \frac{1}{e^{- 0.1 y_{0} - 2} + 1}\right) \left(4 e^{- \frac{y_{0}}{18} - \frac{25}{9}} + \frac{0.1 \left(y_{0} + 25\right)}{1 - 0.0820849986238988 e^{- 0.1 y_{0}}}\right)^{3}}\right) \end{aligned} \]

Parameters

parameter	default
gi	1800
gkc	11
gkv	1700
gl	7
kc	0.183333
rho	0.27
vc	100
vi	100
vk	-75
vl	-40
vm	-50
vn	-30

Usage

import tsdynamics as ts

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

Back to Chem bio systems