#perelson2



#reactions
CD4TproductionV:
 $pool = {1.0}V
  muB * N0 * Tstarstar
VinfectionCD4T:
 {1.0}V + {1.0}T = {1.0}Tstar
  k1 * T * V
activationCD4T:
 {1.0}Tstar = {1.0}Tstarstar
  k2 * Tstar
clearanceV:
 {1.0}V = $pool
  muV * V
deathCD4TactiveInfected:
 {1.0}Ttot + {1.0}Tstarstar = $pool
  muB * Tstarstar
deathCD4Tplain:
 {1.0}T + {1.0}Ttot = $pool
  muT * T
deathCD4TrestingInfected:
 {1.0}Ttot + {1.0}Tstar = $pool
  muT * Tstar
growthCD4T:
 $pool = {1.0}T + {1.0}Ttot
  r * T * (1 - Ttot / Tmax)
sourceCD4T:
 $pool = {1.0}T + {1.0}Ttot
  s


#parameters
N0 = 1000.0
Tmax = 1500.0
k1 = 2.4e-05
k2 = 0.003
muB = 0.24
muT = 0.02
muV = 2.4
r = 0.03
s = 10.0
theta = 1.0
default = 1.0


#initial values
T = 1000.0
Tstar = 0.0
Tstarstar = 0.0
Ttot = 1000.0
V = 0.001


#assignment rules