(* 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], AC9[t], AC9X[t], AC9\[LetterSpace]star[t], AC9\[LetterSpace]starX[t], C3[t], C3\[LetterSpace]star[t], C3\[LetterSpace]starX[t], C9[t], C9X[t], C9\[LetterSpace]star[t], C9\[LetterSpace]starX[t], X[t]
};

initialValues = {
A[0] == 20.0, AC9[0] == 0.0, AC9X[0] == 0.0, AC9\[LetterSpace]star[0] == 0.0, AC9\[LetterSpace]starX[0] == 0.0, C3[0] == 200.0, C3\[LetterSpace]star[0] == 0.0, C3\[LetterSpace]starX[0] == 0.0, C9[0] == 20.0, C9X[0] == 0.0, C9\[LetterSpace]star[0] == 0.0, C9\[LetterSpace]starX[0] == 0.0, X[0] == 40.0
};

rates = {
v1, v10, v11, v12, v13, v14, v15, v16, v17, v18, v19, v2, v20, v21, v22, v23, v24, v25, v26, v27, v28, v3, v4, v5, v6, v7, v8, v9
};

rateEquations = {
v1 -> cytosol*(-(kb1*AC9[t]) + k1*A[t]*C9[t]), v10 -> cytosol*(-(k10b*AC9X[t]) + k10*AC9[t]*X[t]), v11 -> cytosol*(-(k11b*C9\[LetterSpace]starX[t]) + k11*C9\[LetterSpace]star[t]*X[t]), v12 -> cytosol*(-(k12b*AC9\[LetterSpace]starX[t]) + k12*AC9\[LetterSpace]star[t]*X[t]), v13 -> cytosol*(-(k13b*AC9X[t]) + k13*A[t]*C9X[t]), v14 -> cytosol*(-(k14b*AC9\[LetterSpace]starX[t]) + k14*A[t]*C9\[LetterSpace]starX[t]), v15 -> cytosol*(-(k15b*C3\[LetterSpace]starX[t]) + k15*C3\[LetterSpace]star[t]*X[t]), v16 -> cytosol*(k16prod - k16*A[t]), v17 -> cytosol*(k17prod - k17*C9[t]), v18 -> cytosol*(k18prod - k18*X[t]), v19 -> cytosol*k19*C9X[t], v2 -> cytosol*k2*C3[t]*C9[t], v20 -> cytosol*k20*AC9X[t], v21 -> cytosol*k21*AC9[t], v22 -> cytosol*(k22prod - k22*C3[t]), v23 -> cytosol*k23*C3\[LetterSpace]star[t], v24 -> cytosol*k24*C3\[LetterSpace]starX[t], v25 -> cytosol*k25*C9\[LetterSpace]starX[t], v26 -> cytosol*k26*C9\[LetterSpace]star[t], v27 -> cytosol*k27*AC9\[LetterSpace]star[t], v28 -> cytosol*k28*AC9\[LetterSpace]starX[t], v3 -> cytosol*k3*AC9[t]*C3[t], v4 -> cytosol*k4*C3\[LetterSpace]star[t]*C9[t], v5 -> cytosol*k5*AC9[t]*C3\[LetterSpace]star[t], v6 -> cytosol*k6*C3[t]*C9\[LetterSpace]star[t], v7 -> cytosol*k7*AC9\[LetterSpace]star[t]*C3[t], v8 -> cytosol*(-(k8b*AC9\[LetterSpace]star[t]) + k8*A[t]*C9\[LetterSpace]star[t]), v9 -> cytosol*(-(k9b*C9X[t]) + k9*C9[t]*X[t])
};

parameters = {
k1 -> 0.002, k10 -> 0.001, k10b -> 0.001, k11 -> 0.001, k11b -> 0.001, k12 -> 0.001, k12b -> 0.001, k13 -> 0.002, k13b -> 0.1, k14 -> 0.002, k14b -> 0.1, k15 -> 0.003, k15b -> 0.001, k16 -> 0.001, k16prod -> 0.02, k17 -> 0.001, k17prod -> 0.02, k18 -> 0.001, k18prod -> 0.04, k19 -> 0.001, k2 -> 5*^-06, k20 -> 0.001, k21 -> 0.001, k22 -> 0.001, k22prod -> 0.2, k23 -> 0.001, k24 -> 0.001, k25 -> 0.001, k26 -> 0.001, k27 -> 0.001, k28 -> 0.001, k3 -> 0.00035, k4 -> 0.0002, k5 -> 0.0002, k6 -> 5*^-05, k7 -> 0.0035, k8 -> 0.002, k8b -> 0.1, k9 -> 0.001, k9b -> 0.001, kb1 -> 0.1, cytosol -> 1.0
};

assignments = {

};

events = {

};

speciesAnnotations = {

};

reactionAnnotations = {
v4->"http://identifiers.org/ec-code/3.4.22.56", v4->"http://identifiers.org/go/GO:0030693", v5->"http://identifiers.org/ec-code/3.4.22.56", v5->"http://identifiers.org/go/GO:0030693", v6->"http://identifiers.org/ec-code/3.4.22.62", v6->"http://identifiers.org/go/GO:0030693", v7->"http://identifiers.org/ec-code/3.4.22.62", v7->"http://identifiers.org/go/GO:0030693"
};

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

(* Time evolution *)
odes = {
A'[t] == 1.0*v16 -1.0*v1 -1.0*v13 -1.0*v8 -1.0*v14, AC9'[t] == 1.0*v1 +1.0*v3 -1.0*v3 -1.0*v10 -1.0*v5 -1.0*v21, AC9X'[t] == 1.0*v10 +1.0*v13 -1.0*v20, AC9\[LetterSpace]star'[t] == 1.0*v5 +1.0*v8 +1.0*v7 -1.0*v12 -1.0*v7 -1.0*v27, AC9\[LetterSpace]starX'[t] == 1.0*v12 +1.0*v14 -1.0*v28, C3'[t] == 1.0*v22 -1.0*v2 -1.0*v3 -1.0*v6 -1.0*v7, C3\[LetterSpace]star'[t] == 1.0*v2 +1.0*v3 +1.0*v4 +1.0*v5 +1.0*v6 +1.0*v7 -1.0*v15 -1.0*v4 -1.0*v5 -1.0*v23, C3\[LetterSpace]starX'[t] == 1.0*v15 -1.0*v24, C9'[t] == 1.0*v2 +1.0*v17 -1.0*v1 -1.0*v2 -1.0*v9 -1.0*v4, C9X'[t] == 1.0*v9 -1.0*v13 -1.0*v19, C9\[LetterSpace]star'[t] == 1.0*v4 +1.0*v6 -1.0*v8 -1.0*v11 -1.0*v6 -1.0*v26, C9\[LetterSpace]starX'[t] == 1.0*v11 -1.0*v14 -1.0*v25, X'[t] == 1.0*v18 -1.0*v9 -1.0*v10 -1.0*v15 -1.0*v11 -1.0*v12
};
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]}]