(* 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 = {
species\[LetterSpace]1[t], species\[LetterSpace]10[t], species\[LetterSpace]11[t], species\[LetterSpace]2[t], species\[LetterSpace]3[t], species\[LetterSpace]4[t], species\[LetterSpace]5[t], species\[LetterSpace]6[t], species\[LetterSpace]7[t], species\[LetterSpace]8[t], species\[LetterSpace]9[t]
};

initialValues = {
species\[LetterSpace]1[0] == 300.0, species\[LetterSpace]10[0] == 200.0, species\[LetterSpace]11[0] == 20.0, species\[LetterSpace]2[0] == 0.0, species\[LetterSpace]3[0] == 1199.99994221325, species\[LetterSpace]4[0] == 0.0, species\[LetterSpace]5[0] == 0.0, species\[LetterSpace]6[0] == 1199.99994221325, species\[LetterSpace]7[0] == 0.0, species\[LetterSpace]8[0] == 0.0, species\[LetterSpace]9[0] == 100.0
};

rates = {
reaction\[LetterSpace]1, reaction\[LetterSpace]10, reaction\[LetterSpace]2, reaction\[LetterSpace]3, reaction\[LetterSpace]4, reaction\[LetterSpace]5, reaction\[LetterSpace]6, reaction\[LetterSpace]7, reaction\[LetterSpace]8, reaction\[LetterSpace]9
};

rateEquations = {
reaction\[LetterSpace]1 -> compartment\[LetterSpace]1*function\[LetterSpace]1[species\[LetterSpace]1[t], parameter\[LetterSpace]1, reaction\[LetterSpace]1\[LetterSpace]k1, species\[LetterSpace]11[t]], reaction\[LetterSpace]10 -> compartment\[LetterSpace]1*function\[LetterSpace]10[species\[LetterSpace]7[t], species\[LetterSpace]5[t], species\[LetterSpace]4[t], species\[LetterSpace]8[t], species\[LetterSpace]3[t], species\[LetterSpace]6[t], species\[LetterSpace]10[t], parameter\[LetterSpace]10, parameter\[LetterSpace]13, parameter\[LetterSpace]14, parameter\[LetterSpace]12, parameter\[LetterSpace]9, reaction\[LetterSpace]10\[LetterSpace]k10b], reaction\[LetterSpace]2 -> compartment\[LetterSpace]1*function\[LetterSpace]2[species\[LetterSpace]2[t], species\[LetterSpace]1[t], species\[LetterSpace]9[t], species\[LetterSpace]5[t], species\[LetterSpace]4[t], species\[LetterSpace]3[t], parameter\[LetterSpace]11, parameter\[LetterSpace]2, parameter\[LetterSpace]5, parameter\[LetterSpace]6, reaction\[LetterSpace]2\[LetterSpace]k2a], reaction\[LetterSpace]3 -> compartment\[LetterSpace]1*function\[LetterSpace]3[reaction\[LetterSpace]3\[LetterSpace]k3, species\[LetterSpace]2[t], species\[LetterSpace]3[t], parameter\[LetterSpace]3, species\[LetterSpace]4[t], parameter\[LetterSpace]4], reaction\[LetterSpace]4 -> compartment\[LetterSpace]1*function\[LetterSpace]4[reaction\[LetterSpace]4\[LetterSpace]k4, species\[LetterSpace]2[t], species\[LetterSpace]4[t], parameter\[LetterSpace]4, species\[LetterSpace]3[t], parameter\[LetterSpace]3], reaction\[LetterSpace]5 -> compartment\[LetterSpace]1*function\[LetterSpace]5[species\[LetterSpace]5[t], species\[LetterSpace]4[t], species\[LetterSpace]7[t], species\[LetterSpace]8[t], species\[LetterSpace]9[t], species\[LetterSpace]3[t], species\[LetterSpace]6[t], reaction\[LetterSpace]5\[LetterSpace]k5a, parameter\[LetterSpace]5, species\[LetterSpace]1[t], species\[LetterSpace]2[t], reaction\[LetterSpace]5\[LetterSpace]k5b, species\[LetterSpace]10[t], parameter\[LetterSpace]13, parameter\[LetterSpace]6, parameter\[LetterSpace]11, parameter\[LetterSpace]2, parameter\[LetterSpace]14, parameter\[LetterSpace]12, parameter\[LetterSpace]10, parameter\[LetterSpace]9], reaction\[LetterSpace]6 -> compartment\[LetterSpace]1*function\[LetterSpace]6[species\[LetterSpace]9[t], species\[LetterSpace]4[t], species\[LetterSpace]5[t], species\[LetterSpace]7[t], species\[LetterSpace]8[t], species\[LetterSpace]3[t], species\[LetterSpace]6[t], reaction\[LetterSpace]6\[LetterSpace]k6a, parameter\[LetterSpace]6, species\[LetterSpace]1[t], species\[LetterSpace]2[t], species\[LetterSpace]10[t], parameter\[LetterSpace]14, parameter\[LetterSpace]2, parameter\[LetterSpace]11, parameter\[LetterSpace]5, reaction\[LetterSpace]6\[LetterSpace]k6b, parameter\[LetterSpace]13, parameter\[LetterSpace]12, parameter\[LetterSpace]10, parameter\[LetterSpace]9], reaction\[LetterSpace]7 -> compartment\[LetterSpace]1*function\[LetterSpace]7[reaction\[LetterSpace]7\[LetterSpace]k7, species\[LetterSpace]5[t], species\[LetterSpace]6[t], parameter\[LetterSpace]7, species\[LetterSpace]7[t], parameter\[LetterSpace]8], reaction\[LetterSpace]8 -> compartment\[LetterSpace]1*function\[LetterSpace]8[reaction\[LetterSpace]8\[LetterSpace]k7, species\[LetterSpace]5[t], species\[LetterSpace]7[t], parameter\[LetterSpace]8, species\[LetterSpace]6[t], parameter\[LetterSpace]7], reaction\[LetterSpace]9 -> compartment\[LetterSpace]1*function\[LetterSpace]9[species\[LetterSpace]8[t], species\[LetterSpace]9[t], species\[LetterSpace]4[t], species\[LetterSpace]7[t], species\[LetterSpace]3[t], species\[LetterSpace]6[t], species\[LetterSpace]10[t], parameter\[LetterSpace]9, parameter\[LetterSpace]13, parameter\[LetterSpace]14, parameter\[LetterSpace]12, parameter\[LetterSpace]10, reaction\[LetterSpace]9\[LetterSpace]k9b]
};

parameters = {
parameter\[LetterSpace]1 -> 100.0, parameter\[LetterSpace]10 -> 108.6, parameter\[LetterSpace]11 -> 3*^+51, parameter\[LetterSpace]12 -> 3*^+51, parameter\[LetterSpace]13 -> 24.3, parameter\[LetterSpace]14 -> 108.6, parameter\[LetterSpace]2 -> 54.3, parameter\[LetterSpace]3 -> 50.5, parameter\[LetterSpace]4 -> 500.0, parameter\[LetterSpace]5 -> 24.3, parameter\[LetterSpace]6 -> 108.6, parameter\[LetterSpace]7 -> 50.5, parameter\[LetterSpace]8 -> 500.0, parameter\[LetterSpace]9 -> 24.3, reaction\[LetterSpace]1\[LetterSpace]k1 -> 1.0, reaction\[LetterSpace]2\[LetterSpace]k2a -> 0.086, reaction\[LetterSpace]3\[LetterSpace]k3 -> 0.01, reaction\[LetterSpace]4\[LetterSpace]k4 -> 15.0, reaction\[LetterSpace]5\[LetterSpace]k5a -> 0.092, reaction\[LetterSpace]5\[LetterSpace]k5b -> 0.092, reaction\[LetterSpace]6\[LetterSpace]k6a -> 0.086, reaction\[LetterSpace]6\[LetterSpace]k6b -> 0.086, reaction\[LetterSpace]7\[LetterSpace]k7 -> 0.01, reaction\[LetterSpace]8\[LetterSpace]k7 -> 15.0, reaction\[LetterSpace]9\[LetterSpace]k9b -> 0.092, reaction\[LetterSpace]10\[LetterSpace]k10b -> 0.086, compartment\[LetterSpace]1 -> 1.0, compartment\[LetterSpace]2 -> 1.0
};

assignments = {
function\[LetterSpace]10[MK\[LetterSpace]P_,MKK\[LetterSpace]PP_,MKK\[LetterSpace]P_,MK\[LetterSpace]PP_,MKK_,MK_,P2_,K10b_,K5b_,K6b_,Kse2_,K9b_,k10b_] -> (k10b*MK\[LetterSpace]P*P2)/(K10b*(1 + MK/Kse2 + MKK/Kse2 + MKK\[LetterSpace]P/K6b + MKK\[LetterSpace]PP/K5b + MK\[LetterSpace]P/K10b + MK\[LetterSpace]PP/K9b)), function\[LetterSpace]8[k7_,MKK\[LetterSpace]PP_,MK\[LetterSpace]P_,K8_,MK_,K7_] -> (k7*MKK\[LetterSpace]PP*MK\[LetterSpace]P)/(K8*(1 + MK/K7 + MK\[LetterSpace]P/K8)), function\[LetterSpace]9[MK\[LetterSpace]PP_,MKK\[LetterSpace]PP_,MKK\[LetterSpace]P_,MK\[LetterSpace]P_,MKK_,MK_,P2_,K9b_,K5b_,K6b_,Kse2_,K10b_,k9b_] -> (k9b*MK\[LetterSpace]PP*P2)/(K9b*(1 + MK/Kse2 + MKK/Kse2 + MKK\[LetterSpace]P/K6b + MKK\[LetterSpace]PP/K5b + MK\[LetterSpace]P/K10b + MK\[LetterSpace]PP/K9b)), function\[LetterSpace]7[k7_,MKK\[LetterSpace]PP_,MK_,K7_,MK\[LetterSpace]P_,K8_] -> (k7*MK*MKK\[LetterSpace]PP)/(K7*(1 + MK/K7 + MK\[LetterSpace]P/K8)), function\[LetterSpace]5[MKK\[LetterSpace]PP_,MKK\[LetterSpace]P_,MK\[LetterSpace]P_,MK\[LetterSpace]PP_,P1_,MKK_,MK_,k5a_,K5a_,MKKK_,MKKK\[LetterSpace]P_,k5b_,P2_,K5b_,K6a_,Kse1_,K2a_,K6b_,Kse2_,K10b_,K9b_] -> (k5a*MKK\[LetterSpace]PP*P1)/(K5a*(1 + MKK/Kse1 + MKKK/Kse1 + MKKK\[LetterSpace]P/K2a + MKK\[LetterSpace]P/K6a + MKK\[LetterSpace]PP/K5a)) + (k5b*MKK\[LetterSpace]PP*P2)/(K5b*(1 + MK/Kse2 + MKK/Kse2 + MKK\[LetterSpace]P/K6b + MKK\[LetterSpace]PP/K5b + MK\[LetterSpace]P/K10b + MK\[LetterSpace]PP/K9b)), function\[LetterSpace]4[k4_,MKKK\[LetterSpace]P_,MKK\[LetterSpace]P_,K4_,MKK_,K3_] -> (k4*MKKK\[LetterSpace]P*MKK\[LetterSpace]P)/(K4*(1 + MKK/K3 + MKK\[LetterSpace]P/K4)), function\[LetterSpace]1[MKKK_,K1_,k1_,Sig_] -> (k1*MKKK*Sig)/(K1 + MKKK), function\[LetterSpace]6[P1_,MKK\[LetterSpace]P_,MKK\[LetterSpace]PP_,MK\[LetterSpace]P_,MK\[LetterSpace]PP_,MKK_,MK_,k6a_,K6a_,MKKK_,MKKK\[LetterSpace]P_,P2_,K6b_,K2a_,Kse1_,K5a_,k6b_,K5b_,Kse2_,K10b_,K9b_] -> (k6a*MKK\[LetterSpace]P*P1)/(K6a*(1 + MKK/Kse1 + MKKK/Kse1 + MKKK\[LetterSpace]P/K2a + MKK\[LetterSpace]P/K6a + MKK\[LetterSpace]PP/K5a)) + (k6b*MKK\[LetterSpace]P*P2)/(K6b*(1 + MK/Kse2 + MKK/Kse2 + MKK\[LetterSpace]P/K6b + MKK\[LetterSpace]PP/K5b + MK\[LetterSpace]P/K10b + MK\[LetterSpace]PP/K9b)), function\[LetterSpace]3[k3_,MKKK\[LetterSpace]P_,MKK_,K3_,MKK\[LetterSpace]P_,K4_] -> (k3*MKK*MKKK\[LetterSpace]P)/(K3*(1 + MKK/K3 + MKK\[LetterSpace]P/K4)), function\[LetterSpace]2[MKKK\[LetterSpace]P_,MKKK_,P1_,MKK\[LetterSpace]PP_,MKK\[LetterSpace]P_,MKK_,Kse1_,K2a_,K5a_,K6a_,k2a_] -> (k2a*MKKK\[LetterSpace]P*P1)/(K2a*(1 + MKK/Kse1 + MKKK/Kse1 + MKKK\[LetterSpace]P/K2a + MKK\[LetterSpace]P/K6a + MKK\[LetterSpace]PP/K5a))
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
species\[LetterSpace]1'[t] == 1.0*reaction\[LetterSpace]2 -1.0*reaction\[LetterSpace]1, species\[LetterSpace]10'[t] == 0.0 , species\[LetterSpace]11'[t] == 0.0 , species\[LetterSpace]2'[t] == 1.0*reaction\[LetterSpace]1 -1.0*reaction\[LetterSpace]2, species\[LetterSpace]3'[t] == 1.0*reaction\[LetterSpace]6 -1.0*reaction\[LetterSpace]3, species\[LetterSpace]4'[t] == 1.0*reaction\[LetterSpace]3 +1.0*reaction\[LetterSpace]5 -1.0*reaction\[LetterSpace]4 -1.0*reaction\[LetterSpace]6, species\[LetterSpace]5'[t] == 1.0*reaction\[LetterSpace]4 -1.0*reaction\[LetterSpace]5, species\[LetterSpace]6'[t] == 1.0*reaction\[LetterSpace]10 -1.0*reaction\[LetterSpace]7, species\[LetterSpace]7'[t] == 1.0*reaction\[LetterSpace]7 +1.0*reaction\[LetterSpace]9 -1.0*reaction\[LetterSpace]8 -1.0*reaction\[LetterSpace]10, species\[LetterSpace]8'[t] == 1.0*reaction\[LetterSpace]8 -1.0*reaction\[LetterSpace]9, species\[LetterSpace]9'[t] == 0.0 
};
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]}]