(* 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 = {
ACALD[t], BPG[t], F16P[t], F6P[t], G6P[t], GLCi[t], NADH[t], P2G[t], P3G[t], PEP[t], PYR[t], TRIO[t]
};

initialValues = {
ACALD[0] == 0.04, BPG[0] == 1*^-05, F16P[0] == 16.39, F6P[0] == 0.78, G6P[0] == 4.45, GLCi[0] == 0.1, NADH[0] == 0.29, P2G[0] == 0.13, P3G[0] == 1.0, PEP[0] == 0.12, PYR[0] == 3.9, TRIO[0] == 1.0
};

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]2, v\[LetterSpace]3, v\[LetterSpace]4, v\[LetterSpace]5, v\[LetterSpace]6, v\[LetterSpace]7, v\[LetterSpace]8, v\[LetterSpace]9
};

rateEquations = {
v\[LetterSpace]1 -> (VMAXGLT*(GLCo/KMGLTGLCo - GLCi[t]/(KEQGLT*KMGLTGLCo)))/(1 + GLCo/KMGLTGLCo + GLCi[t]/KMGLTGLCi + (0.91*GLCo*GLCi[t])/(KMGLTGLCi*KMGLTGLCo)), v\[LetterSpace]10 -> (VMAXPGK*((ADP*KEQPGK*BPG[t])/(KMPGKATP*KMPGKP3G) - (ATP*P3G[t])/(KMPGKATP*KMPGKP3G)))/((1 + ADP/KMPGKADP + ATP/KMPGKATP)*(1 + BPG[t]/KMPGKBPG + P3G[t]/KMPGKP3G)), v\[LetterSpace]11 -> (VMAXPGM*(-(P2G[t]/(KEQPGM*KMPGMP3G)) + P3G[t]/KMPGMP3G))/(1 + P2G[t]/KMPGMP2G + P3G[t]/KMPGMP3G), v\[LetterSpace]12 -> (VMAXENO*(P2G[t]/KMENOP2G - PEP[t]/(KEQENO*KMENOP2G)))/(1 + P2G[t]/KMENOP2G + PEP[t]/KMENOPEP), v\[LetterSpace]13 -> (ADP*VMAXPYK*PEP[t]*(1 + PEP[t]/KMPYKPEP)^(-1 + n10))/((ADP + KMPYKADP)*KMPYKPEP*(L10*((1 + ATP/KMPYKATP)/(1 + F16P[t]/KMPYKF16P))^n10 + (1 + PEP[t]/KMPYKPEP)^n10)), v\[LetterSpace]14 -> (VMAXPDC*PYR[t]^NHPDC)/(KMPDCPYR^NHPDC*(1 + PYR[t]^NHPDC/KMPDCPYR^NHPDC)), v\[LetterSpace]15 -> KSUC, v\[LetterSpace]16 -> -((VMAXADH*((ETOH*(NADt - NADH[t]))/(KIADHNAD*KMADHETOH) - (ACALD[t]*NADH[t])/(KEQADH*KIADHNAD*KMADHETOH)))/(1 + (ETOH*KMADHNAD)/(KIADHNAD*KMADHETOH) + (KMADHNADH*ACALD[t])/(KIADHNADH*KMADHACALD) + (NADt - NADH[t])/KIADHNAD + (ETOH*(NADt - NADH[t]))/(KIADHNAD*KMADHETOH) + (ETOH*ACALD[t]*(NADt - NADH[t]))/(KIADHACALD*KIADHNAD*KMADHETOH) + (KMADHNADH*ACALD[t]*(NADt - NADH[t]))/(KIADHNAD*KIADHNADH*KMADHACALD) + NADH[t]/KIADHNADH + (ETOH*KMADHNAD*NADH[t])/(KIADHNAD*KIADHNADH*KMADHETOH) + (ACALD[t]*NADH[t])/(KIADHNADH*KMADHACALD) + (ETOH*ACALD[t]*NADH[t])/(KIADHETOH*KIADHNADH*KMADHACALD))), v\[LetterSpace]17 -> KACE*ACALD[t], v\[LetterSpace]2 -> (VMAXHK*(-((ADP*G6P[t])/(KEQHK*KMHKATP*KMHKGLCi)) + (ATP*GLCi[t])/(KMHKATP*KMHKGLCi)))/((1 + ADP/KMHKADP + ATP/KMHKATP)*(1 + T6P/KIHKT6P + G6P[t]/KMHKG6P + GLCi[t]/KMHKGLCi)), v\[LetterSpace]3 -> KTRE1, v\[LetterSpace]4 -> (VMAXPGI*(-(F6P[t]/(KEQPGI*KMPGIG6P)) + G6P[t]/KMPGIG6P))/(1 + F6P[t]/KMPGIF6P + G6P[t]/KMPGIG6P), v\[LetterSpace]5 -> KTRE2, v\[LetterSpace]6 -> (ATP*gR*VMAXPFK*F6P[t]*(1 + (ATP*F6P[t])/(KMPFKATP*KMPFKF6P) + (ATP*gR*F6P[t])/(KMPFKATP*KMPFKF6P)))/(KMPFKATP*KMPFKF6P*(((1 + (ATP*CiPFKATP)/KiPFKATP)^2*(1 + (ATP*CPFKATP)/KMPFKATP)^2*(1 + (AMP*CPFKAMP)/KPFKAMP)^2*L0*(1 + (CPFKF26BP*F26BP)/KPFKF26BP + (CPFKF16BP*F16P[t])/KPFKF16BP)^2)/((1 + ATP/KiPFKATP)^2*(1 + AMP/KPFKAMP)^2*(1 + F26BP/KPFKF26BP + F16P[t]/KPFKF16BP)^2) + (1 + (ATP*F6P[t])/(KMPFKATP*KMPFKF6P) + (ATP*gR*F6P[t])/(KMPFKATP*KMPFKF6P))^2)), v\[LetterSpace]7 -> (VMAXALD*(F16P[t]/KMALDF16P - (KEQTPI*TRIO[t]^2)/(KEQALD*(1 + KEQTPI)^2*KMALDF16P)))/(1 + F16P[t]/KMALDF16P + TRIO[t]/((1 + KEQTPI)*KMALDDHAP) + (KEQTPI*TRIO[t])/((1 + KEQTPI)*KMALDGAP) + (KEQTPI*F16P[t]*TRIO[t])/((1 + KEQTPI)*KMALDF16P*KMALDGAPi) + (KEQTPI*TRIO[t]^2)/((1 + KEQTPI)^2*KMALDDHAP*KMALDGAP)), v\[LetterSpace]8 -> KGLY, v\[LetterSpace]9 -> (CGAPDH*(-((VMAXGAPDHr*BPG[t]*NADH[t])/(KMGAPDHBPG*KMGAPDHNADH)) + (KEQTPI*VMAXGAPDHf*(NADt - NADH[t])*TRIO[t])/((1 + KEQTPI)*KMGAPDHGAP*KMGAPDHNAD)))/((1 + (NADt - NADH[t])/KMGAPDHNAD + NADH[t]/KMGAPDHNADH)*(1 + BPG[t]/KMGAPDHBPG + (KEQTPI*TRIO[t])/((1 + KEQTPI)*KMGAPDHGAP)))
};

parameters = {
ACE -> 10.0, ADP -> 1.09, AMP -> 0.37, ATP -> 4.7, CGAPDH -> 1.0, CO2 -> 1.0, CPFKAMP -> 0.0845, CPFKATP -> 3.0, CPFKF16BP -> 0.397, CPFKF26BP -> 0.0174, CiPFKATP -> 100.0, ETOH -> 25.0, EXTERNAL -> 1.0, F26BP -> 0.006, GLCo -> 50.0, GLY -> 10.0, KACE -> 0.5, KEQADH -> 6.9*^-05, KEQALD -> 0.069, KEQENO -> 6.7, KEQGLT -> 1.0, KEQHK -> 3800.0, KEQPGI -> 0.314, KEQPGK -> 3200.0, KEQPGM -> 0.19, KEQTPI -> 0.045, KGLY -> 21.5, KIADHACALD -> 1.1, KIADHETOH -> 90.0, KIADHNAD -> 0.92, KIADHNADH -> 0.031, KIHKT6P -> 0.04, KMADHACALD -> 1.11, KMADHETOH -> 17.0, KMADHNAD -> 0.17, KMADHNADH -> 0.11, KMALDDHAP -> 2.4, KMALDF16P -> 0.3, KMALDGAP -> 2.0, KMALDGAPi -> 10.0, KMENOP2G -> 0.04, KMENOPEP -> 0.5, KMGAPDHBPG -> 0.51, KMGAPDHGAP -> 0.39, KMGAPDHNAD -> 2.85, KMGAPDHNADH -> 0.007, KMGLTGLCi -> 7.0, KMGLTGLCo -> 7.0, KMHKADP -> 0.23, KMHKATP -> 0.15, KMHKG6P -> 30.0, KMHKGLCi -> 0.08, KMPDCPYR -> 6.36, KMPFKATP -> 0.71, KMPFKF6P -> 0.1, KMPGIF6P -> 0.3, KMPGIG6P -> 1.4, KMPGKADP -> 0.2, KMPGKATP -> 0.3, KMPGKBPG -> 0.003, KMPGKP3G -> 0.53, KMPGMP2G -> 0.08, KMPGMP3G -> 1.2, KMPYKADP -> 0.3, KMPYKATP -> 9.3, KMPYKF16P -> 0.2, KMPYKPEP -> 0.19, KPFKAMP -> 0.0995, KPFKF16BP -> 0.111, KPFKF26BP -> 0.000682, KSUC -> 0.2, KTRE1 -> 1.9, KTRE2 -> 0.0, KiPFKATP -> 0.65, L0 -> 0.66, L10 -> 60000.0, NADt -> 1.59, NHPDC -> 1.9, T6P -> 0.59, TREH -> 2.0, VMAXADH -> 744.0, VMAXALD -> 161.0, VMAXENO -> 272.0, VMAXGAPDHf -> 1450.0, VMAXGAPDHr -> 840.0, VMAXGLT -> 95.0, VMAXHK -> 227.0, VMAXPDC -> 172.0, VMAXPFK -> 93.0, VMAXPGI -> 856.0, VMAXPGK -> 1962.0, VMAXPGM -> 403.0, VMAXPYK -> 480.0, gR -> 5.12, n10 -> 4.0, default\[LetterSpace]compartment -> 1.0
};

assignments = {

};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
ACALD'[t] == 1.0*v\[LetterSpace]14 -1.0*v\[LetterSpace]16 -1.0*v\[LetterSpace]17, BPG'[t] == 1.0*v\[LetterSpace]9 -1.0*v\[LetterSpace]10, F16P'[t] == 1.0*v\[LetterSpace]6 -1.0*v\[LetterSpace]7, F6P'[t] == 1.0*v\[LetterSpace]4 -1.0*v\[LetterSpace]6, G6P'[t] == 1.0*v\[LetterSpace]2 -2.0*v\[LetterSpace]5 -1.0*v\[LetterSpace]4, GLCi'[t] == 1.0*v\[LetterSpace]1 +2.0*v\[LetterSpace]3 -1.0*v\[LetterSpace]2, NADH'[t] == 3.0*v\[LetterSpace]15 +1.0*v\[LetterSpace]17 +1.0*v\[LetterSpace]9 -1.0*v\[LetterSpace]16 -1.0*v\[LetterSpace]8, P2G'[t] == 1.0*v\[LetterSpace]11 -1.0*v\[LetterSpace]12, P3G'[t] == 1.0*v\[LetterSpace]10 -1.0*v\[LetterSpace]11, PEP'[t] == 1.0*v\[LetterSpace]12 -1.0*v\[LetterSpace]13, PYR'[t] == 1.0*v\[LetterSpace]13 -1.0*v\[LetterSpace]14 -2.0*v\[LetterSpace]15, TRIO'[t] == 2.0*v\[LetterSpace]7 -1.0*v\[LetterSpace]8 -1.0*v\[LetterSpace]9
};
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]}]