(* 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 = {
php[t], pser[t]
};

initialValues = {
php[0] == 0.6, pser[0] == 0.09
};

rates = {
PDH, PSA, PSP
};

rateEquations = {
PDH -> (cell*p3g*PDH\[LetterSpace]kcatA*serA)/(PDH\[LetterSpace]KAp3g*(1 + ser/PDH\[LetterSpace]KiAser)*(1 + p3g/PDH\[LetterSpace]KAp3g + php[t]/PDH\[LetterSpace]KAphp)), PSA -> (cell*PSA\[LetterSpace]kcatC*serC*php[t])/(PSA\[LetterSpace]KCphp*(1 + php[t]/PSA\[LetterSpace]KCphp + pser[t]/PSA\[LetterSpace]KCpser)), PSP -> (cell*PSP\[LetterSpace]kcatB*serB*pser[t])/(PSP\[LetterSpace]KBpser*(1 + ser/PSP\[LetterSpace]KBser + pser[t]/PSP\[LetterSpace]KBpser))
};

parameters = {
p3g -> 2.36, ser -> 4.9, serA -> 1.15, serB -> 0.25, serC -> 0.1, PDH\[LetterSpace]kcatA -> 0.55, PDH\[LetterSpace]KAp3g -> 1.2, PDH\[LetterSpace]KAphp -> 0.0032, PDH\[LetterSpace]KiAser -> 0.0038, PSA\[LetterSpace]kcatC -> 1.75, PSA\[LetterSpace]KCphp -> 0.0015, PSA\[LetterSpace]KCpser -> 0.0017, PSP\[LetterSpace]kcatB -> 1.43, PSP\[LetterSpace]KBpser -> 0.0015, PSP\[LetterSpace]KBser -> 0.15, cell -> 1.0
};

assignments = {

};

events = {

};

speciesAnnotations = {
p3g[t]->"http://identifiers.org/chebi/CHEBI:58272", php[t]->"http://identifiers.org/chebi/CHEBI:18110", pser[t]->"http://identifiers.org/chebi/CHEBI:57524", ser[t]->"http://identifiers.org/chebi/CHEBI:17115", serA[t]->"http://identifiers.org/uniprot/P0A9T0", serB[t]->"http://identifiers.org/uniprot/P0AGB0", serC[t]->"http://identifiers.org/uniprot/P23721"
};

reactionAnnotations = {
PDH->"http://identifiers.org/ec-code/1.1.1.95", PSA->"http://identifiers.org/ec-code/2.6.1.52", PSP->"http://identifiers.org/ec-code/3.1.3.3"
};

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

(* Time evolution *)
odes = {
php'[t] == 1.0*PDH -1.0*PSA, pser'[t] == 1.0*PSA -1.0*PSP
};
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]}]