#murray1



#reactions
v1:
 $pool = {1.0}V
  NN * II
v10:
 {1.0}L = $pool
  muL * L
v2:
 $pool = {1.0}V
  Nm * M
v3:
 {1.0}V = $pool
  c * V
v4:
 $pool = {1.0}II
  (1 - eff) * (1 - effi) * kM * M
v5:
 $pool = {1.0}II
  (1 - eff) * (1 - effi) * kL * L
v6:
 {1.0}II = $pool
  delta * II
v7:
 $pool = {1.0}M
  (1 - eff) * (1 - effi) * k2 * V
v8:
 {1.0}M = $pool
  muM * M
v9:
 $pool = {1.0}L
  (1 - eff) * (1 - effi) * k2L * V


#parameters
NN = 300.0
Nm = 0.0
c = 3.0
delta = 0.69
eff = 0.985
effi = 0.82
k2 = 0.016
k2L = 0.00016
kL = 0.011
kM = 0.022
muL = 0.0035
muM = 0.05
default = 1.0


#initial values
II = 316.0
L = 1000.0
M = 10000.0
V = 32000.0


#assignment rules
!F HIV_DNA = M + L