#lavrentovich1



#reactions
R1:
 $pool = {1.0}X
  compartment * vin
R2:
 {1.0}X = $pool
  compartment * kout * X
R3:
 {1.0}Y = {1.0}X
  ER * 4 * vM3 * pow(k_CaA, n) * (pow(X, n) / ((pow(X, n) + pow(k_CaA, n)) * (pow(X, n) + pow(k_CaI, n)))) * (pow(Z, m) / (pow(Z, m) + pow(kip3, m))) * (Y - X)
R4:
 {1.0}X = {1.0}Y
  compartment * vM2 * pow(X, 2) / (pow(X, 2) + pow(k2, 2))
R5:
 {1.0}Y = {1.0}X
  ER * kf * (Y - X)
R6:
 $pool = {1.0}Z
  compartment * vp * pow(X, 2) / (pow(X, 2) + pow(kp, 2))
R7:
 {1.0}Z = $pool
  compartment * kdeg * Z


#parameters
k2 = 0.1
k_CaA = 0.15
k_CaI = 0.15
kdeg = 0.08
kf = 0.5
kip3 = 0.1
kout = 0.5
kp = 0.3
m = 2.2
n = 2.02
vM2 = 15.0
vM3 = 40.0
vin = 0.05
vp = 0.05
ER = 1.0
compartment = 1.0


#initial values
X = 0.1
Y = 1.5
Z = 0.1


#assignment rules