(* 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 = {
C3[t], C3\[LetterSpace]star[t], C8[t], C8\[LetterSpace]star[t], FL[t], FS[t], IKK[t], L[t], L\[LetterSpace]RF[t], L\[LetterSpace]RF\[LetterSpace]C8[t], L\[LetterSpace]RF\[LetterSpace]C8\[LetterSpace]FS[t], L\[LetterSpace]RF\[LetterSpace]FL[t], L\[LetterSpace]RF\[LetterSpace]FL\[LetterSpace]FL[t], L\[LetterSpace]RF\[LetterSpace]FL\[LetterSpace]FS[t], L\[LetterSpace]RF\[LetterSpace]FS[t], L\[LetterSpace]RF\[LetterSpace]FS\[LetterSpace]FS[t], NF\[LetterSpace]kB\[LetterSpace]IkB[t], NF\[LetterSpace]kB\[LetterSpace]IkB\[LetterSpace]P[t], NF\[LetterSpace]kB\[LetterSpace]star[t], RF[t], p43\[LetterSpace]FLIP[t], p43\[LetterSpace]FLIP\[LetterSpace]IKK\[LetterSpace]star[t], p43\[LetterSpace]p41[t]
};

initialValues = {
C3[0] == 1.443404, C3\[LetterSpace]star[0] == 0.0, C8[0] == 64.47652, C8\[LetterSpace]star[0] == 0.0, FL[0] == 7.398562, FS[0] == 5.083923, IKK[0] == 5.772825, L[0] == 113.22, L\[LetterSpace]RF[0] == 0.0, L\[LetterSpace]RF\[LetterSpace]C8[0] == 0.0, L\[LetterSpace]RF\[LetterSpace]C8\[LetterSpace]FS[0] == 0.0, L\[LetterSpace]RF\[LetterSpace]FL[0] == 0.0, L\[LetterSpace]RF\[LetterSpace]FL\[LetterSpace]FL[0] == 0.0, L\[LetterSpace]RF\[LetterSpace]FL\[LetterSpace]FS[0] == 0.0, L\[LetterSpace]RF\[LetterSpace]FS[0] == 0.0, L\[LetterSpace]RF\[LetterSpace]FS\[LetterSpace]FS[0] == 0.0, NF\[LetterSpace]kB\[LetterSpace]IkB[0] == 4.739546, NF\[LetterSpace]kB\[LetterSpace]IkB\[LetterSpace]P[0] == 0.0, NF\[LetterSpace]kB\[LetterSpace]star[0] == 0.0, RF[0] == 91.26592, p43\[LetterSpace]FLIP[0] == 0.0, p43\[LetterSpace]FLIP\[LetterSpace]IKK\[LetterSpace]star[0] == 0.0, p43\[LetterSpace]p41[0] == 0.0
};

rates = {
reaction\[LetterSpace]1, reaction\[LetterSpace]10, reaction\[LetterSpace]11, reaction\[LetterSpace]12, reaction\[LetterSpace]13, reaction\[LetterSpace]14, reaction\[LetterSpace]15, reaction\[LetterSpace]16, reaction\[LetterSpace]17, reaction\[LetterSpace]18, reaction\[LetterSpace]19, reaction\[LetterSpace]2, reaction\[LetterSpace]20, reaction\[LetterSpace]21, reaction\[LetterSpace]22, reaction\[LetterSpace]23, 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 -> default*k1*L[t]*RF[t], reaction\[LetterSpace]10 -> default*k7*FS[t]*L\[LetterSpace]RF\[LetterSpace]FL[t], reaction\[LetterSpace]11 -> default*k5*C8[t]*L\[LetterSpace]RF\[LetterSpace]FS[t], reaction\[LetterSpace]12 -> default*k6*FL[t]*L\[LetterSpace]RF\[LetterSpace]FS[t], reaction\[LetterSpace]13 -> default*k7*FS[t]*L\[LetterSpace]RF\[LetterSpace]FS[t], reaction\[LetterSpace]14 -> default*k8*p43\[LetterSpace]p41[t]^2, reaction\[LetterSpace]15 -> default*k9*C3[t]*C8\[LetterSpace]star[t], reaction\[LetterSpace]16 -> default*k10*C3\[LetterSpace]star[t]*C8[t], reaction\[LetterSpace]17 -> default*k11*C8\[LetterSpace]star[t], reaction\[LetterSpace]18 -> default*k12*C3\[LetterSpace]star[t], reaction\[LetterSpace]19 -> default*k13*IKK[t]*p43\[LetterSpace]FLIP[t], reaction\[LetterSpace]2 -> default*k2*C8[t]*L\[LetterSpace]RF[t], reaction\[LetterSpace]20 -> default*k14*NF\[LetterSpace]kB\[LetterSpace]IkB[t]*p43\[LetterSpace]FLIP\[LetterSpace]IKK\[LetterSpace]star[t], reaction\[LetterSpace]21 -> default*k15*NF\[LetterSpace]kB\[LetterSpace]IkB\[LetterSpace]P[t], reaction\[LetterSpace]22 -> default*k16*p43\[LetterSpace]FLIP\[LetterSpace]IKK\[LetterSpace]star[t], reaction\[LetterSpace]23 -> default*k17*NF\[LetterSpace]kB\[LetterSpace]star[t], reaction\[LetterSpace]3 -> default*k3*FL[t]*L\[LetterSpace]RF[t], reaction\[LetterSpace]4 -> default*k4*FS[t]*L\[LetterSpace]RF[t], reaction\[LetterSpace]5 -> default*k5*C8[t]*L\[LetterSpace]RF\[LetterSpace]C8[t], reaction\[LetterSpace]6 -> default*k6*FL[t]*L\[LetterSpace]RF\[LetterSpace]C8[t], reaction\[LetterSpace]7 -> default*k7*FS[t]*L\[LetterSpace]RF\[LetterSpace]C8[t], reaction\[LetterSpace]8 -> default*k5*C8[t]*L\[LetterSpace]RF\[LetterSpace]FL[t], reaction\[LetterSpace]9 -> default*k6*FL[t]*L\[LetterSpace]RF\[LetterSpace]FL[t]
};

parameters = {
k1 -> 1.0, k10 -> 0.1205258, k11 -> 0.02891451, k12 -> 0.1502914, k13 -> 0.0007204261, k14 -> 0.3588224, k15 -> 3.684162, k16 -> 0.02229912, k17 -> 0.0064182, k2 -> 0.0001277248, k3 -> 0.6693316, k4 -> 1*^-05, k5 -> 0.0005946569, k6 -> 0.9999999, k7 -> 0.8875063, k8 -> 0.0008044378, k9 -> 0.002249759, default -> 1.0
};

assignments = {

};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {

};

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

(* Time evolution *)
odes = {
C3'[t] ==  -1.0*reaction\[LetterSpace]15, C3\[LetterSpace]star'[t] == 1.0*reaction\[LetterSpace]15 +1.0*reaction\[LetterSpace]16 -1.0*reaction\[LetterSpace]16 -1.0*reaction\[LetterSpace]18, C8'[t] ==  -1.0*reaction\[LetterSpace]2 -1.0*reaction\[LetterSpace]5 -1.0*reaction\[LetterSpace]8 -1.0*reaction\[LetterSpace]11 -1.0*reaction\[LetterSpace]16, C8\[LetterSpace]star'[t] == 1.0*reaction\[LetterSpace]14 +1.0*reaction\[LetterSpace]15 -1.0*reaction\[LetterSpace]15 -1.0*reaction\[LetterSpace]17, FL'[t] ==  -1.0*reaction\[LetterSpace]3 -1.0*reaction\[LetterSpace]6 -1.0*reaction\[LetterSpace]9 -1.0*reaction\[LetterSpace]12, FS'[t] ==  -1.0*reaction\[LetterSpace]4 -1.0*reaction\[LetterSpace]7 -1.0*reaction\[LetterSpace]10 -1.0*reaction\[LetterSpace]13, IKK'[t] ==  -1.0*reaction\[LetterSpace]19, L'[t] ==  -1.0*reaction\[LetterSpace]1, L\[LetterSpace]RF'[t] == 1.0*reaction\[LetterSpace]1 -1.0*reaction\[LetterSpace]2 -1.0*reaction\[LetterSpace]3 -1.0*reaction\[LetterSpace]4, L\[LetterSpace]RF\[LetterSpace]C8'[t] == 1.0*reaction\[LetterSpace]2 -1.0*reaction\[LetterSpace]5 -1.0*reaction\[LetterSpace]6 -1.0*reaction\[LetterSpace]7, L\[LetterSpace]RF\[LetterSpace]C8\[LetterSpace]FS'[t] == 1.0*reaction\[LetterSpace]7 +1.0*reaction\[LetterSpace]11 , L\[LetterSpace]RF\[LetterSpace]FL'[t] == 1.0*reaction\[LetterSpace]3 -1.0*reaction\[LetterSpace]8 -1.0*reaction\[LetterSpace]9 -1.0*reaction\[LetterSpace]10, L\[LetterSpace]RF\[LetterSpace]FL\[LetterSpace]FL'[t] == 1.0*reaction\[LetterSpace]9 , L\[LetterSpace]RF\[LetterSpace]FL\[LetterSpace]FS'[t] == 1.0*reaction\[LetterSpace]10 +1.0*reaction\[LetterSpace]12 , L\[LetterSpace]RF\[LetterSpace]FS'[t] == 1.0*reaction\[LetterSpace]4 -1.0*reaction\[LetterSpace]11 -1.0*reaction\[LetterSpace]12 -1.0*reaction\[LetterSpace]13, L\[LetterSpace]RF\[LetterSpace]FS\[LetterSpace]FS'[t] == 1.0*reaction\[LetterSpace]13 , NF\[LetterSpace]kB\[LetterSpace]IkB'[t] ==  -1.0*reaction\[LetterSpace]20, NF\[LetterSpace]kB\[LetterSpace]IkB\[LetterSpace]P'[t] == 1.0*reaction\[LetterSpace]20 -1.0*reaction\[LetterSpace]21, NF\[LetterSpace]kB\[LetterSpace]star'[t] == 1.0*reaction\[LetterSpace]21 -1.0*reaction\[LetterSpace]23, RF'[t] ==  -1.0*reaction\[LetterSpace]1, p43\[LetterSpace]FLIP'[t] == 1.0*reaction\[LetterSpace]6 +1.0*reaction\[LetterSpace]8 -1.0*reaction\[LetterSpace]19, p43\[LetterSpace]FLIP\[LetterSpace]IKK\[LetterSpace]star'[t] == 1.0*reaction\[LetterSpace]19 +1.0*reaction\[LetterSpace]20 -1.0*reaction\[LetterSpace]20 -1.0*reaction\[LetterSpace]22, p43\[LetterSpace]p41'[t] == 1.0*reaction\[LetterSpace]5 +1.0*reaction\[LetterSpace]5 +1.0*reaction\[LetterSpace]16 -1.0*reaction\[LetterSpace]14 -1.0*reaction\[LetterSpace]14
};
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]}]