#nyabadza3

FIX: EXT

#reactions
v_1:
 {1.0}EXT = {1.0}Sn
  Pi_n
v_10:
 {1.0}Ins = {1.0}An
  Rho_1 * Ins
v_11:
 {1.0}Ss = {1.0}Is
  p * (1 - Psi) * Lambda * Ss
v_12:
 {1.0}Is = {1.0}EXT
  Mu * Is
v_13:
 {1.0}Is = {1.0}It
  Sigma * Is
v_14:
 {1.0}Is = {1.0}As
  Rho_2 * Is
v_15:
 {1.0}It = {1.0}EXT
  Mu * It
v_16:
 {1.0}It = {1.0}As
  Rho_3 * It
v_17:
 {1.0}An = {1.0}EXT
  Mu * An
v_18:
 {1.0}An = {1.0}EXT
  Delta_1 * An
v_19:
 {1.0}As = {1.0}EXT
  Mu * As
v_2:
 {1.0}Sn = {1.0}Ins
  Lambda * Sn
v_20:
 {1.0}As = {1.0}EXT
  Delta_2 * As
v_3:
 {1.0}Sn = {1.0}EXT
  Mu * Sn
v_4:
 {1.0}Sn = {1.0}Ss
  k * Sn
v_5:
 {1.0}EXT = {1.0}Ss
  Pi_s
v_6:
 {1.0}Ss = {1.0}Ins
  q * (1 - Psi) * Lambda * Ss
v_7:
 {1.0}Ss = {1.0}EXT
  Mu * Ss
v_8:
 {1.0}Ins = {1.0}EXT
  Mu * Ins
v_9:
 {1.0}Ins = {1.0}Is
  k * Ins


#parameters
Beta = 0.4584
Delta_1 = 0.293
Delta_2 = 0.3379
Eta_1 = 0.7203
Eta_2 = 1.3785
Mu = 0.032
Phi_1 = 0.1
Phi_2 = 7.3471e-05
Pi_n = 920540.0
Pi_s = 750000.0
Psi = 0.2915
Rho_1 = 0.2372
Rho_2 = 0.1828
Rho_3 = 0.13
Sigma = 0.0655
Theta = 0.1605
c = 1.5892
k = 0.0388
m = 45.1202
p = 0.1
q = 0.9
EXT = 1.0
default = 1.0


#initial values
An = 20000.0
As = 0.0
Ins = 145000.0
Is = 0.0
It = 0.0
Sn = 18373000.0
Ss = 0.0


#assignment rules
!F Lambda = c*(1.0 - Theta)*Beta*exp ((-m*(Delta_1*An + Delta_2*As))/(Sn + Ss + Ins + Is + It + An + As))*((Ins + Eta_1*(Is + Phi_1*It) + Eta_2*(An + Phi_2*As))/(Sn + Ss + Ins + Is + It + An + As))
!F people_with_HIV = Ins + Is +It