begin name
perelson2
end name

begin reactions
CD4TproductionV {1.0}$EXTERNAL = {1.0}V
VinfectionCD4T {1.0}V + {1.0}T = {1.0}Tstar
activationCD4T {1.0}Tstar = {1.0}Tstarstar
clearanceV {1.0}V = {1.0}$EXTERNAL
deathCD4TactiveInfected {1.0}Ttot + {1.0}Tstarstar = {1.0}$EXTERNAL
deathCD4Tplain {1.0}T + {1.0}Ttot = {1.0}$EXTERNAL
deathCD4TrestingInfected {1.0}Ttot + {1.0}Tstar = {1.0}$EXTERNAL
growthCD4T {1.0}$EXTERNAL = {1.0}T + {1.0}Ttot
sourceCD4T {1.0}$EXTERNAL = {1.0}T + {1.0}Ttot
end reactions

begin rate equations
CD4TproductionV = muB*N0*Tstarstar[t]
VinfectionCD4T = k1*T[t]*V[t]
activationCD4T = k2*Tstar[t]
clearanceV = muV*V[t]
deathCD4TactiveInfected = muB*Tstarstar[t]
deathCD4Tplain = muT*T[t]
deathCD4TrestingInfected = muT*Tstar[t]
growthCD4T = r*T[t]*(1 - Ttot[t]/Tmax)
sourceCD4T = s
end rate equations

begin 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
end parameters

begin initial conditions
T[0] = Ti
Tstar[0] = Tstari
Tstarstar[0] = Tstarstari
Ttot[0] = Ttoti
V[0] = Vi
end initial conditions

begin initial values
Ti = 1000.0
Tstari = 0.0
Tstarstari = 0.0
Ttoti = 1000.0
Vi = 0.001
end initial values

begin assignment rules
end assignment rules

begin function definitions
end function definitions

begin events
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units