#goldbeter1Demo2Cell

FIX: s

#reactions
v1:
 {1.0}s = {1.0}alpha
  k1
v12:
 {1.0}s = {1.0}alpha2
  k1
v2:
 {1.0}alpha = {1.0}gamma
  k2 * (alpha * es * (1 + alpha * es) * pow(1 + gamma, 2) + el * ve * alpha * c * ep * (1 + alpha * c * ep)) / (el * pow(1 + alpha * c * ep, 2) + pow(1 + gamma, 2) * pow(1 + alpha * es, 2))
v22:
 {1.0}alpha2 = {1.0}gamma2
  k2 * (alpha2 * es * (1 + alpha2 * es) * pow(1 + gamma2, 2) + el * ve * alpha2 * c * ep * (1 + alpha2 * c * ep)) / (el * pow(1 + alpha2 * c * ep, 2) + pow(1 + gamma2, 2) * pow(1 + alpha2 * es, 2))
v3:
 {1.0}gamma = {1.0}p
  k3 * (gamma - p)
v32:
 {1.0}gamma2 = {1.0}p
  k3 * (gamma2 - p)
v4:
 {1.0}p = $pool
  k4 * p


#parameters
EXTERNAL = 0.0
c = 1e-05
el = 1000000.0
ep = 0.9090909
es = 0.9090909
k1 = 0.7
k2 = 3.98
k3 = 0.07
k4 = 10.0
ve = 1.0
s = 0.0
default_compartment = 1.0


#initial values
alpha = 27.0
alpha2 = 84.0
gamma = 10.0
gamma2 = 2.0
p = 0.0


#assignment rules