begin name
vernoux2
end name

begin reactions
reacDIA {1.0}A + {1.0}I = {1.0}DIA
reacDIAxA {1.0}DIA = {1.0}A
reacDII {1.0}I + {1.0}I = {1.0}DII
reacDIIxI {1.0}DII = {1.0}I
reacdegrA {1.0}A = {1.0}$EXTERNAL
reacdegrDIA {1.0}DIA = {1.0}$EXTERNAL
reacdegrDII {1.0}DII = {1.0}$EXTERNAL
reacdegrI {1.0}I = {1.0}$EXTERNAL
reacdegrR {1.0}R = {1.0}$EXTERNAL
reacprodA {1.0}$EXTERNAL = {1.0}A
reacprodI {1.0}$EXTERNAL = {1.0}I
reacprodR {1.0}$EXTERNAL = {1.0}R
end reactions

begin rate equations
reacDIA = -(kprimeIA*DIA[t]) + kIA*A[t]*I[t]
reacDIAxA = (aux*dI*gammaI*Kaux*DIA[t])/(1 + aux*Kaux)
reacDII = -(kprimeII*DII[t]) + kII*I[t]^2
reacDIIxI = (aux*dI*gammaI*Kaux*DII[t])/(1 + aux*Kaux)
reacdegrA = dA*A[t]
reacdegrDIA = dIA*DIA[t]
reacdegrDII = dII*DII[t]
reacdegrI = (aux*dI*gammaI*Kaux*I[t])/(1 + aux*Kaux)
reacdegrR = dr*R[t]
reacprodA = piA
reacprodI = piI*R[t]
reacprodR = (1 + (fc*A[t]*(1 + (fA*wA*A[t])/Bd))/Bd)/(1 + kAm + (A[t]*(1 + (wA*A[t])/Bd))/Bd + (wD*DIA[t])/Bd + (wI*A[t]*I[t])/(Bd*KIA))
end rate equations

begin parameters
Bd = 100.0
KIA = 10.0
KII = 10.0
Kaux = 1.0
dA = 0.003
dI = 0.05
dIA = 0.003
dII = 0.003
dr = 0.007
fA = 10.0
fc = 10.0
gammaI = 10.0
kAm = 10.0
kIA = 1.0
kII = 1.0
piA = 1.0
piI = 1.0
wA = 10.0
wD = 10.0
wI = 10.0
compartment0000001 = 1.0
end parameters

begin initial conditions
A[0] = Ai
DIA[0] = DIAi
DII[0] = DIIi
I[0] = Ii
R[0] = Ri
end initial conditions

begin initial values
Ai = 10.0
DIAi = 10.0
DIIi = 10.0
Ii = 10.0
Ri = 1.0
end initial values

begin assignment rules
auxbasal := 1/((-1 + gammaI)*Kaux)
kprimeII := kII*KII
kprimeIA := kIA*KIA
aux := 1.11 + 1.11*Sin[0.007854*(-200 + 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