begin name
li2
end name

begin reactions
v1 {1.0}$EXTERNAL = {1.0}H
v10 {1.0}$EXTERNAL = {1.0}E
v11 {1.0}$EXTERNAL = {1.0}E
v2 {1.0}H = {1.0}$EXTERNAL
v3 {1.0}H = {1.0}I
v4 {1.0}I = {1.0}$EXTERNAL
v5 {1.0}I = {1.0}$EXTERNAL
v6 {1.0}$EXTERNAL = {1.0}M
v7 {1.0}M = {1.0}$EXTERNAL
v8 {1.0}M = {1.0}$EXTERNAL
v9 {1.0}E = {1.0}$EXTERNAL
end reactions

begin rate equations
v1 = lambda
v10 = (k1*E[t]*I[t])/(1 + beta*I[t])
v11 = (k2*E[t]*M[t])/(1 + gamma*M[t])
v2 = d1*H[t]
v3 = alfa*H[t]*M[t]
v4 = delta*I[t]
v5 = (p1*E[t]*I[t])/(1 + beta*I[t])
v6 = r*I[t]
v7 = mu*M[t]
v8 = (p2*E[t]*M[t])/(1 + gamma*M[t])
v9 = d2*E[t]
end rate equations

begin parameters
alfa = "9e-07"
beta = 0.0005
d1 = 0.0083
d2 = 0.04
delta = 1.0
gamma = 0.000667
k1 = "2.5e-05"
k2 = "1.03e-05"
lambda = 41500.0
mu = 48.0
p1 = "1e-08"
p2 = "1e-08"
r = 12.0
default = 1.0
end parameters

begin initial conditions
E[0] = Ei
H[0] = Hi
I[0] = Ii
M[0] = Mi
end initial conditions

begin initial values
Ei = 0.0001
Hi = 5000000.0
Ii = 0.0
Mi = 10000.0
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