#vilar1

FIX: EmptySet

#reactions
Reaction1:
 {1.0}A + {1.0}R = {1.0}C
  A * R * gammaC
Reaction10:
 {1.0}MA = {1.0}A + {1.0}MA
  MA * betaA
Reaction11:
 {1.0}A + {1.0}DR = {1.0}DRp
  A * DR * gammaR
Reaction12:
 {1.0}DRp = {1.0}A + {1.0}DR
  DRp * thetaR
Reaction13:
 {1.0}DR = {1.0}DR + {1.0}MR
  DR * alphaR
Reaction14:
 {1.0}DRp = {1.0}DRp + {1.0}MR
  DRp * alphaRp
Reaction15:
 {1.0}MR = {1.0}EmptySet
  MR * deltaMR
Reaction16:
 {1.0}MR = {1.0}MR + {1.0}R
  MR * betaR
Reaction2:
 {1.0}A = {1.0}EmptySet
  A * deltaA
Reaction3:
 {1.0}C = {1.0}R
  C * deltaA
Reaction4:
 {1.0}R = {1.0}EmptySet
  R * deltaR
Reaction5:
 {1.0}A + {1.0}DA = {1.0}DAp
  A * DA * gammaA
Reaction6:
 {1.0}DAp = {1.0}A + {1.0}DA
  DAp * thetaA
Reaction7:
 {1.0}DA = {1.0}DA + {1.0}MA
  DA * alphaA
Reaction8:
 {1.0}DAp = {1.0}DAp + {1.0}MA
  DAp * alphaAp
Reaction9:
 {1.0}MA = {1.0}EmptySet
  MA * deltaMA


#parameters
EmptySet = 0.0
Reaction1_gammaC = 2.0
Reaction2_deltaA = 1.0
Reaction3_deltaA = 1.0
Reaction4_deltaR = 0.2
Reaction5_gammaA = 1.0
Reaction6_thetaA = 50.0
Reaction7_alphaA = 50.0
Reaction8_alphaAp = 500.0
Reaction9_deltaMA = 10.0
Reaction10_betaA = 50.0
Reaction11_gammaR = 1.0
Reaction12_thetaR = 100.0
Reaction13_alphaR = 0.01
Reaction14_alphaRp = 50.0
Reaction15_deltaMR = 0.5
Reaction16_betaR = 5.0
deterministicOscillator = 1.0


#initial values
A = 0.0
C = 0.0
DA = 1.0
DAp = 0.0
DR = 1.0
DRp = 0.0
MA = 0.0
MR = 0.0
R = 0.0


#assignment rules