(* Generated by JWS Online *)

(* This is an experimental feature of JWS Online. Please report any mistakes.*)

(* Note that the following notable SBML entities or features are not supported in notebook outputyet: *)
    (* Events *)
    (* Constraints *)
    (* Units and UnitDefinitions *)
    (* AlgebraicRules *)
    (* conversionFactors *)

variables = {
adpPF[t], adpRBCi[t], ampRBCi[t], atpPF[t], atpRBCi[t], b13pgPF[t], b13pgRBCi[t], b23pgRBCi[t], bpgspRBCi[t], bpgspb13pgRBCi[t], bpgspb23pgRBCi[t], bpgsppRBCi[t], bpgsppp2gRBCi[t], bpgsppp3gRBCi[t], bpgsppphosRBCi[t], dhapPF[t], dhapRBCi[t], ery4pRBCi[t], f16bpPF[t], f16p2RBCi[t], f6pPF[t], f6pRBCi[t], g16p2RBCi[t], g3pPF[t], g6pPF[t], g6pRBCi[t], gapPF[t], gapRBCi[t], glcPF[t], glcRBCi[t], gshRBCi[t], gssgRBCi[t], hbRBCi[t], hbadpRBCi[t], hbatpRBCi[t], hbb13pgRBCi[t], hbb23pgRBCi[t], hbmgatpRBCi[t], lacPF[t], lacRBCi[t], mgRBCi[t], mgadpRBCi[t], mgatpRBCi[t], mgb13pgRBCi[t], mgb23pgRBCi[t], mgf16p2RBCi[t], mgg16p2RBCi[t], mgphosRBCi[t], nadPF[t], nadRBCi[t], nadhPF[t], nadhRBCi[t], nadpRBCi[t], nadphRBCi[t], p2gPF[t], p2gRBCi[t], p3gPF[t], p3gRBCi[t], p6gRBCi[t], p6glRBCi[t], pepPF[t], pepRBCi[t], phosRBCi[t], pyrPF[t], pyrRBCi[t], rib5pRBCi[t], ru5pRBCi[t], sed7pRBCi[t], tkRBCi[t], tkgRBCi[t], tkgery4pRBCi[t], tkgrib5pRBCi[t], tkxu5pRBCi[t], xu5pRBCi[t]
};

initialValues = {
adpPF[0] == 0.032, adpRBCi[0] == 0.0279, ampRBCi[0] == 27/10000, atpPF[0] == 0.064, atpRBCi[0] == 0.189, b13pgPF[0] == 0.0032, b13pgRBCi[0] == 63/1000000, b23pgRBCi[0] == 0.603, bpgspRBCi[0] == 0.00034199999999999996, bpgspb13pgRBCi[0] == 0, bpgspb23pgRBCi[0] == 0, bpgsppRBCi[0] == 0, bpgsppp2gRBCi[0] == 0, bpgsppp3gRBCi[0] == 0, bpgsppphosRBCi[0] == 0, dhapPF[0] == 0.0032, dhapRBCi[0] == 0.0017100000000000001, ery4pRBCi[0] == 9/10000, f16bpPF[0] == 0.0032, f16p2RBCi[0] == 0.00024300000000000002, f6pPF[0] == 0.0032, f6pRBCi[0] == 0.00117, g16p2RBCi[0] == 0.01098, g3pPF[0] == 0.0032, g6pPF[0] == 0.0032, g6pRBCi[0] == 9/2500, gapPF[0] == 0.0032, gapRBCi[0] == 0.000513, glcPF[0] == 0.0032, glcRBCi[0] == 0.09, gshRBCi[0] == 0.28800000000000003, gssgRBCi[0] == 81/10000000, hbRBCi[0] == 0.63, hbadpRBCi[0] == 0, hbatpRBCi[0] == 0, hbb13pgRBCi[0] == 0, hbb23pgRBCi[0] == 0, hbmgatpRBCi[0] == 0, lacPF[0] == 0.0032, lacRBCi[0] == 0.126, mgRBCi[0] == 0.27, mgadpRBCi[0] == 0, mgatpRBCi[0] == 0, mgb13pgRBCi[0] == 0, mgb23pgRBCi[0] == 0, mgf16p2RBCi[0] == 0, mgg16p2RBCi[0] == 0, mgphosRBCi[0] == 0, nadPF[0] == 0.08, nadRBCi[0] == 27/5000, nadhPF[0] == 0.016, nadhRBCi[0] == 0.000012599999999999998, nadpRBCi[0] == 0.000011249999999999999, nadphRBCi[0] == 0.00576, p2gPF[0] == 0.0032, p2gRBCi[0] == 9/10000, p3gPF[0] == 0.0032, p3gRBCi[0] == 0.00576, p6gRBCi[0] == 0.000012599999999999998, p6glRBCi[0] == 1.2599999999999999*^-8, pepPF[0] == 0.0032, pepRBCi[0] == 0.00207, phosRBCi[0] == 0.09, pyrPF[0] == 0.0032, pyrRBCi[0] == 27/5000, rib5pRBCi[0] == 9/10000, ru5pRBCi[0] == 9/10000, sed7pRBCi[0] == 9/10000, tkRBCi[0] == 0.0000297, tkgRBCi[0] == 0, tkgery4pRBCi[0] == 0, tkgrib5pRBCi[0] == 0, tkxu5pRBCi[0] == 0, xu5pRBCi[0] == 9/100000
};

rates = {
vPFvALD, vPFvATPASE, vPFvENO, vPFvG3PDH, vPFvGAPDH, vPFvGLCtr, vPFvGLYtr, vPFvHK, vPFvLACtr, vPFvLDH, vPFvPFK, vPFvPGI, vPFvPGK, vPFvPGM, vPFvPK, vPFvPYRtr, vPFvTPI, vRBCivAK, vRBCivALD, vRBCivATPASE, vRBCivBPGSP1, vRBCivBPGSP2, vRBCivBPGSP3, vRBCivBPGSP4, vRBCivBPGSP5, vRBCivBPGSP6, vRBCivBPGSP7, vRBCivBPGSP8, vRBCivBPGSP9, vRBCivENO, vRBCivG6PDH, vRBCivGAPDH, vRBCivGLCTRANSPORT, vRBCivGSSGR, vRBCivHBADP, vRBCivHBATP, vRBCivHBB13PG, vRBCivHBB23PG, vRBCivHBMGATP, vRBCivHK, vRBCivLACTRANSPORT, vRBCivLDH, vRBCivLDHP, vRBCivMGADP, vRBCivMGATP, vRBCivMGB13PG, vRBCivMGB23PG, vRBCivMGF16P2, vRBCivMGG16P2, vRBCivMGPHOS, vRBCivOX, vRBCivOXNADH, vRBCivP6GDH, vRBCivPFK, vRBCivPGI, vRBCivPGK, vRBCivPGLHYDROLYSIS, vRBCivPGM, vRBCivPHOSTRANSPORT, vRBCivPK, vRBCivPYRTRANSPORT, vRBCivR5PI, vRBCivRu5PE, vRBCivTA, vRBCivTIM, vRBCivTK1, vRBCivTK2, vRBCivTK3, vRBCivTK4, vRBCivTK5, vRBCivTK6
};

rateEquations = {
vPFvALD -> (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 -> (VPFvATPASE*atpPF[t]^5)/(KmPFvATPASE^5*Vpf^4*(1 + atpPF[t]^5/(KmPFvATPASE^5*Vpf^5))), vPFvENO -> (VPFvENO*p2gPF[t]*(1 - pepPF[t]/(KeqPFvENO*p2gPF[t])))/(Kp2gPFvENO*(1 + p2gPF[t]/(Kp2gPFvENO*Vpf) + pepPF[t]/(KpepPFvENO*Vpf))), vPFvG3PDH -> (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 -> (-((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 -> (-((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*Vpf*(-(glyEXT/vBld) + g3pPF[t]/Vpf), vPFvHK -> (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 -> 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 -> (-((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 -> (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 -> (VPFvPGI*(1 - f6pPF[t]/(KeqPFvPGI*g6pPF[t]))*g6pPF[t])/(Kg6pPFvPGI*(1 + f6pPF[t]/(Kf6pPFvPGI*Vpf) + g6pPF[t]/(Kg6pPFvPGI*Vpf))), vPFvPGK -> (KadpPFvPGK*Kb13pgPFvPGK*KeqPFvPGK*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 -> (VPFvPGM*(1 - p2gPF[t]/(KeqPFvPGM*p3gPF[t]))*p3gPF[t])/(Kp3gPFvPGM*(1 + p2gPF[t]/(Kp2gPFvPGM*Vpf) + p3gPF[t]/(Kp3gPFvPGM*Vpf))), vPFvPK -> (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 -> 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 -> (VPFvTPI*dhapPF[t]*(1 - gapPF[t]/(KeqPFvTPI*dhapPF[t])))/(KdhapPFvTPI*(1 + dhapPF[t]/(KdhapPFvTPI*Vpf) + gapPF[t]/(KgapPFvTPI*Vpf) + pepPF[t]/(KipepPFvTPI*Vpf))), vRBCivAK -> TsRBC*Vrbci*((K1appRBCvAK*adpRBCi[t]*mgadpRBCi[t])/Vrbci^2 - (K2appRBCvAK*ampRBCi[t]*mgatpRBCi[t])/Vrbci^2), vRBCivALD -> (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*TsRBC*mgatpRBCi[t], vRBCivBPGSP1 -> TsRBC*Vrbci*(-((K2RBCvBPGSP1*bpgspb13pgRBCi[t])/Vrbci) + (K1appRBCvBPGSP1*b13pgRBCi[t]*bpgspRBCi[t])/Vrbci^2), vRBCivBPGSP2 -> K3appRBCvBPGSP2*TsRBC*bpgspb13pgRBCi[t], vRBCivBPGSP3 -> TsRBC*Vrbci*(-((K5RBCvBPGSP3*bpgsppp3gRBCi[t])/Vrbci) + (K4appRBCvBPGSP3*bpgsppRBCi[t]*p3gRBCi[t])/Vrbci^2), vRBCivBPGSP4 -> TsRBC*Vrbci*(-((K7RBCvBPGSP4*bpgsppp2gRBCi[t])/Vrbci) + (K6appRBCvBPGSP4*bpgsppRBCi[t]*p2gRBCi[t])/Vrbci^2), vRBCivBPGSP5 -> TsRBC*Vrbci*(-((K9RBCvBPGSP5*bpgspb23pgRBCi[t])/Vrbci) + (K8RBCvBPGSP5*bpgsppp3gRBCi[t])/Vrbci), vRBCivBPGSP6 -> TsRBC*Vrbci*(-((K11RBCvBPGSP6*bpgspb23pgRBCi[t])/Vrbci) + (K10RBCvBPGSP6*bpgsppp2gRBCi[t])/Vrbci), vRBCivBPGSP7 -> TsRBC*Vrbci*((K12RBCvBPGSP7*bpgspb23pgRBCi[t])/Vrbci - (K13appRBCvBPGSP7*b23pgRBCi[t]*bpgspRBCi[t])/Vrbci^2), vRBCivBPGSP8 -> TsRBC*Vrbci*(-((K15RBCvBPGSP8*bpgsppphosRBCi[t])/Vrbci) + (K14RBCvBPGSP8*bpgsppRBCi[t]*phosRBCi[t])/Vrbci^2), vRBCivBPGSP9 -> K16RBCvBPGSP9*TsRBC*bpgsppphosRBCi[t], vRBCivENO -> (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 -> (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 -> (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 -> 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 -> (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 -> TsRBC*Vrbci*(-((KdRBCvHBADP*hbadpRBCi[t])/Vrbci) + (KaappRBCvHBADP*adpRBCi[t]*hbRBCi[t])/Vrbci^2), vRBCivHBATP -> TsRBC*Vrbci*(-((KdRBCvHBATP*hbatpRBCi[t])/Vrbci) + (KaappRBCvHBATP*atpRBCi[t]*hbRBCi[t])/Vrbci^2), vRBCivHBB13PG -> TsRBC*Vrbci*(-((KdRBCvHBB13PG*hbb13pgRBCi[t])/Vrbci) + (KaappRBCvHBB13PG*b13pgRBCi[t]*hbRBCi[t])/Vrbci^2), vRBCivHBB23PG -> TsRBC*Vrbci*(-((KdRBCvHBBPG*hbb23pgRBCi[t])/Vrbci) + (KaappRBCvHBBPG*b23pgRBCi[t]*hbRBCi[t])/Vrbci^2), vRBCivHBMGATP -> TsRBC*Vrbci*(-((KdRBCvHBMGATP*hbmgatpRBCi[t])/Vrbci) + (KaappRBCvHBMGATP*hbRBCi[t]*mgatpRBCi[t])/Vrbci^2), vRBCivHK -> (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 -> TsRBC*vBld*(-((KiRBCvLACTRANSPORT*lacEXT)/vBld) + (KoRBCvLACTRANSPORT*lacRBCi[t])/Vrbci), vRBCivLDH -> (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 -> (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 -> TsRBC*Vrbci*(-((KdRBCvMGADP*mgadpRBCi[t])/Vrbci) + (KaappRBCvMGADP*adpRBCi[t]*mgRBCi[t])/Vrbci^2), vRBCivMGATP -> TsRBC*Vrbci*(-((KdRBCvMGATP*mgatpRBCi[t])/Vrbci) + (KaappRBCvMGATP*atpRBCi[t]*mgRBCi[t])/Vrbci^2), vRBCivMGB13PG -> TsRBC*Vrbci*(-((KdRBCvMGB13PG*mgb13pgRBCi[t])/Vrbci) + (KaappRBCvMGB13PG*b13pgRBCi[t]*mgRBCi[t])/Vrbci^2), vRBCivMGB23PG -> TsRBC*Vrbci*(-((KdRBCvMGB23PG*mgb23pgRBCi[t])/Vrbci) + (KaappRBCvMGB23PG*b23pgRBCi[t]*mgRBCi[t])/Vrbci^2), vRBCivMGF16P2 -> TsRBC*Vrbci*(-((KdRBCvMGF16P2*mgf16p2RBCi[t])/Vrbci) + (KaappRBCvMGF16P2*f16p2RBCi[t]*mgRBCi[t])/Vrbci^2), vRBCivMGG16P2 -> TsRBC*Vrbci*(-((KdRBCvMGG16P2*mgg16p2RBCi[t])/Vrbci) + (KaappRBCvMGG16P2*g16p2RBCi[t]*mgRBCi[t])/Vrbci^2), vRBCivMGPHOS -> TsRBC*Vrbci*(-((KdRBCvMGPHOS*mgphosRBCi[t])/Vrbci) + (KaappRBCvMGPHOS*mgRBCi[t]*phosRBCi[t])/Vrbci^2), vRBCivOX -> KRBCvOX*TsRBC*gshRBCi[t], vRBCivOXNADH -> KRBCvOXNADH*TsRBC*nadhRBCi[t], vRBCivP6GDH -> (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 -> (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*TsRBC*Vrbci*(-((KcatrRBCvPGI*f6pRBCi[t])/(Kmf6pRBCvPGI*Vrbci)) + (KcatfRBCvPGI*g6pRBCi[t])/(Kmg6pRBCvPGI*Vrbci)))/(1 + f6pRBCi[t]/(Kmf6pRBCvPGI*Vrbci) + g6pRBCi[t]/(Kmg6pRBCvPGI*Vrbci)), vRBCivPGK -> (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 -> TsRBC*Vrbci*((KhydrolRBCvPGLHYDROLYSIS*p6glRBCi[t])/Vrbci + (Klactonase3RBCvPGLHYDROLYSIS*VRBCvPGLHYDROLYSIS*p6glRBCi[t])/(Vrbci*((Klactonase2RBCvPGLHYDROLYSIS + Klactonase3RBCvPGLHYDROLYSIS)/Klactonase1RBCvPGLHYDROLYSIS + p6glRBCi[t]/Vrbci))), vRBCivPGM -> (TsRBC*Vrbci*VRBCvPGM*(-((KcatrRBCvPGM*p2gRBCi[t])/(Kmp2gRBCvPGM*Vrbci)) + (KcatfRBCvPGM*p3gRBCi[t])/(Kmp3gRBCvPGM*Vrbci)))/(1 + p2gRBCi[t]/(Kmp2gRBCvPGM*Vrbci) + p3gRBCi[t]/(Kmp3gRBCvPGM*Vrbci)), vRBCivPHOSTRANSPORT -> TsRBC*vBld*(-((KiRBCvPHOSTRANSPORT*phosEXT)/vBld) + (KoRBCvPHOSTRANSPORT*phosRBCi[t])/Vrbci), vRBCivPK -> (TsRBC*Vrbci*VRBCvPK*((KcatfRBCvPK*mgadpRBCi[t]*pepRBCi[t])/(KrmgadpRBCvPK*KrpepRBCvPK*Vrbci^2) - (KcatrRBCvPK*mgatpRBCi[t]*pyrRBCi[t])/(KrmgatpRBCvPK*KrpyrRBCvPK*Vrbci^2)))/((1 + LRBCvPK)*(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 -> TsRBC*vBld*(-((KiRBCvPYRTRANSPORT*pyrEXT)/vBld) + (KoRBCvPYRTRANSPORT*pyrRBCi[t])/Vrbci), vRBCivR5PI -> (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 -> (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 -> (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 -> (TsRBC*Vrbci*VRBCvTIM*(-((KcatrRBCvTIM*dhapRBCi[t])/(KmdhapRBCvTIM*Vrbci)) + (KcatfRBCvTIM*gapRBCi[t])/(KmgapRBCvTIM*Vrbci)))/(1 + dhapRBCi[t]/(KmdhapRBCvTIM*Vrbci) + gapRBCi[t]/(KmgapRBCvTIM*Vrbci)), vRBCivTK1 -> TsRBC*Vrbci*(-((K2RBCvTK1*tkxu5pRBCi[t])/Vrbci) + (K1RBCvTK1*tkRBCi[t]*xu5pRBCi[t])/Vrbci^2), vRBCivTK2 -> TsRBC*Vrbci*(-((K4RBCvTK2*gapRBCi[t]*tkgRBCi[t])/Vrbci^2) + (K3RBCvTK2*tkxu5pRBCi[t])/Vrbci), vRBCivTK3 -> TsRBC*Vrbci*((K5RBCvTK3*rib5pRBCi[t]*tkgRBCi[t])/Vrbci^2 - (K6RBCvTK3*tkgrib5pRBCi[t])/Vrbci), vRBCivTK4 -> TsRBC*Vrbci*((K7RBCvTK4*tkgrib5pRBCi[t])/Vrbci - (K8RBCvTK4*sed7pRBCi[t]*tkRBCi[t])/Vrbci^2), vRBCivTK5 -> TsRBC*Vrbci*(-((K10RBCvTK5*tkgery4pRBCi[t])/Vrbci) + (K9RBCvTK5*ery4pRBCi[t]*tkgRBCi[t])/Vrbci^2), vRBCivTK6 -> TsRBC*Vrbci*((K11RBCvTK6*tkgery4pRBCi[t])/Vrbci - (K12RBCvTK6*f6pRBCi[t]*tkRBCi[t])/Vrbci^2)
};

parameters = {
CRBCvMGADP -> 0.00062, CRBCvMGATP -> 8.4*^-05, CRBCvMGB13PG -> 0.0032, CRBCvMGB23PG -> 0.0032, CRBCvMGF16P2 -> 0.0083, CRBCvMGG16P2 -> 0.0083, CRBCvMGPHOS -> 40800.0, ConcADPRBC -> 0.00031, ConcAMPRBC -> 3*^-05, ConcATPRBC -> 0.0021, ConcAdpPF -> 0.001, ConcAtpPF -> 0.002, ConcB13PGRBC -> 7*^-07, ConcB13pgPF -> 0.0001, ConcB23PGRBC -> 0.0067, ConcBPGSPB13PGRBC -> 0.0, ConcBPGSPB23PGRBC -> 0.0, ConcBPGSPPP2GRBC -> 0.0, ConcBPGSPPP3GRBC -> 0.0, ConcBPGSPPRBC -> 0.0, ConcBPGSPRBC -> 3.8*^-06, ConcBpgsppphosRBC -> 0.0, ConcCo2RBC -> 0.0012, ConcDhapPF -> 0.0001, ConcDhapRBC -> 1.9*^-05, ConcEry4PRBC -> 1*^-05, ConcF16P2RBC -> 2.7*^-06, ConcF16bpPF -> 0.0001, ConcF6PRBC -> 1.3*^-05, ConcF6pPF -> 0.0001, ConcG16P2RBC -> 0.000122, ConcG6PRBC -> 4*^-05, ConcG6pPF -> 0.0001, ConcGSHRBC -> 0.0032, ConcGapPF -> 0.0001, ConcGapRBC -> 5.7*^-06, ConcGlcEXT -> 0.005, ConcGlcPF -> 0.0001, ConcGlcRBC -> 0.001, ConcGly3pPF -> 0.0001, ConcGlyEXT -> 0.0, ConcGssgRBC -> 9*^-08, ConcHbADPRBC -> 0.0, ConcHbATPRBC -> 0.0, ConcHbB13PGRBC -> 0.0, ConcHbB23PGRBC -> 0.0, ConcHbRBC -> 0.007, ConcHbmgatpRBC -> 0.0, ConcLacEXT -> 0.0, 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.25*^-07, ConcNADRBC -> 6*^-05, ConcNadPF -> 0.0025, ConcNadhPF -> 0.0005, ConcNadhRBC -> 1.4*^-07, ConcNadphRBC -> 6.4*^-05, ConcP2GARBC -> 1*^-05, ConcP2gPF -> 0.0001, ConcP3GARBC -> 6.4*^-05, ConcP3gPF -> 0.0001, ConcP6GLRBC -> 1.4*^-10, ConcP6GRBC -> 1.4*^-07, ConcPEPRBC -> 2.3*^-05, ConcPepPF -> 0.0001, ConcPhosEXT -> 0.00564, ConcPhosRBC -> 0.001, ConcPyrEXT -> 0.0, ConcPyrPF -> 0.0001, ConcPyrRBC -> 6*^-05, ConcRib5PRBC -> 1*^-05, ConcRu5PRBC -> 1*^-05, ConcSed7PRBC -> 1*^-05, ConcTKGEry4PRBC -> 0.0, ConcTKGRBC -> 0.0, ConcTKGRib5PRBC -> 0.0, ConcTKRBC -> 3.3*^-07, ConcTkxu5pRBC -> 0.0, ConcXu5PRBC -> 1*^-06, ERBCvHK -> 2.5*^-08, ERBCvPGI -> 2.18*^-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.4*^-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.91251*^-08, KadpPFvHK -> 0.000846735, KadpPFvPFK -> 0.00074176, KadpPFvPGK -> 0.00015, KadpPFvPK -> 0.000317, KahbRBC -> 2.51189*^-07, KatpPFvHK -> 0.00069647, KatpPFvPFK -> 0.0007862, KatpPFvPGK -> 0.00077, Kb13pgPFvGAPDH -> 2.359*^-05, Kb13pgPFvPGK -> 1.34*^-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 -> 3*^-05, Kdig6pRBCvHK -> 1*^-05, KdigshRBCvHK -> 0.003, KeqPFvALD -> 9*^-05, KeqPFvENO -> 4.6, KeqPFvG3PDH -> 32600.0, KeqPFvHK -> 1310.0, KeqPFvPGI -> 0.33, KeqPFvPGK -> 3200.0, KeqPFvPGM -> 0.19, KeqPFvTPI -> 0.04545, Kf16bpPFvALD -> 6.84*^-05, Kf16bpPFvPFK -> 0.003626, Kf6pPFvPFK -> 0.00109454, Kf6pPFvPGI -> 9.6651*^-05, KfRBCvLDHP -> 0.00346, Kg3pPFvG3PDH -> 0.00398, Kg6pPFvHK -> 4.3*^-05, Kg6pPFvPGI -> 0.00100774, KgapPFvALD -> 4.6*^-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.55544*^-21, Kib13pgRBCvPGK -> 1.6*^-06, Kib23pgRBCvALD -> 0.0015, Kidb13pgRBCvGAPDH -> 1*^-06, KidgapRBCvGAPDH -> 3.1*^-05, KidhapRBCvALD -> 1.1*^-05, KidpyrRBCvLDH -> 0.000101, Kif16p2RBCvALD -> 1.98*^-05, Kig6pRBCvHK -> 4.7*^-05, KigapRBCvGAPDH -> 4.76523*^-19, KiglcRBCvHK -> 4.7*^-05, KilacPFvPYRtr -> 0.000358, KilacRBCvLDH -> 0.00733, KimgRBCvENO -> 4.6*^-05, KimgadpRBCvHK -> 0.001, KimgadpRBCvPGK -> 8*^-05, KimgatpRBCvHK -> 0.001, KimgatpRBCvPGK -> 0.000186, KinadRBCvGAPDH -> 4.5*^-05, KinadRBCvLDH -> 0.000503, KinadhRBCvGAPDH -> 1*^-05, KinadhRBCvLDH -> 2.45*^-06, Kip2gRBCvENO -> 0.00014, Kip3gRBCvPGK -> 0.000205, KipepPFvPK -> 0.00292, KipepPFvTPI -> 1.59*^-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.71*^-07, Kmb13pgRBCvPGK -> 2*^-06, KmdhapRBCvALD -> 3.5*^-05, KmdhapRBCvTIM -> 0.0001624, Kmf16p2RBCvALD -> 7.1*^-06, Kmf16p2RBCvPFK -> 0.0005, Kmf6pRBCvPFK -> 7.5*^-05, Kmf6pRBCvPGI -> 7.1*^-05, Kmg6pRBCvPGI -> 0.000181, KmgapRBCvALD -> 0.000189, KmgapRBCvGAPDH -> 9.5*^-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.6*^-05, KmmgadpRBCvHK -> 0.001, KmmgadpRBCvPFK -> 0.00054, KmmgadpRBCvPGK -> 0.0001, KmmgatpRBCvHK -> 0.001, KmmgatpRBCvPFK -> 6.8*^-05, KmmgatpRBCvPGK -> 0.001, KmnadRBCvGAPDH -> 4.5*^-05, KmnadRBCvLDH -> 0.000107, KmnadhRBCvGAPDH -> 3.3*^-06, KmnadhRBCvLDH -> 8.44*^-06, Kmp2gRBCvENO -> 0.00014, Kmp2gRBCvPGM -> 2.56*^-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 -> 9*^-05, KnadhPFvGAPDH -> 2.87*^-05, KnadhPFvLDH -> 4.6*^-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.43*^-07, KrampRBCvPFK -> 0.0003, Krf16p2RBCvPK -> 5*^-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.7*^-07, VRBCvENO -> 2.2*^-07, VRBCvG6PDH -> 9.3*^-08, VRBCvGAPDH -> 7.66*^-06, VRBCvGLCTRANSPORT -> 0.047, VRBCvGSSGR -> 1.25*^-07, VRBCvLDH -> 3.43*^-06, VRBCvP6GDH -> 2.1*^-06, VRBCvPFK -> 1.1*^-07, VRBCvPGK -> 2.74*^-06, VRBCvPGLHYDROLYSIS -> 1.4*^-05, VRBCvPGM -> 4.1*^-07, VRBCvPK -> 8.7*^-08, VRBCvR5PI -> 1.42*^-05, VRBCvRu5PE -> 4.22*^-06, VRBCvTA -> 6.9*^-07, VRBCvTIM -> 1.14*^-06, VfPFvG3PDH -> 0.012, VfPFvLDH -> 0.542612419668, VrPFvGAPDH -> 0.243199143455, VrPFvLDH -> 0.290792291203, VrPFvPGK -> 0.353533190557, alphaPFvGLCtr -> 0.91, bpgxRBC -> 3.8*^-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.86746*^-08, ki2 -> 2.08*^-05, nRBCvPFK -> 5.0, nadpxRBC -> 6.4125*^-05, pHConversionFactor -> 1.0, phPF -> 7.2, phRBC -> 7.2, pt -> 0.07, tkxRBC -> 3.3*^-07, xatpRBC -> -0.00456, xg16p2RBC -> 0.000122, xphosRBC -> 0.000438, Vpf -> 32.0, Vrbci -> 90.0, vBld -> 4050.0
};

assignments = {
LRBCvPK -> (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), KiapppyrRBCvLDH -> 0.284747248950801*(1 + 1.5848931924611098*^-7*10^phRBC)*KipyrRBCvLDH, KiapplacRBCvLDH -> 0.715252751049199*(1 + 6.309573444801929*^6/10^phRBC)*KilacRBCvLDH, KcatfappRBCvGAPDH -> KcatfRBCvGAPDH/(1 + 10^(-10 + phRBC) + 3.16227766016838*^7/10^phRBC), KaappRBCvMGF16P2 -> (CRBCvMGF16P2*KaRBCvMGF16P2*(hRBC*KhfRBCvMGF16P2*KmghfRBCvMGF16P2 + KRBCvMGF16P2))/(1 + hRBC*KhfRBCvMGF16P2 + hRBC^2*Kh2fRBCvMGF16P2*KhfRBCvMGF16P2 + KkfRBCvMGF16P2*kRBC + hRBC*KhfRBCvMGF16P2*KkhfRBCvMGF16P2*kRBC), KaappRBCvMGG16P2 -> (CRBCvMGG16P2*KaRBCvMGG16P2*(hRBC*KhfRBCvMGG16P2*KmghfRBCvMGG16P2 + KRBCvMGG16P2))/(1 + hRBC*KhfRBCvMGG16P2 + hRBC^2*Kh2fRBCvMGG16P2*KhfRBCvMGG16P2 + KkfRBCvMGG16P2*kRBC + hRBC*KhfRBCvMGG16P2*KkhfRBCvMGG16P2*kRBC), KeqRBCvPHOSTRANSPORT -> (1 + 10^(-6.75 + phRBC))/(10^(-6.75 + phRBC)/RtvRBC^2 + RtvRBC^(-1)), K6appRBCvBPGSP4 -> (1.75857757502918*K6RBCvBPGSP4)/(1 + 4.78630092322638*^28/10^(4*phRBC)), KiappgapRBCvGAPDH -> KigapRBCvGAPDH/(1 + 10^(-10 + phRBC) + 3.16227766016838*^7/10^phRBC), glyEXT -> ConcGlyEXT*vBld, KaappRBCvHBB13PG -> HbpHRBC*KaRBCvHBB13PG, hPF -> pHConversionFactor/10^phPF, KaappRBCvHBATP -> HbpHRBC*KaRBCvHBATP, KaappRBCvHBMGATP -> HbpHRBC*KaRBCvHBMGATP, KiappnadhRBCvGAPDH -> 6.30957344480193*^-8*10^phRBC*KinadhRBCvGAPDH, trophToRBC -> Vpf/Vrbci, lacEXT -> ConcLacEXT*vBld, KaappRBCvHBADP -> HbpHRBC*KaRBCvHBADP, KaappRBCvHBBPG -> HbpHRBC*KaRBCvHBBPG, KcatrappRBCvHK -> KcatrRBCvHK/(1 + 1.0471285480509*^7/10^phRBC + 2.8183829312644505*^-10*10^phRBC), K4appRBCvBPGSP3 -> (1.75857757502918*K4RBCvBPGSP3)/(1 + 4.78630092322638*^28/10^(4*phRBC)), KaappRBCvMGADP -> (CRBCvMGADP*KaRBCvMGADP*(hRBC*KhadpRBCvMGADP*KmghadpRBCvMGADP + KRBCvMGADP))/(1 + hRBC*KhadpRBCvMGADP + KkadpRBCvMGADP*kRBC), K2appRBCvAK -> K2RBCvAK*(1 + hRBC*KhampRBCvAK + KkampRBCvAK*kRBC), KmapplacRBCvLDH -> 0.715252751049199*(1 + 6.309573444801929*^6/10^phRBC)*KmlacRBCvLDH, Kiappb13pgRBCvGAPDH -> Kib13pgRBCvGAPDH/(1 + 10^(-10 + phRBC) + 3.16227766016838*^7/10^phRBC), KoRBCvPHOSTRANSPORT -> KiRBCvPHOSTRANSPORT/KeqRBCvPHOSTRANSPORT, co2RBCi -> ConcCo2RBC*Vrbci, phosEXT -> ConcPhosEXT*vBld, KoRBCvPYRTRANSPORT -> KiRBCvPYRTRANSPORT/RtvRBC, K13appRBCvBPGSP7 -> (1.75857757502918*K13RBCvBPGSP7)/(1 + 4.78630092322638*^28/10^(4*phRBC)), KoRBCvLACTRANSPORT -> KiRBCvLACTRANSPORT/KeqRBCvLACTRANSPORT, KmapppyrRBCvLDH -> 0.284747248950801*(1 + 1.5848931924611098*^-7*10^phRBC)*KmpyrRBCvLDH, hRBC -> pHConversionFactor/10^phRBC, glcEXT -> ConcGlcEXT*vBld, pyrEXT -> ConcPyrEXT*vBld, KmappnadhRBCvGAPDH -> 6.30957344480193*^-8*10^phRBC*KmnadhRBCvGAPDH, KeqRBCvLACTRANSPORT -> (1 + 10^(-3.73 + phRBC))/(1 + 10^(-3.73 + phRBC)/RtvRBC), K1appRBCvBPGSP1 -> (3.51188643150958*K1RBCvBPGSP1)/(1 + 1.5848931924611098*^-7*10^phRBC), KcatfappRBCvHK -> KcatfRBCvHK/(1 + 1.0471285480509*^7/10^phRBC + 2.8183829312644505*^-10*10^phRBC), KcatrappRBCvGAPDH -> KcatrRBCvGAPDH/(1 + 10^(-10 + phRBC) + 3.16227766016838*^7/10^phRBC), 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), K3appRBCvBPGSP2 -> (1.75857757502918*K3RBCvBPGSP2)/(1 + 4.78630092322638*^28/10^(4*phRBC)), KaappRBCvMGATP -> (CRBCvMGATP*KaRBCvMGATP*(hRBC*KhatpRBCvMGATP*KmghatpRBCvMGATP + KRBCvMGATP))/(1 + KkatpRBCvMGATP*kRBC + KhatpRBCvMGATP/(10^phRBC*pHConversionFactor)), K1appRBCvAK -> K1RBCvAK*(1 + KkadpRBCvAK*kRBC + (KhadpRBCvAK*pHConversionFactor)/10^phRBC), HbpHRBC -> (1 + (2.51188643150958*^14*KahbRBC^2)/pHConversionFactor^2 + (3.16978638492223*^7*KahbRBC)/pHConversionFactor)/(1 + (2*KahbRBC)/hRBC + KahbRBC^2/hRBC^2), KaappRBCvMGB13PG -> (CRBCvMGB13PG*KaRBCvMGB13PG*(hRBC*KhbpgRBCvMGB13PG*KmghbpgRBCvMGB13PG + KRBCvMGB13PG))/(1 + hRBC*KhbpgRBCvMGB13PG + hRBC^2*Kh2bpgRBCvMGB13PG*KhbpgRBCvMGB13PG + KkbpgRBCvMGB13PG*kRBC + hRBC*KhbpgRBCvMGB13PG*KkhbpgRBCvMGB13PG*kRBC)
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

units = {
{"time" -> "min", "metabolite" -> "fmol", "extent" -> "fmol"}
};

(* Time evolution *)
odes = {
adpPF'[t] == 1.0*vPFvATPASE +1.0*vPFvHK +1.0*vPFvPFK -1.0*vPFvPGK -1.0*vPFvPK, adpRBCi'[t] ==  -1.0*vRBCivAK -1.0*vRBCivHBADP -1.0*vRBCivMGADP, ampRBCi'[t] == 1.0*vRBCivAK , atpPF'[t] == 1.0*vPFvPGK +1.0*vPFvPK -1.0*vPFvATPASE -1.0*vPFvHK -1.0*vPFvPFK, atpRBCi'[t] ==  -1.0*vRBCivHBATP -1.0*vRBCivMGATP, b13pgPF'[t] == 1.0*vPFvGAPDH -1.0*vPFvPGK, b13pgRBCi'[t] == 1.0*vRBCivGAPDH -1.0*vRBCivBPGSP1 -1.0*vRBCivHBB13PG -1.0*vRBCivMGB13PG -1.0*vRBCivPGK, b23pgRBCi'[t] == 1.0*vRBCivBPGSP7 -1.0*vRBCivHBB23PG -1.0*vRBCivMGB23PG, bpgspRBCi'[t] == 1.0*vRBCivBPGSP7 +1.0*vRBCivBPGSP9 -1.0*vRBCivBPGSP1, bpgspb13pgRBCi'[t] == 1.0*vRBCivBPGSP1 -1.0*vRBCivBPGSP2, bpgspb23pgRBCi'[t] == 1.0*vRBCivBPGSP5 +1.0*vRBCivBPGSP6 -1.0*vRBCivBPGSP7, bpgsppRBCi'[t] == 1.0*vRBCivBPGSP2 -1.0*vRBCivBPGSP3 -1.0*vRBCivBPGSP4 -1.0*vRBCivBPGSP8, bpgsppp2gRBCi'[t] == 1.0*vRBCivBPGSP4 -1.0*vRBCivBPGSP6, bpgsppp3gRBCi'[t] == 1.0*vRBCivBPGSP3 -1.0*vRBCivBPGSP5, bpgsppphosRBCi'[t] == 1.0*vRBCivBPGSP8 -1.0*vRBCivBPGSP9, dhapPF'[t] == 1.0*vPFvALD -1.0*vPFvG3PDH -1.0*vPFvTPI, dhapRBCi'[t] == 1.0*vRBCivALD +1.0*vRBCivTIM , ery4pRBCi'[t] == 1.0*vRBCivTA -1.0*vRBCivTK5, f16bpPF'[t] == 1.0*vPFvPFK -1.0*vPFvALD, f16p2RBCi'[t] == 1.0*vRBCivPFK -1.0*vRBCivALD -1.0*vRBCivMGF16P2, f6pPF'[t] == 1.0*vPFvPGI -1.0*vPFvPFK, f6pRBCi'[t] == 1.0*vRBCivPGI +1.0*vRBCivTA +1.0*vRBCivTK6 -1.0*vRBCivPFK, g16p2RBCi'[t] ==  -1.0*vRBCivMGG16P2, g3pPF'[t] == 1.0*vPFvG3PDH -1.0*vPFvGLYtr, g6pPF'[t] == 1.0*vPFvHK -1.0*vPFvPGI, g6pRBCi'[t] == 1.0*vRBCivHK -1.0*vRBCivG6PDH -1.0*vRBCivPGI, gapPF'[t] == 1.0*vPFvALD +1.0*vPFvTPI -1.0*vPFvGAPDH, gapRBCi'[t] == 1.0*vRBCivALD +1.0*vRBCivTK2 -1.0*vRBCivGAPDH -1.0*vRBCivTA -1.0*vRBCivTIM, glcPF'[t] == 1.0*vPFvGLCtr -1.0*vPFvHK, glcRBCi'[t] == 1.0*vRBCivGLCTRANSPORT -1.0*vPFvGLCtr -1.0*vRBCivHK, gshRBCi'[t] == 2.0*vRBCivGSSGR -2.0*vRBCivOX, gssgRBCi'[t] == 1.0*vRBCivOX -1.0*vRBCivGSSGR, hbRBCi'[t] ==  -1.0*vRBCivHBADP -1.0*vRBCivHBATP -1.0*vRBCivHBB13PG -1.0*vRBCivHBB23PG -1.0*vRBCivHBMGATP, hbadpRBCi'[t] == 1.0*vRBCivHBADP , hbatpRBCi'[t] == 1.0*vRBCivHBATP , hbb13pgRBCi'[t] == 1.0*vRBCivHBB13PG , hbb23pgRBCi'[t] == 1.0*vRBCivHBB23PG , hbmgatpRBCi'[t] == 1.0*vRBCivHBMGATP , lacPF'[t] == 1.0*vPFvLDH -1.0*vPFvLACtr, lacRBCi'[t] == 1.0*vPFvLACtr +1.0*vRBCivLDH +1.0*vRBCivLDHP -1.0*vRBCivLACTRANSPORT, mgRBCi'[t] ==  -1.0*vRBCivMGADP -1.0*vRBCivMGATP -1.0*vRBCivMGB13PG -1.0*vRBCivMGB23PG -1.0*vRBCivMGF16P2 -1.0*vRBCivMGG16P2 -1.0*vRBCivMGPHOS, mgadpRBCi'[t] == 1.0*vRBCivATPASE +1.0*vRBCivHK +1.0*vRBCivMGADP +1.0*vRBCivPFK -1.0*vRBCivAK -1.0*vRBCivPGK -1.0*vRBCivPK, mgatpRBCi'[t] == 1.0*vRBCivAK +1.0*vRBCivMGATP +1.0*vRBCivPGK +1.0*vRBCivPK -1.0*vRBCivATPASE -1.0*vRBCivHBMGATP -1.0*vRBCivHK -1.0*vRBCivPFK, mgb13pgRBCi'[t] == 1.0*vRBCivMGB13PG , mgb23pgRBCi'[t] == 1.0*vRBCivMGB23PG , mgf16p2RBCi'[t] == 1.0*vRBCivMGF16P2 , mgg16p2RBCi'[t] == 1.0*vRBCivMGG16P2 , mgphosRBCi'[t] == 1.0*vRBCivMGPHOS , nadPF'[t] == 1.0*vPFvG3PDH +1.0*vPFvLDH -1.0*vPFvGAPDH, nadRBCi'[t] == 1.0*vRBCivLDH +1.0*vRBCivOXNADH -1.0*vRBCivGAPDH, nadhPF'[t] == 1.0*vPFvGAPDH -1.0*vPFvG3PDH -1.0*vPFvLDH, nadhRBCi'[t] == 1.0*vRBCivGAPDH -1.0*vRBCivLDH -1.0*vRBCivOXNADH, nadpRBCi'[t] == 1.0*vRBCivGSSGR +1.0*vRBCivLDHP -1.0*vRBCivG6PDH -1.0*vRBCivP6GDH, nadphRBCi'[t] == 1.0*vRBCivG6PDH +1.0*vRBCivP6GDH -1.0*vRBCivGSSGR -1.0*vRBCivLDHP, p2gPF'[t] == 1.0*vPFvPGM -1.0*vPFvENO, p2gRBCi'[t] == 1.0*vRBCivPGM -1.0*vRBCivBPGSP4 -1.0*vRBCivENO, p3gPF'[t] == 1.0*vPFvPGK -1.0*vPFvPGM, p3gRBCi'[t] == 1.0*vRBCivBPGSP2 +1.0*vRBCivPGK -1.0*vRBCivBPGSP3 -1.0*vRBCivPGM, p6gRBCi'[t] == 1.0*vRBCivPGLHYDROLYSIS -1.0*vRBCivP6GDH, p6glRBCi'[t] == 1.0*vRBCivG6PDH -1.0*vRBCivPGLHYDROLYSIS, pepPF'[t] == 1.0*vPFvENO -1.0*vPFvPK, pepRBCi'[t] == 1.0*vRBCivENO -1.0*vRBCivPK, phosRBCi'[t] == 1.0*vRBCivATPASE +2.0*vRBCivBPGSP9 -1.0*vRBCivBPGSP8 -1.0*vRBCivGAPDH -1.0*vRBCivMGPHOS -1.0*vRBCivPHOSTRANSPORT, pyrPF'[t] == 1.0*vPFvPK -1.0*vPFvLDH -1.0*vPFvPYRtr, pyrRBCi'[t] == 1.0*vPFvPYRtr +1.0*vRBCivPK -1.0*vRBCivLDH -1.0*vRBCivLDHP -1.0*vRBCivPYRTRANSPORT, rib5pRBCi'[t] == 1.0*vRBCivR5PI -1.0*vRBCivTK3, ru5pRBCi'[t] == 1.0*vRBCivP6GDH -1.0*vRBCivR5PI -1.0*vRBCivRu5PE, sed7pRBCi'[t] == 1.0*vRBCivTK4 -1.0*vRBCivTA, tkRBCi'[t] == 1.0*vRBCivTK4 +1.0*vRBCivTK6 -1.0*vRBCivTK1, tkgRBCi'[t] == 1.0*vRBCivTK2 -1.0*vRBCivTK3 -1.0*vRBCivTK5, tkgery4pRBCi'[t] == 1.0*vRBCivTK5 -1.0*vRBCivTK6, tkgrib5pRBCi'[t] == 1.0*vRBCivTK3 -1.0*vRBCivTK4, tkxu5pRBCi'[t] == 1.0*vRBCivTK1 -1.0*vRBCivTK2, xu5pRBCi'[t] == 1.0*vRBCivRu5PE -1.0*vRBCivTK1
};
timeCourse = NDSolve[Join[odes, initialValues]//.rateEquations//.assignments//.parameters, variables, {t, 0, 100}];

(* Steady-state solution initialized with result of time evolution *)
findRootEquations = odes /.D[_[t],t]->0;
findRootVariables = Partition[Flatten[{#, #/.timeCourse/.t->100} &/@variables],2];
steadyStateVariables = FindRoot[findRootEquations//.rateEquations//.assignments//.parameters, findRootVariables, MaxIterations->100]
fluxes = #//.assignments//.parameters/.steadyStateVariables&/@rateEquations

(* Plot the time evolution of the variables *)
plotTable=Table[Plot[variables[[i]]/.parameters/.timeCourse,{t,0,100},PlotLegends->variables[[i]],PlotRange->Full],{i,Length[variables]}]