(* 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 = {
X1[t], X10[t], X11[t], X12[t], X13[t], X14[t], X15[t], X16[t], X17[t], X18[t], X19[t], X2[t], X20[t], X21[t], X22[t], X23[t], X24[t], X25[t], X26[t], X27[t], X28[t], X3[t], X4[t], X5[t], X6[t], X7[t], X8[t], X9[t]
};

initialValues = {
X1[0] == 100.0, X10[0] == 0.9881, X11[0] == 42.7224, X12[0] == 0.0008, X13[0] == 6.1879, X14[0] == 8.8121, X15[0] == 3.4392, X16[0] == 200.0, X17[0] == 0.0, X18[0] == 112.5585, X19[0] == 6.486, X2[0] == 0.0, X20[0] == 296.1137, X21[0] == 3.8863, X22[0] == 297.8897, X23[0] == 2.1103, X24[0] == 180.9595, X25[0] == 418.1788, X26[0] == 0.8619, X27[0] == 10.263, X28[0] == 0.85544, X3[0] == 0.0153, X4[0] == 0.0076, X5[0] == 49.1372, X6[0] == 0.0015, X7[0] == 96.6019, X8[0] == 0.002, X9[0] == 0.002
};

rates = {
R1, R10, R11, R12, R13, R14, R15, R16, R17, R18, R19, R2, R20, R21, R22, R23, R24, R25, R26, R27, R28, R29, R3, R30, R31, R32, R33, R4, R5, R6, R7, R8, R9
};

rateEquations = {
R1 -> cytoplasm*k1*W*X1[t], R10 -> cytoplasm*k10*X9[t], R11 -> cytoplasm*k11*X10[t], R12 -> cytoplasm*V12, R13 -> k13*nucleus*X11[t], R14 -> nucleus*(k14 + k21*(X11[t] + X14[t])), R15 -> cytoplasm*k15*X12[t], R16 -> nucleus*(k\[LetterSpace]plus16*X11[t]*X13[t] - k\[LetterSpace]minus16*X14[t]), R17 -> cytoplasm*(-(k\[LetterSpace]minus17*X15[t]) + k\[LetterSpace]plus17*X11[t]*X7[t]), R18 -> (cytoplasm*Ki*Vmax1*W*X16[t])/((Km1 + X16[t])*(Ki + X23[t])), R19 -> (cytoplasm*Vmax2*X17[t])/(Km2 + X17[t]), R2 -> cytoplasm*k2*X2[t], R20 -> (cytoplasm*kcat1*X17[t]*X18[t])/(Km3 + X18[t]), R21 -> (cytoplasm*Vmax3*X19[t])/(Km4 + X19[t]), R22 -> (cytoplasm*kcat2*X19[t]*X20[t])/(Km5 + X20[t]), R23 -> (cytoplasm*Vmax4*X21[t])/(Km6 + X21[t]), R24 -> (cytoplasm*kcat3*X21[t]*X22[t])/(Km7 + X22[t]), R25 -> (cytoplasm*Vmax5*X23[t])/(Km8 + X23[t]), R26 -> (cytoplasm*kcat4*X23[t]*X24[t])/(Km9 + X24[t]), R27 -> cytoplasm*(-(k19*X24[t]) + k18*X18[t]*X25[t]), R28 -> (cytoplasm*Vmax6*X26[t])/(Km10 + X26[t]), R29 -> (cytoplasm*kcat5*X14[t]^n1)/(Km11^n1 + X14[t]^n1), R3 -> cytoplasm*k3*X2[t]*X4[t], R30 -> cytoplasm*k20*X27[t], R31 -> (cytoplasm*kcat6*X18[t]*X27[t])/(Km12 + X18[t]), R32 -> (cytoplasm*kcat7*X23[t]*X5[t])/(Km13 + X5[t]), R33 -> (cytoplasm*Vmax7*X28[t])/(Km14 + X28[t]), R4 -> cytoplasm*k4*X4[t], R5 -> cytoplasm*k5*X3[t], R6 -> cytoplasm*(-(k\[LetterSpace]minus6*X4[t]) + k\[LetterSpace]plus6*X5[t]*X6[t]), R7 -> cytoplasm*(-(k\[LetterSpace]minus7*X6[t]) + k\[LetterSpace]plus7*X12[t]*X7[t]), R8 -> cytoplasm*(k\[LetterSpace]plus8*X11[t]*X3[t] - k\[LetterSpace]minus8*X8[t]), R9 -> cytoplasm*k9*X8[t]
};

parameters = {
Ki -> 9.0, Km1 -> 10.0, Km10 -> 12.0, Km11 -> 15.0, Km12 -> 15.0, Km13 -> 15.0, Km14 -> 15.0, Km2 -> 8.0, Km3 -> 15.0, Km4 -> 15.0, Km5 -> 15.0, Km6 -> 15.0, Km7 -> 15.0, Km8 -> 15.0, Km9 -> 9.0, V12 -> 0.423, Vmax1 -> 150.0, Vmax2 -> 15.0, Vmax3 -> 45.0, Vmax4 -> 45.0, Vmax5 -> 45.0, Vmax6 -> 45.0, Vmax7 -> 45.0, W -> 0.0, flag\[LetterSpace]for\[LetterSpace]wnt\[LetterSpace]signal -> 0.0, k1 -> 0.182, k10 -> 206.0, k11 -> 0.417, k13 -> 0.000257, k14 -> 8.22*^-05, k15 -> 0.167, k18 -> 0.15, k19 -> 39.0, k2 -> 0.0182, k20 -> 0.015, k21 -> 1*^-06, k3 -> 0.05, k4 -> 0.267, k5 -> 0.133, k9 -> 206.0, k\[LetterSpace]minus16 -> 30.0, k\[LetterSpace]minus17 -> 1200.0, k\[LetterSpace]minus6 -> 0.909, k\[LetterSpace]minus7 -> 50.0, k\[LetterSpace]minus8 -> 120.0, k\[LetterSpace]plus16 -> 1.0, k\[LetterSpace]plus17 -> 1.0, k\[LetterSpace]plus6 -> 0.0909, k\[LetterSpace]plus7 -> 1.0, k\[LetterSpace]plus8 -> 1.0, kcat1 -> 1.5, kcat2 -> 1.5, kcat3 -> 1.5, kcat4 -> 1.5, kcat5 -> 0.6, kcat6 -> 1.5, kcat7 -> 1.5, n1 -> 2.0, cytoplasm -> 1.0, nucleus -> 1.0
};

assignments = {
X13X14 -> X13[t] + X14[t]
};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {
R18->"http://identifiers.org/go/GO:0032856", R22->"http://identifiers.org/ec-code/2.7.11.1", R22->"http://identifiers.org/go/GO:0004709", R22->"http://identifiers.org/go/GO:0000186", R32->"http://identifiers.org/ec-code/2.7.11.24", R32->"http://identifiers.org/go/GO:0004707", R32->"http://identifiers.org/go/GO:0016310", R8->"http://identifiers.org/go/GO:0008013", R9->"http://identifiers.org/ec-code/2.7.11.26", R9->"http://identifiers.org/go/GO:0004674"
};

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

(* Time evolution *)
odes = {
X1'[t] == 1.0*R2 -1.0*R1, X10'[t] == 1.0*R10 -1.0*R11, X11'[t] == 1.0*R12 -1.0*R8 -1.0*R13 -1.0*R16 -1.0*R17, X12'[t] == 1.0*R14 -1.0*R7 -1.0*R15, X13'[t] ==  -1.0*R16, X14'[t] == 1.0*R16 , X15'[t] == 1.0*R17 , X16'[t] == 1.0*R19 -1.0*R18, X17'[t] == 1.0*R18 -1.0*R19, X18'[t] == 1.0*R21 +1.0*R26 -1.0*R20 -1.0*R27 -1.0*R31, X19'[t] == 1.0*R20 +1.0*R31 -1.0*R21, X2'[t] == 1.0*R1 -1.0*R2, X20'[t] == 1.0*R23 -1.0*R22, X21'[t] == 1.0*R22 -1.0*R23, X22'[t] == 1.0*R25 -1.0*R24, X23'[t] == 1.0*R24 -1.0*R25, X24'[t] == 1.0*R27 -1.0*R26, X25'[t] == 1.0*R28 -1.0*R27, X26'[t] == 1.0*R26 -1.0*R28, X27'[t] == 1.0*R29 -1.0*R30, X28'[t] == 1.0*R32 -1.0*R33, X3'[t] == 1.0*R4 +1.0*R10 -1.0*R5 -1.0*R8, X4'[t] == 1.0*R5 +1.0*R6 -1.0*R3 -1.0*R4, X5'[t] == 1.0*R3 +1.0*R33 -1.0*R6 -1.0*R32, X6'[t] == 1.0*R3 +1.0*R7 -1.0*R6, X7'[t] ==  -1.0*R7 -1.0*R17, X8'[t] == 1.0*R8 -1.0*R9, X9'[t] == 1.0*R9 -1.0*R10
};
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]}]