#bachar1

FIX: EXT

#reactions
v_1:
 {1.0}EXT = {1.0}S
  b*(S+I1+I2)
v_2:
 {1.0}S = {1.0}I1
  c*Beta*(I1+a*I2)*S/(S+I1+I2)
v_3:
 {1.0}S = {1.0}EXT
  d*S
v_4:
 {1.0}I1 = {1.0}EXT
  d*I1
v_5:
 {1.0}I1 = {1.0}I2
  k1*I1
v_6:
 {1.0}I2 = {1.0}I1
  Alpha*I2
v_7:
 {1.0}I2 = {1.0}EXT
  d*I2
v_8:
 {1.0}I2 = {1.0}EXT
  k2*I2


#parameters
Alpha = 0.2
Beta = 0.5
a = 0.6
b = 0.001
c = 0.3
d = 0.0006
k1 = 0.019
k2 = 0.0159
EXT = 1.0
default = 1.0


#initial values
I1 = 50.0
I2 = 0.0
S = 4950.0


#assignment rules
!F N = S + I1 +I2