begin name
venturelli4
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.9118
aBHBO = -0.2145
aBHBT = -0.2718
aBHBU = -0.2275
aBHBV = -0.1294
aBHCA = -0.3058
aBHCH = -0.3478
aBHDP = -0.9002
aBHEL = 0.1764
aBHER = -0.5455
aBHFP = -0.2307
aBHPC = -0.5286
aBOBH = 0.0
aBOBO = -0.7339
aBOBT = -0.6235
aBOBU = -0.6317
aBOBV = -0.515
aBOCA = 0.0
aBOCH = -0.5069
aBODP = 0.0
aBOEL = 1.7566
aBOER = 0.0
aBOFP = -0.2087
aBOPC = 0.0
aBTBH = 0.0
aBTBO = -0.8193
aBTBT = -0.9067
aBTBU = -0.7538
aBTBV = -0.7552
aBTCA = 0.0
aBTCH = -0.0864
aBTDP = 0.0
aBTEL = 2.2719
aBTER = -0.7375
aBTFP = -0.7033
aBTPC = 0.0
aBUBH = 0.0
aBUBO = -0.9208
aBUBT = -0.8155
aBUBU = -0.8804
aBUBV = -0.5498
aBUCA = 0.0
aBUCH = 0.0656
aBUDP = 0.0
aBUEL = 3.3782
aBUER = 0.0
aBUFP = -0.7822
aBUPC = 0.0
aBVBH = 0.1366
aBVBO = -0.5556
aBVBT = -0.6416
aBVBU = -0.5837
aBVBV = -0.6597
aBVCA = -0.6566
aBVCH = -0.0468
aBVDP = -0.1062
aBVEL = 1.3026
aBVER = 0.0
aBVFP = -0.6379
aBVPC = 0.0
aCABH = 0.4526
aCABO = -0.2776
aCABT = -0.2736
aCABU = -0.2614
aCABV = -0.1679
aCACA = -0.829
aCACH = 0.3118
aCADP = 0.0
aCAEL = -0.4475
aCAER = 0.0
aCAFP = -1.121
aCAPC = -0.6709
aCHBH = 0.0
aCHBO = -0.4645
aCHBT = -0.6319
aCHBU = -0.1511
aCHBV = 0.039
aCHCA = -0.2413
aCHCH = -1.4543
aCHDP = -2.157
aCHEL = 0.0239
aCHER = -0.4394
aCHFP = -0.5069
aCHPC = -0.771
aDPBH = 0.0
aDPBO = -0.2028
aDPBT = -0.1999
aDPBU = -0.1763
aDPBV = -0.026
aDPCA = 0.042
aDPCH = -0.151
aDPDP = -1.2535
aDPEL = 0.1756
aDPER = 0.0
aDPFP = 0.0
aDPPC = -0.4333
aELBH = 0.6924
aELBO = -0.1141
aELBT = -0.1688
aELBU = -0.1241
aELBV = -0.0493
aELCA = -1.098
aELCH = 0.0
aELDP = -0.4084
aELEL = -2.4418
aELER = 0.0
aELFP = -0.1531
aELPC = -1.0774
aERBH = 1.3433
aERBO = -0.0203
aERBT = 0.0
aERBU = -0.0613
aERBV = -0.0296
aERCA = 0.0
aERCH = 1.0831
aERDP = 0.0
aEREL = -0.1385
aERER = -1.2705
aERFP = -0.1676
aERPC = 0.0
aFPBH = 0.9613
aFPBO = -0.0993
aFPBT = -0.0685
aFPBU = 0.2313
aFPBV = 0.7585
aFPCA = 0.0
aFPCH = 0.4481
aFPDP = 1.0147
aFPEL = -0.768
aFPER = 0.0
aFPFP = -1.0382
aFPPC = -0.402
aPCBH = 0.0
aPCBO = -0.2646
aPCBT = -0.3033
aPCBU = -0.3237
aPCBV = -0.2017
aPCCA = -0.5596
aPCCH = 0.265
aPCDP = -0.9771
aPCEL = -0.9041
aPCER = -0.8171
aPCFP = -0.4053
aPCPC = -0.6217
uBH = 0.2453
uBO = 0.4777
uBT = 0.5979
uBU = 0.5841
uBV = 0.4573
uCA = 0.2464
uCH = 0.5025
uDP = 0.2321
uEL = 0.4021
uER = 0.1558
uFP = 0.2192
uPC = 0.2375
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
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])
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])
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])
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
event24 = (t > 24, {xER[t] -> 0.05*xER[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], 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]}, 0)
event48 = (t > 48, {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