begin name
bachar1
end name

begin reactions
v1 {1.0}$EXT = {1.0}S
v2 {1.0}S = {1.0}I1
v3 {1.0}S = {1.0}$EXT
v4 {1.0}I1 = {1.0}$EXT
v5 {1.0}I1 = {1.0}I2
v6 {1.0}I2 = {1.0}I1
v7 {1.0}I2 = {1.0}$EXT
v8 {1.0}I2 = {1.0}$EXT
end reactions

begin rate equations
v1 = b*(I1[t] + I2[t] + S[t])
v2 = (Beta*c*(I1[t] + a*I2[t])*S[t])/(I1[t] + I2[t] + S[t])
v3 = d*S[t]
v4 = d*I1[t]
v5 = k1*I1[t]
v6 = Alpha*I2[t]
v7 = d*I2[t]
v8 = k2*I2[t]
end rate equations

begin parameters
Alpha = 0.2
Beta = 0.5
a = 0.6
b = 0.001
c = 0.3
d = 0.0006
k1 = 0.019
k2 = 0.0159
EXT = 1.0
default = 1.0
end parameters

begin initial conditions
I1[0] = I1i
I2[0] = I2i
S[0] = Si
end initial conditions

begin initial values
I1i = 50.0
I2i = 0.0
Si = 4950.0
end initial values

begin assignment rules
N := I1[t] + I2[t] + S[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