begin name
dutoit3
end name

begin reactions
vPFvALD {1.0}f16bpPF = {1.0}gapPF + {1.0}dhapPF
vPFvATPASE {1.0}atpPF = {1.0}adpPF
vPFvENO {1.0}p2gPF = {1.0}pepPF
vPFvG3PDH {1.0}nadhPF + {1.0}dhapPF = {1.0}nadPF + {1.0}g3pPF
vPFvGAPDH {1.0}nadPF + {1.0}gapPF = {1.0}nadhPF + {1.0}b13pgPF
vPFvGLCtr {1.0}glcRBCi = {1.0}glcPF
vPFvGLYtr {1.0}g3pPF = {1.0}$EXTERNAL
vPFvHK {1.0}glcPF + {1.0}atpPF = {1.0}g6pPF + {1.0}adpPF
vPFvLACtr {1.0}lacPF = {1.0}lacRBCi
vPFvLDH {1.0}pyrPF + {1.0}nadhPF = {1.0}nadPF + {1.0}lacPF
vPFvPFK {1.0}f6pPF + {1.0}atpPF = {1.0}f16bpPF + {1.0}adpPF
vPFvPGI {1.0}g6pPF = {1.0}f6pPF
vPFvPGK {1.0}b13pgPF + {1.0}adpPF = {1.0}p3gPF + {1.0}atpPF
vPFvPGM {1.0}p3gPF = {1.0}p2gPF
vPFvPK {1.0}pepPF + {1.0}adpPF = {1.0}pyrPF + {1.0}atpPF
vPFvPYRtr {1.0}pyrPF = {1.0}pyrRBCi
vPFvTPI {1.0}dhapPF = {1.0}gapPF
vRBCivAK {1.0}mgadpRBCi + {1.0}adpRBCi = {1.0}mgatpRBCi + {1.0}ampRBCi
vRBCivALD {1.0}f16p2RBCi = {1.0}gapRBCi + {1.0}dhapRBCi
vRBCivATPASE {1.0}mgatpRBCi = {1.0}phosRBCi + {1.0}mgadpRBCi
vRBCivBPGSP1 {1.0}bpgspRBCi + {1.0}b13pgRBCi = {1.0}bpgspb13pgRBCi
vRBCivBPGSP2 {1.0}bpgspb13pgRBCi = {1.0}p3gRBCi + {1.0}bpgsppRBCi
vRBCivBPGSP3 {1.0}p3gRBCi + {1.0}bpgsppRBCi = {1.0}bpgsppp3gRBCi
vRBCivBPGSP4 {1.0}p2gRBCi + {1.0}bpgsppRBCi = {1.0}bpgsppp2gRBCi
vRBCivBPGSP5 {1.0}bpgsppp3gRBCi = {1.0}bpgspb23pgRBCi
vRBCivBPGSP6 {1.0}bpgsppp2gRBCi = {1.0}bpgspb23pgRBCi
vRBCivBPGSP7 {1.0}bpgspb23pgRBCi = {1.0}bpgspRBCi + {1.0}b23pgRBCi
vRBCivBPGSP8 {1.0}phosRBCi + {1.0}bpgsppRBCi = {1.0}bpgsppphosRBCi
vRBCivBPGSP9 {1.0}bpgsppphosRBCi = {2.0}phosRBCi + {1.0}bpgspRBCi
vRBCivENO {1.0}p2gRBCi = {1.0}pepRBCi
vRBCivG6PDH {1.0}nadpRBCi + {1.0}g6pRBCi = {1.0}p6glRBCi + {1.0}nadphRBCi
vRBCivGAPDH {1.0}phosRBCi + {1.0}nadRBCi + {1.0}gapRBCi = {1.0}nadhRBCi + {1.0}b13pgRBCi
vRBCivGLCTRANSPORT {1.0}$EXTERNAL = {1.0}glcRBCi
vRBCivGSSGR {1.0}nadphRBCi + {1.0}gssgRBCi = {1.0}nadpRBCi + {2.0}gshRBCi
vRBCivHBADP {1.0}hbRBCi + {1.0}adpRBCi = {1.0}hbadpRBCi
vRBCivHBATP {1.0}hbRBCi + {1.0}atpRBCi = {1.0}hbatpRBCi
vRBCivHBB13PG {1.0}hbRBCi + {1.0}b13pgRBCi = {1.0}hbb13pgRBCi
vRBCivHBB23PG {1.0}hbRBCi + {1.0}b23pgRBCi = {1.0}hbb23pgRBCi
vRBCivHBMGATP {1.0}mgatpRBCi + {1.0}hbRBCi = {1.0}hbmgatpRBCi
vRBCivHK {1.0}mgatpRBCi + {1.0}glcRBCi = {1.0}mgadpRBCi + {1.0}g6pRBCi
vRBCivLACTRANSPORT {1.0}lacRBCi = {1.0}$EXTERNAL
vRBCivLDH {1.0}pyrRBCi + {1.0}nadhRBCi = {1.0}nadRBCi + {1.0}lacRBCi
vRBCivLDHP {1.0}pyrRBCi + {1.0}nadphRBCi = {1.0}nadpRBCi + {1.0}lacRBCi
vRBCivMGADP {1.0}mgRBCi + {1.0}adpRBCi = {1.0}mgadpRBCi
vRBCivMGATP {1.0}atpRBCi + {1.0}mgRBCi = {1.0}mgatpRBCi
vRBCivMGB13PG {1.0}mgRBCi + {1.0}b13pgRBCi = {1.0}mgb13pgRBCi
vRBCivMGB23PG {1.0}mgRBCi + {1.0}b23pgRBCi = {1.0}mgb23pgRBCi
vRBCivMGF16P2 {1.0}mgRBCi + {1.0}f16p2RBCi = {1.0}mgf16p2RBCi
vRBCivMGG16P2 {1.0}mgRBCi + {1.0}g16p2RBCi = {1.0}mgg16p2RBCi
vRBCivMGPHOS {1.0}phosRBCi + {1.0}mgRBCi = {1.0}mgphosRBCi
vRBCivOX {2.0}gshRBCi = {1.0}gssgRBCi
vRBCivOXNADH {1.0}nadhRBCi = {1.0}nadRBCi
vRBCivP6GDH {1.0}p6gRBCi + {1.0}nadpRBCi = {1.0}ru5pRBCi + {1.0}nadphRBCi
vRBCivPFK {1.0}mgatpRBCi + {1.0}f6pRBCi = {1.0}mgadpRBCi + {1.0}f16p2RBCi
vRBCivPGI {1.0}g6pRBCi = {1.0}f6pRBCi
vRBCivPGK {1.0}mgadpRBCi + {1.0}b13pgRBCi = {1.0}p3gRBCi + {1.0}mgatpRBCi
vRBCivPGLHYDROLYSIS {1.0}p6glRBCi = {1.0}p6gRBCi
vRBCivPGM {1.0}p3gRBCi = {1.0}p2gRBCi
vRBCivPHOSTRANSPORT {1.0}phosRBCi = {1.0}$EXTERNAL
vRBCivPK {1.0}pepRBCi + {1.0}mgadpRBCi = {1.0}pyrRBCi + {1.0}mgatpRBCi
vRBCivPYRTRANSPORT {1.0}pyrRBCi = {1.0}$EXTERNAL
vRBCivR5PI {1.0}ru5pRBCi = {1.0}rib5pRBCi
vRBCivRu5PE {1.0}ru5pRBCi = {1.0}xu5pRBCi
vRBCivTA {1.0}sed7pRBCi + {1.0}gapRBCi = {1.0}f6pRBCi + {1.0}ery4pRBCi
vRBCivTIM {1.0}gapRBCi = {1.0}dhapRBCi
vRBCivTK1 {1.0}xu5pRBCi + {1.0}tkRBCi = {1.0}tkxu5pRBCi
vRBCivTK2 {1.0}tkxu5pRBCi = {1.0}tkgRBCi + {1.0}gapRBCi
vRBCivTK3 {1.0}tkgRBCi + {1.0}rib5pRBCi = {1.0}tkgrib5pRBCi
vRBCivTK4 {1.0}tkgrib5pRBCi = {1.0}tkRBCi + {1.0}sed7pRBCi
vRBCivTK5 {1.0}tkgRBCi + {1.0}ery4pRBCi = {1.0}tkgery4pRBCi
vRBCivTK6 {1.0}tkgery4pRBCi = {1.0}tkRBCi + {1.0}f6pRBCi
vRBCvAK {1.0}mgadpRBC + {1.0}adpRBC = {1.0}mgatpRBC + {1.0}ampRBC
vRBCvALD {1.0}f16p2RBC = {1.0}gapRBC + {1.0}dhapRBC
vRBCvATPASE {1.0}mgatpRBC = {1.0}phosRBC + {1.0}mgadpRBC
vRBCvBPGSP1 {1.0}bpgspRBC + {1.0}b13pgRBC = {1.0}bpgspb13pgRBC
vRBCvBPGSP2 {1.0}bpgspb13pgRBC = {1.0}p3gRBC + {1.0}bpgsppRBC
vRBCvBPGSP3 {1.0}p3gRBC + {1.0}bpgsppRBC = {1.0}bpgsppp3gRBC
vRBCvBPGSP4 {1.0}p2gRBC + {1.0}bpgsppRBC = {1.0}bpgsppp2gRBC
vRBCvBPGSP5 {1.0}bpgsppp3gRBC = {1.0}bpgspb23pgRBC
vRBCvBPGSP6 {1.0}bpgsppp2gRBC = {1.0}bpgspb23pgRBC
vRBCvBPGSP7 {1.0}bpgspb23pgRBC = {1.0}bpgspRBC + {1.0}b23pgRBC
vRBCvBPGSP8 {1.0}phosRBC + {1.0}bpgsppRBC = {1.0}bpgsppphosRBC
vRBCvBPGSP9 {1.0}bpgsppphosRBC = {2.0}phosRBC + {1.0}bpgspRBC
vRBCvENO {1.0}p2gRBC = {1.0}pepRBC
vRBCvG6PDH {1.0}nadpRBC + {1.0}g6pRBC = {1.0}p6glRBC + {1.0}nadphRBC
vRBCvGAPDH {1.0}phosRBC + {1.0}nadRBC + {1.0}gapRBC = {1.0}nadhRBC + {1.0}b13pgRBC
vRBCvGLCTRANSPORT {1.0}$EXTERNAL = {1.0}glcRBC
vRBCvGSSGR {1.0}nadphRBC + {1.0}gssgRBC = {1.0}nadpRBC + {2.0}gshRBC
vRBCvHBADP {1.0}hbRBC + {1.0}adpRBC = {1.0}hbadpRBC
vRBCvHBATP {1.0}hbRBC + {1.0}atpRBC = {1.0}hbatpRBC
vRBCvHBB13PG {1.0}hbRBC + {1.0}b13pgRBC = {1.0}hbb13pgRBC
vRBCvHBB23PG {1.0}hbRBC + {1.0}b23pgRBC = {1.0}hbb23pgRBC
vRBCvHBMGATP {1.0}mgatpRBC + {1.0}hbRBC = {1.0}hbmgatpRBC
vRBCvHK {1.0}mgatpRBC + {1.0}glcRBC = {1.0}mgadpRBC + {1.0}g6pRBC
vRBCvLACTRANSPORT {1.0}lacRBC = {1.0}$EXTERNAL
vRBCvLDH {1.0}pyrRBC + {1.0}nadhRBC = {1.0}nadRBC + {1.0}lacRBC
vRBCvLDHP {1.0}pyrRBC + {1.0}nadphRBC = {1.0}nadpRBC + {1.0}lacRBC
vRBCvMGADP {1.0}mgRBC + {1.0}adpRBC = {1.0}mgadpRBC
vRBCvMGATP {1.0}mgRBC + {1.0}atpRBC = {1.0}mgatpRBC
vRBCvMGB13PG {1.0}mgRBC + {1.0}b13pgRBC = {1.0}mgb13pgRBC
vRBCvMGB23PG {1.0}mgRBC + {1.0}b23pgRBC = {1.0}mgb23pgRBC
vRBCvMGF16P2 {1.0}mgRBC + {1.0}f16p2RBC = {1.0}mgf16p2RBC
vRBCvMGG16P2 {1.0}mgRBC + {1.0}g16p2RBC = {1.0}mgg16p2RBC
vRBCvMGPHOS {1.0}phosRBC + {1.0}mgRBC = {1.0}mgphosRBC
vRBCvOX {2.0}gshRBC = {1.0}gssgRBC
vRBCvOXNADH {1.0}nadhRBC = {1.0}nadRBC
vRBCvP6GDH {1.0}p6gRBC + {1.0}nadpRBC = {1.0}ru5pRBC + {1.0}nadphRBC
vRBCvPFK {1.0}mgatpRBC + {1.0}f6pRBC = {1.0}mgadpRBC + {1.0}f16p2RBC
vRBCvPGI {1.0}g6pRBC = {1.0}f6pRBC
vRBCvPGK {1.0}mgadpRBC + {1.0}b13pgRBC = {1.0}p3gRBC + {1.0}mgatpRBC
vRBCvPGLHYDROLYSIS {1.0}p6glRBC = {1.0}p6gRBC
vRBCvPGM {1.0}p3gRBC = {1.0}p2gRBC
vRBCvPHOSTRANSPORT {1.0}phosRBC = {1.0}$EXTERNAL
vRBCvPK {1.0}pepRBC + {1.0}mgadpRBC = {1.0}pyrRBC + {1.0}mgatpRBC
vRBCvPYRTRANSPORT {1.0}pyrRBC = {1.0}$EXTERNAL
vRBCvR5PI {1.0}ru5pRBC = {1.0}rib5pRBC
vRBCvRu5PE {1.0}ru5pRBC = {1.0}xu5pRBC
vRBCvTA {1.0}sed7pRBC + {1.0}gapRBC = {1.0}f6pRBC + {1.0}ery4pRBC
vRBCvTIM {1.0}gapRBC = {1.0}dhapRBC
vRBCvTK1 {1.0}xu5pRBC + {1.0}tkRBC = {1.0}tkxu5pRBC
vRBCvTK2 {1.0}tkxu5pRBC = {1.0}tkgRBC + {1.0}gapRBC
vRBCvTK3 {1.0}tkgRBC + {1.0}rib5pRBC = {1.0}tkgrib5pRBC
vRBCvTK4 {1.0}tkgrib5pRBC = {1.0}tkRBC + {1.0}sed7pRBC
vRBCvTK5 {1.0}tkgRBC + {1.0}ery4pRBC = {1.0}tkgery4pRBC
vRBCvTK6 {1.0}tkgery4pRBC = {1.0}tkRBC + {1.0}f6pRBC
end reactions

begin rate equations
vPFvALD = (parMulti*VPFvALD*f16bpPF[t]*(1 - (dhapPF[t]*gapPF[t])/(KeqPFvALD*Vpf*f16bpPF[t])))/(Kf16bpPFvALD*(1 + dhapPF[t]/(KdhapPFvALD*Vpf) + f16bpPF[t]/(Kf16bpPFvALD*Vpf) + gapPF[t]/(KgapPFvALD*Vpf) + (dhapPF[t]*gapPF[t])/(KdhapPFvALD*KgapPFvALD*Vpf^2)))
vPFvATPASE = (parMulti*VPFvATPASE*atpPF[t]^5)/(KmPFvATPASE^5*Vpf^4*(1 + atpPF[t]^5/(KmPFvATPASE^5*Vpf^5)))
vPFvENO = (parMulti*VPFvENO*p2gPF[t]*(1 - pepPF[t]/(KeqPFvENO*p2gPF[t])))/(Kp2gPFvENO*(1 + p2gPF[t]/(Kp2gPFvENO*Vpf) + pepPF[t]/(KpepPFvENO*Vpf)))
vPFvG3PDH = (parMulti*VfPFvG3PDH*dhapPF[t]*nadhPF[t]*(1 - (g3pPF[t]*nadPF[t])/(KeqPFvG3PDH*dhapPF[t]*nadhPF[t])))/(KdhapPFvG3PDH*KnadhPFvG3PDH*Vpf*(1 + dhapPF[t]/(KdhapPFvG3PDH*Vpf) + g3pPF[t]/(Kg3pPFvG3PDH*Vpf))*(1 + nadhPF[t]/(KnadhPFvG3PDH*Vpf) + nadPF[t]/(KnadPFvG3PDH*Vpf)))
vPFvGAPDH = (parMulti*(-((VrPFvGAPDH*b13pgPF[t]*nadhPF[t])/(Kb13pgPFvGAPDH*KnadhPFvGAPDH*Vpf)) + (VPFvGAPDH*gapPF[t]*nadPF[t])/(KgapPFvGAPDH*KnadPFvGAPDH*Vpf)))/((1 + b13pgPF[t]/(Kb13pgPFvGAPDH*Vpf) + gapPF[t]/(KgapPFvGAPDH*Vpf))*(1 + nadhPF[t]/(KnadhPFvGAPDH*Vpf) + nadPF[t]/(KnadPFvGAPDH*Vpf)))
vPFvGLCtr = (parMulti*(-((VPFvGLCtr*glcPF[t])/((1 + cyt/ki1)*KmPFvGLCtr)) + (Vpf*VPFvGLCtr*glcRBCi[t])/((1 + cyt/ki1)*KmPFvGLCtr*Vrbci)))/(1 + ((1 + cyt/ki2)*glcPF[t])/((1 + cyt/ki1)*KmPFvGLCtr*Vpf) + ((1 + cyt/ki2)*glcRBCi[t])/((1 + cyt/ki1)*KmPFvGLCtr*Vrbci) + (alphaPFvGLCtr*glcPF[t]*glcRBCi[t])/((1 + cyt/ki1)*KmPFvGLCtr^2*Vpf*Vrbci))
vPFvGLYtr = kPFvGLYtr*parMulti*Vpf*(-(glyEXT/vBld) + g3pPF[t]/Vpf)
vPFvHK = (parMulti*VPFvHK*atpPF[t]*(1 - (adpPF[t]*g6pPF[t])/(KeqPFvHK*atpPF[t]*glcPF[t]))*glcPF[t])/(KatpPFvHK*KglcPFvHK*Vpf*(1 + adpPF[t]/(KadpPFvHK*Vpf) + atpPF[t]/(KatpPFvHK*Vpf))*(1 + atpPF[t]/(KiatpPFvHK*Vpf))*(1 + g6pPF[t]/(Kg6pPFvHK*Vpf) + glcPF[t]/(KglcPFvHK*Vpf)))
vPFvLACtr = parMulti*(KdPFvLACtr*lacPF[t]*(1 - (hRBC*Vpf*lacRBCi[t])/(hPF*Vrbci*lacPF[t])) + (VPFvLACtr*lacPF[t]*(1 - (hRBC*Vpf*lacRBCi[t])/(hPF*Vrbci*lacPF[t])))/(KlacPFvLACtr*(1 + lacPF[t]/(KlacPFvLACtr*Vpf) + lacRBCi[t]/(KlacPFvLACtr*Vrbci) + pyrPF[t]/(KipyrPFvLACtr*Vpf) + pyrRBCi[t]/(KipyrPFvLACtr*Vrbci))))
vPFvLDH = (parMulti*(-((VrPFvLDH*lacPF[t]*nadPF[t])/(KlacPFvLDH*KnadPFvLDH*Vpf)) + (VfPFvLDH*nadhPF[t]*pyrPF[t])/(KnadhPFvLDH*KpyrPFvLDH*Vpf)))/((1 + nadhPF[t]/(KnadhPFvLDH*Vpf) + nadPF[t]/(KnadPFvLDH*Vpf))*(1 + lacPF[t]/(KlacPFvLDH*Vpf) + pyrPF[t]/(KpyrPFvLDH*Vpf)))
vPFvPFK = (parMulti*VPFvPFK*atpPF[t]*f6pPF[t])/(KatpPFvPFK*Kf6pPFvPFK*Vpf*(1 + atpPF[t]/(KatpPFvPFK*Vpf))*(1 + adpPF[t]/(KadpPFvPFK*Vpf) + atpPF[t]/(KatpPFvPFK*Vpf))*(1 + f16bpPF[t]/(Kf16bpPFvPFK*Vpf) + f6pPF[t]/(Kf6pPFvPFK*Vpf)))
vPFvPGI = (parMulti*VPFvPGI*(1 - f6pPF[t]/(KeqPFvPGI*g6pPF[t]))*g6pPF[t])/(Kg6pPFvPGI*(1 + f6pPF[t]/(Kf6pPFvPGI*Vpf) + g6pPF[t]/(Kg6pPFvPGI*Vpf)))
vPFvPGK = (KadpPFvPGK*Kb13pgPFvPGK*KeqPFvPGK*parMulti*VrPFvPGK*atpPF[t]*(-1 + (KeqPFvPGK*adpPF[t]*b13pgPF[t])/(atpPF[t]*p3gPF[t]))*p3gPF[t])/(KatpPFvPGK^2*Kp3gPFvPGK^2*Vpf*(1 + adpPF[t]/(KadpPFvPGK*Vpf) + atpPF[t]/(KatpPFvPGK*Vpf))*(1 + b13pgPF[t]/(Kb13pgPFvPGK*Vpf) + p3gPF[t]/(Kp3gPFvPGK*Vpf)))
vPFvPGM = (parMulti*VPFvPGM*(1 - p2gPF[t]/(KeqPFvPGM*p3gPF[t]))*p3gPF[t])/(Kp3gPFvPGM*(1 + p2gPF[t]/(Kp2gPFvPGM*Vpf) + p3gPF[t]/(Kp3gPFvPGM*Vpf)))
vPFvPK = (parMulti*VPFvPK*adpPF[t]*pepPF[t])/(KadpPFvPK*KpepPFvPK*Vpf*(1 + adpPF[t]/(KadpPFvPK*Vpf))*(1 + pepPF[t]/(KpepPFvPK*Vpf))*(1 + LPFvPK*(1 + (adpPF[t]/(KiadpPFvPK*Vpf) + atpPF[t]/(KiatpPFvPK*Vpf))^h)*(1 + (pepPF[t]/(KipepPFvPK*Vpf) + pyrPF[t]/(KipyrPFvPK*Vpf))^h)))
vPFvPYRtr = parMulti*(KdPFvPYRtr*pyrPF[t]*(1 - (hRBC*Vpf*pyrRBCi[t])/(hPF*Vrbci*pyrPF[t])) + (VPFvPYRtr*pyrPF[t]*(1 - (hRBC*Vpf*pyrRBCi[t])/(hPF*Vrbci*pyrPF[t])))/(KpyrPFvPYRtr*(1 + lacPF[t]/(KilacPFvPYRtr*Vpf) + lacRBCi[t]/(KilacPFvPYRtr*Vrbci) + pyrPF[t]/(KpyrPFvPYRtr*Vpf) + pyrRBCi[t]/(KpyrPFvPYRtr*Vrbci))))
vPFvTPI = (parMulti*VPFvTPI*dhapPF[t]*(1 - gapPF[t]/(KeqPFvTPI*dhapPF[t])))/(KdhapPFvTPI*(1 + dhapPF[t]/(KdhapPFvTPI*Vpf) + gapPF[t]/(KgapPFvTPI*Vpf) + pepPF[t]/(KipepPFvTPI*Vpf)))
vRBCivAK = parMulti*TsRBC*Vrbci*((K1appRBCvAK*adpRBCi[t]*mgadpRBCi[t])/Vrbci^2 - (K2appRBCvAK*ampRBCi[t]*mgatpRBCi[t])/Vrbci^2)
vRBCivALD = (parMulti*TsRBC*Vrbci*VRBCvALD*((KcatfRBCvALD*f16p2RBCi[t])/(Kmf16p2RBCvALD*Vrbci) - (KcatrRBCvALD*dhapRBCi[t]*gapRBCi[t])/(KidhapRBCvALD*KmgapRBCvALD*Vrbci^2)))/(1 + dhapRBCi[t]/(KidhapRBCvALD*Vrbci) + f16p2RBCi[t]/(Kmf16p2RBCvALD*Vrbci) + (dhapRBCi[t]*gapRBCi[t])/(KidhapRBCvALD*KmgapRBCvALD*Vrbci^2) + (KmdhapRBCvALD*f16p2RBCi[t]*gapRBCi[t])/(KidhapRBCvALD*Kif16p2RBCvALD*KmgapRBCvALD*Vrbci^2) + (b23pgRBCi[t]/Vrbci + mgb23pgRBCi[t]/Vrbci)/Kib23pgRBCvALD + (KmdhapRBCvALD*gapRBCi[t]*(1 + (b23pgRBCi[t]/Vrbci + mgb23pgRBCi[t]/Vrbci)/Kib23pgRBCvALD))/(KidhapRBCvALD*KmgapRBCvALD*Vrbci))
vRBCivATPASE = KRBCvATPASE*parMulti*TsRBC*mgatpRBCi[t]
vRBCivBPGSP1 = parMulti*TsRBC*Vrbci*(-((K2RBCvBPGSP1*bpgspb13pgRBCi[t])/Vrbci) + (K1appRBCvBPGSP1*b13pgRBCi[t]*bpgspRBCi[t])/Vrbci^2)
vRBCivBPGSP2 = K3appRBCvBPGSP2*parMulti*TsRBC*bpgspb13pgRBCi[t]
vRBCivBPGSP3 = parMulti*TsRBC*Vrbci*(-((K5RBCvBPGSP3*bpgsppp3gRBCi[t])/Vrbci) + (K4appRBCvBPGSP3*bpgsppRBCi[t]*p3gRBCi[t])/Vrbci^2)
vRBCivBPGSP4 = parMulti*TsRBC*Vrbci*(-((K7RBCvBPGSP4*bpgsppp2gRBCi[t])/Vrbci) + (K6appRBCvBPGSP4*bpgsppRBCi[t]*p2gRBCi[t])/Vrbci^2)
vRBCivBPGSP5 = parMulti*TsRBC*Vrbci*(-((K9RBCvBPGSP5*bpgspb23pgRBCi[t])/Vrbci) + (K8RBCvBPGSP5*bpgsppp3gRBCi[t])/Vrbci)
vRBCivBPGSP6 = parMulti*TsRBC*Vrbci*(-((K11RBCvBPGSP6*bpgspb23pgRBCi[t])/Vrbci) + (K10RBCvBPGSP6*bpgsppp2gRBCi[t])/Vrbci)
vRBCivBPGSP7 = parMulti*TsRBC*Vrbci*((K12RBCvBPGSP7*bpgspb23pgRBCi[t])/Vrbci - (K13appRBCvBPGSP7*b23pgRBCi[t]*bpgspRBCi[t])/Vrbci^2)
vRBCivBPGSP8 = parMulti*TsRBC*Vrbci*(-((K15RBCvBPGSP8*bpgsppphosRBCi[t])/Vrbci) + (K14RBCvBPGSP8*bpgsppRBCi[t]*phosRBCi[t])/Vrbci^2)
vRBCivBPGSP9 = K16RBCvBPGSP9*parMulti*TsRBC*bpgsppphosRBCi[t]
vRBCivENO = (parMulti*TsRBC*Vrbci*VRBCvENO*((KcatfRBCvENO*mgRBCi[t]*p2gRBCi[t])/(KimgRBCvENO*Kmp2gRBCvENO*Vrbci^2) - (KcatrRBCvENO*mgRBCi[t]*pepRBCi[t])/(KipepRBCvENO*KmmgRBCvENO*Vrbci^2)))/(1 + mgRBCi[t]/(KimgRBCvENO*Vrbci) + p2gRBCi[t]/(Kip2gRBCvENO*Vrbci) + (mgRBCi[t]*p2gRBCi[t])/(KimgRBCvENO*Kmp2gRBCvENO*Vrbci^2) + pepRBCi[t]/(KipepRBCvENO*Vrbci) + (mgRBCi[t]*pepRBCi[t])/(KipepRBCvENO*KmmgRBCvENO*Vrbci^2))
vRBCivG6PDH = (parMulti*TsRBC*Vrbci*VRBCvG6PDH*((K1RBCvG6PDH*K3RBCvG6PDH*K5RBCvG6PDH*K7RBCvG6PDH*K9RBCvG6PDH*g6pRBCi[t]*nadpRBCi[t])/Vrbci^2 - (K10RBCvG6PDH*K2RBCvG6PDH*K4RBCvG6PDH*K6RBCvG6PDH*K8RBCvG6PDH*nadphRBCi[t]*p6glRBCi[t])/Vrbci^2))/(K2RBCvG6PDH*(K4RBCvG6PDH*K6RBCvG6PDH + K4RBCvG6PDH*K7RBCvG6PDH + K5RBCvG6PDH*K7RBCvG6PDH)*K9RBCvG6PDH + (K3RBCvG6PDH*K5RBCvG6PDH*K7RBCvG6PDH*K9RBCvG6PDH*g6pRBCi[t])/Vrbci + (K10RBCvG6PDH*K2RBCvG6PDH*(K4RBCvG6PDH*K6RBCvG6PDH + K5RBCvG6PDH*K6RBCvG6PDH + K5RBCvG6PDH*K7RBCvG6PDH)*nadphRBCi[t])/Vrbci + (K10RBCvG6PDH*K3RBCvG6PDH*K5RBCvG6PDH*K7RBCvG6PDH*g6pRBCi[t]*nadphRBCi[t])/Vrbci^2 + (K1RBCvG6PDH*(K4RBCvG6PDH*K6RBCvG6PDH + K4RBCvG6PDH*K7RBCvG6PDH + K5RBCvG6PDH*K7RBCvG6PDH)*K9RBCvG6PDH*nadpRBCi[t])/Vrbci + (K1RBCvG6PDH*K3RBCvG6PDH*(K5RBCvG6PDH*K7RBCvG6PDH + K5RBCvG6PDH*K9RBCvG6PDH + K6RBCvG6PDH*K9RBCvG6PDH + K7RBCvG6PDH*K9RBCvG6PDH)*g6pRBCi[t]*nadpRBCi[t])/Vrbci^2 + (K2RBCvG6PDH*K4RBCvG6PDH*K6RBCvG6PDH*K8RBCvG6PDH*p6glRBCi[t])/Vrbci + (K10RBCvG6PDH*(K2RBCvG6PDH*K4RBCvG6PDH + K2RBCvG6PDH*K5RBCvG6PDH + K2RBCvG6PDH*K6RBCvG6PDH + K4RBCvG6PDH*K6RBCvG6PDH)*K8RBCvG6PDH*nadphRBCi[t]*p6glRBCi[t])/Vrbci^2 + (K10RBCvG6PDH*K3RBCvG6PDH*(K5RBCvG6PDH + K6RBCvG6PDH)*K8RBCvG6PDH*g6pRBCi[t]*nadphRBCi[t]*p6glRBCi[t])/Vrbci^3 + (K1RBCvG6PDH*K4RBCvG6PDH*K6RBCvG6PDH*K8RBCvG6PDH*nadpRBCi[t]*p6glRBCi[t])/Vrbci^2 + (K1RBCvG6PDH*K3RBCvG6PDH*(K5RBCvG6PDH + K6RBCvG6PDH)*K8RBCvG6PDH*g6pRBCi[t]*nadpRBCi[t]*p6glRBCi[t])/Vrbci^3)
vRBCivGAPDH = (parMulti*TsRBC*Vrbci*VRBCvGAPDH*(-((KcatrappRBCvGAPDH*b13pgRBCi[t]*nadhRBCi[t])/(Kiappb13pgRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbci^2)) + (KcatfappRBCvGAPDH*gapRBCi[t]*nadRBCi[t]*phosRBCi[t])/(KiappgapRBCvGAPDH*KiphosRBCvGAPDH*KmnadRBCvGAPDH*Vrbci^3)))/((b13pgRBCi[t]*(1 + gapRBCi[t]/(KidgapRBCvGAPDH*Vrbci)))/(Kiappb13pgRBCvGAPDH*Vrbci) + (gapRBCi[t]*(1 + gapRBCi[t]/(KidgapRBCvGAPDH*Vrbci)))/(KiappgapRBCvGAPDH*Vrbci) + (Kmb13pgRBCvGAPDH*nadhRBCi[t])/(Kiappb13pgRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbci) + (b13pgRBCi[t]*nadhRBCi[t])/(Kiappb13pgRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbci^2) + (gapRBCi[t]*nadhRBCi[t])/(KiappgapRBCvGAPDH*KiappnadhRBCvGAPDH*Vrbci^2) + (b13pgRBCi[t]*nadRBCi[t])/(Kiappb13pgRBCvGAPDH*KinadRBCvGAPDH*Vrbci^2) + (gapRBCi[t]*nadRBCi[t])/(KiappgapRBCvGAPDH*KinadRBCvGAPDH*Vrbci^2) + (gapRBCi[t]*(1 + gapRBCi[t]/(KidgapRBCvGAPDH*Vrbci))*phosRBCi[t])/(KiappgapRBCvGAPDH*KiphosRBCvGAPDH*Vrbci^2) + (Kmb13pgRBCvGAPDH*nadhRBCi[t]*phosRBCi[t])/(Kiappb13pgRBCvGAPDH*KiphosRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbci^2) + (Kmb13pgRBCvGAPDH*b13pgRBCi[t]*nadhRBCi[t]*phosRBCi[t])/(Kiappb13pgRBCvGAPDH*Kidb13pgRBCvGAPDH*KiphosRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbci^3) + (gapRBCi[t]*nadhRBCi[t]*phosRBCi[t])/(KiappgapRBCvGAPDH*KiappnadhRBCvGAPDH*KiphosRBCvGAPDH*Vrbci^3) + (KmgapRBCvGAPDH*nadRBCi[t]*phosRBCi[t])/(KiappgapRBCvGAPDH*KiphosRBCvGAPDH*KmnadRBCvGAPDH*Vrbci^2) + (KmgapRBCvGAPDH*b13pgRBCi[t]*nadRBCi[t]*phosRBCi[t])/(KiappgapRBCvGAPDH*Kidb13pgRBCvGAPDH*KiphosRBCvGAPDH*KmnadRBCvGAPDH*Vrbci^3) + (gapRBCi[t]*nadRBCi[t]*phosRBCi[t])/(KiappgapRBCvGAPDH*KiphosRBCvGAPDH*KmnadRBCvGAPDH*Vrbci^3))
vRBCivGLCTRANSPORT = parMulti*TsRBC*Vrbci*((glcAlpha*glcEXT*VRBCvGLCTRANSPORT)/((kAlpha + (glcAlpha*glcEXT)/vBld + ((1 - glcAlpha)*glcEXT*kAlpha)/(kBeta*vBld))*vBld) + ((1 - glcAlpha)*glcEXT*VRBCvGLCTRANSPORT)/((kBeta + ((1 - glcAlpha)*glcEXT)/vBld + (glcAlpha*glcEXT*kBeta)/(kAlpha*vBld))*vBld) - (glcAlpha*VRBCvGLCTRANSPORT*glcRBCi[t])/(Vrbci*(kAlpha + (glcAlpha*glcRBCi[t])/Vrbci + ((1 - glcAlpha)*kAlpha*glcRBCi[t])/(kBeta*Vrbci))) - ((1 - glcAlpha)*VRBCvGLCTRANSPORT*glcRBCi[t])/(Vrbci*(kBeta + ((1 - glcAlpha)*glcRBCi[t])/Vrbci + (glcAlpha*kBeta*glcRBCi[t])/(kAlpha*Vrbci))))
vRBCivGSSGR = (parMulti*TsRBC*Vrbci*VRBCvGSSGR*((K11RBCvGSSGR*K1RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR*gssgRBCi[t]*nadphRBCi[t])/Vrbci^2 - (K10RBCvGSSGR*K12RBCvGSSGR*K2RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBCi[t]^2*nadpRBCi[t])/Vrbci^3))/(K11RBCvGSSGR*K2RBCvGSSGR*(K4RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K7RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR)*K9RBCvGSSGR + (K11RBCvGSSGR*K2RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBCi[t])/Vrbci + (K10RBCvGSSGR*K2RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBCi[t]^2)/Vrbci^2 + (K11RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR*gssgRBCi[t])/Vrbci + (K11RBCvGSSGR*K1RBCvGSSGR*(K4RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K7RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR)*K9RBCvGSSGR*nadphRBCi[t])/Vrbci + (K11RBCvGSSGR*K1RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBCi[t]*nadphRBCi[t])/Vrbci^2 + (K10RBCvGSSGR*K1RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBCi[t]^2*nadphRBCi[t])/Vrbci^3 + (K1RBCvGSSGR*K3RBCvGSSGR*(K11RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR + K11RBCvGSSGR*K5RBCvGSSGR*K9RBCvGSSGR + K11RBCvGSSGR*K6RBCvGSSGR*K9RBCvGSSGR + K11RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR)*gssgRBCi[t]*nadphRBCi[t])/Vrbci^2 + (K10RBCvGSSGR*K1RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*gshRBCi[t]*gssgRBCi[t]*nadphRBCi[t])/Vrbci^3 + (K11RBCvGSSGR*K1RBCvGSSGR*K3RBCvGSSGR*(K5RBCvGSSGR + K6RBCvGSSGR)*K8RBCvGSSGR*gshRBCi[t]*gssgRBCi[t]*nadphRBCi[t])/Vrbci^3 + (K10RBCvGSSGR*K1RBCvGSSGR*K3RBCvGSSGR*(K5RBCvGSSGR + K6RBCvGSSGR)*K8RBCvGSSGR*gshRBCi[t]^2*gssgRBCi[t]*nadphRBCi[t])/Vrbci^4 + (K12RBCvGSSGR*K2RBCvGSSGR*(K4RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K7RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR)*K9RBCvGSSGR*nadpRBCi[t])/Vrbci + (K10RBCvGSSGR*K12RBCvGSSGR*K2RBCvGSSGR*(K4RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K7RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR)*gshRBCi[t]*nadpRBCi[t])/Vrbci^2 + (K12RBCvGSSGR*K2RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBCi[t]*nadpRBCi[t])/Vrbci^2 + (K10RBCvGSSGR*K12RBCvGSSGR*(K2RBCvGSSGR*K4RBCvGSSGR + K2RBCvGSSGR*K5RBCvGSSGR + K2RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K6RBCvGSSGR)*K8RBCvGSSGR*gshRBCi[t]^2*nadpRBCi[t])/Vrbci^3 + (K12RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR*gssgRBCi[t]*nadpRBCi[t])/Vrbci^2 + (K10RBCvGSSGR*K12RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*gshRBCi[t]*gssgRBCi[t]*nadpRBCi[t])/Vrbci^3 + (K10RBCvGSSGR*K12RBCvGSSGR*K3RBCvGSSGR*(K5RBCvGSSGR + K6RBCvGSSGR)*K8RBCvGSSGR*gshRBCi[t]^2*gssgRBCi[t]*nadpRBCi[t])/Vrbci^4)
vRBCivHBADP = parMulti*TsRBC*Vrbci*(-((KdRBCvHBADP*hbadpRBCi[t])/Vrbci) + (KaappRBCvHBADP*adpRBCi[t]*hbRBCi[t])/Vrbci^2)
vRBCivHBATP = parMulti*TsRBC*Vrbci*(-((KdRBCvHBATP*hbatpRBCi[t])/Vrbci) + (KaappRBCvHBATP*atpRBCi[t]*hbRBCi[t])/Vrbci^2)
vRBCivHBB13PG = parMulti*TsRBC*Vrbci*(-((KdRBCvHBB13PG*hbb13pgRBCi[t])/Vrbci) + (KaappRBCvHBB13PG*b13pgRBCi[t]*hbRBCi[t])/Vrbci^2)
vRBCivHBB23PG = parMulti*TsRBC*Vrbci*(-((KdRBCvHBBPG*hbb23pgRBCi[t])/Vrbci) + (KaappRBCvHBBPG*b23pgRBCi[t]*hbRBCi[t])/Vrbci^2)
vRBCivHBMGATP = parMulti*TsRBC*Vrbci*(-((KdRBCvHBMGATP*hbmgatpRBCi[t])/Vrbci) + (KaappRBCvHBMGATP*hbRBCi[t]*mgatpRBCi[t])/Vrbci^2)
vRBCivHK = (parMulti*TsRBC*Vrbci*(-((ERBCvHK*KcatrappRBCvHK*g6pRBCi[t]*mgadpRBCi[t])/(Kig6pRBCvHK*KmmgadpRBCvHK*Vrbci^2)) + (ERBCvHK*KcatfappRBCvHK*glcRBCi[t]*mgatpRBCi[t])/(KiglcRBCvHK*KmmgatpRBCvHK*Vrbci^2)))/(1 + g6pRBCi[t]/(Kig6pRBCvHK*Vrbci) + glcRBCi[t]/(KiglcRBCvHK*Vrbci) + (b23pgRBCi[t]*glcRBCi[t])/(KdibpgRBCvHK*KiglcRBCvHK*Vrbci^2) + (g16p2RBCi[t]*glcRBCi[t])/(Kdig16p2RBCvHK*KiglcRBCvHK*Vrbci^2) + (g6pRBCi[t]*glcRBCi[t])/(Kdig6pRBCvHK*KiglcRBCvHK*Vrbci^2) + (glcRBCi[t]*gshRBCi[t])/(KdigshRBCvHK*KiglcRBCvHK*Vrbci^2) + mgadpRBCi[t]/(KimgadpRBCvHK*Vrbci) + (g6pRBCi[t]*mgadpRBCi[t])/(Kig6pRBCvHK*KmmgadpRBCvHK*Vrbci^2) + mgatpRBCi[t]/(KimgatpRBCvHK*Vrbci) + (glcRBCi[t]*mgatpRBCi[t])/(KiglcRBCvHK*KmmgatpRBCvHK*Vrbci^2))
vRBCivLACTRANSPORT = parMulti*TsRBC*vBld*(-((KiRBCvLACTRANSPORT*lacEXT)/vBld) + (KoRBCvLACTRANSPORT*lacRBCi[t])/Vrbci)
vRBCivLDH = (parMulti*TsRBC*Vrbci*VRBCvLDH*(-((KcatrRBCvLDH*lacRBCi[t]*nadRBCi[t])/(KinadRBCvLDH*KmapplacRBCvLDH*Vrbci^2)) + (KcatfRBCvLDH*nadhRBCi[t]*pyrRBCi[t])/(KinadhRBCvLDH*KmapppyrRBCvLDH*Vrbci^2)))/(nadhRBCi[t]/(KinadhRBCvLDH*Vrbci) + (KmnadRBCvLDH*lacRBCi[t]*nadhRBCi[t])/(KinadhRBCvLDH*KinadRBCvLDH*KmapplacRBCvLDH*Vrbci^2) + nadRBCi[t]/(KinadRBCvLDH*Vrbci) + (lacRBCi[t]*nadRBCi[t])/(KinadRBCvLDH*KmapplacRBCvLDH*Vrbci^2) + (nadhRBCi[t]*pyrRBCi[t])/(KinadhRBCvLDH*KmapppyrRBCvLDH*Vrbci^2) + (lacRBCi[t]*nadhRBCi[t]*pyrRBCi[t])/(KiapplacRBCvLDH*KinadhRBCvLDH*KmapppyrRBCvLDH*Vrbci^3) + (KmnadhRBCvLDH*nadRBCi[t]*pyrRBCi[t])/(KinadhRBCvLDH*KinadRBCvLDH*KmapppyrRBCvLDH*Vrbci^2) + (lacRBCi[t]*nadRBCi[t]*pyrRBCi[t])/(KiapppyrRBCvLDH*KinadRBCvLDH*KmapplacRBCvLDH*Vrbci^3) + (1 + pyrRBCi[t]/(KidpyrRBCvLDH*Vrbci))*(1 + (KmnadRBCvLDH*lacRBCi[t])/(KinadRBCvLDH*KmapplacRBCvLDH*Vrbci) + (KmnadhRBCvLDH*pyrRBCi[t])/(KinadhRBCvLDH*KmapppyrRBCvLDH*Vrbci)))
vRBCivLDHP = (parMulti*TsRBC*Vrbci*(-((KrRBCvLDHP*lacRBCi[t]*nadpRBCi[t])/(KmlacRBCvLDHP*Vrbci^2)) + (KfRBCvLDHP*nadphRBCi[t]*pyrRBCi[t])/(KmpyrRBCvLDHP*Vrbci^2)))/(1 + lacRBCi[t]/(KmlacRBCvLDHP*Vrbci) + pyrRBCi[t]/(KmpyrRBCvLDHP*Vrbci))
vRBCivMGADP = parMulti*TsRBC*Vrbci*(-((KdRBCvMGADP*mgadpRBCi[t])/Vrbci) + (KaappRBCvMGADP*adpRBCi[t]*mgRBCi[t])/Vrbci^2)
vRBCivMGATP = parMulti*TsRBC*Vrbci*(-((KdRBCvMGATP*mgatpRBCi[t])/Vrbci) + (KaappRBCvMGATP*atpRBCi[t]*mgRBCi[t])/Vrbci^2)
vRBCivMGB13PG = parMulti*TsRBC*Vrbci*(-((KdRBCvMGB13PG*mgb13pgRBCi[t])/Vrbci) + (KaappRBCvMGB13PG*b13pgRBCi[t]*mgRBCi[t])/Vrbci^2)
vRBCivMGB23PG = parMulti*TsRBC*Vrbci*(-((KdRBCvMGB23PG*mgb23pgRBCi[t])/Vrbci) + (KaappRBCvMGB23PG*b23pgRBCi[t]*mgRBCi[t])/Vrbci^2)
vRBCivMGF16P2 = parMulti*TsRBC*Vrbci*(-((KdRBCvMGF16P2*mgf16p2RBCi[t])/Vrbci) + (KaappRBCvMGF16P2*f16p2RBCi[t]*mgRBCi[t])/Vrbci^2)
vRBCivMGG16P2 = parMulti*TsRBC*Vrbci*(-((KdRBCvMGG16P2*mgg16p2RBCi[t])/Vrbci) + (KaappRBCvMGG16P2*g16p2RBCi[t]*mgRBCi[t])/Vrbci^2)
vRBCivMGPHOS = parMulti*TsRBC*Vrbci*(-((KdRBCvMGPHOS*mgphosRBCi[t])/Vrbci) + (KaappRBCvMGPHOS*mgRBCi[t]*phosRBCi[t])/Vrbci^2)
vRBCivOX = KRBCvOX*parMulti*TsRBC*gshRBCi[t]
vRBCivOXNADH = KRBCvOXNADH*parMulti*TsRBC*nadhRBCi[t]
vRBCivP6GDH = (parMulti*TsRBC*Vrbci*VRBCvP6GDH*((K11RBCvP6GDH*K1RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH*nadpRBCi[t]*p6gRBCi[t])/Vrbci^2 - (co2RBCi*K10RBCvP6GDH*K12RBCvP6GDH*K2RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*nadphRBCi[t]*ru5pRBCi[t])/Vrbci^3))/(K11RBCvP6GDH*K2RBCvP6GDH*(K4RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K7RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH)*K9RBCvP6GDH + (co2RBCi*K11RBCvP6GDH*K2RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH)/Vrbci + (co2RBCi*K12RBCvP6GDH*K2RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*nadphRBCi[t])/Vrbci^2 + (K12RBCvP6GDH*K2RBCvP6GDH*(K4RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K7RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH)*K9RBCvP6GDH*nadphRBCi[t])/Vrbci + (co2RBCi*K11RBCvP6GDH*K1RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*nadpRBCi[t])/Vrbci^2 + (K11RBCvP6GDH*K1RBCvP6GDH*(K4RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K7RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH)*K9RBCvP6GDH*nadpRBCi[t])/Vrbci + (K11RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH*p6gRBCi[t])/Vrbci + (K12RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH*nadphRBCi[t]*p6gRBCi[t])/Vrbci^2 + (co2RBCi*K11RBCvP6GDH*K1RBCvP6GDH*K3RBCvP6GDH*(K5RBCvP6GDH + K6RBCvP6GDH)*K8RBCvP6GDH*nadpRBCi[t]*p6gRBCi[t])/Vrbci^3 + (K1RBCvP6GDH*K3RBCvP6GDH*(K11RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH + K11RBCvP6GDH*K5RBCvP6GDH*K9RBCvP6GDH + K11RBCvP6GDH*K6RBCvP6GDH*K9RBCvP6GDH + K11RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH)*nadpRBCi[t]*p6gRBCi[t])/Vrbci^2 + (co2RBCi*K10RBCvP6GDH*K2RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*ru5pRBCi[t])/Vrbci^2 + (co2RBCi*K10RBCvP6GDH*K12RBCvP6GDH*(K2RBCvP6GDH*K4RBCvP6GDH + K2RBCvP6GDH*K5RBCvP6GDH + K2RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K6RBCvP6GDH)*K8RBCvP6GDH*nadphRBCi[t]*ru5pRBCi[t])/Vrbci^3 + (K10RBCvP6GDH*K12RBCvP6GDH*K2RBCvP6GDH*(K4RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K7RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH)*nadphRBCi[t]*ru5pRBCi[t])/Vrbci^2 + (co2RBCi*K10RBCvP6GDH*K1RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*nadpRBCi[t]*ru5pRBCi[t])/Vrbci^3 + (co2RBCi*K10RBCvP6GDH*K12RBCvP6GDH*K3RBCvP6GDH*(K5RBCvP6GDH + K6RBCvP6GDH)*K8RBCvP6GDH*nadphRBCi[t]*p6gRBCi[t]*ru5pRBCi[t])/Vrbci^4 + (K10RBCvP6GDH*K12RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*nadphRBCi[t]*p6gRBCi[t]*ru5pRBCi[t])/Vrbci^3 + (co2RBCi*K10RBCvP6GDH*K1RBCvP6GDH*K3RBCvP6GDH*(K5RBCvP6GDH + K6RBCvP6GDH)*K8RBCvP6GDH*nadpRBCi[t]*p6gRBCi[t]*ru5pRBCi[t])/Vrbci^4 + (K10RBCvP6GDH*K1RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*nadpRBCi[t]*p6gRBCi[t]*ru5pRBCi[t])/Vrbci^3)
vRBCivPFK = (parMulti*TsRBC*Vrbci*VRBCvPFK*(-((KcatfRBCvPFK*f16p2RBCi[t]*mgadpRBCi[t])/(Kmf16p2RBCvPFK*KmmgadpRBCvPFK*Vrbci^2)) + (KcatfRBCvPFK*f6pRBCi[t]*mgatpRBCi[t])/(Kmf6pRBCvPFK*KmmgatpRBCvPFK*Vrbci^2)))/((1 + f16p2RBCi[t]/(Kmf16p2RBCvPFK*Vrbci) + f6pRBCi[t]/(Kmf6pRBCvPFK*Vrbci) + mgadpRBCi[t]/(KmmgadpRBCvPFK*Vrbci) + (f16p2RBCi[t]*mgadpRBCi[t])/(Kmf16p2RBCvPFK*KmmgadpRBCvPFK*Vrbci^2) + mgatpRBCi[t]/(KmmgatpRBCvPFK*Vrbci) + (f6pRBCi[t]*mgatpRBCi[t])/(Kmf6pRBCvPFK*KmmgatpRBCvPFK*Vrbci^2))*(1 + ((hRBC/KaRBCvPFK)^nRBCvPFK*(1 + atpRBCi[t]/(KtatpRBCvPFK*Vrbci))^4*(1 + b23pgRBCi[t]/(Ktb23pgRBCvPFK*Vrbci))^4*(1 + mgRBCi[t]/(KtmgRBCvPFK*Vrbci))^4)/((1 + ampRBCi[t]/(KrampRBCvPFK*Vrbci))^4*(1 + f16p2RBCi[t]/(Kmf16p2RBCvPFK*Vrbci) + f6pRBCi[t]/(Kmf6pRBCvPFK*Vrbci))^4*(1 + g16p2RBCi[t]/(Krg16p2RBCvPFK*Vrbci))^4*(1 + phosRBCi[t]/(KrphosRBCvPFK*Vrbci))^4)))
vRBCivPGI = (ERBCvPGI*parMulti*TsRBC*Vrbci*(-((KcatrRBCvPGI*f6pRBCi[t])/(Kmf6pRBCvPGI*Vrbci)) + (KcatfRBCvPGI*g6pRBCi[t])/(Kmg6pRBCvPGI*Vrbci)))/(1 + f6pRBCi[t]/(Kmf6pRBCvPGI*Vrbci) + g6pRBCi[t]/(Kmg6pRBCvPGI*Vrbci))
vRBCivPGK = (parMulti*TsRBC*Vrbci*VRBCvPGK*((KcatfRBCvPGK*b13pgRBCi[t]*mgadpRBCi[t])/(KimgadpRBCvPGK*Kmb13pgRBCvPGK*Vrbci^2) - (KcatrRBCvPGK*mgatpRBCi[t]*p3gRBCi[t])/(KimgatpRBCvPGK*Kmp3gRBCvPGK*Vrbci^2)))/(1 + b13pgRBCi[t]/(Kib13pgRBCvPGK*Vrbci) + mgadpRBCi[t]/(KimgadpRBCvPGK*Vrbci) + (b13pgRBCi[t]*mgadpRBCi[t])/(KimgadpRBCvPGK*Kmb13pgRBCvPGK*Vrbci^2) + mgatpRBCi[t]/(KimgatpRBCvPGK*Vrbci) + p3gRBCi[t]/(Kip3gRBCvPGK*Vrbci) + (mgatpRBCi[t]*p3gRBCi[t])/(KimgatpRBCvPGK*Kmp3gRBCvPGK*Vrbci^2))
vRBCivPGLHYDROLYSIS = parMulti*TsRBC*Vrbci*((KhydrolRBCvPGLHYDROLYSIS*p6glRBCi[t])/Vrbci + (Klactonase3RBCvPGLHYDROLYSIS*VRBCvPGLHYDROLYSIS*p6glRBCi[t])/(Vrbci*((Klactonase2RBCvPGLHYDROLYSIS + Klactonase3RBCvPGLHYDROLYSIS)/Klactonase1RBCvPGLHYDROLYSIS + p6glRBCi[t]/Vrbci)))
vRBCivPGM = (parMulti*TsRBC*Vrbci*VRBCvPGM*(-((KcatrRBCvPGM*p2gRBCi[t])/(Kmp2gRBCvPGM*Vrbci)) + (KcatfRBCvPGM*p3gRBCi[t])/(Kmp3gRBCvPGM*Vrbci)))/(1 + p2gRBCi[t]/(Kmp2gRBCvPGM*Vrbci) + p3gRBCi[t]/(Kmp3gRBCvPGM*Vrbci))
vRBCivPHOSTRANSPORT = parMulti*TsRBC*vBld*(-((KiRBCvPHOSTRANSPORT*phosEXT)/vBld) + (KoRBCvPHOSTRANSPORT*phosRBCi[t])/Vrbci)
vRBCivPK = (parMulti*TsRBC*Vrbci*VRBCvPK*((KcatfRBCvPK*mgadpRBCi[t]*pepRBCi[t])/(KrmgadpRBCvPK*KrpepRBCvPK*Vrbci^2) - (KcatrRBCvPK*mgatpRBCi[t]*pyrRBCi[t])/(KrmgatpRBCvPK*KrpyrRBCvPK*Vrbci^2)))/((1 + LRBCivPK)*(1 + mgadpRBCi[t]/(KrmgadpRBCvPK*Vrbci) + mgatpRBCi[t]/(KrmgatpRBCvPK*Vrbci) + pepRBCi[t]/(KrpepRBCvPK*Vrbci) + (mgadpRBCi[t]*pepRBCi[t])/(KrmgadpRBCvPK*KrpepRBCvPK*Vrbci^2) + pyrRBCi[t]/(KrpyrRBCvPK*Vrbci) + (mgatpRBCi[t]*pyrRBCi[t])/(KrmgatpRBCvPK*KrpyrRBCvPK*Vrbci^2)))
vRBCivPYRTRANSPORT = parMulti*TsRBC*vBld*(-((KiRBCvPYRTRANSPORT*pyrEXT)/vBld) + (KoRBCvPYRTRANSPORT*pyrRBCi[t])/Vrbci)
vRBCivR5PI = (parMulti*TsRBC*Vrbci*VRBCvR5PI*(-((K2RBCvR5PI*K4RBCvR5PI*rib5pRBCi[t])/((K2RBCvR5PI + K3RBCvR5PI)*Vrbci)) + (K1RBCvR5PI*K3RBCvR5PI*ru5pRBCi[t])/((K2RBCvR5PI + K3RBCvR5PI)*Vrbci)))/(1 + (K1RBCvR5PI*ru5pRBCi[t])/((K2RBCvR5PI + K3RBCvR5PI)*Vrbci) + (K4RBCvR5PI*xu5pRBCi[t])/((K2RBCvR5PI + K3RBCvR5PI)*Vrbci))
vRBCivRu5PE = (parMulti*TsRBC*Vrbci*VRBCvRu5PE*((K1RBCvRu5PE*K3RBCvRu5PE*ru5pRBCi[t])/((K2RBCvRu5PE + K3RBCvRu5PE)*Vrbci) - (K2RBCvRu5PE*K4RBCvRu5PE*xu5pRBCi[t])/((K2RBCvRu5PE + K3RBCvRu5PE)*Vrbci)))/(1 + (K1RBCvRu5PE*ru5pRBCi[t])/((K2RBCvRu5PE + K3RBCvRu5PE)*Vrbci) + (K4RBCvRu5PE*xu5pRBCi[t])/((K2RBCvRu5PE + K3RBCvRu5PE)*Vrbci))
vRBCivTA = (parMulti*TsRBC*Vrbci*VRBCvTA*(-((K2RBCvTA*K4RBCvTA*K6RBCvTA*K8RBCvTA*ery4pRBCi[t]*f6pRBCi[t])/Vrbci^2) + (K1RBCvTA*K3RBCvTA*K5RBCvTA*K7RBCvTA*gapRBCi[t]*sed7pRBCi[t])/Vrbci^2))/((K2RBCvTA*K4RBCvTA*(K6RBCvTA + K7RBCvTA)*ery4pRBCi[t])/Vrbci + ((K2RBCvTA + K3RBCvTA)*K6RBCvTA*K8RBCvTA*f6pRBCi[t])/Vrbci + (K4RBCvTA*(K2RBCvTA + K6RBCvTA)*K8RBCvTA*ery4pRBCi[t]*f6pRBCi[t])/Vrbci^2 + ((K2RBCvTA + K3RBCvTA)*K5RBCvTA*K7RBCvTA*gapRBCi[t])/Vrbci + ((K2RBCvTA + K3RBCvTA)*K5RBCvTA*K8RBCvTA*f6pRBCi[t]*gapRBCi[t])/Vrbci^2 + (K1RBCvTA*K3RBCvTA*(K6RBCvTA + K7RBCvTA)*sed7pRBCi[t])/Vrbci + (K1RBCvTA*K4RBCvTA*(K6RBCvTA + K7RBCvTA)*ery4pRBCi[t]*sed7pRBCi[t])/Vrbci^2 + (K1RBCvTA*K5RBCvTA*(K3RBCvTA + K7RBCvTA)*gapRBCi[t]*sed7pRBCi[t])/Vrbci^2)
vRBCivTIM = (parMulti*TsRBC*Vrbci*VRBCvTIM*(-((KcatrRBCvTIM*dhapRBCi[t])/(KmdhapRBCvTIM*Vrbci)) + (KcatfRBCvTIM*gapRBCi[t])/(KmgapRBCvTIM*Vrbci)))/(1 + dhapRBCi[t]/(KmdhapRBCvTIM*Vrbci) + gapRBCi[t]/(KmgapRBCvTIM*Vrbci))
vRBCivTK1 = parMulti*TsRBC*Vrbci*(-((K2RBCvTK1*tkxu5pRBCi[t])/Vrbci) + (K1RBCvTK1*tkRBCi[t]*xu5pRBCi[t])/Vrbci^2)
vRBCivTK2 = parMulti*TsRBC*Vrbci*(-((K4RBCvTK2*gapRBCi[t]*tkgRBCi[t])/Vrbci^2) + (K3RBCvTK2*tkxu5pRBCi[t])/Vrbci)
vRBCivTK3 = parMulti*TsRBC*Vrbci*((K5RBCvTK3*rib5pRBCi[t]*tkgRBCi[t])/Vrbci^2 - (K6RBCvTK3*tkgrib5pRBCi[t])/Vrbci)
vRBCivTK4 = parMulti*TsRBC*Vrbci*((K7RBCvTK4*tkgrib5pRBCi[t])/Vrbci - (K8RBCvTK4*sed7pRBCi[t]*tkRBCi[t])/Vrbci^2)
vRBCivTK5 = parMulti*TsRBC*Vrbci*(-((K10RBCvTK5*tkgery4pRBCi[t])/Vrbci) + (K9RBCvTK5*ery4pRBCi[t]*tkgRBCi[t])/Vrbci^2)
vRBCivTK6 = parMulti*TsRBC*Vrbci*((K11RBCvTK6*tkgery4pRBCi[t])/Vrbci - (K12RBCvTK6*f6pRBCi[t]*tkRBCi[t])/Vrbci^2)
vRBCvAK = parMultiUninf*TsRBC*Vrbcu*((K1appRBCvAK*adpRBC[t]*mgadpRBC[t])/Vrbcu^2 - (K2appRBCvAK*ampRBC[t]*mgatpRBC[t])/Vrbcu^2)
vRBCvALD = (parMultiUninf*TsRBC*Vrbcu*VRBCvALD*((KcatfRBCvALD*f16p2RBC[t])/(Kmf16p2RBCvALD*Vrbcu) - (KcatrRBCvALD*dhapRBC[t]*gapRBC[t])/(KidhapRBCvALD*KmgapRBCvALD*Vrbcu^2)))/(1 + dhapRBC[t]/(KidhapRBCvALD*Vrbcu) + f16p2RBC[t]/(Kmf16p2RBCvALD*Vrbcu) + (dhapRBC[t]*gapRBC[t])/(KidhapRBCvALD*KmgapRBCvALD*Vrbcu^2) + (KmdhapRBCvALD*f16p2RBC[t]*gapRBC[t])/(KidhapRBCvALD*Kif16p2RBCvALD*KmgapRBCvALD*Vrbcu^2) + (b23pgRBC[t]/Vrbcu + mgb23pgRBC[t]/Vrbcu)/Kib23pgRBCvALD + (KmdhapRBCvALD*gapRBC[t]*(1 + (b23pgRBC[t]/Vrbcu + mgb23pgRBC[t]/Vrbcu)/Kib23pgRBCvALD))/(KidhapRBCvALD*KmgapRBCvALD*Vrbcu))
vRBCvATPASE = KRBCvATPASE*parMultiUninf*TsRBC*mgatpRBC[t]
vRBCvBPGSP1 = parMultiUninf*TsRBC*Vrbcu*(-((K2RBCvBPGSP1*bpgspb13pgRBC[t])/Vrbcu) + (K1appRBCvBPGSP1*b13pgRBC[t]*bpgspRBC[t])/Vrbcu^2)
vRBCvBPGSP2 = K3appRBCvBPGSP2*parMultiUninf*TsRBC*bpgspb13pgRBC[t]
vRBCvBPGSP3 = parMultiUninf*TsRBC*Vrbcu*(-((K5RBCvBPGSP3*bpgsppp3gRBC[t])/Vrbcu) + (K4appRBCvBPGSP3*bpgsppRBC[t]*p3gRBC[t])/Vrbcu^2)
vRBCvBPGSP4 = parMultiUninf*TsRBC*Vrbcu*(-((K7RBCvBPGSP4*bpgsppp2gRBC[t])/Vrbcu) + (K6appRBCvBPGSP4*bpgsppRBC[t]*p2gRBC[t])/Vrbcu^2)
vRBCvBPGSP5 = parMultiUninf*TsRBC*Vrbcu*(-((K9RBCvBPGSP5*bpgspb23pgRBC[t])/Vrbcu) + (K8RBCvBPGSP5*bpgsppp3gRBC[t])/Vrbcu)
vRBCvBPGSP6 = parMultiUninf*TsRBC*Vrbcu*(-((K11RBCvBPGSP6*bpgspb23pgRBC[t])/Vrbcu) + (K10RBCvBPGSP6*bpgsppp2gRBC[t])/Vrbcu)
vRBCvBPGSP7 = parMultiUninf*TsRBC*Vrbcu*((K12RBCvBPGSP7*bpgspb23pgRBC[t])/Vrbcu - (K13appRBCvBPGSP7*b23pgRBC[t]*bpgspRBC[t])/Vrbcu^2)
vRBCvBPGSP8 = parMultiUninf*TsRBC*Vrbcu*(-((K15RBCvBPGSP8*bpgsppphosRBC[t])/Vrbcu) + (K14RBCvBPGSP8*bpgsppRBC[t]*phosRBC[t])/Vrbcu^2)
vRBCvBPGSP9 = K16RBCvBPGSP9*parMultiUninf*TsRBC*bpgsppphosRBC[t]
vRBCvENO = (parMultiUninf*TsRBC*Vrbcu*VRBCvENO*((KcatfRBCvENO*mgRBC[t]*p2gRBC[t])/(KimgRBCvENO*Kmp2gRBCvENO*Vrbcu^2) - (KcatrRBCvENO*mgRBC[t]*pepRBC[t])/(KipepRBCvENO*KmmgRBCvENO*Vrbcu^2)))/(1 + mgRBC[t]/(KimgRBCvENO*Vrbcu) + p2gRBC[t]/(Kip2gRBCvENO*Vrbcu) + (mgRBC[t]*p2gRBC[t])/(KimgRBCvENO*Kmp2gRBCvENO*Vrbcu^2) + pepRBC[t]/(KipepRBCvENO*Vrbcu) + (mgRBC[t]*pepRBC[t])/(KipepRBCvENO*KmmgRBCvENO*Vrbcu^2))
vRBCvG6PDH = (parMultiUninf*TsRBC*Vrbcu*VRBCvG6PDH*((K1RBCvG6PDH*K3RBCvG6PDH*K5RBCvG6PDH*K7RBCvG6PDH*K9RBCvG6PDH*g6pRBC[t]*nadpRBC[t])/Vrbcu^2 - (K10RBCvG6PDH*K2RBCvG6PDH*K4RBCvG6PDH*K6RBCvG6PDH*K8RBCvG6PDH*nadphRBC[t]*p6glRBC[t])/Vrbcu^2))/(K2RBCvG6PDH*(K4RBCvG6PDH*K6RBCvG6PDH + K4RBCvG6PDH*K7RBCvG6PDH + K5RBCvG6PDH*K7RBCvG6PDH)*K9RBCvG6PDH + (K3RBCvG6PDH*K5RBCvG6PDH*K7RBCvG6PDH*K9RBCvG6PDH*g6pRBC[t])/Vrbcu + (K10RBCvG6PDH*K2RBCvG6PDH*(K4RBCvG6PDH*K6RBCvG6PDH + K5RBCvG6PDH*K6RBCvG6PDH + K5RBCvG6PDH*K7RBCvG6PDH)*nadphRBC[t])/Vrbcu + (K10RBCvG6PDH*K3RBCvG6PDH*K5RBCvG6PDH*K7RBCvG6PDH*g6pRBC[t]*nadphRBC[t])/Vrbcu^2 + (K1RBCvG6PDH*(K4RBCvG6PDH*K6RBCvG6PDH + K4RBCvG6PDH*K7RBCvG6PDH + K5RBCvG6PDH*K7RBCvG6PDH)*K9RBCvG6PDH*nadpRBC[t])/Vrbcu + (K1RBCvG6PDH*K3RBCvG6PDH*(K5RBCvG6PDH*K7RBCvG6PDH + K5RBCvG6PDH*K9RBCvG6PDH + K6RBCvG6PDH*K9RBCvG6PDH + K7RBCvG6PDH*K9RBCvG6PDH)*g6pRBC[t]*nadpRBC[t])/Vrbcu^2 + (K2RBCvG6PDH*K4RBCvG6PDH*K6RBCvG6PDH*K8RBCvG6PDH*p6glRBC[t])/Vrbcu + (K10RBCvG6PDH*(K2RBCvG6PDH*K4RBCvG6PDH + K2RBCvG6PDH*K5RBCvG6PDH + K2RBCvG6PDH*K6RBCvG6PDH + K4RBCvG6PDH*K6RBCvG6PDH)*K8RBCvG6PDH*nadphRBC[t]*p6glRBC[t])/Vrbcu^2 + (K10RBCvG6PDH*K3RBCvG6PDH*(K5RBCvG6PDH + K6RBCvG6PDH)*K8RBCvG6PDH*g6pRBC[t]*nadphRBC[t]*p6glRBC[t])/Vrbcu^3 + (K1RBCvG6PDH*K4RBCvG6PDH*K6RBCvG6PDH*K8RBCvG6PDH*nadpRBC[t]*p6glRBC[t])/Vrbcu^2 + (K1RBCvG6PDH*K3RBCvG6PDH*(K5RBCvG6PDH + K6RBCvG6PDH)*K8RBCvG6PDH*g6pRBC[t]*nadpRBC[t]*p6glRBC[t])/Vrbcu^3)
vRBCvGAPDH = (parMultiUninf*TsRBC*Vrbcu*VRBCvGAPDH*(-((KcatrappRBCvGAPDH*b13pgRBC[t]*nadhRBC[t])/(Kiappb13pgRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbcu^2)) + (KcatfappRBCvGAPDH*gapRBC[t]*nadRBC[t]*phosRBC[t])/(KiappgapRBCvGAPDH*KiphosRBCvGAPDH*KmnadRBCvGAPDH*Vrbcu^3)))/((b13pgRBC[t]*(1 + gapRBC[t]/(KidgapRBCvGAPDH*Vrbcu)))/(Kiappb13pgRBCvGAPDH*Vrbcu) + (gapRBC[t]*(1 + gapRBC[t]/(KidgapRBCvGAPDH*Vrbcu)))/(KiappgapRBCvGAPDH*Vrbcu) + (Kmb13pgRBCvGAPDH*nadhRBC[t])/(Kiappb13pgRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbcu) + (b13pgRBC[t]*nadhRBC[t])/(Kiappb13pgRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbcu^2) + (gapRBC[t]*nadhRBC[t])/(KiappgapRBCvGAPDH*KiappnadhRBCvGAPDH*Vrbcu^2) + (b13pgRBC[t]*nadRBC[t])/(Kiappb13pgRBCvGAPDH*KinadRBCvGAPDH*Vrbcu^2) + (gapRBC[t]*nadRBC[t])/(KiappgapRBCvGAPDH*KinadRBCvGAPDH*Vrbcu^2) + (gapRBC[t]*(1 + gapRBC[t]/(KidgapRBCvGAPDH*Vrbcu))*phosRBC[t])/(KiappgapRBCvGAPDH*KiphosRBCvGAPDH*Vrbcu^2) + (Kmb13pgRBCvGAPDH*nadhRBC[t]*phosRBC[t])/(Kiappb13pgRBCvGAPDH*KiphosRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbcu^2) + (Kmb13pgRBCvGAPDH*b13pgRBC[t]*nadhRBC[t]*phosRBC[t])/(Kiappb13pgRBCvGAPDH*Kidb13pgRBCvGAPDH*KiphosRBCvGAPDH*KmappnadhRBCvGAPDH*Vrbcu^3) + (gapRBC[t]*nadhRBC[t]*phosRBC[t])/(KiappgapRBCvGAPDH*KiappnadhRBCvGAPDH*KiphosRBCvGAPDH*Vrbcu^3) + (KmgapRBCvGAPDH*nadRBC[t]*phosRBC[t])/(KiappgapRBCvGAPDH*KiphosRBCvGAPDH*KmnadRBCvGAPDH*Vrbcu^2) + (KmgapRBCvGAPDH*b13pgRBC[t]*nadRBC[t]*phosRBC[t])/(KiappgapRBCvGAPDH*Kidb13pgRBCvGAPDH*KiphosRBCvGAPDH*KmnadRBCvGAPDH*Vrbcu^3) + (gapRBC[t]*nadRBC[t]*phosRBC[t])/(KiappgapRBCvGAPDH*KiphosRBCvGAPDH*KmnadRBCvGAPDH*Vrbcu^3))
vRBCvGLCTRANSPORT = parMultiUninf*TsRBC*Vrbcu*((glcAlpha*glcEXT*VRBCvGLCTRANSPORT)/((kAlpha + (glcAlpha*glcEXT)/vBld + ((1 - glcAlpha)*glcEXT*kAlpha)/(kBeta*vBld))*vBld) + ((1 - glcAlpha)*glcEXT*VRBCvGLCTRANSPORT)/((kBeta + ((1 - glcAlpha)*glcEXT)/vBld + (glcAlpha*glcEXT*kBeta)/(kAlpha*vBld))*vBld) - (glcAlpha*VRBCvGLCTRANSPORT*glcRBC[t])/(Vrbcu*(kAlpha + (glcAlpha*glcRBC[t])/Vrbcu + ((1 - glcAlpha)*kAlpha*glcRBC[t])/(kBeta*Vrbcu))) - ((1 - glcAlpha)*VRBCvGLCTRANSPORT*glcRBC[t])/(Vrbcu*(kBeta + ((1 - glcAlpha)*glcRBC[t])/Vrbcu + (glcAlpha*kBeta*glcRBC[t])/(kAlpha*Vrbcu))))
vRBCvGSSGR = (parMultiUninf*TsRBC*Vrbcu*VRBCvGSSGR*((K11RBCvGSSGR*K1RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR*gssgRBC[t]*nadphRBC[t])/Vrbcu^2 - (K10RBCvGSSGR*K12RBCvGSSGR*K2RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBC[t]^2*nadpRBC[t])/Vrbcu^3))/(K11RBCvGSSGR*K2RBCvGSSGR*(K4RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K7RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR)*K9RBCvGSSGR + (K11RBCvGSSGR*K2RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBC[t])/Vrbcu + (K10RBCvGSSGR*K2RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBC[t]^2)/Vrbcu^2 + (K11RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR*gssgRBC[t])/Vrbcu + (K11RBCvGSSGR*K1RBCvGSSGR*(K4RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K7RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR)*K9RBCvGSSGR*nadphRBC[t])/Vrbcu + (K11RBCvGSSGR*K1RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBC[t]*nadphRBC[t])/Vrbcu^2 + (K10RBCvGSSGR*K1RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBC[t]^2*nadphRBC[t])/Vrbcu^3 + (K1RBCvGSSGR*K3RBCvGSSGR*(K11RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR + K11RBCvGSSGR*K5RBCvGSSGR*K9RBCvGSSGR + K11RBCvGSSGR*K6RBCvGSSGR*K9RBCvGSSGR + K11RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR)*gssgRBC[t]*nadphRBC[t])/Vrbcu^2 + (K10RBCvGSSGR*K1RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*gshRBC[t]*gssgRBC[t]*nadphRBC[t])/Vrbcu^3 + (K11RBCvGSSGR*K1RBCvGSSGR*K3RBCvGSSGR*(K5RBCvGSSGR + K6RBCvGSSGR)*K8RBCvGSSGR*gshRBC[t]*gssgRBC[t]*nadphRBC[t])/Vrbcu^3 + (K10RBCvGSSGR*K1RBCvGSSGR*K3RBCvGSSGR*(K5RBCvGSSGR + K6RBCvGSSGR)*K8RBCvGSSGR*gshRBC[t]^2*gssgRBC[t]*nadphRBC[t])/Vrbcu^4 + (K12RBCvGSSGR*K2RBCvGSSGR*(K4RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K7RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR)*K9RBCvGSSGR*nadpRBC[t])/Vrbcu + (K10RBCvGSSGR*K12RBCvGSSGR*K2RBCvGSSGR*(K4RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K7RBCvGSSGR + K5RBCvGSSGR*K7RBCvGSSGR)*gshRBC[t]*nadpRBC[t])/Vrbcu^2 + (K12RBCvGSSGR*K2RBCvGSSGR*K4RBCvGSSGR*K6RBCvGSSGR*K8RBCvGSSGR*gshRBC[t]*nadpRBC[t])/Vrbcu^2 + (K10RBCvGSSGR*K12RBCvGSSGR*(K2RBCvGSSGR*K4RBCvGSSGR + K2RBCvGSSGR*K5RBCvGSSGR + K2RBCvGSSGR*K6RBCvGSSGR + K4RBCvGSSGR*K6RBCvGSSGR)*K8RBCvGSSGR*gshRBC[t]^2*nadpRBC[t])/Vrbcu^3 + (K12RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*K9RBCvGSSGR*gssgRBC[t]*nadpRBC[t])/Vrbcu^2 + (K10RBCvGSSGR*K12RBCvGSSGR*K3RBCvGSSGR*K5RBCvGSSGR*K7RBCvGSSGR*gshRBC[t]*gssgRBC[t]*nadpRBC[t])/Vrbcu^3 + (K10RBCvGSSGR*K12RBCvGSSGR*K3RBCvGSSGR*(K5RBCvGSSGR + K6RBCvGSSGR)*K8RBCvGSSGR*gshRBC[t]^2*gssgRBC[t]*nadpRBC[t])/Vrbcu^4)
vRBCvHBADP = parMultiUninf*TsRBC*Vrbcu*(-((KdRBCvHBADP*hbadpRBC[t])/Vrbcu) + (KaappRBCvHBADP*adpRBC[t]*hbRBC[t])/Vrbcu^2)
vRBCvHBATP = parMultiUninf*TsRBC*Vrbcu*(-((KdRBCvHBATP*hbatpRBC[t])/Vrbcu) + (KaappRBCvHBATP*atpRBC[t]*hbRBC[t])/Vrbcu^2)
vRBCvHBB13PG = parMultiUninf*TsRBC*Vrbcu*(-((KdRBCvHBB13PG*hbb13pgRBC[t])/Vrbcu) + (KaappRBCvHBB13PG*b13pgRBC[t]*hbRBC[t])/Vrbcu^2)
vRBCvHBB23PG = parMultiUninf*TsRBC*Vrbcu*(-((KdRBCvHBBPG*hbb23pgRBC[t])/Vrbcu) + (KaappRBCvHBBPG*b23pgRBC[t]*hbRBC[t])/Vrbcu^2)
vRBCvHBMGATP = parMultiUninf*TsRBC*Vrbcu*(-((KdRBCvHBMGATP*hbmgatpRBC[t])/Vrbcu) + (KaappRBCvHBMGATP*hbRBC[t]*mgatpRBC[t])/Vrbcu^2)
vRBCvHK = (parMultiUninf*TsRBC*Vrbcu*(-((ERBCvHK*KcatrappRBCvHK*g6pRBC[t]*mgadpRBC[t])/(Kig6pRBCvHK*KmmgadpRBCvHK*Vrbcu^2)) + (ERBCvHK*KcatfappRBCvHK*glcRBC[t]*mgatpRBC[t])/(KiglcRBCvHK*KmmgatpRBCvHK*Vrbcu^2)))/(1 + g6pRBC[t]/(Kig6pRBCvHK*Vrbcu) + glcRBC[t]/(KiglcRBCvHK*Vrbcu) + (b23pgRBC[t]*glcRBC[t])/(KdibpgRBCvHK*KiglcRBCvHK*Vrbcu^2) + (g16p2RBC[t]*glcRBC[t])/(Kdig16p2RBCvHK*KiglcRBCvHK*Vrbcu^2) + (g6pRBC[t]*glcRBC[t])/(Kdig6pRBCvHK*KiglcRBCvHK*Vrbcu^2) + (glcRBC[t]*gshRBC[t])/(KdigshRBCvHK*KiglcRBCvHK*Vrbcu^2) + mgadpRBC[t]/(KimgadpRBCvHK*Vrbcu) + (g6pRBC[t]*mgadpRBC[t])/(Kig6pRBCvHK*KmmgadpRBCvHK*Vrbcu^2) + mgatpRBC[t]/(KimgatpRBCvHK*Vrbcu) + (glcRBC[t]*mgatpRBC[t])/(KiglcRBCvHK*KmmgatpRBCvHK*Vrbcu^2))
vRBCvLACTRANSPORT = parMultiUninf*TsRBC*vBld*(-((KiRBCvLACTRANSPORT*lacEXT)/vBld) + (KoRBCvLACTRANSPORT*lacRBC[t])/Vrbcu)
vRBCvLDH = (parMultiUninf*TsRBC*Vrbcu*VRBCvLDH*(-((KcatrRBCvLDH*lacRBC[t]*nadRBC[t])/(KinadRBCvLDH*KmapplacRBCvLDH*Vrbcu^2)) + (KcatfRBCvLDH*nadhRBC[t]*pyrRBC[t])/(KinadhRBCvLDH*KmapppyrRBCvLDH*Vrbcu^2)))/(nadhRBC[t]/(KinadhRBCvLDH*Vrbcu) + (KmnadRBCvLDH*lacRBC[t]*nadhRBC[t])/(KinadhRBCvLDH*KinadRBCvLDH*KmapplacRBCvLDH*Vrbcu^2) + nadRBC[t]/(KinadRBCvLDH*Vrbcu) + (lacRBC[t]*nadRBC[t])/(KinadRBCvLDH*KmapplacRBCvLDH*Vrbcu^2) + (nadhRBC[t]*pyrRBC[t])/(KinadhRBCvLDH*KmapppyrRBCvLDH*Vrbcu^2) + (lacRBC[t]*nadhRBC[t]*pyrRBC[t])/(KiapplacRBCvLDH*KinadhRBCvLDH*KmapppyrRBCvLDH*Vrbcu^3) + (KmnadhRBCvLDH*nadRBC[t]*pyrRBC[t])/(KinadhRBCvLDH*KinadRBCvLDH*KmapppyrRBCvLDH*Vrbcu^2) + (lacRBC[t]*nadRBC[t]*pyrRBC[t])/(KiapppyrRBCvLDH*KinadRBCvLDH*KmapplacRBCvLDH*Vrbcu^3) + (1 + pyrRBC[t]/(KidpyrRBCvLDH*Vrbcu))*(1 + (KmnadRBCvLDH*lacRBC[t])/(KinadRBCvLDH*KmapplacRBCvLDH*Vrbcu) + (KmnadhRBCvLDH*pyrRBC[t])/(KinadhRBCvLDH*KmapppyrRBCvLDH*Vrbcu)))
vRBCvLDHP = (parMultiUninf*TsRBC*Vrbcu*(-((KrRBCvLDHP*lacRBC[t]*nadpRBC[t])/(KmlacRBCvLDHP*Vrbcu^2)) + (KfRBCvLDHP*nadphRBC[t]*pyrRBC[t])/(KmpyrRBCvLDHP*Vrbcu^2)))/(1 + lacRBC[t]/(KmlacRBCvLDHP*Vrbcu) + pyrRBC[t]/(KmpyrRBCvLDHP*Vrbcu))
vRBCvMGADP = parMultiUninf*TsRBC*Vrbcu*(-((KdRBCvMGADP*mgadpRBC[t])/Vrbcu) + (KaappRBCvMGADP*adpRBC[t]*mgRBC[t])/Vrbcu^2)
vRBCvMGATP = parMultiUninf*TsRBC*Vrbcu*(-((KdRBCvMGATP*mgatpRBC[t])/Vrbcu) + (KaappRBCvMGATP*atpRBC[t]*mgRBC[t])/Vrbcu^2)
vRBCvMGB13PG = parMultiUninf*TsRBC*Vrbcu*(-((KdRBCvMGB13PG*mgb13pgRBC[t])/Vrbcu) + (KaappRBCvMGB13PG*b13pgRBC[t]*mgRBC[t])/Vrbcu^2)
vRBCvMGB23PG = parMultiUninf*TsRBC*Vrbcu*(-((KdRBCvMGB23PG*mgb23pgRBC[t])/Vrbcu) + (KaappRBCvMGB23PG*b23pgRBC[t]*mgRBC[t])/Vrbcu^2)
vRBCvMGF16P2 = parMultiUninf*TsRBC*Vrbcu*(-((KdRBCvMGF16P2*mgf16p2RBC[t])/Vrbcu) + (KaappRBCvMGF16P2*f16p2RBC[t]*mgRBC[t])/Vrbcu^2)
vRBCvMGG16P2 = parMultiUninf*TsRBC*Vrbcu*(-((KdRBCvMGG16P2*mgg16p2RBC[t])/Vrbcu) + (KaappRBCvMGG16P2*g16p2RBC[t]*mgRBC[t])/Vrbcu^2)
vRBCvMGPHOS = parMultiUninf*TsRBC*Vrbcu*(-((KdRBCvMGPHOS*mgphosRBC[t])/Vrbcu) + (KaappRBCvMGPHOS*mgRBC[t]*phosRBC[t])/Vrbcu^2)
vRBCvOX = KRBCvOX*parMultiUninf*TsRBC*gshRBC[t]
vRBCvOXNADH = KRBCvOXNADH*parMultiUninf*TsRBC*nadhRBC[t]
vRBCvP6GDH = (parMultiUninf*TsRBC*Vrbcu*VRBCvP6GDH*((K11RBCvP6GDH*K1RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH*nadpRBC[t]*p6gRBC[t])/Vrbcu^2 - (co2RBC*K10RBCvP6GDH*K12RBCvP6GDH*K2RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*nadphRBC[t]*ru5pRBC[t])/Vrbcu^3))/(K11RBCvP6GDH*K2RBCvP6GDH*(K4RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K7RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH)*K9RBCvP6GDH + (co2RBC*K11RBCvP6GDH*K2RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH)/Vrbcu + (co2RBC*K12RBCvP6GDH*K2RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*nadphRBC[t])/Vrbcu^2 + (K12RBCvP6GDH*K2RBCvP6GDH*(K4RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K7RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH)*K9RBCvP6GDH*nadphRBC[t])/Vrbcu + (co2RBC*K11RBCvP6GDH*K1RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*nadpRBC[t])/Vrbcu^2 + (K11RBCvP6GDH*K1RBCvP6GDH*(K4RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K7RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH)*K9RBCvP6GDH*nadpRBC[t])/Vrbcu + (K11RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH*p6gRBC[t])/Vrbcu + (K12RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH*nadphRBC[t]*p6gRBC[t])/Vrbcu^2 + (co2RBC*K11RBCvP6GDH*K1RBCvP6GDH*K3RBCvP6GDH*(K5RBCvP6GDH + K6RBCvP6GDH)*K8RBCvP6GDH*nadpRBC[t]*p6gRBC[t])/Vrbcu^3 + (K1RBCvP6GDH*K3RBCvP6GDH*(K11RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH + K11RBCvP6GDH*K5RBCvP6GDH*K9RBCvP6GDH + K11RBCvP6GDH*K6RBCvP6GDH*K9RBCvP6GDH + K11RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH*K9RBCvP6GDH)*nadpRBC[t]*p6gRBC[t])/Vrbcu^2 + (co2RBC*K10RBCvP6GDH*K2RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*ru5pRBC[t])/Vrbcu^2 + (co2RBC*K10RBCvP6GDH*K12RBCvP6GDH*(K2RBCvP6GDH*K4RBCvP6GDH + K2RBCvP6GDH*K5RBCvP6GDH + K2RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K6RBCvP6GDH)*K8RBCvP6GDH*nadphRBC[t]*ru5pRBC[t])/Vrbcu^3 + (K10RBCvP6GDH*K12RBCvP6GDH*K2RBCvP6GDH*(K4RBCvP6GDH*K6RBCvP6GDH + K4RBCvP6GDH*K7RBCvP6GDH + K5RBCvP6GDH*K7RBCvP6GDH)*nadphRBC[t]*ru5pRBC[t])/Vrbcu^2 + (co2RBC*K10RBCvP6GDH*K1RBCvP6GDH*K4RBCvP6GDH*K6RBCvP6GDH*K8RBCvP6GDH*nadpRBC[t]*ru5pRBC[t])/Vrbcu^3 + (co2RBC*K10RBCvP6GDH*K12RBCvP6GDH*K3RBCvP6GDH*(K5RBCvP6GDH + K6RBCvP6GDH)*K8RBCvP6GDH*nadphRBC[t]*p6gRBC[t]*ru5pRBC[t])/Vrbcu^4 + (K10RBCvP6GDH*K12RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*nadphRBC[t]*p6gRBC[t]*ru5pRBC[t])/Vrbcu^3 + (co2RBC*K10RBCvP6GDH*K1RBCvP6GDH*K3RBCvP6GDH*(K5RBCvP6GDH + K6RBCvP6GDH)*K8RBCvP6GDH*nadpRBC[t]*p6gRBC[t]*ru5pRBC[t])/Vrbcu^4 + (K10RBCvP6GDH*K1RBCvP6GDH*K3RBCvP6GDH*K5RBCvP6GDH*K7RBCvP6GDH*nadpRBC[t]*p6gRBC[t]*ru5pRBC[t])/Vrbcu^3)
vRBCvPFK = (parMultiUninf*TsRBC*Vrbcu*VRBCvPFK*(-((KcatfRBCvPFK*f16p2RBC[t]*mgadpRBC[t])/(Kmf16p2RBCvPFK*KmmgadpRBCvPFK*Vrbcu^2)) + (KcatfRBCvPFK*f6pRBC[t]*mgatpRBC[t])/(Kmf6pRBCvPFK*KmmgatpRBCvPFK*Vrbcu^2)))/((1 + f16p2RBC[t]/(Kmf16p2RBCvPFK*Vrbcu) + f6pRBC[t]/(Kmf6pRBCvPFK*Vrbcu) + mgadpRBC[t]/(KmmgadpRBCvPFK*Vrbcu) + (f16p2RBC[t]*mgadpRBC[t])/(Kmf16p2RBCvPFK*KmmgadpRBCvPFK*Vrbcu^2) + mgatpRBC[t]/(KmmgatpRBCvPFK*Vrbcu) + (f6pRBC[t]*mgatpRBC[t])/(Kmf6pRBCvPFK*KmmgatpRBCvPFK*Vrbcu^2))*(1 + ((hRBC/KaRBCvPFK)^nRBCvPFK*(1 + atpRBC[t]/(KtatpRBCvPFK*Vrbcu))^4*(1 + b23pgRBC[t]/(Ktb23pgRBCvPFK*Vrbcu))^4*(1 + mgRBC[t]/(KtmgRBCvPFK*Vrbcu))^4)/((1 + ampRBC[t]/(KrampRBCvPFK*Vrbcu))^4*(1 + f16p2RBC[t]/(Kmf16p2RBCvPFK*Vrbcu) + f6pRBC[t]/(Kmf6pRBCvPFK*Vrbcu))^4*(1 + g16p2RBC[t]/(Krg16p2RBCvPFK*Vrbcu))^4*(1 + phosRBC[t]/(KrphosRBCvPFK*Vrbcu))^4)))
vRBCvPGI = (ERBCvPGI*parMultiUninf*TsRBC*Vrbcu*(-((KcatrRBCvPGI*f6pRBC[t])/(Kmf6pRBCvPGI*Vrbcu)) + (KcatfRBCvPGI*g6pRBC[t])/(Kmg6pRBCvPGI*Vrbcu)))/(1 + f6pRBC[t]/(Kmf6pRBCvPGI*Vrbcu) + g6pRBC[t]/(Kmg6pRBCvPGI*Vrbcu))
vRBCvPGK = (parMultiUninf*TsRBC*Vrbcu*VRBCvPGK*((KcatfRBCvPGK*b13pgRBC[t]*mgadpRBC[t])/(KimgadpRBCvPGK*Kmb13pgRBCvPGK*Vrbcu^2) - (KcatrRBCvPGK*mgatpRBC[t]*p3gRBC[t])/(KimgatpRBCvPGK*Kmp3gRBCvPGK*Vrbcu^2)))/(1 + b13pgRBC[t]/(Kib13pgRBCvPGK*Vrbcu) + mgadpRBC[t]/(KimgadpRBCvPGK*Vrbcu) + (b13pgRBC[t]*mgadpRBC[t])/(KimgadpRBCvPGK*Kmb13pgRBCvPGK*Vrbcu^2) + mgatpRBC[t]/(KimgatpRBCvPGK*Vrbcu) + p3gRBC[t]/(Kip3gRBCvPGK*Vrbcu) + (mgatpRBC[t]*p3gRBC[t])/(KimgatpRBCvPGK*Kmp3gRBCvPGK*Vrbcu^2))
vRBCvPGLHYDROLYSIS = parMultiUninf*TsRBC*Vrbcu*((KhydrolRBCvPGLHYDROLYSIS*p6glRBC[t])/Vrbcu + (Klactonase3RBCvPGLHYDROLYSIS*VRBCvPGLHYDROLYSIS*p6glRBC[t])/(Vrbcu*((Klactonase2RBCvPGLHYDROLYSIS + Klactonase3RBCvPGLHYDROLYSIS)/Klactonase1RBCvPGLHYDROLYSIS + p6glRBC[t]/Vrbcu)))
vRBCvPGM = (parMultiUninf*TsRBC*Vrbcu*VRBCvPGM*(-((KcatrRBCvPGM*p2gRBC[t])/(Kmp2gRBCvPGM*Vrbcu)) + (KcatfRBCvPGM*p3gRBC[t])/(Kmp3gRBCvPGM*Vrbcu)))/(1 + p2gRBC[t]/(Kmp2gRBCvPGM*Vrbcu) + p3gRBC[t]/(Kmp3gRBCvPGM*Vrbcu))
vRBCvPHOSTRANSPORT = parMultiUninf*TsRBC*vBld*(-((KiRBCvPHOSTRANSPORT*phosEXT)/vBld) + (KoRBCvPHOSTRANSPORT*phosRBC[t])/Vrbcu)
vRBCvPK = (parMultiUninf*TsRBC*Vrbcu*VRBCvPK*((KcatfRBCvPK*mgadpRBC[t]*pepRBC[t])/(KrmgadpRBCvPK*KrpepRBCvPK*Vrbcu^2) - (KcatrRBCvPK*mgatpRBC[t]*pyrRBC[t])/(KrmgatpRBCvPK*KrpyrRBCvPK*Vrbcu^2)))/((1 + LRBCvPK)*(1 + mgadpRBC[t]/(KrmgadpRBCvPK*Vrbcu) + mgatpRBC[t]/(KrmgatpRBCvPK*Vrbcu) + pepRBC[t]/(KrpepRBCvPK*Vrbcu) + (mgadpRBC[t]*pepRBC[t])/(KrmgadpRBCvPK*KrpepRBCvPK*Vrbcu^2) + pyrRBC[t]/(KrpyrRBCvPK*Vrbcu) + (mgatpRBC[t]*pyrRBC[t])/(KrmgatpRBCvPK*KrpyrRBCvPK*Vrbcu^2)))
vRBCvPYRTRANSPORT = parMultiUninf*TsRBC*vBld*(-((KiRBCvPYRTRANSPORT*pyrEXT)/vBld) + (KoRBCvPYRTRANSPORT*pyrRBC[t])/Vrbcu)
vRBCvR5PI = (parMultiUninf*TsRBC*Vrbcu*VRBCvR5PI*(-((K2RBCvR5PI*K4RBCvR5PI*rib5pRBC[t])/((K2RBCvR5PI + K3RBCvR5PI)*Vrbcu)) + (K1RBCvR5PI*K3RBCvR5PI*ru5pRBC[t])/((K2RBCvR5PI + K3RBCvR5PI)*Vrbcu)))/(1 + (K1RBCvR5PI*ru5pRBC[t])/((K2RBCvR5PI + K3RBCvR5PI)*Vrbcu) + (K4RBCvR5PI*xu5pRBC[t])/((K2RBCvR5PI + K3RBCvR5PI)*Vrbcu))
vRBCvRu5PE = (parMultiUninf*TsRBC*Vrbcu*VRBCvRu5PE*((K1RBCvRu5PE*K3RBCvRu5PE*ru5pRBC[t])/((K2RBCvRu5PE + K3RBCvRu5PE)*Vrbcu) - (K2RBCvRu5PE*K4RBCvRu5PE*xu5pRBC[t])/((K2RBCvRu5PE + K3RBCvRu5PE)*Vrbcu)))/(1 + (K1RBCvRu5PE*ru5pRBC[t])/((K2RBCvRu5PE + K3RBCvRu5PE)*Vrbcu) + (K4RBCvRu5PE*xu5pRBC[t])/((K2RBCvRu5PE + K3RBCvRu5PE)*Vrbcu))
vRBCvTA = (parMultiUninf*TsRBC*Vrbcu*VRBCvTA*(-((K2RBCvTA*K4RBCvTA*K6RBCvTA*K8RBCvTA*ery4pRBC[t]*f6pRBC[t])/Vrbcu^2) + (K1RBCvTA*K3RBCvTA*K5RBCvTA*K7RBCvTA*gapRBC[t]*sed7pRBC[t])/Vrbcu^2))/((K2RBCvTA*K4RBCvTA*(K6RBCvTA + K7RBCvTA)*ery4pRBC[t])/Vrbcu + ((K2RBCvTA + K3RBCvTA)*K6RBCvTA*K8RBCvTA*f6pRBC[t])/Vrbcu + (K4RBCvTA*(K2RBCvTA + K6RBCvTA)*K8RBCvTA*ery4pRBC[t]*f6pRBC[t])/Vrbcu^2 + ((K2RBCvTA + K3RBCvTA)*K5RBCvTA*K7RBCvTA*gapRBC[t])/Vrbcu + ((K2RBCvTA + K3RBCvTA)*K5RBCvTA*K8RBCvTA*f6pRBC[t]*gapRBC[t])/Vrbcu^2 + (K1RBCvTA*K3RBCvTA*(K6RBCvTA + K7RBCvTA)*sed7pRBC[t])/Vrbcu + (K1RBCvTA*K4RBCvTA*(K6RBCvTA + K7RBCvTA)*ery4pRBC[t]*sed7pRBC[t])/Vrbcu^2 + (K1RBCvTA*K5RBCvTA*(K3RBCvTA + K7RBCvTA)*gapRBC[t]*sed7pRBC[t])/Vrbcu^2)
vRBCvTIM = (parMultiUninf*TsRBC*Vrbcu*VRBCvTIM*(-((KcatrRBCvTIM*dhapRBC[t])/(KmdhapRBCvTIM*Vrbcu)) + (KcatfRBCvTIM*gapRBC[t])/(KmgapRBCvTIM*Vrbcu)))/(1 + dhapRBC[t]/(KmdhapRBCvTIM*Vrbcu) + gapRBC[t]/(KmgapRBCvTIM*Vrbcu))
vRBCvTK1 = parMultiUninf*TsRBC*Vrbcu*(-((K2RBCvTK1*tkxu5pRBC[t])/Vrbcu) + (K1RBCvTK1*tkRBC[t]*xu5pRBC[t])/Vrbcu^2)
vRBCvTK2 = parMultiUninf*TsRBC*Vrbcu*(-((K4RBCvTK2*gapRBC[t]*tkgRBC[t])/Vrbcu^2) + (K3RBCvTK2*tkxu5pRBC[t])/Vrbcu)
vRBCvTK3 = parMultiUninf*TsRBC*Vrbcu*((K5RBCvTK3*rib5pRBC[t]*tkgRBC[t])/Vrbcu^2 - (K6RBCvTK3*tkgrib5pRBC[t])/Vrbcu)
vRBCvTK4 = parMultiUninf*TsRBC*Vrbcu*((K7RBCvTK4*tkgrib5pRBC[t])/Vrbcu - (K8RBCvTK4*sed7pRBC[t]*tkRBC[t])/Vrbcu^2)
vRBCvTK5 = parMultiUninf*TsRBC*Vrbcu*(-((K10RBCvTK5*tkgery4pRBC[t])/Vrbcu) + (K9RBCvTK5*ery4pRBC[t]*tkgRBC[t])/Vrbcu^2)
vRBCvTK6 = parMultiUninf*TsRBC*Vrbcu*((K11RBCvTK6*tkgery4pRBC[t])/Vrbcu - (K12RBCvTK6*f6pRBC[t]*tkRBC[t])/Vrbcu^2)
end rate equations

begin parameters
CRBCvMGADP = 0.00062
CRBCvMGATP = "8.4e-05"
CRBCvMGB13PG = 0.0032
CRBCvMGB23PG = 0.0032
CRBCvMGF16P2 = 0.0083
CRBCvMGG16P2 = 0.0083
CRBCvMGPHOS = 40800.0
ConcADPRBC = 0.00031
ConcAMPRBC = "3e-05"
ConcATPRBC = 0.0021
ConcAdpPF = 0.001
ConcAtpPF = 0.002
ConcB13PGRBC = "7e-07"
ConcB13pgPF = 0.0001
ConcB23PGRBC = 0.0067
ConcBPGSPB13PGRBC = 0.0
ConcBPGSPB23PGRBC = 0.0
ConcBPGSPPP2GRBC = 0.0
ConcBPGSPPP3GRBC = 0.0
ConcBPGSPPRBC = 0.0
ConcBPGSPRBC = "3.8e-06"
ConcBpgsppphosRBC = 0.0
ConcCo2RBC = 0.0012
ConcDhapPF = 0.0001
ConcDhapRBC = "1.9e-05"
ConcEry4PRBC = "1e-05"
ConcF16P2RBC = "2.7e-06"
ConcF16bpPF = 0.0001
ConcF6PRBC = "1.3e-05"
ConcF6pPF = 0.0001
ConcG16P2RBC = 0.000122
ConcG6PRBC = "4e-05"
ConcG6pPF = 0.0001
ConcGSHRBC = 0.0032
ConcGapPF = 0.0001
ConcGapRBC = "5.7e-06"
ConcGlcEXT = 0.005
ConcGlcPF = 0.0001
ConcGlcRBC = 0.001
ConcGly3pPF = 0.0001
ConcGlyEXT = 0.0
ConcGssgRBC = "9e-08"
ConcHbADPRBC = 0.0
ConcHbATPRBC = 0.0
ConcHbB13PGRBC = 0.0
ConcHbB23PGRBC = 0.0
ConcHbRBC = 0.007
ConcHbmgatpRBC = 0.0
ConcLacEXT = 0.00182
ConcLacPF = 0.0001
ConcLacRBC = 0.0014
ConcMgADPRBC = 0.0
ConcMgB13PGRBC = 0.0
ConcMgB23PGRBC = 0.0
ConcMgF16P2RBC = 0.0
ConcMgRBC = 0.003
ConcMgatpRBC = 0.0
ConcMgg16p2RBC = 0.0
ConcMgphosRBC = 0.0
ConcNADPRBC = "1.25e-07"
ConcNADRBC = "6e-05"
ConcNadPF = 0.0025
ConcNadhPF = 0.0005
ConcNadhRBC = "1.4e-07"
ConcNadphRBC = "6.4e-05"
ConcP2GARBC = "1e-05"
ConcP2gPF = 0.0001
ConcP3GARBC = "6.4e-05"
ConcP3gPF = 0.0001
ConcP6GLRBC = "1.4e-10"
ConcP6GRBC = "1.4e-07"
ConcPEPRBC = "2.3e-05"
ConcPepPF = 0.0001
ConcPhosEXT = 0.00192
ConcPhosRBC = 0.001
ConcPyrEXT = "8.5e-05"
ConcPyrPF = 0.0001
ConcPyrRBC = "6e-05"
ConcRib5PRBC = "1e-05"
ConcRu5PRBC = "1e-05"
ConcSed7PRBC = "1e-05"
ConcTKGEry4PRBC = 0.0
ConcTKGRBC = 0.0
ConcTKGRib5PRBC = 0.0
ConcTKRBC = "3.3e-07"
ConcTkxu5pRBC = 0.0
ConcXu5PRBC = "1e-06"
ERBCvHK = "2.5e-08"
ERBCvPGI = "2.18e-07"
K10RBCvBPGSP6 = 1979.0
K10RBCvG6PDH = 1400000000.0
K10RBCvGSSGR = 50000000.0
K10RBCvP6GDH = 225000.0
K10RBCvTK5 = 175.0
K11RBCvBPGSP6 = 0.01
K11RBCvGSSGR = 7000.0
K11RBCvP6GDH = 300.0
K11RBCvTK6 = 40.0
K12RBCvBPGSP7 = 1000.0
K12RBCvGSSGR = 100000000.0
K12RBCvP6GDH = 4950000.0
K12RBCvTK6 = 21300.0
K13RBCvBPGSP7 = 1800000.0
K14RBCvBPGSP8 = 1000000000.0
K15RBCvBPGSP8 = 610000.0
K16RBCvBPGSP9 = 0.19
K1RBCvAK = 4300.0
K1RBCvBPGSP1 = 80000000.0
K1RBCvG6PDH = 110000000.0
K1RBCvGSSGR = 85000000.0
K1RBCvP6GDH = 2400000.0
K1RBCvR5PI = 60900.0
K1RBCvRu5PE = 3910000.0
K1RBCvTA = 580000.0
K1RBCvTK1 = 216000.0
K2RBCvAK = 1400.0
K2RBCvBPGSP1 = 400.0
K2RBCvG6PDH = 870.0
K2RBCvGSSGR = 510.0
K2RBCvP6GDH = 410.0
K2RBCvR5PI = 33.3
K2RBCvRu5PE = 438.0
K2RBCvTA = 45.3
K2RBCvTK1 = 38.0
K3RBCvBPGSP2 = 9.9
K3RBCvG6PDH = 26000000.0
K3RBCvGSSGR = 1000000000.0
K3RBCvP6GDH = 2000000000.0
K3RBCvR5PI = 14.2
K3RBCvRu5PE = 305.0
K3RBCvTA = 16.3
K3RBCvTK2 = 34.0
K4RBCvBPGSP3 = 185000000.0
K4RBCvG6PDH = 300.0
K4RBCvGSSGR = 72000.0
K4RBCvP6GDH = 26000.0
K4RBCvR5PI = 21600.0
K4RBCvRu5PE = 1490000.0
K4RBCvTA = 1010000.0
K4RBCvTK2 = 156000.0
K5RBCvBPGSP3 = 100000000.0
K5RBCvG6PDH = 750.0
K5RBCvGSSGR = 810.0
K5RBCvP6GDH = 48.0
K5RBCvTA = 490000.0
K5RBCvTK3 = 329000.0
K6RBCvBPGSP4 = 1000.0
K6RBCvG6PDH = 2000.0
K6RBCvGSSGR = 1000.0
K6RBCvP6GDH = 30.0
K6RBCvTA = 60.0
K6RBCvTK3 = 175.0
K7RBCvBPGSP4 = 1000.0
K7RBCvG6PDH = 220000.0
K7RBCvGSSGR = 1000000.0
K7RBCvP6GDH = 630.0
K7RBCvTA = 17.0
K7RBCvTK4 = 40.0
K8RBCvBPGSP5 = 10000.0
K8RBCvG6PDH = 1100000000.0
K8RBCvGSSGR = 50000000.0
K8RBCvP6GDH = 36000.0
K8RBCvTA = 79000.0
K8RBCvTK4 = 44800.0
K9RBCvBPGSP5 = 0.55
K9RBCvG6PDH = 10000.0
K9RBCvGSSGR = 1000000.0
K9RBCvP6GDH = 800.0
K9RBCvTK5 = 2240000.0
KRBCvATPASE = 0.000585
KRBCvMGADP = 3290.0
KRBCvMGATP = 43200.0
KRBCvMGB13PG = 7410.0
KRBCvMGB23PG = 7410.0
KRBCvMGF16P2 = 363.0
KRBCvMGG16P2 = 363.0
KRBCvOX = "3.4e-05"
KRBCvOXNADH = 0.0163
KaRBCvHBADP = 300000.0
KaRBCvHBATP = 432000.0
KaRBCvHBB13PG = 380000.0
KaRBCvHBBPG = 300000.0
KaRBCvHBMGATP = 46800.0
KaRBCvMGADP = 2760000.0
KaRBCvMGATP = 31200000.0
KaRBCvMGB13PG = 228000.0
KaRBCvMGB23PG = 804000.0
KaRBCvMGF16P2 = 480000.0
KaRBCvMGG16P2 = 480000.0
KaRBCvMGPHOS = 40800.0
KaRBCvPFK = "8.91251e-08"
KadpPFvHK = 0.000846735
KadpPFvPFK = 0.00074176
KadpPFvPGK = 0.00015
KadpPFvPK = 0.000317
KahbRBC = "2.51189e-07"
KatpPFvHK = 0.00069647
KatpPFvPFK = 0.0007862
KatpPFvPGK = 0.00077
Kb13pgPFvGAPDH = "2.359e-05"
Kb13pgPFvPGK = "1.34e-05"
KcatfRBCvALD = 68.0
KcatfRBCvENO = 190.0
KcatfRBCvGAPDH = 695.304
KcatfRBCvHK = 299.16
KcatfRBCvLDH = 458.0
KcatfRBCvPFK = 822.0
KcatfRBCvPGI = 1470.0
KcatfRBCvPGK = 2290.0
KcatfRBCvPGM = 795.0
KcatfRBCvPK = 1386.0
KcatfRBCvTIM = 14560.0
KcatrRBCvALD = 234.0
KcatrRBCvENO = 50.0
KcatrRBCvGAPDH = 512.487
KcatrRBCvHK = 1.92792
KcatrRBCvLDH = 115.0
KcatrRBCvPFK = 36.0
KcatrRBCvPGI = 1760.0
KcatrRBCvPGK = 917.0
KcatrRBCvPGM = 714.0
KcatrRBCvPK = 3.26
KcatrRBCvTIM = 1280.0
KdPFvLACtr = 0.005
KdPFvPYRtr = 0.0007
KdRBCvHBADP = 1200.0
KdRBCvHBATP = 1200.0
KdRBCvHBB13PG = 1200.0
KdRBCvHBBPG = 1200.0
KdRBCvHBMGATP = 1200.0
KdRBCvMGADP = 1200.0
KdRBCvMGATP = 1200.0
KdRBCvMGB13PG = 1200.0
KdRBCvMGB23PG = 1200.0
KdRBCvMGF16P2 = 1200.0
KdRBCvMGG16P2 = 1200.0
KdRBCvMGPHOS = 1200.0
KdhapPFvALD = 0.00011
KdhapPFvG3PDH = 0.00034
KdhapPFvTPI = 0.001954
KdibpgRBCvHK = 0.004
Kdig16p2RBCvHK = "3e-05"
Kdig6pRBCvHK = "1e-05"
KdigshRBCvHK = 0.003
KeqPFvALD = "9e-05"
KeqPFvENO = 4.6
KeqPFvG3PDH = 32600.0
KeqPFvHK = 1310.0
KeqPFvPGI = 0.33
KeqPFvPGK = 3200.0
KeqPFvPGM = 0.19
KeqPFvTPI = 0.04545
Kf16bpPFvALD = "6.84e-05"
Kf16bpPFvPFK = 0.003626
Kf6pPFvPFK = 0.00109454
Kf6pPFvPGI = "9.6651e-05"
KfRBCvLDHP = 0.00346
Kg3pPFvG3PDH = 0.00398
Kg6pPFvHK = "4.3e-05"
Kg6pPFvPGI = 0.00100774
KgapPFvALD = "4.6e-05"
KgapPFvGAPDH = 0.000917
KgapPFvTPI = 0.000337
KglcPFvHK = 0.000168613
Kh2bpgRBCvMGB13PG = 4270000.0
Kh2bpgRBCvMGB23PG = 4270000.0
Kh2fRBCvMGF16P2 = 1120000.0
Kh2fRBCvMGG16P2 = 1120000.0
KhadpRBCvAK = 5420000.0
KhadpRBCvMGADP = 5420000.0
KhampRBCvAK = 3090000.0
KhatpRBCvMGATP = 9070000.0
KhbpgRBCvMGB13PG = 162000000.0
KhbpgRBCvMGB23PG = 162000000.0
KhfRBCvMGF16P2 = 7560000.0
KhfRBCvMGG16P2 = 7560000.0
KhphosRBCvMGPHOS = 5680000.0
KhydrolRBCvPGLHYDROLYSIS = 0.00071
KiRBCvLACTRANSPORT = 0.0036
KiRBCvPHOSTRANSPORT = 0.00056
KiRBCvPYRTRANSPORT = 0.018
KiadpPFvPK = 0.002
KiatpPFvHK = 0.026
KiatpPFvPK = 0.00182
Kib13pgRBCvGAPDH = "4.55544e-21"
Kib13pgRBCvPGK = "1.6e-06"
Kib23pgRBCvALD = 0.0015
Kidb13pgRBCvGAPDH = "1e-06"
KidgapRBCvGAPDH = "3.1e-05"
KidhapRBCvALD = "1.1e-05"
KidpyrRBCvLDH = 0.000101
Kif16p2RBCvALD = "1.98e-05"
Kig6pRBCvHK = "4.7e-05"
KigapRBCvGAPDH = "4.76523e-19"
KiglcRBCvHK = "4.7e-05"
KilacPFvPYRtr = 0.000358
KilacRBCvLDH = 0.00733
KimgRBCvENO = "4.6e-05"
KimgadpRBCvHK = 0.001
KimgadpRBCvPGK = "8e-05"
KimgatpRBCvHK = 0.001
KimgatpRBCvPGK = 0.000186
KinadRBCvGAPDH = "4.5e-05"
KinadRBCvLDH = 0.000503
KinadhRBCvGAPDH = "1e-05"
KinadhRBCvLDH = "2.45e-06"
Kip2gRBCvENO = 0.00014
Kip3gRBCvPGK = 0.000205
KipepPFvPK = 0.00292
KipepPFvTPI = "1.59e-05"
KipepRBCvENO = 0.0001105
KiphosRBCvGAPDH = 0.00316
KipyrPFvLACtr = 0.00163
KipyrPFvPK = 0.105
KipyrRBCvLDH = 0.000228
KkadpRBCvAK = 4.8
KkadpRBCvMGADP = 4.8
KkampRBCvAK = 1.8
KkatpRBCvMGATP = 14.0
KkbpgRBCvMGB13PG = 85.1
KkbpgRBCvMGB23PG = 85.1
KkfRBCvMGF16P2 = 10.7
KkfRBCvMGG16P2 = 10.7
KkhbpgRBCvMGB13PG = 8.9
KkhbpgRBCvMGB23PG = 8.9
KkhfRBCvMGF16P2 = 3.3
KkhfRBCvMGG16P2 = 3.3
KkphosRBCvMGPHOS = 3.0
KlacPFvLACtr = 0.0038
KlacPFvLDH = 0.003611
Klactonase1RBCvPGLHYDROLYSIS = 13000000.0
Klactonase2RBCvPGLHYDROLYSIS = 1000.0
Klactonase3RBCvPGLHYDROLYSIS = 29.0
KmPFvATPASE = 0.0045
KmPFvGLCtr = 0.000213
Kmb13pgRBCvGAPDH = "6.71e-07"
Kmb13pgRBCvPGK = "2e-06"
KmdhapRBCvALD = "3.5e-05"
KmdhapRBCvTIM = 0.0001624
Kmf16p2RBCvALD = "7.1e-06"
Kmf16p2RBCvPFK = 0.0005
Kmf6pRBCvPFK = "7.5e-05"
Kmf6pRBCvPGI = "7.1e-05"
Kmg6pRBCvPGI = 0.000181
KmgapRBCvALD = 0.000189
KmgapRBCvGAPDH = "9.5e-05"
KmgapRBCvTIM = 0.000446
KmghadpRBCvMGADP = 107.0
KmghatpRBCvMGATP = 748.0
KmghbpgRBCvMGB13PG = 513.0
KmghbpgRBCvMGB23PG = 513.0
KmghfRBCvMGF16P2 = 89.0
KmghfRBCvMGG16P2 = 89.0
KmlacRBCvLDH = 0.00107
KmlacRBCvLDHP = 0.000414
KmmgRBCvENO = "4.6e-05"
KmmgadpRBCvHK = 0.001
KmmgadpRBCvPFK = 0.00054
KmmgadpRBCvPGK = 0.0001
KmmgatpRBCvHK = 0.001
KmmgatpRBCvPFK = "6.8e-05"
KmmgatpRBCvPGK = 0.001
KmnadRBCvGAPDH = "4.5e-05"
KmnadRBCvLDH = 0.000107
KmnadhRBCvGAPDH = "3.3e-06"
KmnadhRBCvLDH = "8.44e-06"
Kmp2gRBCvENO = 0.00014
Kmp2gRBCvPGM = "2.56e-05"
Kmp3gRBCvPGK = 0.0011
Kmp3gRBCvPGM = 0.000168
KmpepRBCvENO = 0.0001105
KmphosRBCvGAPDH = 0.00316
KmpyrRBCvLDH = 0.000137
KmpyrRBCvLDHP = 0.000414
KnadPFvG3PDH = 0.00051
KnadPFvGAPDH = 0.0005662
KnadPFvLDH = 0.000234
KnadhPFvG3PDH = "9e-05"
KnadhPFvGAPDH = "2.87e-05"
KnadhPFvLDH = "4.6e-05"
Kp2gPFvENO = 0.000521
Kp2gPFvPGM = 0.000318
Kp3gPFvPGK = 0.000267
Kp3gPFvPGM = 0.00173
KpepPFvENO = 0.00129
KpepPFvPK = 0.000406
KpyrPFvLDH = 0.000133
KpyrPFvPYRtr = 0.0157
KrRBCvLDHP = "5.43e-07"
KrampRBCvPFK = 0.0003
Krf16p2RBCvPK = "5e-06"
Krg16p2RBCvPFK = 0.01
Krg16p2RBCvPK = 0.0001
KrmgadpRBCvPK = 0.000474
KrmgatpRBCvPK = 0.003
KrpepRBCvPK = 0.000225
KrphosRBCvPFK = 0.03
KrpyrRBCvPK = 0.002
KtatpRBCvPFK = 0.0001
KtatpRBCvPK = 0.00339
Ktb23pgRBCvPFK = 0.005
KtmgRBCvPFK = 0.004
LPFvPK = 0.255
RtvRBC = 0.69
TsRBC = 60.0
VPFvALD = 0.0569593147718
VPFvATPASE = 0.345
VPFvENO = 0.286295503195
VPFvGAPDH = 0.656263383259
VPFvGLCtr = 0.0469
VPFvHK = 0.083
VPFvLACtr = 0.598
VPFvPFK = 0.41
VPFvPGI = 0.800428265477
VPFvPGM = 0.186937901488
VPFvPK = 0.762
VPFvPYRtr = 0.216
VPFvTPI = 1.50599571726
VRBCvALD = "3.7e-07"
VRBCvENO = "2.2e-07"
VRBCvG6PDH = "9.3e-08"
VRBCvGAPDH = "7.66e-06"
VRBCvGLCTRANSPORT = 0.047
VRBCvGSSGR = "1.25e-07"
VRBCvLDH = "3.43e-06"
VRBCvP6GDH = "2.1e-06"
VRBCvPFK = "1.1e-07"
VRBCvPGK = "2.74e-06"
VRBCvPGLHYDROLYSIS = "1.4e-05"
VRBCvPGM = "4.1e-07"
VRBCvPK = "8.7e-08"
VRBCvR5PI = "1.42e-05"
VRBCvRu5PE = "4.22e-06"
VRBCvTA = "6.9e-07"
VRBCvTIM = "1.14e-06"
VfPFvG3PDH = 0.012
VfPFvLDH = 0.542612419668
VrPFvGAPDH = 0.243199143455
VrPFvLDH = 0.290792291203
VrPFvPGK = 0.353533190557
alphaPFvGLCtr = 0.91
bpgxRBC = "3.8e-06"
cyt = 0.0
glcAlpha = 0.36
gssgxRBC = 0.00160009
h = 4.0
hbxRBC = 0.007
kAlpha = 0.0083
kBeta = 0.007
kPFvGLYtr = 2.0
kRBC = 0.15
ki1 = "8.86746e-08"
ki2 = "2.08e-05"
nRBCvPFK = 5.0
nadpxRBC = "6.4125e-05"
pHConversionFactor = 1.0
par = 1.0
phPF = 7.2
phRBC = 7.2
tkxRBC = "3.3e-07"
trophToRBC = 0.355555555556
xatpRBC = -0.00456
xg16p2RBC = 0.000122
xphosRBC = 0.000438
Vrbc = 90.0
vBld = 128.5
end parameters

begin initial conditions
adpPF[0] = adpPFi
adpRBC[0] = adpRBCi
adpRBCi[0] = adpRBCii
ampRBC[0] = ampRBCi
ampRBCi[0] = ampRBCii
atpPF[0] = atpPFi
atpRBC[0] = atpRBCi
atpRBCi[0] = atpRBCii
b13pgPF[0] = b13pgPFi
b13pgRBC[0] = b13pgRBCi
b13pgRBCi[0] = b13pgRBCii
b23pgRBC[0] = b23pgRBCi
b23pgRBCi[0] = b23pgRBCii
bpgspRBC[0] = bpgspRBCi
bpgspRBCi[0] = bpgspRBCii
bpgspb13pgRBC[0] = bpgspb13pgRBCi
bpgspb13pgRBCi[0] = bpgspb13pgRBCii
bpgspb23pgRBC[0] = bpgspb23pgRBCi
bpgspb23pgRBCi[0] = bpgspb23pgRBCii
bpgsppRBC[0] = bpgsppRBCi
bpgsppRBCi[0] = bpgsppRBCii
bpgsppp2gRBC[0] = bpgsppp2gRBCi
bpgsppp2gRBCi[0] = bpgsppp2gRBCii
bpgsppp3gRBC[0] = bpgsppp3gRBCi
bpgsppp3gRBCi[0] = bpgsppp3gRBCii
bpgsppphosRBC[0] = bpgsppphosRBCi
bpgsppphosRBCi[0] = bpgsppphosRBCii
dhapPF[0] = dhapPFi
dhapRBC[0] = dhapRBCi
dhapRBCi[0] = dhapRBCii
ery4pRBC[0] = ery4pRBCi
ery4pRBCi[0] = ery4pRBCii
f16bpPF[0] = f16bpPFi
f16p2RBC[0] = f16p2RBCi
f16p2RBCi[0] = f16p2RBCii
f6pPF[0] = f6pPFi
f6pRBC[0] = f6pRBCi
f6pRBCi[0] = f6pRBCii
g16p2RBC[0] = g16p2RBCi
g16p2RBCi[0] = g16p2RBCii
g3pPF[0] = g3pPFi
g6pPF[0] = g6pPFi
g6pRBC[0] = g6pRBCi
g6pRBCi[0] = g6pRBCii
gapPF[0] = gapPFi
gapRBC[0] = gapRBCi
gapRBCi[0] = gapRBCii
glcPF[0] = glcPFi
glcRBC[0] = glcRBCi
glcRBCi[0] = glcRBCii
gshRBC[0] = gshRBCi
gshRBCi[0] = gshRBCii
gssgRBC[0] = gssgRBCi
gssgRBCi[0] = gssgRBCii
hbRBC[0] = hbRBCi
hbRBCi[0] = hbRBCii
hbadpRBC[0] = hbadpRBCi
hbadpRBCi[0] = hbadpRBCii
hbatpRBC[0] = hbatpRBCi
hbatpRBCi[0] = hbatpRBCii
hbb13pgRBC[0] = hbb13pgRBCi
hbb13pgRBCi[0] = hbb13pgRBCii
hbb23pgRBC[0] = hbb23pgRBCi
hbb23pgRBCi[0] = hbb23pgRBCii
hbmgatpRBC[0] = hbmgatpRBCi
hbmgatpRBCi[0] = hbmgatpRBCii
lacPF[0] = lacPFi
lacRBC[0] = lacRBCi
lacRBCi[0] = lacRBCii
mgRBC[0] = mgRBCi
mgRBCi[0] = mgRBCii
mgadpRBC[0] = mgadpRBCi
mgadpRBCi[0] = mgadpRBCii
mgatpRBC[0] = mgatpRBCi
mgatpRBCi[0] = mgatpRBCii
mgb13pgRBC[0] = mgb13pgRBCi
mgb13pgRBCi[0] = mgb13pgRBCii
mgb23pgRBC[0] = mgb23pgRBCi
mgb23pgRBCi[0] = mgb23pgRBCii
mgf16p2RBC[0] = mgf16p2RBCi
mgf16p2RBCi[0] = mgf16p2RBCii
mgg16p2RBC[0] = mgg16p2RBCi
mgg16p2RBCi[0] = mgg16p2RBCii
mgphosRBC[0] = mgphosRBCi
mgphosRBCi[0] = mgphosRBCii
nadPF[0] = nadPFi
nadRBC[0] = nadRBCi
nadRBCi[0] = nadRBCii
nadhPF[0] = nadhPFi
nadhRBC[0] = nadhRBCi
nadhRBCi[0] = nadhRBCii
nadpRBC[0] = nadpRBCi
nadpRBCi[0] = nadpRBCii
nadphRBC[0] = nadphRBCi
nadphRBCi[0] = nadphRBCii
p2gPF[0] = p2gPFi
p2gRBC[0] = p2gRBCi
p2gRBCi[0] = p2gRBCii
p3gPF[0] = p3gPFi
p3gRBC[0] = p3gRBCi
p3gRBCi[0] = p3gRBCii
p6gRBC[0] = p6gRBCi
p6gRBCi[0] = p6gRBCii
p6glRBC[0] = p6glRBCi
p6glRBCi[0] = p6glRBCii
pepPF[0] = pepPFi
pepRBC[0] = pepRBCi
pepRBCi[0] = pepRBCii
phosRBC[0] = phosRBCi
phosRBCi[0] = phosRBCii
pyrPF[0] = pyrPFi
pyrRBC[0] = pyrRBCi
pyrRBCi[0] = pyrRBCii
rib5pRBC[0] = rib5pRBCi
rib5pRBCi[0] = rib5pRBCii
ru5pRBC[0] = ru5pRBCi
ru5pRBCi[0] = ru5pRBCii
sed7pRBC[0] = sed7pRBCi
sed7pRBCi[0] = sed7pRBCii
tkRBC[0] = tkRBCi
tkRBCi[0] = tkRBCii
tkgRBC[0] = tkgRBCi
tkgRBCi[0] = tkgRBCii
tkgery4pRBC[0] = tkgery4pRBCi
tkgery4pRBCi[0] = tkgery4pRBCii
tkgrib5pRBC[0] = tkgrib5pRBCi
tkgrib5pRBCi[0] = tkgrib5pRBCii
tkxu5pRBC[0] = tkxu5pRBCi
tkxu5pRBCi[0] = tkxu5pRBCii
xu5pRBC[0] = xu5pRBCi
xu5pRBCi[0] = xu5pRBCii
end initial conditions

begin initial values
adpPFi = 0.032000000000039996
adpRBCi = 0.00031
adpRBCii = 0.0279
ampRBCi = 3/100000
ampRBCii = 27/10000
atpPFi = 0.06400000000007999
atpRBCi = 0.0021
atpRBCii = 0.189
b13pgPFi = 0.003200000000004
b13pgRBCi = 7/10000000
b13pgRBCii = 63/1000000
b23pgRBCi = 0.0067
b23pgRBCii = 0.603
bpgspRBCi = 3.7999999999999996*^-6
bpgspRBCii = 0.00034199999999999996
bpgspb13pgRBCi = 0
bpgspb13pgRBCii = 0
bpgspb23pgRBCi = 0
bpgspb23pgRBCii = 0
bpgsppRBCi = 0
bpgsppRBCii = 0
bpgsppp2gRBCi = 0
bpgsppp2gRBCii = 0
bpgsppp3gRBCi = 0
bpgsppp3gRBCii = 0
bpgsppphosRBCi = 0
bpgsppphosRBCii = 0
dhapPFi = 0.003200000000004
dhapRBCi = 0.000019
dhapRBCii = 0.0017100000000000001
ery4pRBCi = 1/100000
ery4pRBCii = 9/10000
f16bpPFi = 0.003200000000004
f16p2RBCi = 2.7*^-6
f16p2RBCii = 0.00024300000000000002
f6pPFi = 0.003200000000004
f6pRBCi = 0.000013000000000000001
f6pRBCii = 0.00117
g16p2RBCi = 0.000122
g16p2RBCii = 0.01098
g3pPFi = 0.003200000000004
g6pPFi = 0.003200000000004
g6pRBCi = 1/25000
g6pRBCii = 9/2500
gapPFi = 0.003200000000004
gapRBCi = 5.7*^-6
gapRBCii = 0.000513
glcPFi = 0.003200000000004
glcRBCi = 0.001
glcRBCii = 0.09
gshRBCi = 0.0032
gshRBCii = 0.28800000000000003
gssgRBCi = 9/100000000
gssgRBCii = 81/10000000
hbRBCi = 0.007
hbRBCii = 0.63
hbadpRBCi = 0
hbadpRBCii = 0
hbatpRBCi = 0
hbatpRBCii = 0
hbb13pgRBCi = 0
hbb13pgRBCii = 0
hbb23pgRBCi = 0
hbb23pgRBCii = 0
hbmgatpRBCi = 0
hbmgatpRBCii = 0
lacPFi = 0.003200000000004
lacRBCi = 0.0014
lacRBCii = 0.126
mgRBCi = 0.003
mgRBCii = 0.27
mgadpRBCi = 0
mgadpRBCii = 0
mgatpRBCi = 0
mgatpRBCii = 0
mgb13pgRBCi = 0
mgb13pgRBCii = 0
mgb23pgRBCi = 0
mgb23pgRBCii = 0
mgf16p2RBCi = 0
mgf16p2RBCii = 0
mgg16p2RBCi = 0
mgg16p2RBCii = 0
mgphosRBCi = 0
mgphosRBCii = 0
nadPFi = 0.08000000000009999
nadRBCi = 3/50000
nadRBCii = 27/5000
nadhPFi = 0.016000000000019998
nadhRBCi = 1.3999999999999998*^-7
nadhRBCii = 0.000012599999999999998
nadpRBCi = 1.25*^-7
nadpRBCii = 0.000011249999999999999
nadphRBCi = 0.00006400000000000001
nadphRBCii = 0.00576
p2gPFi = 0.003200000000004
p2gRBCi = 1/100000
p2gRBCii = 9/10000
p3gPFi = 0.003200000000004
p3gRBCi = 0.00006400000000000001
p3gRBCii = 0.00576
p6gRBCi = 1.3999999999999998*^-7
p6gRBCii = 0.000012599999999999998
p6glRBCi = 1.4*^-10
p6glRBCii = 1.2599999999999999*^-8
pepPFi = 0.003200000000004
pepRBCi = 0.000023
pepRBCii = 0.00207
phosRBCi = 0.001
phosRBCii = 0.09
pyrPFi = 0.003200000000004
pyrRBCi = 3/50000
pyrRBCii = 27/5000
rib5pRBCi = 1/100000
rib5pRBCii = 9/10000
ru5pRBCi = 1/100000
ru5pRBCii = 9/10000
sed7pRBCi = 1/100000
sed7pRBCii = 9/10000
tkRBCi = 3.2999999999999996*^-7
tkRBCii = 0.0000297
tkgRBCi = 0
tkgRBCii = 0
tkgery4pRBCi = 0
tkgery4pRBCii = 0
tkgrib5pRBCi = 0
tkgrib5pRBCii = 0
tkxu5pRBCi = 0
tkxu5pRBCii = 0
xu5pRBCi = 1/1000000
xu5pRBCii = 9/100000
end initial values

begin assignment rules
KiapplacRBCvLDH := 0.715252751049199*(1 + 6.309573444801929*^6/10^phRBC)*KilacRBCvLDH
K1appRBCvBPGSP1 := (3.51188643150958*K1RBCvBPGSP1)/(1 + 1.5848931924611098*^-7*10^phRBC)
KiapppyrRBCvLDH := 0.284747248950801*(1 + 1.5848931924611098*^-7*10^phRBC)*KipyrRBCvLDH
KcatfappRBCvGAPDH := KcatfRBCvGAPDH/(1 + 10^(-10 + phRBC) + 3.16227766016838*^7/10^phRBC)
KaappRBCvMGG16P2 := (CRBCvMGG16P2*KaRBCvMGG16P2*(hRBC*KhfRBCvMGG16P2*KmghfRBCvMGG16P2 + KRBCvMGG16P2))/(1 + hRBC*KhfRBCvMGG16P2 + hRBC^2*Kh2fRBCvMGG16P2*KhfRBCvMGG16P2 + KkfRBCvMGG16P2*kRBC + hRBC*KhfRBCvMGG16P2*KkhfRBCvMGG16P2*kRBC)
K2appRBCvAK := K2RBCvAK*(1 + hRBC*KhampRBCvAK + KkampRBCvAK*kRBC)
K1appRBCvAK := K1RBCvAK*(1 + KkadpRBCvAK*kRBC + (KhadpRBCvAK*pHConversionFactor)/10^phRBC)
Vpf := Piecewise[{{trophToRBC*Vrbci, par > 0}}, 1]
KcatrappRBCvGAPDH := KcatrRBCvGAPDH/(1 + 10^(-10 + phRBC) + 3.16227766016838*^7/10^phRBC)
parMultiUninf := Piecewise[{{1, par < 1}}, 0]
parMulti := Piecewise[{{1, par > 0}}, 0]
hRBC := pHConversionFactor/10^phRBC
glcEXT := ConcGlcEXT*vBld
glyEXT := ConcGlyEXT*vBld
phosEXT := ConcPhosEXT*vBld
pyrEXT := ConcPyrEXT*vBld
lacEXT := ConcLacEXT*vBld
KmapppyrRBCvLDH := 0.284747248950801*(1 + 1.5848931924611098*^-7*10^phRBC)*KmpyrRBCvLDH
KaappRBCvMGADP := (CRBCvMGADP*KaRBCvMGADP*(hRBC*KhadpRBCvMGADP*KmghadpRBCvMGADP + KRBCvMGADP))/(1 + hRBC*KhadpRBCvMGADP + KkadpRBCvMGADP*kRBC)
KcatfappRBCvHK := KcatfRBCvHK/(1 + 1.0471285480509*^7/10^phRBC + 2.8183829312644505*^-10*10^phRBC)
K13appRBCvBPGSP7 := (1.75857757502918*K13RBCvBPGSP7)/(1 + 4.78630092322638*^28/10^(4*phRBC))
K6appRBCvBPGSP4 := (1.75857757502918*K6RBCvBPGSP4)/(1 + 4.78630092322638*^28/10^(4*phRBC))
K3appRBCvBPGSP2 := (1.75857757502918*K3RBCvBPGSP2)/(1 + 4.78630092322638*^28/10^(4*phRBC))
K4appRBCvBPGSP3 := (1.75857757502918*K4RBCvBPGSP3)/(1 + 4.78630092322638*^28/10^(4*phRBC))
co2RBCi := ConcCo2RBC*Vrbci
KaappRBCvHBATP := HbpHRBC*KaRBCvHBATP
KoRBCvPHOSTRANSPORT := KiRBCvPHOSTRANSPORT/KeqRBCvPHOSTRANSPORT
KoRBCvLACTRANSPORT := KiRBCvLACTRANSPORT/KeqRBCvLACTRANSPORT
KaappRBCvHBMGATP := HbpHRBC*KaRBCvHBMGATP
co2RBC := ConcCo2RBC*Vrbcu
KaappRBCvHBB13PG := HbpHRBC*KaRBCvHBB13PG
hPF := pHConversionFactor/10^phPF
KiappnadhRBCvGAPDH := 6.30957344480193*^-8*10^phRBC*KinadhRBCvGAPDH
Vrbci := Piecewise[{{par*Vrbc, par > 0}}, 1]
KmapplacRBCvLDH := 0.715252751049199*(1 + 6.309573444801929*^6/10^phRBC)*KmlacRBCvLDH
KmappnadhRBCvGAPDH := 6.30957344480193*^-8*10^phRBC*KmnadhRBCvGAPDH
KaappRBCvHBADP := HbpHRBC*KaRBCvHBADP
Vrbcu := Piecewise[{{Vrbc - parMulti*Vrbci, par < 1}}, 1]
Kiappb13pgRBCvGAPDH := Kib13pgRBCvGAPDH/(1 + 10^(-10 + phRBC) + 3.16227766016838*^7/10^phRBC)
KaappRBCvHBBPG := HbpHRBC*KaRBCvHBBPG
KoRBCvPYRTRANSPORT := KiRBCvPYRTRANSPORT/RtvRBC
KeqRBCvPHOSTRANSPORT := (1 + 10^(-6.75 + phRBC))/(10^(-6.75 + phRBC)/RtvRBC^2 + RtvRBC^(-1))
KcatrappRBCvHK := KcatrRBCvHK/(1 + 1.0471285480509*^7/10^phRBC + 2.8183829312644505*^-10*10^phRBC)
KiappgapRBCvGAPDH := KigapRBCvGAPDH/(1 + 10^(-10 + phRBC) + 3.16227766016838*^7/10^phRBC)
KaappRBCvMGATP := (CRBCvMGATP*KaRBCvMGATP*(hRBC*KhatpRBCvMGATP*KmghatpRBCvMGATP + KRBCvMGATP))/(1 + KkatpRBCvMGATP*kRBC + KhatpRBCvMGATP/(10^phRBC*pHConversionFactor))
KaappRBCvMGB23PG := (CRBCvMGB23PG*KaRBCvMGB23PG*(hRBC*KhbpgRBCvMGB23PG*KmghbpgRBCvMGB23PG + KRBCvMGB23PG))/(1 + hRBC*KhbpgRBCvMGB23PG + hRBC^2*Kh2bpgRBCvMGB23PG*KhbpgRBCvMGB23PG + KkbpgRBCvMGB23PG*kRBC + hRBC*KhbpgRBCvMGB23PG*KkhbpgRBCvMGB23PG*kRBC)
KaappRBCvMGPHOS := (KaRBCvMGPHOS*(1 + KkphosRBCvMGPHOS*kRBC + 6.30957344480193*^-8*KhphosRBCvMGPHOS*pHConversionFactor))/(1 + hRBC*KhphosRBCvMGPHOS + KkphosRBCvMGPHOS*kRBC)
LRBCvPK := (1.5848931924611098*^-7*pHConversionFactor*(1 + atpRBC[t]/(KtatpRBCvPK*Vrbcu))^4)/(hRBC*(1 + f16p2RBC[t]/(Krf16p2RBCvPK*Vrbcu) + g16p2RBC[t]/(Krg16p2RBCvPK*Vrbcu))^4*(1 + pepRBC[t]/(KrpepRBCvPK*Vrbcu) + pyrRBC[t]/(KrpyrRBCvPK*Vrbcu))^4)
HbpHRBC := (1 + (2.51188643150958*^14*KahbRBC^2)/pHConversionFactor^2 + (3.16978638492223*^7*KahbRBC)/pHConversionFactor)/(1 + (2*KahbRBC)/hRBC + KahbRBC^2/hRBC^2)
KeqRBCvLACTRANSPORT := (1 + 10^(-3.73 + phRBC))/(1 + 10^(-3.73 + phRBC)/RtvRBC)
KaappRBCvMGF16P2 := (CRBCvMGF16P2*KaRBCvMGF16P2*(hRBC*KhfRBCvMGF16P2*KmghfRBCvMGF16P2 + KRBCvMGF16P2))/(1 + hRBC*KhfRBCvMGF16P2 + hRBC^2*Kh2fRBCvMGF16P2*KhfRBCvMGF16P2 + KkfRBCvMGF16P2*kRBC + hRBC*KhfRBCvMGF16P2*KkhfRBCvMGF16P2*kRBC)
LRBCivPK := (1.5848931924611098*^-7*pHConversionFactor*(1 + atpRBCi[t]/(KtatpRBCvPK*Vrbci))^4)/(hRBC*(1 + f16p2RBCi[t]/(Krf16p2RBCvPK*Vrbci) + g16p2RBCi[t]/(Krg16p2RBCvPK*Vrbci))^4*(1 + pepRBCi[t]/(KrpepRBCvPK*Vrbci) + pyrRBCi[t]/(KrpyrRBCvPK*Vrbci))^4)
KaappRBCvMGB13PG := (CRBCvMGB13PG*KaRBCvMGB13PG*(hRBC*KhbpgRBCvMGB13PG*KmghbpgRBCvMGB13PG + KRBCvMGB13PG))/(1 + hRBC*KhbpgRBCvMGB13PG + hRBC^2*Kh2bpgRBCvMGB13PG*KhbpgRBCvMGB13PG + KkbpgRBCvMGB13PG*kRBC + hRBC*KhbpgRBCvMGB13PG*KkhbpgRBCvMGB13PG*kRBC)
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