(* 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 = {
a[t], b[t], c[t], d[t], e[t], p[t], w[t]
};

initialValues = {
a[0] == 61.4274855570532, b[0] == 25.0461652186299, c[0] == 1.0, d[0] == 1.0, e[0] == 1.0, p[0] == 17.0663407928578, w[0] == 1.76742036253853
};

rates = {
v\[LetterSpace]1, v\[LetterSpace]10, 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 -> (aext*e1*kcat1)/(kmaext1*(1 + aext/kmaext1 + a[t]/kma1)), v\[LetterSpace]10 -> (e10*kcat10*w[t])/(kmw10*(1 + wext/kmwext10 + w[t]/kmw10)), v\[LetterSpace]2 -> (e2*kcat2*a[t])/(kma2*(1 + a[t]/kma2 + b[t]/kmb2)*(1 + w[t]/kmw2)), v\[LetterSpace]3 -> (e3*kcat3*a[t])/(kma3*(1 + sext/kmsext3)^2*(1 + a[t]/kma3 + c[t]/kmc3)), v\[LetterSpace]4 -> (e4*kcat4*a[t])/(kma4*(1 + a[t]/kma4 + d[t]/kmd4)), v\[LetterSpace]5 -> (e5*kcat5*c[t])/(kmc5*(1 + b[t]/kmb5 + c[t]/kmc5)), v\[LetterSpace]6 -> (e6*kcat6*b[t])/(kmb6*(1 + b[t]/kmb6 + p[t]/kmp6)), v\[LetterSpace]7 -> (e7*kcat7*c[t]*d[t])/(kmc7*kmd7*(1 + uext/kmuext7 + d[t]/kmd7)*(1 + c[t]/kmc7 + e[t]/kme7)), v\[LetterSpace]8 -> (e8*kcat8*e[t])/(kme8*(1 + e[t]/kme8 + p[t]/kmp8)*(1 + w[t]/kmw8)), v\[LetterSpace]9 -> (e9*kcat9*p[t])/(kmp9*(1 + fext/kmfext9 + p[t]/kmp9)*(1 + w[t]/kmw9))
};

parameters = {
EXTERNAL -> 1.0, aext -> 100.0, e1 -> 1.17533712011841, e10 -> 4.08468285189287, e2 -> 1.47281664811495, e3 -> 0.0, e4 -> 0.0, e5 -> 0.0, e6 -> 0.764727360489035, e7 -> 0.0, e8 -> 0.0, e9 -> 2.50243601938474, fext -> 2.0, kcat1 -> 1.0, kcat10 -> 1.0, kcat2 -> 2.0, kcat3 -> 6.0, kcat4 -> 1.0, kcat5 -> 15.0, kcat6 -> 2.0, kcat7 -> 1.0, kcat8 -> 1.0, kcat9 -> 1.0, kma1 -> 1.0, kma2 -> 1.0, kma3 -> 0.1, kma4 -> 1.0, kmaext1 -> 1.0, kmb2 -> 1.0, kmb5 -> 10.0, kmb6 -> 1.0, kmc3 -> 10.0, kmc5 -> 0.1, kmc7 -> 1.0, kmd4 -> 1.0, kmd7 -> 1.0, kme7 -> 1.0, kme8 -> 1.0, kmfext9 -> 1.0, kmp6 -> 1.0, kmp8 -> 1.0, kmp9 -> 1.0, kmsext3 -> 10.0, kmuext7 -> 1.0, kmw10 -> 1.0, kmw2 -> 1.0, kmw8 -> 1.0, kmw9 -> 1.0, kmwext10 -> 1.0, sext -> 0.0, uext -> 2.0, wext -> 2.0, default\[LetterSpace]compartment -> 1.0
};

assignments = {

};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
a'[t] == 1.0*v\[LetterSpace]1 -1.0*v\[LetterSpace]3 -1.0*v\[LetterSpace]2 -1.0*v\[LetterSpace]4, b'[t] == 1.0*v\[LetterSpace]5 +1.0*v\[LetterSpace]2 -1.0*v\[LetterSpace]6, c'[t] == 1.0*v\[LetterSpace]3 -1.0*v\[LetterSpace]5 -1.0*v\[LetterSpace]7, d'[t] == 1.0*v\[LetterSpace]4 -1.0*v\[LetterSpace]7, e'[t] == 1.0*v\[LetterSpace]7 -1.0*v\[LetterSpace]8, p'[t] == 1.0*v\[LetterSpace]6 +1.0*v\[LetterSpace]8 -1.0*v\[LetterSpace]9, w'[t] == 1.0*v\[LetterSpace]2 +1.0*v\[LetterSpace]9 +2.0*v\[LetterSpace]8 -1.0*v\[LetterSpace]10
};
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]}]