(* 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 = {
species\[LetterSpace]1[t], species\[LetterSpace]10[t], species\[LetterSpace]11[t], species\[LetterSpace]12[t], species\[LetterSpace]13[t], species\[LetterSpace]14[t], species\[LetterSpace]15[t], species\[LetterSpace]16[t], species\[LetterSpace]17[t], species\[LetterSpace]18[t], species\[LetterSpace]19[t], species\[LetterSpace]2[t], species\[LetterSpace]20[t], species\[LetterSpace]21[t], species\[LetterSpace]22[t], species\[LetterSpace]23[t], species\[LetterSpace]24[t], species\[LetterSpace]3[t], species\[LetterSpace]4[t], species\[LetterSpace]5[t], species\[LetterSpace]6[t], species\[LetterSpace]7[t], species\[LetterSpace]8[t], species\[LetterSpace]9[t]
};

initialValues = {
species\[LetterSpace]1[0] == 10.0, species\[LetterSpace]10[0] == 0.0, species\[LetterSpace]11[0] == 0.240501857508, species\[LetterSpace]12[0] == 1.51899628498, species\[LetterSpace]13[0] == 4.24050185751, species\[LetterSpace]14[0] == 0.5, species\[LetterSpace]15[0] == 0.5, species\[LetterSpace]16[0] == 10.0, species\[LetterSpace]17[0] == 8.47911460193, species\[LetterSpace]18[0] == 2.5, species\[LetterSpace]19[0] == 2.0, species\[LetterSpace]2[0] == 1.31652277625, species\[LetterSpace]20[0] == 2.0, species\[LetterSpace]21[0] == 0.5, species\[LetterSpace]22[0] == 10.5208853981, species\[LetterSpace]23[0] == 0.0, species\[LetterSpace]24[0] == 0.0, species\[LetterSpace]3[0] == 0.341738611875, species\[LetterSpace]4[0] == 0.0, species\[LetterSpace]5[0] == 0.0, species\[LetterSpace]6[0] == 2.24173861188, species\[LetterSpace]7[0] == 0.0, species\[LetterSpace]8[0] == 2.23134594788, species\[LetterSpace]9[0] == 2.76865405212
};

rates = {
vAKc, vAKg, vALD, vANTI, vAU, vENO, vGAPDH, vGK, vGPDH, vGPO, vGT, vGlyT, vHK, vPFK, vPGI, vPGK, vPGM, vPGT, vPK, vPT, vTPI
};

rateEquations = {
vAKc -> compartment\[LetterSpace]1*function\[LetterSpace]22[vAKc\[LetterSpace]k, species\[LetterSpace]3[t], species\[LetterSpace]6[t], vAKc\[LetterSpace]keqak, species\[LetterSpace]2[t]], vAKg -> compartment\[LetterSpace]2*function\[LetterSpace]22[vAKg\[LetterSpace]k, species\[LetterSpace]11[t], species\[LetterSpace]13[t], vAKg\[LetterSpace]keqak, species\[LetterSpace]12[t]], vALD -> compartment\[LetterSpace]2*function\[LetterSpace]14[vALD\[LetterSpace]Vmax\[LetterSpace]v5, species\[LetterSpace]16[t], species\[LetterSpace]18[t], species\[LetterSpace]17[t], vALD\[LetterSpace]Keq\[LetterSpace]v5, species\[LetterSpace]11[t], species\[LetterSpace]12[t], species\[LetterSpace]13[t], vALD\[LetterSpace]r\[LetterSpace]v5, vALD\[LetterSpace]KGAP\[LetterSpace]v5, vALD\[LetterSpace]KGAPi\[LetterSpace]v5], vANTI -> vANTI\[LetterSpace]k1*species\[LetterSpace]22[t]*species\[LetterSpace]8[t] - vANTI\[LetterSpace]k2*species\[LetterSpace]17[t]*species\[LetterSpace]9[t], vAU -> compartment\[LetterSpace]1*function\[LetterSpace]20[vAU\[LetterSpace]k, species\[LetterSpace]3[t], species\[LetterSpace]2[t]], vENO -> compartment\[LetterSpace]1*RaENO*function\[LetterSpace]12[species\[LetterSpace]5[t], species\[LetterSpace]4[t], vENO\[LetterSpace]Kms, vENO\[LetterSpace]Kmp, vENO\[LetterSpace]Vf, vENO\[LetterSpace]Vr], vGAPDH -> compartment\[LetterSpace]2*function\[LetterSpace]15[vGAPDH\[LetterSpace]Vmax\[LetterSpace]v7, species\[LetterSpace]18[t], vGAPDH\[LetterSpace]KGAP\[LetterSpace]v7, species\[LetterSpace]19[t], vGAPDH\[LetterSpace]KNAD\[LetterSpace]v7, vGAPDH\[LetterSpace]r\[LetterSpace]v7, species\[LetterSpace]21[t], vGAPDH\[LetterSpace]KBPGA13\[LetterSpace]v7, species\[LetterSpace]20[t], vGAPDH\[LetterSpace]KNADH\[LetterSpace]v7], vGK -> compartment\[LetterSpace]2*function\[LetterSpace]21[vGK\[LetterSpace]Vmax\[LetterSpace]v14, species\[LetterSpace]22[t], vGK\[LetterSpace]KGly3Pg\[LetterSpace]v14, species\[LetterSpace]12[t], vGK\[LetterSpace]KADPg\[LetterSpace]v14, vGK\[LetterSpace]r\[LetterSpace]v14, species\[LetterSpace]24[t], vGK\[LetterSpace]KGlycerol\[LetterSpace]v14, species\[LetterSpace]11[t], vGK\[LetterSpace]KATPg\[LetterSpace]v14], vGPDH -> compartment\[LetterSpace]2*function\[LetterSpace]16[vGPDH\[LetterSpace]Vmax\[LetterSpace]v8, species\[LetterSpace]17[t], vGPDH\[LetterSpace]KDHAPg\[LetterSpace]v8, species\[LetterSpace]20[t], vGPDH\[LetterSpace]KNADH\[LetterSpace]v8, vGPDH\[LetterSpace]r\[LetterSpace]v8, species\[LetterSpace]19[t], vGPDH\[LetterSpace]KNAD\[LetterSpace]v8, species\[LetterSpace]22[t], vGPDH\[LetterSpace]KGly3Pg\[LetterSpace]v8], vGPO -> compartment\[LetterSpace]1*function\[LetterSpace]17[species\[LetterSpace]9[t], vGPO\[LetterSpace]Km, vGPO\[LetterSpace]V], vGT -> function\[LetterSpace]10[vGT\[LetterSpace]Vmax\[LetterSpace]v1, species\[LetterSpace]25, species\[LetterSpace]10[t], vGT\[LetterSpace]KGlc, vGT\[LetterSpace]Alpha\[LetterSpace]v1], vGlyT -> -(species\[LetterSpace]27*vGlyT\[LetterSpace]k2) + vGlyT\[LetterSpace]k1*species\[LetterSpace]24[t], vHK -> compartment\[LetterSpace]2*RaHXK*function\[LetterSpace]11[vHK\[LetterSpace]Vmax\[LetterSpace]v2, species\[LetterSpace]10[t], vHK\[LetterSpace]KGlcInt\[LetterSpace]v2, species\[LetterSpace]11[t], vHK\[LetterSpace]KATPg\[LetterSpace]v2, species\[LetterSpace]12[t], vHK\[LetterSpace]KADPg\[LetterSpace]v2, species\[LetterSpace]14[t], vHK\[LetterSpace]KGlc6P\[LetterSpace]v2], vPFK -> compartment\[LetterSpace]2*RaPFK*function\[LetterSpace]13[vPFK\[LetterSpace]Vmax\[LetterSpace]v4, vPFK\[LetterSpace]Ki1Fru16BP\[LetterSpace]v4, species\[LetterSpace]16[t], species\[LetterSpace]15[t], vPFK\[LetterSpace]KFru6P\[LetterSpace]v4, species\[LetterSpace]11[t], vPFK\[LetterSpace]KATPg\[LetterSpace]v4, vPFK\[LetterSpace]Ki2Fru16BP\[LetterSpace]v4], vPGI -> compartment\[LetterSpace]2*function\[LetterSpace]12[species\[LetterSpace]14[t], species\[LetterSpace]15[t], vPGI\[LetterSpace]Kms, vPGI\[LetterSpace]Kmp, vPGI\[LetterSpace]Vf, vPGI\[LetterSpace]Vr], vPGK -> compartment\[LetterSpace]2*function\[LetterSpace]18[vPGK\[LetterSpace]Vmax\[LetterSpace]v11, species\[LetterSpace]21[t], vPGK\[LetterSpace]KBPGA13\[LetterSpace]v11, species\[LetterSpace]12[t], vPGK\[LetterSpace]KADPg\[LetterSpace]v11, vPGK\[LetterSpace]r\[LetterSpace]v11, species\[LetterSpace]23[t], vPGK\[LetterSpace]KPGA3\[LetterSpace]v11, species\[LetterSpace]11[t], vPGK\[LetterSpace]KATPg\[LetterSpace]v11], vPGM -> compartment\[LetterSpace]1*RaPGAM*function\[LetterSpace]12[species\[LetterSpace]7[t], species\[LetterSpace]5[t], vPGM\[LetterSpace]Kms, vPGM\[LetterSpace]Kmp, vPGM\[LetterSpace]Vf, vPGM\[LetterSpace]Vr], vPGT -> vPGT\[LetterSpace]k1*species\[LetterSpace]23[t] - vPGT\[LetterSpace]k2*species\[LetterSpace]7[t], vPK -> compartment\[LetterSpace]1*RaPYK*function\[LetterSpace]19[vPK\[LetterSpace]Vmax\[LetterSpace]v12, species\[LetterSpace]4[t], vPK\[LetterSpace]PK\[LetterSpace]n, species\[LetterSpace]2[t], vPK\[LetterSpace]KADP\[LetterSpace]v12, species\[LetterSpace]3[t]], vPT -> function\[LetterSpace]17[species\[LetterSpace]1[t], vPT\[LetterSpace]Km, vPT\[LetterSpace]V], vTPI -> compartment\[LetterSpace]2*function\[LetterSpace]12[species\[LetterSpace]17[t], species\[LetterSpace]18[t], vTPI\[LetterSpace]Kms, vTPI\[LetterSpace]Kmp, vTPI\[LetterSpace]Vf, vTPI\[LetterSpace]Vr]
};

parameters = {
RaENO -> 1.0, RaHXK -> 1.0, RaPFK -> 1.0, RaPGAM -> 1.0, RaPYK -> 1.0, species\[LetterSpace]25 -> 5.0, species\[LetterSpace]26 -> 0.0, species\[LetterSpace]27 -> 0.0, vHK\[LetterSpace]KATPg\[LetterSpace]v2 -> 0.116, vGT\[LetterSpace]Vmax\[LetterSpace]v1 -> 108.9, vGT\[LetterSpace]KGlc -> 1.0, vGT\[LetterSpace]Alpha\[LetterSpace]v1 -> 0.75, vHK\[LetterSpace]Vmax\[LetterSpace]v2 -> 1929.0, vHK\[LetterSpace]KGlcInt\[LetterSpace]v2 -> 0.1, vHK\[LetterSpace]KADPg\[LetterSpace]v2 -> 0.126, vHK\[LetterSpace]KGlc6P\[LetterSpace]v2 -> 12.0, vPGI\[LetterSpace]Kms -> 0.4, vPGI\[LetterSpace]Kmp -> 0.12, vPGI\[LetterSpace]Vf -> 1305.0, vPGI\[LetterSpace]Vr -> 1305.0, vPFK\[LetterSpace]Vmax\[LetterSpace]v4 -> 1708.0, vPFK\[LetterSpace]Ki1Fru16BP\[LetterSpace]v4 -> 15.8, vPFK\[LetterSpace]KFru6P\[LetterSpace]v4 -> 0.82, vPFK\[LetterSpace]KATPg\[LetterSpace]v4 -> 0.026, vPFK\[LetterSpace]Ki2Fru16BP\[LetterSpace]v4 -> 10.7, vALD\[LetterSpace]Vmax\[LetterSpace]v5 -> 560.0, vALD\[LetterSpace]Keq\[LetterSpace]v5 -> 0.069, vALD\[LetterSpace]r\[LetterSpace]v5 -> 1.19, vALD\[LetterSpace]KGAP\[LetterSpace]v5 -> 0.067, vALD\[LetterSpace]KGAPi\[LetterSpace]v5 -> 0.098, vTPI\[LetterSpace]Kms -> 1.2, vTPI\[LetterSpace]Kmp -> 0.25, vTPI\[LetterSpace]Vf -> 999.3, vTPI\[LetterSpace]Vr -> 5696.01, vGAPDH\[LetterSpace]Vmax\[LetterSpace]v7 -> 720.9, vGAPDH\[LetterSpace]KGAP\[LetterSpace]v7 -> 0.15, vGAPDH\[LetterSpace]KNAD\[LetterSpace]v7 -> 0.45, vGAPDH\[LetterSpace]r\[LetterSpace]v7 -> 0.67, vGAPDH\[LetterSpace]KBPGA13\[LetterSpace]v7 -> 0.1, vGAPDH\[LetterSpace]KNADH\[LetterSpace]v7 -> 0.02, vGPDH\[LetterSpace]Vmax\[LetterSpace]v8 -> 465.0, vGPDH\[LetterSpace]KDHAPg\[LetterSpace]v8 -> 0.1, vGPDH\[LetterSpace]KNADH\[LetterSpace]v8 -> 0.01, vGPDH\[LetterSpace]r\[LetterSpace]v8 -> 0.28, vGPDH\[LetterSpace]KNAD\[LetterSpace]v8 -> 0.4, vGPDH\[LetterSpace]KGly3Pg\[LetterSpace]v8 -> 2.0, vGPO\[LetterSpace]Km -> 1.7, vGPO\[LetterSpace]V -> 368.0, vPT\[LetterSpace]Km -> 1.96, vPT\[LetterSpace]V -> 200.0, vPGK\[LetterSpace]Vmax\[LetterSpace]v11 -> 2862.0, vPGK\[LetterSpace]KBPGA13\[LetterSpace]v11 -> 0.003, vPGK\[LetterSpace]KADPg\[LetterSpace]v11 -> 0.1, vPGK\[LetterSpace]r\[LetterSpace]v11 -> 0.47, vPGK\[LetterSpace]KPGA3\[LetterSpace]v11 -> 1.62, vPGK\[LetterSpace]KATPg\[LetterSpace]v11 -> 0.29, vPK\[LetterSpace]Vmax\[LetterSpace]v12 -> 1020.0, vPK\[LetterSpace]PK\[LetterSpace]n -> 2.5, vPK\[LetterSpace]KADP\[LetterSpace]v12 -> 0.114, vAU\[LetterSpace]k -> 50.0, vGK\[LetterSpace]Vmax\[LetterSpace]v14 -> 200.0, vGK\[LetterSpace]KGly3Pg\[LetterSpace]v14 -> 3.83, vGK\[LetterSpace]KADPg\[LetterSpace]v14 -> 0.56, vGK\[LetterSpace]r\[LetterSpace]v14 -> 60.86, vGK\[LetterSpace]KGlycerol\[LetterSpace]v14 -> 0.44, vGK\[LetterSpace]KATPg\[LetterSpace]v14 -> 0.24, vPGM\[LetterSpace]Kms -> 0.27, vPGM\[LetterSpace]Kmp -> 0.11, vPGM\[LetterSpace]Vf -> 225.0, vPGM\[LetterSpace]Vr -> 495.0, vENO\[LetterSpace]Kms -> 0.054, vENO\[LetterSpace]Kmp -> 0.24, vENO\[LetterSpace]Vf -> 598.0, vENO\[LetterSpace]Vr -> 394.68, vAKc\[LetterSpace]k -> 1000000.0, vAKc\[LetterSpace]keqak -> 0.442, vAKg\[LetterSpace]k -> 1000000.0, vAKg\[LetterSpace]keqak -> 0.442, vPGT\[LetterSpace]k1 -> 1000000.0, vPGT\[LetterSpace]k2 -> 1000000.0, vANTI\[LetterSpace]k1 -> 1000000.0, vANTI\[LetterSpace]k2 -> 1000000.0, vGlyT\[LetterSpace]k1 -> 1000000.0, vGlyT\[LetterSpace]k2 -> 1000000.0, compartment\[LetterSpace]1 -> 1.0, compartment\[LetterSpace]2 -> 1.0, compartment\[LetterSpace]3 -> 1.0
};

assignments = {
function\[LetterSpace]17[substrate_,Km_,V_] -> (substrate*V)/(Km + substrate), function\[LetterSpace]12[substrate_,product_,Kms_,Kmp_,Vf_,Vr_] -> ((substrate*Vf)/Kms - (product*Vr)/Kmp)/(1 + product/Kmp + substrate/Kms), function\[LetterSpace]20[k_,atpc_,adpc_] -> (atpc*k)/adpc, function\[LetterSpace]18[Vmax\[LetterSpace]v11_,BPGA13_,KBPGA13\[LetterSpace]v11_,ADPg_,KADPg\[LetterSpace]v11_,r\[LetterSpace]v11_,PGA3_,KPGA3\[LetterSpace]v11_,ATPg_,KATPg\[LetterSpace]v11_] -> (((ADPg*BPGA13)/(KADPg\[LetterSpace]v11*KBPGA13\[LetterSpace]v11) - (ATPg*PGA3*r\[LetterSpace]v11)/(KATPg\[LetterSpace]v11*KPGA3\[LetterSpace]v11))*Vmax\[LetterSpace]v11)/((1 + ADPg/KADPg\[LetterSpace]v11 + ATPg/KATPg\[LetterSpace]v11)*(1 + BPGA13/KBPGA13\[LetterSpace]v11 + PGA3/KPGA3\[LetterSpace]v11)), function\[LetterSpace]11[Vmax\[LetterSpace]v2_,GlucoseInt_,KGlcInt\[LetterSpace]v2_,ATPg_,KATPg\[LetterSpace]v2_,ADPg_,KADPg\[LetterSpace]v2_,Glc6P_,KGlc6P\[LetterSpace]v2_] -> (ATPg*GlucoseInt*Vmax\[LetterSpace]v2)/((1 + ADPg/KADPg\[LetterSpace]v2 + ATPg/KATPg\[LetterSpace]v2)*KATPg\[LetterSpace]v2*(1 + Glc6P/KGlc6P\[LetterSpace]v2 + GlucoseInt/KGlcInt\[LetterSpace]v2)*KGlcInt\[LetterSpace]v2), function\[LetterSpace]22[k_,atp_,amp_,keqak_,adp_] -> k*(amp*atp - adp^2*keqak), function\[LetterSpace]10[Vmax\[LetterSpace]v1_,GlucoseExt_,GlucoseInt_,KGlc_,Alpha\[LetterSpace]v1_] -> ((GlucoseExt - GlucoseInt)*Vmax\[LetterSpace]v1)/(GlucoseExt + GlucoseInt + (Alpha\[LetterSpace]v1*GlucoseExt*GlucoseInt)/KGlc + KGlc), function\[LetterSpace]14[Vmax\[LetterSpace]v5_,Fru16BP_,GAP_,DHAPg_,Keq\[LetterSpace]v5_,ATPg_,ADPg_,AMPg_,r\[LetterSpace]v5_,KGAP\[LetterSpace]v5_,KGAPi\[LetterSpace]v5_] -> ((Fru16BP - (DHAPg*GAP)/Keq\[LetterSpace]v5)*Vmax\[LetterSpace]v5)/(0.009*(1 + 0.6622516556291391*ADPg + 0.273972602739726*AMPg + 1.4705882352941175*ATPg) + Fru16BP + (Fru16BP*GAP)/KGAPi\[LetterSpace]v5 + (0.015*(1 + 0.6622516556291391*ADPg + 0.273972602739726*AMPg + 1.4705882352941175*ATPg)*GAP)/(Keq\[LetterSpace]v5*r\[LetterSpace]v5) + (DHAPg*GAP)/(Keq\[LetterSpace]v5*r\[LetterSpace]v5) + (DHAPg*KGAP\[LetterSpace]v5)/(Keq\[LetterSpace]v5*r\[LetterSpace]v5)), function\[LetterSpace]21[Vmax\[LetterSpace]v14_,Gly3Pg_,KGly3Pg\[LetterSpace]v14_,ADPg_,KADPg\[LetterSpace]v14_,r\[LetterSpace]v14_,Glycerol_,KGlycerol\[LetterSpace]v14_,ATPg_,KATPg\[LetterSpace]v14_] -> (((ADPg*Gly3Pg)/(KADPg\[LetterSpace]v14*KGly3Pg\[LetterSpace]v14) - (ATPg*Glycerol*r\[LetterSpace]v14)/(KATPg\[LetterSpace]v14*KGlycerol\[LetterSpace]v14))*Vmax\[LetterSpace]v14)/((1 + ADPg/KADPg\[LetterSpace]v14 + ATPg/KATPg\[LetterSpace]v14)*(1 + Gly3Pg/KGly3Pg\[LetterSpace]v14 + Glycerol/KGlycerol\[LetterSpace]v14)), function\[LetterSpace]13[Vmax\[LetterSpace]v4_,Ki1Fru16BP\[LetterSpace]v4_,Fru16BP_,Fru6P_,KFru6P\[LetterSpace]v4_,ATPg_,KATPg\[LetterSpace]v4_,Ki2Fru16BP\[LetterSpace]v4_] -> (ATPg*Fru6P*Ki1Fru16BP\[LetterSpace]v4*Vmax\[LetterSpace]v4)/((1 + ATPg/KATPg\[LetterSpace]v4)*KATPg\[LetterSpace]v4*KFru6P\[LetterSpace]v4*(Fru16BP + Ki1Fru16BP\[LetterSpace]v4)*(1 + Fru6P/KFru6P\[LetterSpace]v4 + Fru16BP/Ki2Fru16BP\[LetterSpace]v4)), function\[LetterSpace]16[Vmax\[LetterSpace]v8_,DHAPg_,KDHAPg\[LetterSpace]v8_,NADH_,KNADH\[LetterSpace]v8_,r\[LetterSpace]v8_,NAD_,KNAD\[LetterSpace]v8_,Gly3Pg_,KGly3Pg\[LetterSpace]v8_] -> (((DHAPg*NADH)/(KDHAPg\[LetterSpace]v8*KNADH\[LetterSpace]v8) - (Gly3Pg*NAD*r\[LetterSpace]v8)/(KGly3Pg\[LetterSpace]v8*KNAD\[LetterSpace]v8))*Vmax\[LetterSpace]v8)/((1 + DHAPg/KDHAPg\[LetterSpace]v8 + Gly3Pg/KGly3Pg\[LetterSpace]v8)*(1 + NAD/KNAD\[LetterSpace]v8 + NADH/KNADH\[LetterSpace]v8)), function\[LetterSpace]19[Vmax\[LetterSpace]v12_,PEP_,PK\[LetterSpace]n_,ADPc_,KADP\[LetterSpace]v12_,ATPc_] -> (2.941176470588235^PK\[LetterSpace]n*ADPc*(PEP/(1 + 1.5625*ADPc + 1.7543859649122808*ATPc))^PK\[LetterSpace]n*Vmax\[LetterSpace]v12)/((1 + ADPc/KADP\[LetterSpace]v12)*KADP\[LetterSpace]v12*(1 + 2.941176470588235^PK\[LetterSpace]n*(PEP/(1 + 1.5625*ADPc + 1.7543859649122808*ATPc))^PK\[LetterSpace]n)), function\[LetterSpace]15[Vmax\[LetterSpace]v7_,GAP_,KGAP\[LetterSpace]v7_,NAD_,KNAD\[LetterSpace]v7_,r\[LetterSpace]v7_,BPGA13_,KBPGA13\[LetterSpace]v7_,NADH_,KNADH\[LetterSpace]v7_] -> (((GAP*NAD)/(KGAP\[LetterSpace]v7*KNAD\[LetterSpace]v7) - (BPGA13*NADH*r\[LetterSpace]v7)/(KBPGA13\[LetterSpace]v7*KNADH\[LetterSpace]v7))*Vmax\[LetterSpace]v7)/((1 + BPGA13/KBPGA13\[LetterSpace]v7 + GAP/KGAP\[LetterSpace]v7)*(1 + NAD/KNAD\[LetterSpace]v7 + NADH/KNADH\[LetterSpace]v7))
};

events = {

};

speciesAnnotations = {
species\[LetterSpace]10[t]->"http://identifiers.org/chebi/CHEBI:4167", species\[LetterSpace]10[t]->"http://identifiers.org/kegg.compound/C00031", species\[LetterSpace]10[t]->"http://identifiers.org/pubchem.substance/3587", species\[LetterSpace]11[t]->"http://identifiers.org/chebi/CHEBI:15422", species\[LetterSpace]11[t]->"http://identifiers.org/kegg.compound/C00002", species\[LetterSpace]11[t]->"http://identifiers.org/pubchem.substance/3304", species\[LetterSpace]11[t]->"http://identifiers.org/cas/56-65-5", species\[LetterSpace]12[t]->"http://identifiers.org/chebi/CHEBI:16761", species\[LetterSpace]12[t]->"http://identifiers.org/kegg.compound/C00008", species\[LetterSpace]12[t]->"http://identifiers.org/pubchem.substance/3310", species\[LetterSpace]12[t]->"http://identifiers.org/cas/20398-34-9", species\[LetterSpace]13[t]->"http://identifiers.org/chebi/CHEBI:16027", species\[LetterSpace]13[t]->"http://identifiers.org/kegg.compound/C00020", species\[LetterSpace]13[t]->"http://identifiers.org/pubchem.substance/3322", species\[LetterSpace]13[t]->"http://identifiers.org/cas/61-19-8", species\[LetterSpace]14[t]->"http://identifiers.org/chebi/CHEBI:4170", species\[LetterSpace]14[t]->"http://identifiers.org/kegg.compound/C00668", species\[LetterSpace]14[t]->"http://identifiers.org/pubchem.substance/3937", species\[LetterSpace]15[t]->"http://identifiers.org/chebi/CHEBI:16084", species\[LetterSpace]15[t]->"http://identifiers.org/kegg.compound/C05345", species\[LetterSpace]15[t]->"http://identifiers.org/pubchem.substance/7723", species\[LetterSpace]16[t]->"http://identifiers.org/chebi/CHEBI:28013", species\[LetterSpace]16[t]->"http://identifiers.org/kegg.compound/C05378", species\[LetterSpace]16[t]->"http://identifiers.org/pubchem.substance/7752", species\[LetterSpace]17[t]->"http://identifiers.org/chebi/CHEBI:16108", species\[LetterSpace]17[t]->"http://identifiers.org/kegg.compound/C00111", species\[LetterSpace]17[t]->"http://identifiers.org/pubchem.substance/3411", species\[LetterSpace]17[t]->"http://identifiers.org/3dmet/B00029", species\[LetterSpace]18[t]->"http://identifiers.org/chebi/CHEBI:29052", species\[LetterSpace]18[t]->"http://identifiers.org/kegg.compound/C00118", species\[LetterSpace]18[t]->"http://identifiers.org/pubchem.substance/3418", species\[LetterSpace]18[t]->"http://identifiers.org/cas/591-57-1", species\[LetterSpace]19[t]->"http://identifiers.org/chebi/CHEBI:15846", species\[LetterSpace]19[t]->"http://identifiers.org/kegg.compound/C00003", species\[LetterSpace]19[t]->"http://identifiers.org/pubchem.substance/3305", species\[LetterSpace]19[t]->"http://identifiers.org/cas/53-84-9", species\[LetterSpace]2[t]->"http://identifiers.org/chebi/CHEBI:16761", species\[LetterSpace]2[t]->"http://identifiers.org/kegg.compound/C00008", species\[LetterSpace]2[t]->"http://identifiers.org/pubchem.substance/3310", species\[LetterSpace]2[t]->"http://identifiers.org/cas/20398-34-9", species\[LetterSpace]20[t]->"http://identifiers.org/chebi/CHEBI:16908", species\[LetterSpace]20[t]->"http://identifiers.org/kegg.compound/C00004", species\[LetterSpace]20[t]->"http://identifiers.org/pubchem.substance/3306", species\[LetterSpace]21[t]->"http://identifiers.org/chebi/CHEBI:16001", species\[LetterSpace]21[t]->"http://identifiers.org/kegg.compound/C00236", species\[LetterSpace]21[t]->"http://identifiers.org/pubchem.substance/3535", species\[LetterSpace]21[t]->"http://identifiers.org/cas/38168-82-0", species\[LetterSpace]22[t]->"http://identifiers.org/chebi/CHEBI:15978", species\[LetterSpace]22[t]->"http://identifiers.org/kegg.compound/C00093", species\[LetterSpace]22[t]->"http://identifiers.org/pubchem.substance/3393", species\[LetterSpace]22[t]->"http://identifiers.org/cas/57-03-4", species\[LetterSpace]23[t]->"http://identifiers.org/chebi/CHEBI:17794", species\[LetterSpace]23[t]->"http://identifiers.org/kegg.compound/C00197", species\[LetterSpace]23[t]->"http://identifiers.org/pubchem.substance/3497", species\[LetterSpace]24[t]->"http://identifiers.org/chebi/CHEBI:17754", species\[LetterSpace]24[t]->"http://identifiers.org/kegg.compound/C00116", species\[LetterSpace]24[t]->"http://identifiers.org/pubchem.substance/3416", species\[LetterSpace]24[t]->"http://identifiers.org/3dmet/B00032", species\[LetterSpace]24[t]->"http://identifiers.org/cas/56-81-5", species\[LetterSpace]25[t]->"http://identifiers.org/chebi/CHEBI:4167", species\[LetterSpace]25[t]->"http://identifiers.org/kegg.compound/C00031", species\[LetterSpace]27[t]->"http://identifiers.org/chebi/CHEBI:17754", species\[LetterSpace]27[t]->"http://identifiers.org/kegg.compound/C00116", species\[LetterSpace]27[t]->"http://identifiers.org/pubchem.substance/3416", species\[LetterSpace]27[t]->"http://identifiers.org/3dmet/B00032", species\[LetterSpace]27[t]->"http://identifiers.org/cas/56-81-5", species\[LetterSpace]3[t]->"http://identifiers.org/chebi/CHEBI:15422", species\[LetterSpace]3[t]->"http://identifiers.org/kegg.compound/C00002", species\[LetterSpace]3[t]->"http://identifiers.org/pubchem.substance/3304", species\[LetterSpace]3[t]->"http://identifiers.org/cas/56-65-5", species\[LetterSpace]4[t]->"http://identifiers.org/kegg.compound/C00074", species\[LetterSpace]4[t]->"http://identifiers.org/pubchem.substance/3374", species\[LetterSpace]4[t]->"http://identifiers.org/3dmet/B00019", species\[LetterSpace]4[t]->"http://identifiers.org/cas/138-08-9", species\[LetterSpace]4[t]->"http://identifiers.org/chebi/CHEBI:18021", species\[LetterSpace]5[t]->"http://identifiers.org/chebi/CHEBI:17835", species\[LetterSpace]5[t]->"http://identifiers.org/kegg.compound/C00631", species\[LetterSpace]5[t]->"http://identifiers.org/pubchem.substance/3904", species\[LetterSpace]6[t]->"http://identifiers.org/chebi/CHEBI:16027", species\[LetterSpace]6[t]->"http://identifiers.org/kegg.compound/C00020", species\[LetterSpace]6[t]->"http://identifiers.org/pubchem.substance/3322", species\[LetterSpace]6[t]->"http://identifiers.org/cas/61-19-8", species\[LetterSpace]7[t]->"http://identifiers.org/chebi/CHEBI:16108", species\[LetterSpace]7[t]->"http://identifiers.org/chebi/CHEBI:17794", species\[LetterSpace]7[t]->"http://identifiers.org/kegg.compound/C00197", species\[LetterSpace]7[t]->"http://identifiers.org/pubchem.substance/3497", species\[LetterSpace]8[t]->"http://identifiers.org/chebi/CHEBI:16108", species\[LetterSpace]8[t]->"http://identifiers.org/kegg.compound/C00111", species\[LetterSpace]8[t]->"http://identifiers.org/pubchem.substance/3411", species\[LetterSpace]8[t]->"http://identifiers.org/3dmet/B00029", species\[LetterSpace]9[t]->"http://identifiers.org/chebi/CHEBI:15978", species\[LetterSpace]9[t]->"http://identifiers.org/kegg.compound/C00093", species\[LetterSpace]9[t]->"http://identifiers.org/pubchem.substance/3393", species\[LetterSpace]9[t]->"http://identifiers.org/cas/57-03-4"
};

reactionAnnotations = {
vAKc->"http://identifiers.org/kegg.reaction/R00127", vAKc->"http://identifiers.org/ec-code/2.7.4.3", vAKc->"http://identifiers.org/go/GO:0004017", vAKc->"http://identifiers.org/reactome/REACT\[LetterSpace]643", vAKg->"http://identifiers.org/kegg.reaction/R00127", vAKg->"http://identifiers.org/ec-code/2.7.4.3", vAKg->"http://identifiers.org/go/GO:0004017", vAKg->"http://identifiers.org/reactome/REACT\[LetterSpace]643", vENO->"http://identifiers.org/kegg.reaction/R00658", vENO->"http://identifiers.org/ec-code/4.2.1.11", vENO->"http://identifiers.org/go/GO:0004634", vENO->"http://identifiers.org/reactome/REACT\[LetterSpace]1400", vPGM->"http://identifiers.org/kegg.reaction/R01518", vPGM->"http://identifiers.org/ec-code/5.4.2.1", vPGM->"http://identifiers.org/go/GO:0004619", vPGM->"http://identifiers.org/reactome/REACT\[LetterSpace]576"
};

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

(* Time evolution *)
odes = {
species\[LetterSpace]1'[t] == 1.0*vPK -1.0*vPT, species\[LetterSpace]10'[t] == 1.0*vGT -1.0*vHK, species\[LetterSpace]11'[t] == 1.0*vPGK +1.0*vGK -1.0*vHK -1.0*vPFK -1.0*vAKg, species\[LetterSpace]12'[t] == 1.0*vHK +1.0*vPFK +2.0*vAKg -1.0*vPGK -1.0*vGK, species\[LetterSpace]13'[t] ==  -1.0*vAKg, species\[LetterSpace]14'[t] == 1.0*vHK -1.0*vPGI, species\[LetterSpace]15'[t] == 1.0*vPGI -1.0*vPFK, species\[LetterSpace]16'[t] == 1.0*vPFK -1.0*vALD, species\[LetterSpace]17'[t] == 1.0*vALD +1.0*vANTI -1.0*vTPI -1.0*vGPDH, species\[LetterSpace]18'[t] == 1.0*vALD +1.0*vTPI -1.0*vGAPDH, species\[LetterSpace]19'[t] == 1.0*vGPDH -1.0*vGAPDH, species\[LetterSpace]2'[t] == 1.0*vAU +2.0*vAKc -1.0*vPK, species\[LetterSpace]20'[t] == 1.0*vGAPDH -1.0*vGPDH, species\[LetterSpace]21'[t] == 1.0*vGAPDH -1.0*vPGK, species\[LetterSpace]22'[t] == 1.0*vGPDH -1.0*vGK -1.0*vANTI, species\[LetterSpace]23'[t] == 1.0*vPGK -1.0*vPGT, species\[LetterSpace]24'[t] == 1.0*vGK -1.0*vGlyT, species\[LetterSpace]3'[t] == 1.0*vPK -1.0*vAU -1.0*vAKc, species\[LetterSpace]4'[t] == 1.0*vENO -1.0*vPK, species\[LetterSpace]5'[t] == 1.0*vPGM -1.0*vENO, species\[LetterSpace]6'[t] ==  -1.0*vAKc, species\[LetterSpace]7'[t] == 1.0*vPGT -1.0*vPGM, species\[LetterSpace]8'[t] == 1.0*vGPO -1.0*vANTI, species\[LetterSpace]9'[t] == 1.0*vANTI -1.0*vGPO
};
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]}]