#naresh2

FIX: EXT

#reactions
v_1:
 {1.0}EXT = {1.0}S
  Q_0
v_10:
 {1.0}I2 = {1.0}A
  Delta_2*I2
v_11:
 {1.0}I2 = {1.0}EXT
  Mu*I2
v_12:
 {1.0}A = {1.0}EXT
  Alpha*A
v_13:
 {1.0}A = {1.0}EXT
  Mu*A
v_2:
 {1.0}S = {1.0}I1
  (Beta_1*S*I1)/(S+I1+I2+A)
v_3:
 {1.0}S = {1.0}I2
  (Epsilon*Beta_2*S*I2)/(S+I1+I2+A)
v_4:
 {1.0}S = {1.0}EXT
  Mu*S
v_5:
 {1.0}S = {1.0}I1
  ((1-Epsilon)*Beta_2*I2*S)/(S+I1+I2+A)
v_6:
 {1.0}I1 = {1.0}I2
  Theta*I1
v_7:
 {1.0}I1 = {1.0}A
  Delta_1*I1
v_8:
 {1.0}I1 = {1.0}EXT
  Mu*I1
v_9:
 {1.0}I1 = {1.0}I2
  k*I1*I2


#parameters
Alpha = 1.0
Beta_1 = 1.5
Beta_2 = 0.5
Delta_1 = 0.2
Delta_2 = 0.1
Epsilon = 0.01
Mu = 0.02
Q_0 = 2000.0
Theta = 0.015
k = 1e-05
EXT = 1.0
default = 1.0


#initial values
A = 200.0
I1 = 5000.0
I2 = 1000.0
S = 3800.0


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