(* 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], fm[t], g[t], p[t]
};

initialValues = {
a[0] == 43.6059113414968, b[0] == 1.0, c[0] == 39.152640553106, d[0] == 875.728067679313, fm[0] == 0.0627769167303887, g[0] == 1.0, p[0] == 1.071877291459
};

rates = {
v\[LetterSpace]1, v\[LetterSpace]10, v\[LetterSpace]11, 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*g[t])/(kmg10*(1 + d[t]/kmd10 + g[t]/kmg10)), v\[LetterSpace]11 -> (e11*kcat11*fm[t]*p[t])/(kmfm11*kmp11*(1 + (ftot - fm[t])/kmfp11 + fm[t]/kmfm11)*(1 + pext/kmpext11 + p[t]/kmp11)), v\[LetterSpace]2 -> (bext*e2*kcat2*a[t]*(ftot - fm[t])*(1 - (b[t]*fm[t])/(bext*keq2*(ftot - fm[t]))))/(kaa2*kmbext2*kmfp2*(1 + a[t]/kaa2)*(1 + bext/kmbext2 + b[t]/kmb2)*(1 + (ftot - fm[t])/kmfp2 + fm[t]/kmfm2)), v\[LetterSpace]3 -> (e3*kcat3*a[t])/(kma3*(1 + a[t]/kma3 + b[t]/kmb3)), v\[LetterSpace]4 -> (e4*kcat4*a[t]*(1 - c[t]/(keq4*a[t])))/(kma4*(1 + a[t]/kma4 + c[t]/kmc4)), v\[LetterSpace]5 -> (e5*kcat5*a[t]*c[t])/(kac5*kma5*(1 + c[t]/kac5)*(1 + a[t]/kma5 + d[t]/kmd5)), v\[LetterSpace]6 -> (e6*kcat6*b[t]*(1 - c[t]/(keq6*b[t])))/(kmb6*(1 + b[t]/kmb6 + c[t]/kmc6)), v\[LetterSpace]7 -> (e7*kcat7*b[t])/(kmb7*(1 + c[t]/kic7)*(1 + b[t]/kmb7 + p[t]/kmp7)), v\[LetterSpace]8 -> (e8*kcat8*c[t]*d[t]*(ftot - fm[t]))/(kmc8*kmd8*kmfp8*(1 + d[t]/kmd8)*(1 + (ftot - fm[t])/kmfp8 + fm[t]/kmfm8)*(1 + c[t]/kmc8 + p[t]/kmp8)), v\[LetterSpace]9 -> (e9*kcat9*d[t])/(kmd9*(1 + d[t]/kmd9 + g[t]/kmg9))
};

parameters = {
EXTERNAL -> 1.0, aext -> 2.32175122459047, bext -> 1.82097585950486, e1 -> 1.97972425219533, e10 -> 1*^-06, e11 -> 0.548430069295763, e2 -> 1*^-06, e3 -> 1*^-06, e4 -> 1.8765818035848, e5 -> 0.154469876349102, e6 -> 1*^-06, e7 -> 1*^-06, e8 -> 12.9070774037041, e9 -> 1*^-06, ftot -> 5.0, kaa2 -> 0.537724121637483, kac5 -> 0.45109181910517, kcat1 -> 4.35128295344753, kcat10 -> 2.04583181726299, kcat11 -> 15.6308547611266, kcat2 -> 2.08470697013596, kcat3 -> 0.659201960419019, kcat4 -> 1.21740513898481, kcat5 -> 7.48882370581632, kcat6 -> 0.26424160469917, kcat7 -> 0.515186023045187, kcat8 -> 0.148031799420018, kcat9 -> 4.00713635555074, keq2 -> 2.0, keq4 -> 2.0, keq6 -> 2.0, kic7 -> 1.01557360659626, kma1 -> 2.65103841264887, kma3 -> 0.183420075929795, kma4 -> 6.89306707926306, kma5 -> 0.329345484841003, kmaext1 -> 0.232175122459047, kmb2 -> 0.190218968479388, kmb3 -> 0.965624239632328, kmb6 -> 0.890368839214066, kmb7 -> 5.5364593839414, kmbext2 -> 18.2097585950486, kmc4 -> 6.56779619616803, kmc6 -> 55.7611490930134, kmc8 -> 4.70771101695726, kmd10 -> 0.417200201130384, kmd5 -> 33.0655370631572, kmd8 -> 1.79688793384447, kmd9 -> 0.104280856663846, kmfm11 -> 0.140510704125535, kmfm2 -> 8.42083191123236, kmfm8 -> 0.660727064408495, kmfp11 -> 7.26391447299304, kmfp2 -> 1.27736476875397, kmfp8 -> 3.23838653900196, kmg10 -> 1.38293464784513, kmg9 -> 0.43294178975732, kmp11 -> 0.571956386391092, kmp7 -> 11.3398543717022, kmp8 -> 6.64006979927262, kmpext11 -> 1.32831286086816, pext -> 0.132831286086816, 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]5 -1.0*v\[LetterSpace]4, b'[t] == 1.0*v\[LetterSpace]3 +1.0*v\[LetterSpace]2 -1.0*v\[LetterSpace]7 -1.0*v\[LetterSpace]6, c'[t] == 1.0*v\[LetterSpace]6 +1.0*v\[LetterSpace]4 -1.0*v\[LetterSpace]8, d'[t] == 1.0*v\[LetterSpace]10 +1.0*v\[LetterSpace]5 -1.0*v\[LetterSpace]8 -1.0*v\[LetterSpace]9, fm'[t] == 1.0*v\[LetterSpace]8 +1.0*v\[LetterSpace]2 -1.0*v\[LetterSpace]11, g'[t] == 1.0*v\[LetterSpace]9 -1.0*v\[LetterSpace]10, p'[t] == 1.0*v\[LetterSpace]7 +1.0*v\[LetterSpace]8 -1.0*v\[LetterSpace]11
};
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]}]