begin name
venturelli3
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.9241
aBHBO = -0.2199
aBHBT = -0.2641
aBHBU = -0.2284
aBHBV = -0.131
aBHCA = -0.312
aBHCH = -0.3284
aBHDP = -0.9182
aBHEL = 0.1821
aBHER = -0.5488
aBHFP = -0.2367
aBHPC = -0.5294
aBOBH = 0.0
aBOBO = -0.7134
aBOBT = -0.5769
aBOBU = -0.5064
aBOBV = -0.5089
aBOCA = 0.0
aBOCH = -0.4668
aBODP = 0.0
aBOEL = 1.8096
aBOER = -0.0391
aBOFP = -0.1758
aBOPC = -0.0001
aBTBH = 0.0582
aBTBO = -0.8237
aBTBT = -0.9597
aBTBU = -0.972
aBTBV = -0.6819
aBTCA = 0.0
aBTCH = -0.0727
aBTDP = 0.0
aBTEL = 1.9668
aBTER = 0.0
aBTFP = -0.4498
aBTPC = 0.0
aBUBH = 0.0
aBUBO = -0.7487
aBUBT = -0.9377
aBUBU = -0.9059
aBUBV = -0.5581
aBUCA = 0.0
aBUCH = 0.0564
aBUDP = 0.0
aBUEL = 3.0816
aBUER = 0.0
aBUFP = -0.8309
aBUPC = 0.0
aBVBH = 0.1337
aBVBO = -0.5768
aBVBT = -0.6133
aBVBU = -0.5773
aBVBV = -0.6578
aBVCA = -0.654
aBVCH = -0.0464
aBVDP = -0.108
aBVEL = 1.2848
aBVER = 0.0
aBVFP = -0.6349
aBVPC = 0.0
aCABH = 0.4301
aCABO = -0.283
aCABT = -0.2737
aCABU = -0.2612
aCABV = -0.1701
aCACA = -0.818
aCACH = 0.3033
aCADP = 0.0
aCAEL = -0.4766
aCAER = 0.0
aCAFP = -1.1212
aCAPC = -0.685
aCHBH = -0.3846
aCHBO = -0.4679
aCHBT = -0.5906
aCHBU = -0.2163
aCHBV = 0.0
aCHCA = -0.3186
aCHCH = -1.242
aCHDP = -2.3006
aCHEL = 0.0
aCHER = -0.5077
aCHFP = -0.6943
aCHPC = -0.7845
aDPBH = 0.0
aDPBO = -0.187
aDPBT = -0.1915
aDPBU = -0.1628
aDPBV = -0.0172
aDPCA = 0.0532
aDPCH = -0.1472
aDPDP = -1.3228
aDPEL = 0.1989
aDPER = 0.0
aDPFP = 0.0
aDPPC = -0.4407
aELBH = 0.8285
aELBO = -0.0932
aELBT = -0.099
aELBU = -0.0776
aELBV = -0.0229
aELCA = -0.9925
aELCH = 0.0
aELDP = -0.0003
aELEL = -2.6202
aELER = 0.0
aELFP = 0.0
aELPC = -1.091
aERBH = 1.4712
aERBO = 0.0
aERBT = -0.0315
aERBU = -0.0285
aERBV = -0.0178
aERCA = 0.0
aERCH = 1.1905
aERDP = 0.0
aEREL = 0.0
aERER = -1.3219
aERFP = -0.0319
aERPC = 0.0
aFPBH = 0.9091
aFPBO = -0.114
aFPBT = -0.0742
aFPBU = 0.1988
aFPBV = 0.7253
aFPCA = 0.0
aFPCH = 0.4306
aFPDP = 0.9762
aFPEL = -0.8005
aFPER = 0.0
aFPFP = -0.9946
aFPPC = -0.423
aPCBH = 0.0
aPCBO = -0.2718
aPCBT = -0.3041
aPCBU = -0.3241
aPCBV = -0.2039
aPCCA = -0.5827
aPCCH = 0.2406
aPCDP = -1.0009
aPCEL = -0.9852
aPCER = -0.9156
aPCFP = -0.4513
aPCPC = -0.6154
uBH = 0.2472
uBO = 0.4624
uBT = 0.6259
uBU = 0.5991
uBV = 0.4562
uCA = 0.2464
uCH = 0.4682
uDP = 0.2361
uEL = 0.4052
uER = 0.1507
uFP = 0.2164
uPC = 0.2381
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
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])
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])
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])
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])
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])
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
event24 = (t > 24, {xER[t] -> 0.05*xER[t], xDP[t] -> 0.05*xDP[t], xCH[t] -> 0.05*xCH[t], xEL[t] -> 0.05*xEL[t], xBT[t] -> 0.05*xBT[t], xBV[t] -> 0.05*xBV[t], xPC[t] -> 0.05*xPC[t], xBU[t] -> 0.05*xBU[t], xBH[t] -> 0.05*xBH[t], xCA[t] -> 0.05*xCA[t], xBO[t] -> 0.05*xBO[t], xFP[t] -> 0.05*xFP[t]}, 0)
event48 = (t > 48, {xER[t] -> 0.05*xER[t], xDP[t] -> 0.05*xDP[t], xFP[t] -> 0.05*xFP[t], xEL[t] -> 0.05*xEL[t], xBT[t] -> 0.05*xBT[t], xBO[t] -> 0.05*xBO[t], xPC[t] -> 0.05*xPC[t], xCA[t] -> 0.05*xCA[t], xBH[t] -> 0.05*xBH[t], xBU[t] -> 0.05*xBU[t], xBV[t] -> 0.05*xBV[t], xCH[t] -> 0.05*xCH[t]}, 0)
end events

begin process annotations
end process annotations

begin species annotations
end species annotations

begin units
end units