#komarova1



#reactions
R1:
 $pool = {1.0}x1
  Ratex1appearance + alpha1*(1 - Fig7AR1alphaOnOff + Fig7AR1alphaOnOff*x1**g11*x2**g21)
R2:
 {1.0}x1 = $pool
  beta1*x1
R3:
 $pool = {1.0}x2
  alpha2*x1**g12*x2**g22
R4:
 {1.0}x2 = $pool
  beta2*x2
R5:
 {1.0}z = $pool
  k1*y1
R6:
 $pool = {1.0}z
  k2*y2


#parameters
Fig7AR1alphaOnOff = 1.0
Ratex1appearanceOnOff = 0.0
Ratex1appearanceOnOffPiecewise = 0.0
Ratex1appearanceduration = 1.0
Ratex1appearancemagnitude = 7.0
Ratex1appearancetimefirst = 10.0
Ratex1appearancetimesecond = 150.0
alpha1 = 7.0
alpha2 = 7.0
beta1basalvalue = 0.2
beta1increaseOnOff = 1.0
beta1increaseOntime = 500.0
beta1increasevalue = 0.0
beta1increasevalue2 = 0.03
beta2 = 0.02
flagformation = 1.0
flagresorption = 1.0
g11 = 1.09
g12 = 1.0
g21 = -0.5
g22 = 0.09
initialvaluex1 = 5.08098
initialvaluex2 = 1397.73
initialvaluez = 90.0
k1 = 0.26
k2 = 0.0008
plotzyaxisdenominator = 113.749
default = 1.0


#initial values
x1 = 5.08098
x2 = 1397.73
z = 90.0


#assignment rules
!F y1 = flagresorption*(-x1bar + x1)
!F x2bar = (beta1/alpha1)^(g12/gamma)*(beta2/alpha2)^((1.0 - g11)/gamma)
!F x1variable = x1
!F Ratex1appearance = Ratex1appearancemagnitude*Ratex1appearanceOnOff*Ratex1appearanceOnOffPiecewise
!F y2 = flagformation*(-x2bar + x2)
!F x2variable = x2
!F gamma = g12*g21 - (1.0 - g11)*(1.0 - g22)
!F beta1 = beta1basalvalue + beta1increaseOnOff*beta1increasevalue
!F plotzyaxis = (100.0*z)/plotzyaxisdenominator
!F x1bar = (beta1/alpha1)^((1.0 - g22)/gamma)*(beta2/alpha2)^(g21/gamma)
!F plotx2yaxis = -initialvaluex2 + y2 + x2