(* 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] == 0.25, F6P[0] == 0.2, G6P[0] == 1.2, GLCi[0] == 0.1, NADH[0] == 0.29, P2G[0] == 0.06, P3G[0] == 0.54, PEP[0] == 0.5, PYR[0] == 0.14, 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*(-(GLCi[t]/(KEQGLT*KMGLTGLCo)) + Piecewise[{{glc1, t < tswitch}}, glc2]/KMGLTGLCo))/(1 + GLCi[t]/KMGLTGLCi + Piecewise[{{glc1, t < tswitch}}, glc2]/KMGLTGLCo + (0.91*GLCi[t]*Piecewise[{{glc1, t < tswitch}}, glc2])/(KMGLTGLCi*KMGLTGLCo)), v\[LetterSpace]10 -> (VMAXPGK*((ADP*KEQPGK*BPG[t])/(KMPGKATP*KMPGKP3G) - (P3G[t]*Piecewise[{{atp1, t < tswitch}}, atp2])/(KMPGKATP*KMPGKP3G)))/((1 + BPG[t]/KMPGKBPG + P3G[t]/KMPGKP3G)*(1 + ADP/KMPGKADP + Piecewise[{{atp1, t < tswitch}}, atp2]/KMPGKATP)), 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*((1 + PEP[t]/KMPYKPEP)^n10 + L10*((1 + Piecewise[{{atp1, t < tswitch}}, atp2]/KMPYKATP)/(1 + F16P[t]/KMPYKF16P))^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)) + (GLCi[t]*Piecewise[{{atp1, t < tswitch}}, atp2])/(KMHKATP*KMHKGLCi)))/((1 + T6P/KIHKT6P + G6P[t]/KMHKG6P + GLCi[t]/KMHKGLCi)*(1 + ADP/KMHKADP + Piecewise[{{atp1, t < tswitch}}, atp2]/KMHKATP)), 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 -> (gR*VMAXPFK*F6P[t]*Piecewise[{{atp1, t < tswitch}}, atp2]*(1 + (F6P[t]*Piecewise[{{atp1, t < tswitch}}, atp2])/(KMPFKATP*KMPFKF6P) + (gR*F6P[t]*Piecewise[{{atp1, t < tswitch}}, atp2])/(KMPFKATP*KMPFKF6P)))/(KMPFKATP*KMPFKF6P*(((1 + (AMP*CPFKAMP)/KPFKAMP)^2*L0*(1 + (CPFKF26BP*F26BP)/KPFKF26BP + (CPFKF16BP*F16P[t])/KPFKF16BP)^2*(1 + (CiPFKATP*Piecewise[{{atp1, t < tswitch}}, atp2])/KiPFKATP)^2*(1 + (CPFKATP*Piecewise[{{atp1, t < tswitch}}, atp2])/KMPFKATP)^2)/((1 + AMP/KPFKAMP)^2*(1 + F26BP/KPFKF26BP + F16P[t]/KPFKF16BP)^2*(1 + Piecewise[{{atp1, t < tswitch}}, atp2]/KiPFKATP)^2) + (1 + (F6P[t]*Piecewise[{{atp1, t < tswitch}}, atp2])/(KMPFKATP*KMPFKF6P) + (gR*F6P[t]*Piecewise[{{atp1, t < tswitch}}, atp2])/(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.0, AMP -> 0.3, 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.014, 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 -> 17.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 -> 1.0, KMGLTGLCo -> 1.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.9, KTRE1 -> 0.0, KTRE2 -> 0.0, KiPFKATP -> 0.65, L0 -> 0.66, L10 -> 60000.0, NADt -> 1.59, NHPDC -> 1.9, T6P -> 0.2, TREH -> 2.0, VMAXADH -> 932.0, VMAXALD -> 310.0, VMAXENO -> 357.0, VMAXGAPDHf -> 1487.0, VMAXGAPDHr -> 853.0, VMAXGLT -> 160.0, VMAXHK -> 213.0, VMAXPDC -> 395.0, VMAXPFK -> 213.0, VMAXPGI -> 787.0, VMAXPGK -> 2512.0, VMAXPGM -> 856.0, VMAXPYK -> 820.0, atp1 -> 3.0, atp2 -> 1.5, gR -> 5.12, glc1 -> 0.2, glc2 -> 5.0, n10 -> 4.0, tswitch -> 10.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]}]