begin name
venturelli1
end name

begin reactions
v1 {1.0}$x1 = {1.0}xBH
v10 {1.0}$x10 = {1.0}xCH
v11 {1.0}$x11 = {1.0}xDP
v12 {1.0}$x12 = {1.0}xER
v2 {1.0}$x2 = {1.0}xCA
v3 {1.0}$x3 = {1.0}xBU
v4 {1.0}$x4 = {1.0}xPC
v5 {1.0}$x5 = {1.0}xBO
v6 {1.0}$x6 = {1.0}xBV
v7 {1.0}$x7 = {1.0}xBT
v8 {1.0}$x8 = {1.0}xEL
v9 {1.0}$x9 = {1.0}xFP
end reactions

begin rate equations
v1 = xBH[t]*(uBH + aBHBH*xBH[t] + aBHBO*xBO[t] + aBHBT*xBT[t] + aBHBU*xBU[t] + aBHBV*xBV[t] + aBHCA*xCA[t] + aBHCH*xCH[t] + aBHDP*xDP[t] + aBHEL*xEL[t] + aBHER*xER[t] + aBHFP*xFP[t] + aBHPC*xPC[t])
v10 = xCH[t]*(uCH + aCHBH*xBH[t] + aCHBO*xBO[t] + aCHBT*xBT[t] + aCHBU*xBU[t] + aCHBV*xBV[t] + aCHCA*xCA[t] + aCHCH*xCH[t] + aCHDP*xDP[t] + aCHEL*xEL[t] + aCHER*xER[t] + aCHFP*xFP[t] + aCHPC*xPC[t])
v11 = xDP[t]*(uDP + aDPBH*xBH[t] + aDPBO*xBO[t] + aDPBT*xBT[t] + aDPBU*xBU[t] + aDPBV*xBV[t] + aDPCA*xCA[t] + aDPCH*xCH[t] + aDPDP*xDP[t] + aDPEL*xEL[t] + aDPER*xER[t] + aDPFP*xFP[t] + aDPPC*xPC[t])
v12 = xER[t]*(uER + aERBH*xBH[t] + aERBO*xBO[t] + aERBT*xBT[t] + aERBU*xBU[t] + aERBV*xBV[t] + aERCA*xCA[t] + aERCH*xCH[t] + aERDP*xDP[t] + aEREL*xEL[t] + aERER*xER[t] + aERFP*xFP[t] + aERPC*xPC[t])
v2 = xCA[t]*(uCA + aCABH*xBH[t] + aCABO*xBO[t] + aCABT*xBT[t] + aCABU*xBU[t] + aCABV*xBV[t] + aCACA*xCA[t] + aCACH*xCH[t] + aCADP*xDP[t] + aCAEL*xEL[t] + aCAER*xER[t] + aCAFP*xFP[t] + aCAPC*xPC[t])
v3 = xBU[t]*(uBU + aBUBH*xBH[t] + aBUBO*xBO[t] + aBUBT*xBT[t] + aBUBU*xBU[t] + aBUBV*xBV[t] + aBUCA*xCA[t] + aBUCH*xCH[t] + aBUDP*xDP[t] + aBUEL*xEL[t] + aBUER*xER[t] + aBUFP*xFP[t] + aBUPC*xPC[t])
v4 = xPC[t]*(uPC + aPCBH*xBH[t] + aPCBO*xBO[t] + aPCBT*xBT[t] + aPCBU*xBU[t] + aPCBV*xBV[t] + aPCCA*xCA[t] + aPCCH*xCH[t] + aPCDP*xDP[t] + aPCEL*xEL[t] + aPCER*xER[t] + aPCFP*xFP[t] + aPCPC*xPC[t])
v5 = xBO[t]*(uBO + aBOBH*xBH[t] + aBOBO*xBO[t] + aBOBT*xBT[t] + aBOBU*xBU[t] + aBOBV*xBV[t] + aBOCA*xCA[t] + aBOCH*xCH[t] + aBODP*xDP[t] + aBOEL*xEL[t] + aBOER*xER[t] + aBOFP*xFP[t] + aBOPC*xPC[t])
v6 = xBV[t]*(uBV + aBVBH*xBH[t] + aBVBO*xBO[t] + aBVBT*xBT[t] + aBVBU*xBU[t] + aBVBV*xBV[t] + aBVCA*xCA[t] + aBVCH*xCH[t] + aBVDP*xDP[t] + aBVEL*xEL[t] + aBVER*xER[t] + aBVFP*xFP[t] + aBVPC*xPC[t])
v7 = xBT[t]*(uBT + aBTBH*xBH[t] + aBTBO*xBO[t] + aBTBT*xBT[t] + aBTBU*xBU[t] + aBTBV*xBV[t] + aBTCA*xCA[t] + aBTCH*xCH[t] + aBTDP*xDP[t] + aBTEL*xEL[t] + aBTER*xER[t] + aBTFP*xFP[t] + aBTPC*xPC[t])
v8 = xEL[t]*(uEL + aELBH*xBH[t] + aELBO*xBO[t] + aELBT*xBT[t] + aELBU*xBU[t] + aELBV*xBV[t] + aELCA*xCA[t] + aELCH*xCH[t] + aELDP*xDP[t] + aELEL*xEL[t] + aELER*xER[t] + aELFP*xFP[t] + aELPC*xPC[t])
v9 = xFP[t]*(uFP + aFPBH*xBH[t] + aFPBO*xBO[t] + aFPBT*xBT[t] + aFPBU*xBU[t] + aFPBV*xBV[t] + aFPCA*xCA[t] + aFPCH*xCH[t] + aFPDP*xDP[t] + aFPEL*xEL[t] + aFPER*xER[t] + aFPFP*xFP[t] + aFPPC*xPC[t])
end rate equations

begin parameters
aBHBH = -0.7905
aBHBO = 0.0
aBHBT = 0.0
aBHBU = 0.0
aBHBV = 0.0
aBHCA = 0.0
aBHCH = 0.0
aBHDP = 0.0
aBHEL = 0.0
aBHER = 0.0
aBHFP = 0.0
aBHPC = 0.0
aBOBH = 0.0
aBOBO = -1.082
aBOBT = 0.0
aBOBU = 0.0
aBOBV = 0.0
aBOCA = 0.0
aBOCH = 0.0
aBODP = 0.0
aBOEL = 0.0
aBOER = 0.0
aBOFP = 0.0
aBOPC = 0.0
aBTBH = 0.0
aBTBO = 0.0
aBTBT = -1.1639
aBTBU = 0.0
aBTBV = 0.0
aBTCA = 0.0
aBTCH = 0.0
aBTDP = 0.0
aBTEL = 0.0
aBTER = 0.0
aBTFP = 0.0
aBTPC = 0.0
aBUBH = 0.0
aBUBO = 0.0
aBUBT = 0.0
aBUBU = -1.1205
aBUBV = 0.0
aBUCA = 0.0
aBUCH = 0.0
aBUDP = 0.0
aBUEL = 0.0
aBUER = 0.0
aBUFP = 0.0
aBUPC = 0.0
aBVBH = 0.0
aBVBO = 0.0
aBVBT = 0.0
aBVBU = 0.0
aBVBV = -0.7535
aBVCA = 0.0
aBVCH = 0.0
aBVDP = 0.0
aBVEL = 0.0
aBVER = 0.0
aBVFP = 0.0
aBVPC = 0.0
aCABH = 0.0
aCABO = 0.0
aCABT = 0.0
aCABU = 0.0
aCABV = 0.0
aCACA = -0.7839
aCACH = 0.0
aCADP = 0.0
aCAEL = 0.0
aCAER = 0.0
aCAFP = 0.0
aCAPC = 0.0
aCHBH = 0.0
aCHBO = 0.0
aCHBT = 0.0
aCHBU = 0.0
aCHBV = 0.0
aCHCA = 0.0
aCHCH = -1.4921
aCHDP = 0.0
aCHEL = 0.0
aCHER = 0.0
aCHFP = 0.0
aCHPC = 0.0
aDPBH = 0.0
aDPBO = 0.0
aDPBT = 0.0
aDPBU = 0.0
aDPBV = 0.0
aDPCA = 0.0
aDPCH = 0.0
aDPDP = -0.6553
aDPEL = 0.0
aDPER = 0.0
aDPFP = 0.0
aDPPC = 0.0
aELBH = 0.0
aELBO = 0.0
aELBT = 0.0
aELBU = 0.0
aELBV = 0.0
aELCA = 0.0
aELCH = 0.0
aELDP = 0.0
aELEL = -0.6457
aELER = 0.0
aELFP = 0.0
aELPC = 0.0
aERBH = 0.0
aERBO = 0.0
aERBT = 0.0
aERBU = 0.0
aERBV = 0.0
aERCA = 0.0
aERCH = 0.0
aERDP = 0.0
aEREL = 0.0
aERER = -0.5581
aERFP = 0.0
aERPC = 0.0
aFPBH = 0.0
aFPBO = 0.0
aFPBT = 0.0
aFPBU = 0.0
aFPBV = 0.0
aFPCA = 0.0
aFPCH = 0.0
aFPDP = 0.0
aFPEL = 0.0
aFPER = 0.0
aFPFP = -0.5218
aFPPC = 0.0
aPCBH = 0.0
aPCBO = 0.0
aPCBT = 0.0
aPCBU = 0.0
aPCBV = 0.0
aPCCA = 0.0
aPCCH = 0.0
aPCDP = 0.0
aPCEL = 0.0
aPCER = 0.0
aPCFP = 0.0
aPCPC = -0.5146
uBH = 0.2972
uBO = 0.7332
uBT = 0.8777
uBU = 0.7689
uBV = 0.5834
uCA = 0.3454
uCH = 0.7103
uDP = 0.2082
uEL = 0.1715
uER = 0.0904
uFP = 0.1619
uPC = 0.1262
x1 = 1.0
x10 = 1.0
x11 = 1.0
x12 = 1.0
x2 = 1.0
x3 = 1.0
x4 = 1.0
x5 = 1.0
x6 = 1.0
x7 = 1.0
x8 = 1.0
x9 = 1.0
default = 1.0
end parameters

begin initial conditions
xBH[0] = xBHi
xBO[0] = xBOi
xBT[0] = xBTi
xBU[0] = xBUi
xBV[0] = xBVi
xCA[0] = xCAi
xCH[0] = xCHi
xDP[0] = xDPi
xEL[0] = xELi
xER[0] = xERi
xFP[0] = xFPi
xPC[0] = xPCi
end initial conditions

begin initial values
xBHi = 0.01
xBOi = 0.01
xBTi = 0.01
xBUi = 0.01
xBVi = 0.01
xCAi = 0.01
xCHi = 0.01
xDPi = 0.01
xELi = 0.01
xERi = 0.01
xFPi = 0.01
xPCi = 0.01
end initial values

begin assignment rules
fBU := xBU[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fBH := xBH[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fBV := xBV[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fPC := xPC[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fER := xER[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fBO := xBO[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fCH := xCH[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fBT := xBT[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fEL := xEL[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fFP := xFP[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fDP := xDP[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
fCA := xCA[t]/(xBH[t] + xBO[t] + xBT[t] + xBU[t] + xBV[t] + xCA[t] + xCH[t] + xDP[t] + xEL[t] + xER[t] + xFP[t] + xPC[t])
end assignment rules

begin function definitions
end function definitions

begin events
event48 = (t > 48, {xBO[t] -> 0.05*xBO[t], xBV[t] -> 0.05*xBV[t], xBT[t] -> 0.05*xBT[t], xEL[t] -> 0.05*xEL[t], xBH[t] -> 0.05*xBH[t], xCA[t] -> 0.05*xCA[t], xBU[t] -> 0.05*xBU[t], xPC[t] -> 0.05*xPC[t], xFP[t] -> 0.05*xFP[t], xCH[t] -> 0.05*xCH[t], xDP[t] -> 0.05*xDP[t], xER[t] -> 0.05*xER[t]}, 0)
event24 = (t > 24, {xBH[t] -> 0.05*xBH[t], xCA[t] -> 0.05*xCA[t], xBU[t] -> 0.05*xBU[t], xPC[t] -> 0.05*xPC[t], xBO[t] -> 0.05*xBO[t], xBV[t] -> 0.05*xBV[t], xBT[t] -> 0.05*xBT[t], xEL[t] -> 0.05*xEL[t], xFP[t] -> 0.05*xFP[t], xCH[t] -> 0.05*xCH[t], xDP[t] -> 0.05*xDP[t], xER[t] -> 0.05*xER[t]}, 0)
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units