begin name
rong3
end name

begin reactions
v1 {1.0}$EXTERNAL = {1.0}T
v10 {1.0}T1 = {1.0}L
v11 {1.0}L = {1.0}T1
v12 {1.0}$EXTERNAL = {1.0}V + {1.0}Vi
v13 {1.0}Vi = {1.0}Vni
v2 {1.0}T = {1.0}$EXTERNAL
v3 {1.0}T = {1.0}T1
v4 {1.0}L = {1.0}$EXTERNAL
v5 {1.0}L = {1.0}T1
v6 {1.0}T1 = {1.0}$EXTERNAL
v7 {1.0}V + {1.0}Vi = {1.0}$EXTERNAL
v8 {1.0}V + {1.0}Vni = {1.0}$EXTERNAL
v9 {1.0}T1 = {1.0}T
end reactions

begin rate equations
v1 = lambda
v10 = alphaL*k*T[t]*Vi[t]
v11 = alphaL*epsilonRT*k*T[t]*Vi[t]
v12 = delta*NN*T1[t]
v13 = delta*epsilonPI*NN*T1[t]
v2 = dT*T[t]
v3 = k*T[t]*Vi[t]
v4 = dL*L[t]
v5 = a*L[t]
v6 = delta*T1[t]
v7 = c*Vi[t]
v8 = c*Vni[t]
v9 = epsilonRT*k*T[t]*Vi[t]
end rate equations

begin parameters
NN = 4000.0
a = 0.1
alphaL = 0.001
c = 23.0
dL = 0.004
dT = 0.01
delta = 1.0
epsilonPI = 0.65
epsilonRT = 0.0
k = "2.4e-08"
lambda = 10000.0
default = 1.0
end parameters

begin initial conditions
L[0] = Li
T[0] = Ti
T1[0] = T1i
V[0] = Vi
Vi[0] = Vii
Vni[0] = Vnii
end initial conditions

begin initial values
Li = 2.0
Ti = 600000.0
T1i = 0.3
Vi = 50000.0
Vii = 50000.0
Vnii = 0.0
end initial values

begin assignment rules
ratio := Log[Vni[t]/Vi[t]]/Log[10]
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