(* 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 = {
ADPf[t], AMPf[t], ATPf[t], DHAP[t], E4P[t], Fru16P2[t], Fru6P[t], GSH[t], GSSG[t], Glc6P[t], GlcA6P[t], Glcin[t], GraP[t], Gri13P2[t], Gri23P2f[t], Gri2P[t], Gri3P[t], Lac[t], MgADP[t], MgAMP[t], MgATP[t], MgGri23P2[t], Mgf[t], NAD[t], NADH[t], NADPHf[t], NADPf[t], P1NADP[t], P1NADPH[t], P1f[t], P2NADP[t], P2NADPH[t], P2f[t], PEP[t], Phi[t], Pyr[t], Rib5P[t], Rul5P[t], Sed7P[t], Xul5P[t]
};

initialValues = {
ADPf[0] == 0.25, AMPf[0] == 0.0, ATPf[0] == 0.25, DHAP[0] == 0.1492, E4P[0] == 0.0063, Fru16P2[0] == 0.0097, Fru6P[0] == 0.0153, GSH[0] == 3.1136, GSSG[0] == 0.0004, Glc6P[0] == 0.0394, GlcA6P[0] == 0.025, Glcin[0] == 4.5663, GraP[0] == 0.0061, Gri13P2[0] == 0.0005, Gri23P2f[0] == 2.0601, Gri2P[0] == 0.0084, Gri3P[0] == 0.0658, Lac[0] == 1.6803, MgADP[0] == 0.1, MgAMP[0] == 0.0, MgATP[0] == 1.4, MgGri23P2[0] == 0.5, Mgf[0] == 0.8, NAD[0] == 0.0653, NADH[0] == 0.0002, NADPHf[0] == 0.004, NADPf[0] == 0.0, P1NADP[0] == 0.0, P1NADPH[0] == 0.024, P1f[0] == 0.0, P2NADP[0] == 0.0, P2NADPH[0] == 0.024, P2f[0] == 0.0, PEP[0] == 0.0109, Phi[0] == 0.9992, Pyr[0] == 0.084, Rib5P[0] == 0.014, Rul5P[0] == 0.0047, Sed7P[0] == 0.0154, Xul5P[0] == 0.0127
};

rates = {
v\[LetterSpace]1, v\[LetterSpace]10, v\[LetterSpace]11, v\[LetterSpace]12, v\[LetterSpace]13, v\[LetterSpace]14, v\[LetterSpace]15, v\[LetterSpace]16, v\[LetterSpace]17, v\[LetterSpace]18, v\[LetterSpace]19, v\[LetterSpace]2, v\[LetterSpace]20, v\[LetterSpace]21, v\[LetterSpace]22, v\[LetterSpace]23, v\[LetterSpace]24, v\[LetterSpace]25, v\[LetterSpace]26, v\[LetterSpace]27, v\[LetterSpace]28, v\[LetterSpace]29, v\[LetterSpace]3, v\[LetterSpace]30, v\[LetterSpace]31, v\[LetterSpace]32, v\[LetterSpace]33, v\[LetterSpace]34, v\[LetterSpace]35, v\[LetterSpace]36, v\[LetterSpace]37, v\[LetterSpace]38, v\[LetterSpace]4, v\[LetterSpace]5, v\[LetterSpace]6, v\[LetterSpace]7, v\[LetterSpace]8, v\[LetterSpace]9
};

rateEquations = {
v\[LetterSpace]1 -> (Vmaxv0*(Glcout - Glcin[t]/Keqv0))/(KMoutv0*(1 + Glcout/KMoutv0 + Glcin[t]/KMinv0 + (alfav0*Glcout*Glcin[t])/(KMinv0*KMoutv0))), v\[LetterSpace]10 -> (Vmaxv9*(Gri23P2f[t] - Gri3P[t]/Keqv9 + MgGri23P2[t]))/(K23P2Gv9 + Gri23P2f[t] + MgGri23P2[t]), v\[LetterSpace]11 -> (Vmaxv10*(-(Gri2P[t]/Keqv10) + Gri3P[t]))/(K3PGv10*(1 + Gri2P[t]/K2PGv10) + Gri3P[t]), v\[LetterSpace]12 -> (Vmaxv11*(Gri2P[t] - PEP[t]/Keqv11))/(Gri2P[t] + K2PGv11*(1 + PEP[t]/KPEPv11)), v\[LetterSpace]13 -> (Vmaxv12*(MgADP[t]*PEP[t] - (MgATP[t]*Pyr[t])/Keqv12))/((KMgADPv12 + MgADP[t])*(KPEPv12 + PEP[t])*(1 + (L0v12*(1 + (ATPf[t] + MgATP[t])/KATPv12)^4)/((1 + Fru16P2[t]/KFru16P2v12)^4*(1 + PEP[t]/KPEPv12)^4))), v\[LetterSpace]14 -> Vmaxv13*(-((Lac[t]*NAD[t])/Keqv13) + NADH[t]*Pyr[t]), v\[LetterSpace]15 -> kLDHv14*(-((Lac[t]*NADPf[t])/Keqv14) + NADPHf[t]*Pyr[t]), v\[LetterSpace]16 -> kATPasev15*MgATP[t], v\[LetterSpace]17 -> (Vmaxv16*(-((ADPf[t]*MgADP[t])/Keqv16) + AMPf[t]*MgATP[t]))/(KAMPv16*KATPv16*((ADPf[t]*MgADP[t])/KADPv16^2 + (ADPf[t] + MgADP[t])/KADPv16 + (1 + AMPf[t]/KAMPv16)*(1 + MgATP[t]/KATPv16))), v\[LetterSpace]18 -> (Vmaxv17*(Glc6P[t]*NADPf[t] - (GlcA6P[t]*NADPHf[t])/Keqv17))/(KG6Pv17*KNADPv17*(1 + (ATPf[t] + MgATP[t])/KATPv17 + (Gri23P2f[t] + MgGri23P2[t])/KPGA23v17 + ((1 + Glc6P[t]/KG6Pv17)*NADPf[t])/KNADPv17 + NADPHf[t]/KNADPHv17)), v\[LetterSpace]19 -> (Vmaxv18*(GlcA6P[t]*NADPf[t] - (NADPHf[t]*Rul5P[t])/Keqv18))/(K6PG1v18*KNADPv18*((ATPf[t] + MgATP[t])/KATPv18 + (1 + GlcA6P[t]/K6PG1v18 + (Gri23P2f[t] + MgGri23P2[t])/KPGA23v18)*(1 + NADPf[t]/KNADPv18) + ((1 + GlcA6P[t]/K6PG2v18)*NADPHf[t])/KNADPHv18)), v\[LetterSpace]2 -> (Inhibv1*Vmax1v1*Glcin[t]*(-((Glc6P[t]*MgADP[t])/Keqv1) + MgATP[t] + (Vmax2v1*MgATP[t]*Mgf[t])/(KMgATPMgv1*Vmax1v1)))/(KMgATPv1*(KMGlcv1 + Glcin[t])*(1 + Mgf[t]/KMgv1 + (MgATP[t]*(1 + Mgf[t]/KMgATPMgv1))/KMgATPv1 + (1.55 + Glc6P[t]/KGlc6Pv1)*(1 + Mgf[t]/KMgv1) + (Gri23P2f[t] + MgGri23P2[t])/K23P2Gv1 + (Mgf[t]*(Gri23P2f[t] + MgGri23P2[t]))/(KMg23P2Gv1*KMgv1))), v\[LetterSpace]20 -> (Vmaxv19*(-((GSH[t]^2*NADPf[t])/(Keqv19*KGSHv19^2*KNADPv19)) + (GSSG[t]*NADPHf[t])/(KGSSGv19*KNADPHv19)))/(1 + ((1 + (GSH[t]*(1 + GSH[t]/KGSHv19))/KGSHv19)*NADPf[t])/KNADPv19 + ((1 + GSSG[t]/KGSSGv19)*NADPHf[t])/KNADPHv19), v\[LetterSpace]21 -> Kv20*GSH[t], v\[LetterSpace]22 -> (Vmaxv21*(Rul5P[t] - Xul5P[t]/Keqv21))/(Rul5P[t] + KRu5Pv21*(1 + Xul5P[t]/KX5Pv21)), v\[LetterSpace]23 -> (Vmaxv22*(-(Rib5P[t]/Keqv22) + Rul5P[t]))/(KRu5Pv22*(1 + Rib5P[t]/KR5Pv22) + Rul5P[t]), v\[LetterSpace]24 -> (Vmaxv23*(-((GraP[t]*Sed7P[t])/Keqv23) + Rib5P[t]*Xul5P[t]))/(K4v23*Sed7P[t] + GraP[t]*(K3v23 + K5v23*Sed7P[t]) + Rib5P[t]*(K2v23 + K6v23*Sed7P[t]) + K7v23*GraP[t]*Xul5P[t] + (K1v23 + Rib5P[t])*Xul5P[t]), v\[LetterSpace]25 -> (Vmaxv24*(-((E4P[t]*Fru6P[t])/Keqv24) + GraP[t]*Sed7P[t]))/(K4v24*Fru6P[t] + E4P[t]*(K3v24 + K5v24*Fru6P[t]) + (K2v24 + K6v24*Fru6P[t])*GraP[t] + K7v24*E4P[t]*Sed7P[t] + (K1v24 + GraP[t])*Sed7P[t]), v\[LetterSpace]26 -> (Vmaxv25*(-((PRPP*MgAMP[t])/Keqv25) + MgATP[t]*Rib5P[t]))/((KATPv25 + MgATP[t])*(KR5Pv25 + Rib5P[t])), v\[LetterSpace]27 -> (Vmaxv26*(-((Fru6P[t]*GraP[t])/Keqv26) + E4P[t]*Xul5P[t]))/(K4v26*Fru6P[t] + E4P[t]*(K2v26 + K6v26*Fru6P[t]) + (K3v26 + K5v26*Fru6P[t])*GraP[t] + (K1v26 + E4P[t])*Xul5P[t] + K7v26*GraP[t]*Xul5P[t]), v\[LetterSpace]28 -> Vmaxv27*(Phiex - Phi[t]/Keqv27), v\[LetterSpace]29 -> Vmaxv28*(Lacex - Lac[t]/Keqv28), v\[LetterSpace]3 -> (Vmaxv2*(-(Fru6P[t]/Keqv2) + Glc6P[t]))/(KGlc6Pv2*(1 + Fru6P[t]/KFru6Pv2) + Glc6P[t]), v\[LetterSpace]30 -> Vmaxv29*(Pyrex - Pyr[t]/Keqv29), v\[LetterSpace]31 -> EqMult*(MgATP[t] - (ATPf[t]*Mgf[t])/KdATP), v\[LetterSpace]32 -> EqMult*(MgADP[t] - (ADPf[t]*Mgf[t])/KdADP), v\[LetterSpace]33 -> EqMult*(MgAMP[t] - (AMPf[t]*Mgf[t])/KdAMP), v\[LetterSpace]34 -> EqMult*(-((Gri23P2f[t]*Mgf[t])/Kd23P2G) + MgGri23P2[t]), v\[LetterSpace]35 -> EqMult*(-((NADPf[t]*P1f[t])/Kd1) + P1NADP[t]), v\[LetterSpace]36 -> EqMult*(-((NADPHf[t]*P1f[t])/Kd3) + P1NADPH[t]), v\[LetterSpace]37 -> EqMult*(-((NADPf[t]*P2f[t])/Kd2) + P2NADP[t]), v\[LetterSpace]38 -> EqMult*(-((NADPHf[t]*P2f[t])/Kd4) + P2NADPH[t]), v\[LetterSpace]4 -> (Vmaxv3*(-((Fru16P2[t]*MgADP[t])/Keqv3) + Fru6P[t]*MgATP[t]))/((KFru6Pv3 + Fru6P[t])*(KMgATPv3 + MgATP[t])*(1 + (L0v3*(1 + ATPf[t]/KATPv3)^4*(1 + Mgf[t]/KMgv3)^4)/((1 + Fru6P[t]/KFru6Pv3)^4*(1 + (AMPf[t] + MgAMP[t])/KAMPv3)^4))), v\[LetterSpace]5 -> (Vmaxv4*(Fru16P2[t] - (DHAP[t]*GraP[t])/Keqv4))/(KFru16P2v4*(1 + Fru16P2[t]/KFru16P2v4 + GraP[t]/KiGraPv4 + (Fru16P2[t]*GraP[t])/(KFru16P2v4*KiiGraPv4) + (DHAP[t]*(KGraPv4 + GraP[t]))/(KDHAPv4*KiGraPv4))), v\[LetterSpace]6 -> (Vmaxv5*(DHAP[t] - GraP[t]/Keqv5))/(DHAP[t] + KDHAPv5*(1 + GraP[t]/KGraPv5)), v\[LetterSpace]7 -> (Vmaxv6*(-((Gri13P2[t]*NADH[t])/Keqv6) + GraP[t]*NAD[t]*Phi[t]))/(KGraPv6*KNADv6*KPv6*(-1 + (1 + Gri13P2[t]/K13P2Gv6)*(1 + NADH[t]/KNADHv6) + (1 + GraP[t]/KGraPv6)*(1 + NAD[t]/KNADv6)*(1 + Phi[t]/KPv6))), v\[LetterSpace]8 -> (Vmaxv7*(Gri13P2[t]*MgADP[t] - (Gri3P[t]*MgATP[t])/Keqv7))/(K13P2Gv7*KMgADPv7*(-1 + (1 + Gri13P2[t]/K13P2Gv7)*(1 + MgADP[t]/KMgADPv7) + (1 + Gri3P[t]/K3PGv7)*(1 + MgATP[t]/KMgATPv7))), v\[LetterSpace]9 -> (kDPGMv8*(Gri13P2[t] - (Gri23P2f[t] + MgGri23P2[t])/Keqv8))/(1 + (Gri23P2f[t] + MgGri23P2[t])/K23P2Gv8)
};

parameters = {
Atot -> 2.0, EqMult -> 1000.0, GStotal -> 3.114, Inhibv1 -> 1.0, K13P2Gv6 -> 0.0035, K13P2Gv7 -> 0.002, K1v23 -> 0.4177, K1v24 -> 0.00823, K1v26 -> 0.00184, K23P2Gv1 -> 2.7, K23P2Gv8 -> 0.04, K23P2Gv9 -> 0.2, K2PGv10 -> 1.0, K2PGv11 -> 1.0, K2v23 -> 0.3055, K2v24 -> 0.04765, K2v26 -> 0.3055, K3PGv10 -> 5.0, K3PGv7 -> 1.2, K3v23 -> 12.432, K3v24 -> 0.1733, K3v26 -> 0.0548, K4v23 -> 0.00496, K4v24 -> 0.006095, K4v26 -> 0.0003, K5v23 -> 0.41139, K5v24 -> 0.8683, K5v26 -> 0.0287, K6PG1v18 -> 0.01, K6PG2v18 -> 0.058, K6v23 -> 0.00774, K6v24 -> 0.4653, K6v26 -> 0.122, K7v23 -> 48.8, K7v24 -> 2.524, K7v26 -> 0.215, KADPv16 -> 0.11, KAMPv16 -> 0.08, KAMPv3 -> 0.033, KATPv12 -> 3.39, KATPv16 -> 0.09, KATPv17 -> 0.749, KATPv18 -> 0.154, KATPv25 -> 0.03, KATPv3 -> 0.01, KDHAPv4 -> 0.0364, KDHAPv5 -> 0.838, KFru16P2v12 -> 0.005, KFru16P2v4 -> 0.0071, KFru6Pv2 -> 0.071, KFru6Pv3 -> 0.1, KG6Pv17 -> 0.0667, KGSHv19 -> 20.0, KGSSGv19 -> 0.0652, KGlc6Pv1 -> 0.0045, KGlc6Pv2 -> 0.182, KGraPv4 -> 0.1906, KGraPv5 -> 0.428, KGraPv6 -> 0.005, KMGlcv1 -> 0.1, KMg23P2Gv1 -> 3.44, KMgADPv12 -> 0.474, KMgADPv7 -> 0.35, KMgATPMgv1 -> 1.14, KMgATPv1 -> 1.44, KMgATPv3 -> 0.068, KMgATPv7 -> 0.48, KMgv1 -> 1.03, KMgv3 -> 0.44, KMinv0 -> 6.9, KMoutv0 -> 1.7, KNADHv6 -> 0.0083, KNADPHv17 -> 0.00312, KNADPHv18 -> 0.0045, KNADPHv19 -> 0.00852, KNADPv17 -> 0.00367, KNADPv18 -> 0.018, KNADPv19 -> 0.07, KNADv6 -> 0.05, KPEPv11 -> 1.0, KPEPv12 -> 0.225, KPGA23v17 -> 2.289, KPGA23v18 -> 0.12, KPv6 -> 3.9, KR5Pv22 -> 2.2, KR5Pv25 -> 0.57, KRu5Pv21 -> 0.19, KRu5Pv22 -> 0.78, KX5Pv21 -> 0.5, Kd1 -> 0.0002, Kd2 -> 1*^-05, Kd23P2G -> 1.667, Kd3 -> 1*^-05, Kd4 -> 0.0002, KdADP -> 0.76, KdAMP -> 16.64, KdATP -> 0.072, Keqv0 -> 1.0, Keqv1 -> 3900.0, Keqv10 -> 0.145, Keqv11 -> 1.7, Keqv12 -> 13790.0, Keqv13 -> 9090.0, Keqv14 -> 14181.8, Keqv16 -> 0.25, Keqv17 -> 2000.0, Keqv18 -> 141.7, Keqv19 -> 1.04, Keqv2 -> 0.3925, Keqv21 -> 2.7, Keqv22 -> 3.0, Keqv23 -> 1.05, Keqv24 -> 1.05, Keqv25 -> 100000.0, Keqv26 -> 1.2, Keqv27 -> 1.0, Keqv28 -> 1.0, Keqv29 -> 1.0, Keqv3 -> 100000.0, Keqv4 -> 0.114, Keqv5 -> 0.0407, Keqv6 -> 0.000192, Keqv7 -> 1455.0, Keqv8 -> 100000.0, Keqv9 -> 100000.0, KiGraPv4 -> 0.0572, KiiGraPv4 -> 0.176, Kv20 -> 0.03, L0v12 -> 19.0, L0v3 -> 0.001072, Mgtot -> 2.8, NADPtot -> 0.052, NADtot -> 0.0655, Vmax1v1 -> 15.8, Vmax2v1 -> 33.2, Vmaxv0 -> 33.6, Vmaxv10 -> 2000.0, Vmaxv11 -> 1500.0, Vmaxv12 -> 570.0, Vmaxv13 -> 2800000.0, Vmaxv16 -> 1380.0, Vmaxv17 -> 162.0, Vmaxv18 -> 1575.0, Vmaxv19 -> 90.0, Vmaxv2 -> 935.0, Vmaxv21 -> 4634.0, Vmaxv22 -> 730.0, Vmaxv23 -> 23.5, Vmaxv24 -> 27.2, Vmaxv25 -> 1.1, Vmaxv26 -> 23.5, Vmaxv27 -> 100.0, Vmaxv28 -> 10000.0, Vmaxv29 -> 10000.0, Vmaxv3 -> 239.0, Vmaxv4 -> 98.91000366, Vmaxv5 -> 5456.600098, Vmaxv6 -> 4300.0, Vmaxv7 -> 5000.0, Vmaxv9 -> 0.53, alfav0 -> 0.54, kATPasev15 -> 1.68, kDPGMv8 -> 76000.0, kLDHv14 -> 243.4, protein1 -> 0.024, protein2 -> 0.024, Glcout -> 5.0, Lacex -> 1.68, PRPP -> 1.0, Phiex -> 1.0, Pyrex -> 0.084, default\[LetterSpace]compartment -> 1.0
};

assignments = {

};

events = {

};

speciesAnnotations = {
ADPf[t]->"http://identifiers.org/kegg.compound/C00008", AMPf[t]->"http://identifiers.org/kegg.compound/C00020", ATPf[t]->"http://identifiers.org/kegg.compound/C00002", DHAP[t]->"http://identifiers.org/kegg.compound/C00111", E4P[t]->"http://identifiers.org/kegg.compound/C00279", Fru16P2[t]->"http://identifiers.org/kegg.compound/C00354", Fru6P[t]->"http://identifiers.org/kegg.compound/C00085", GSH[t]->"http://identifiers.org/kegg.compound/C00051", GSSG[t]->"http://identifiers.org/kegg.compound/C00127", Glc6P[t]->"http://identifiers.org/kegg.compound/C00092", Glcin[t]->"http://identifiers.org/kegg.compound/C00092", Glcout[t]->"http://identifiers.org/kegg.compound/C00293", Gri13P2[t]->"http://identifiers.org/kegg.compound/C00236", Gri23P2f[t]->"http://identifiers.org/kegg.compound/C01159", Gri2P[t]->"http://identifiers.org/kegg.compound/C00631", Gri3P[t]->"http://identifiers.org/kegg.compound/C00197", Lac[t]->"http://identifiers.org/kegg.compound/C01432", Lacex[t]->"http://identifiers.org/kegg.compound/C01432", Mgf[t]->"http://identifiers.org/kegg.compound/C00305", NAD[t]->"http://identifiers.org/kegg.compound/C00003", NADH[t]->"http://identifiers.org/kegg.compound/C00004", NADPHf[t]->"http://identifiers.org/kegg.compound/C00005", NADPf[t]->"http://identifiers.org/kegg.compound/C00006", P1NADP[t]->"http://identifiers.org/kegg.compound/C00006", P1NADPH[t]->"http://identifiers.org/kegg.compound/C00005", P2NADP[t]->"http://identifiers.org/kegg.compound/C00006", P2NADPH[t]->"http://identifiers.org/kegg.compound/C00005", PEP[t]->"http://identifiers.org/kegg.compound/C00074", Pyr[t]->"http://identifiers.org/kegg.compound/C00022", Pyrex[t]->"http://identifiers.org/kegg.compound/C00022", Rib5P[t]->"http://identifiers.org/kegg.compound/C00117", Rul5P[t]->"http://identifiers.org/kegg.compound/C00199", Sed7P[t]->"http://identifiers.org/kegg.compound/C05382", Xul5P[t]->"http://identifiers.org/kegg.compound/C03291"
};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
ADPf'[t] == 1.0*v\[LetterSpace]17 +1.0*v\[LetterSpace]32 , AMPf'[t] == 1.0*v\[LetterSpace]33 -1.0*v\[LetterSpace]17, ATPf'[t] == 1.0*v\[LetterSpace]31 , DHAP'[t] == 1.0*v\[LetterSpace]5 -1.0*v\[LetterSpace]6, E4P'[t] == 1.0*v\[LetterSpace]25 -1.0*v\[LetterSpace]27, Fru16P2'[t] == 1.0*v\[LetterSpace]4 -1.0*v\[LetterSpace]5, Fru6P'[t] == 1.0*v\[LetterSpace]27 +1.0*v\[LetterSpace]3 +1.0*v\[LetterSpace]25 -1.0*v\[LetterSpace]4, GSH'[t] == 2.0*v\[LetterSpace]20 -2.0*v\[LetterSpace]21, GSSG'[t] == 1.0*v\[LetterSpace]21 -1.0*v\[LetterSpace]20, Glc6P'[t] == 1.0*v\[LetterSpace]2 -1.0*v\[LetterSpace]18 -1.0*v\[LetterSpace]3, GlcA6P'[t] == 1.0*v\[LetterSpace]18 -1.0*v\[LetterSpace]19, Glcin'[t] == 1.0*v\[LetterSpace]1 -1.0*v\[LetterSpace]2, GraP'[t] == 1.0*v\[LetterSpace]27 +1.0*v\[LetterSpace]6 +1.0*v\[LetterSpace]24 +1.0*v\[LetterSpace]5 -1.0*v\[LetterSpace]25 -1.0*v\[LetterSpace]7, Gri13P2'[t] == 1.0*v\[LetterSpace]7 -1.0*v\[LetterSpace]8 -1.0*v\[LetterSpace]9, Gri23P2f'[t] == 1.0*v\[LetterSpace]34 +1.0*v\[LetterSpace]9 -1.0*v\[LetterSpace]10, Gri2P'[t] == 1.0*v\[LetterSpace]11 -1.0*v\[LetterSpace]12, Gri3P'[t] == 1.0*v\[LetterSpace]10 +1.0*v\[LetterSpace]8 -1.0*v\[LetterSpace]11, Lac'[t] == 1.0*v\[LetterSpace]15 +1.0*v\[LetterSpace]29 +1.0*v\[LetterSpace]14 , MgADP'[t] == 1.0*v\[LetterSpace]17 +1.0*v\[LetterSpace]4 +1.0*v\[LetterSpace]16 +1.0*v\[LetterSpace]2 -1.0*v\[LetterSpace]8 -1.0*v\[LetterSpace]13 -1.0*v\[LetterSpace]32, MgAMP'[t] == 1.0*v\[LetterSpace]26 -1.0*v\[LetterSpace]33, MgATP'[t] == 1.0*v\[LetterSpace]8 +1.0*v\[LetterSpace]13 -1.0*v\[LetterSpace]17 -1.0*v\[LetterSpace]4 -1.0*v\[LetterSpace]16 -1.0*v\[LetterSpace]31 -1.0*v\[LetterSpace]26 -1.0*v\[LetterSpace]2, MgGri23P2'[t] ==  -1.0*v\[LetterSpace]34, Mgf'[t] == 1.0*v\[LetterSpace]34 +1.0*v\[LetterSpace]33 +1.0*v\[LetterSpace]31 +1.0*v\[LetterSpace]32 , NAD'[t] == 1.0*v\[LetterSpace]14 -1.0*v\[LetterSpace]7, NADH'[t] == 1.0*v\[LetterSpace]7 -1.0*v\[LetterSpace]14, NADPHf'[t] == 1.0*v\[LetterSpace]36 +1.0*v\[LetterSpace]38 +1.0*v\[LetterSpace]18 +1.0*v\[LetterSpace]19 -1.0*v\[LetterSpace]15 -1.0*v\[LetterSpace]20, NADPf'[t] == 1.0*v\[LetterSpace]37 +1.0*v\[LetterSpace]15 +1.0*v\[LetterSpace]20 +1.0*v\[LetterSpace]35 -1.0*v\[LetterSpace]18 -1.0*v\[LetterSpace]19, P1NADP'[t] ==  -1.0*v\[LetterSpace]35, P1NADPH'[t] ==  -1.0*v\[LetterSpace]36, P1f'[t] == 1.0*v\[LetterSpace]36 +1.0*v\[LetterSpace]35 , P2NADP'[t] ==  -1.0*v\[LetterSpace]37, P2NADPH'[t] ==  -1.0*v\[LetterSpace]38, P2f'[t] == 1.0*v\[LetterSpace]37 +1.0*v\[LetterSpace]38 , PEP'[t] == 1.0*v\[LetterSpace]12 -1.0*v\[LetterSpace]13, Phi'[t] == 1.0*v\[LetterSpace]10 +1.0*v\[LetterSpace]16 +1.0*v\[LetterSpace]28 -1.0*v\[LetterSpace]7, Pyr'[t] == 1.0*v\[LetterSpace]13 +1.0*v\[LetterSpace]30 -1.0*v\[LetterSpace]15 -1.0*v\[LetterSpace]14, Rib5P'[t] == 1.0*v\[LetterSpace]23 -1.0*v\[LetterSpace]24 -1.0*v\[LetterSpace]26, Rul5P'[t] == 1.0*v\[LetterSpace]19 -1.0*v\[LetterSpace]23 -1.0*v\[LetterSpace]22, Sed7P'[t] == 1.0*v\[LetterSpace]24 -1.0*v\[LetterSpace]25, Xul5P'[t] == 1.0*v\[LetterSpace]22 -1.0*v\[LetterSpace]27 -1.0*v\[LetterSpace]24
};
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]}]