begin name
nyabadza3
end name

begin reactions
v1 {1.0}$EXT = {1.0}Sn
v10 {1.0}Ins = {1.0}An
v11 {1.0}Ss = {1.0}Is
v12 {1.0}Is = {1.0}$EXT
v13 {1.0}Is = {1.0}It
v14 {1.0}Is = {1.0}As
v15 {1.0}It = {1.0}$EXT
v16 {1.0}It = {1.0}As
v17 {1.0}An = {1.0}$EXT
v18 {1.0}An = {1.0}$EXT
v19 {1.0}As = {1.0}$EXT
v2 {1.0}Sn = {1.0}Ins
v20 {1.0}As = {1.0}$EXT
v3 {1.0}Sn = {1.0}$EXT
v4 {1.0}Sn = {1.0}Ss
v5 {1.0}$EXT = {1.0}Ss
v6 {1.0}Ss = {1.0}Ins
v7 {1.0}Ss = {1.0}$EXT
v8 {1.0}Ins = {1.0}$EXT
v9 {1.0}Ins = {1.0}Is
end reactions

begin rate equations
v1 = Pin
v10 = Rho1*Ins[t]
v11 = Lambda*p*(1 - Psi)*Ss[t]
v12 = Mu*Is[t]
v13 = Sigma*Is[t]
v14 = Rho2*Is[t]
v15 = Mu*It[t]
v16 = Rho3*It[t]
v17 = Mu*An[t]
v18 = Delta1*An[t]
v19 = Mu*As[t]
v2 = Lambda*Sn[t]
v20 = Delta2*As[t]
v3 = Mu*Sn[t]
v4 = k*Sn[t]
v5 = Pis
v6 = Lambda*(1 - Psi)*q*Ss[t]
v7 = Mu*Ss[t]
v8 = Mu*Ins[t]
v9 = k*Ins[t]
end rate equations

begin parameters
Beta = 0.4584
Delta1 = 0.293
Delta2 = 0.3379
Eta1 = 0.7203
Eta2 = 1.3785
Mu = 0.032
Phi1 = 0.1
Phi2 = "7.3471e-05"
Pin = 920540.0
Pis = 750000.0
Psi = 0.2915
Rho1 = 0.2372
Rho2 = 0.1828
Rho3 = 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
end parameters

begin initial conditions
An[0] = Ani
As[0] = Asi
Ins[0] = Insi
Is[0] = Isi
It[0] = Iti
Sn[0] = Sni
Ss[0] = Ssi
end initial conditions

begin initial values
Ani = 20000.0
Asi = 0.0
Insi = 145000.0
Isi = 0.0
Iti = 0.0
Sni = 18373000.0
Ssi = 0.0
end initial values

begin assignment rules
Lambda := (Beta*c*(1 - Theta)*(Eta2*(An[t] + Phi2*As[t]) + Ins[t] + Eta1*(Is[t] + Phi1*It[t])))/(E^((m*(Delta1*An[t] + Delta2*As[t]))/(An[t] + As[t] + Ins[t] + Is[t] + It[t] + Sn[t] + Ss[t]))*(An[t] + As[t] + Ins[t] + Is[t] + It[t] + Sn[t] + Ss[t]))
peoplewithHIV := Ins[t] + Is[t] + It[t]
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