(* 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 = { Bc[t], Bcp[t], Bn[t], Bnp[t], Cc[t], Ccp[t], In[t], Mb[t], Mc[t], Mp[t], Mr[t], PCc[t], PCcp[t], PCn[t], PCnp[t], Pc[t], Pcp[t], Rc[t], Rn[t] }; initialValues = { Bc[0] == 0.0, Bcp[0] == 0.0, Bn[0] == 0.0, Bnp[0] == 0.0, Cc[0] == 0.0, Ccp[0] == 0.0, In[0] == 0.0, Mb[0] == 0.0, Mc[0] == 0.0, Mp[0] == 0.0, Mr[0] == 0.0, PCc[0] == 0.0, PCcp[0] == 0.0, PCn[0] == 0.0, PCnp[0] == 0.0, Pc[0] == 0.0, Pcp[0] == 0.0, Rc[0] == 0.0, Rn[0] == 0.0 }; rates = { reaction\[LetterSpace]0, 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]24, reaction\[LetterSpace]25, reaction\[LetterSpace]26, reaction\[LetterSpace]27, reaction\[LetterSpace]28, reaction\[LetterSpace]29, reaction\[LetterSpace]3, reaction\[LetterSpace]30, reaction\[LetterSpace]31, reaction\[LetterSpace]32, reaction\[LetterSpace]33, reaction\[LetterSpace]34, reaction\[LetterSpace]35, reaction\[LetterSpace]36, reaction\[LetterSpace]37, reaction\[LetterSpace]38, reaction\[LetterSpace]39, reaction\[LetterSpace]4, reaction\[LetterSpace]40, reaction\[LetterSpace]41, reaction\[LetterSpace]42, reaction\[LetterSpace]43, reaction\[LetterSpace]44, reaction\[LetterSpace]45, reaction\[LetterSpace]46, reaction\[LetterSpace]47, reaction\[LetterSpace]48, reaction\[LetterSpace]49, reaction\[LetterSpace]5, reaction\[LetterSpace]50, reaction\[LetterSpace]51, reaction\[LetterSpace]52, reaction\[LetterSpace]53, reaction\[LetterSpace]54, reaction\[LetterSpace]55, reaction\[LetterSpace]56, reaction\[LetterSpace]6, reaction\[LetterSpace]7, reaction\[LetterSpace]8, reaction\[LetterSpace]9 }; rateEquations = { reaction\[LetterSpace]0 -> cell*function\[LetterSpace]0[reaction\[LetterSpace]0\[LetterSpace]vsb, reaction\[LetterSpace]0\[LetterSpace]K, reaction\[LetterSpace]0\[LetterSpace]m, Rn[t]], reaction\[LetterSpace]1 -> cell*function\[LetterSpace]1[reaction\[LetterSpace]1\[LetterSpace]k, Mb[t]], reaction\[LetterSpace]10 -> cell*function\[LetterSpace]1[reaction\[LetterSpace]10\[LetterSpace]k, Mp[t]], reaction\[LetterSpace]11 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]11\[LetterSpace]V, Pcp[t], reaction\[LetterSpace]11\[LetterSpace]Km], reaction\[LetterSpace]12 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]12\[LetterSpace]V, Pc[t], reaction\[LetterSpace]12\[LetterSpace]Km], reaction\[LetterSpace]13 -> cell*(reaction\[LetterSpace]13\[LetterSpace]k1*Cc[t]*Pc[t] - reaction\[LetterSpace]13\[LetterSpace]k2*PCc[t]), reaction\[LetterSpace]14 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]14\[LetterSpace]V, PCc[t], reaction\[LetterSpace]14\[LetterSpace]Km], reaction\[LetterSpace]15 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]15\[LetterSpace]V, PCcp[t], reaction\[LetterSpace]15\[LetterSpace]Km], reaction\[LetterSpace]16 -> cell*(reaction\[LetterSpace]16\[LetterSpace]k1*PCc[t] - reaction\[LetterSpace]16\[LetterSpace]k2*PCn[t]), reaction\[LetterSpace]17 -> cell*reaction\[LetterSpace]17\[LetterSpace]k1*PCnp[t], reaction\[LetterSpace]18 -> cell*reaction\[LetterSpace]18\[LetterSpace]k1*Bcp[t], reaction\[LetterSpace]19 -> cell*reaction\[LetterSpace]19\[LetterSpace]k1*Bnp[t], reaction\[LetterSpace]2 -> cell*reaction\[LetterSpace]2\[LetterSpace]k1*Mb[t], reaction\[LetterSpace]20 -> cell*function\[LetterSpace]3[reaction\[LetterSpace]20\[LetterSpace]Vs, Bn[t], reaction\[LetterSpace]20\[LetterSpace]n, reaction\[LetterSpace]20\[LetterSpace]K], reaction\[LetterSpace]21 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]21\[LetterSpace]V, PCn[t], reaction\[LetterSpace]21\[LetterSpace]Km], reaction\[LetterSpace]22 -> cell*reaction\[LetterSpace]22\[LetterSpace]k1*Mp[t], reaction\[LetterSpace]23 -> cell*(-(reaction\[LetterSpace]23\[LetterSpace]k2*In[t]) + reaction\[LetterSpace]23\[LetterSpace]k1*Bn[t]*PCn[t]), reaction\[LetterSpace]24 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]24\[LetterSpace]V, Mb[t], reaction\[LetterSpace]24\[LetterSpace]Km], reaction\[LetterSpace]25 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]25\[LetterSpace]V, Mc[t], reaction\[LetterSpace]25\[LetterSpace]Km], reaction\[LetterSpace]26 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]26\[LetterSpace]V, Mp[t], reaction\[LetterSpace]26\[LetterSpace]Km], reaction\[LetterSpace]27 -> cell*reaction\[LetterSpace]27\[LetterSpace]k1*Pc[t], reaction\[LetterSpace]28 -> cell*reaction\[LetterSpace]28\[LetterSpace]k1*Cc[t], reaction\[LetterSpace]29 -> cell*reaction\[LetterSpace]29\[LetterSpace]k1*Pcp[t], reaction\[LetterSpace]3 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]3\[LetterSpace]V, Bc[t], reaction\[LetterSpace]3\[LetterSpace]Km], reaction\[LetterSpace]30 -> cell*reaction\[LetterSpace]30\[LetterSpace]k1*Ccp[t], reaction\[LetterSpace]31 -> cell*reaction\[LetterSpace]31\[LetterSpace]k1*PCcp[t], reaction\[LetterSpace]32 -> cell*reaction\[LetterSpace]32\[LetterSpace]k1*PCc[t], reaction\[LetterSpace]33 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]33\[LetterSpace]V, PCnp[t], reaction\[LetterSpace]33\[LetterSpace]Km], reaction\[LetterSpace]34 -> cell*reaction\[LetterSpace]34\[LetterSpace]k1*Bc[t], reaction\[LetterSpace]35 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]35\[LetterSpace]V, Bcp[t], reaction\[LetterSpace]35\[LetterSpace]Km], reaction\[LetterSpace]36 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]36\[LetterSpace]V, Bn[t], reaction\[LetterSpace]36\[LetterSpace]Km], reaction\[LetterSpace]37 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]37\[LetterSpace]V, Bnp[t], reaction\[LetterSpace]37\[LetterSpace]Km], reaction\[LetterSpace]38 -> cell*reaction\[LetterSpace]38\[LetterSpace]k1*In[t], reaction\[LetterSpace]39 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]39\[LetterSpace]V, In[t], reaction\[LetterSpace]39\[LetterSpace]Km], reaction\[LetterSpace]4 -> cell*(reaction\[LetterSpace]4\[LetterSpace]k1*Bc[t] - reaction\[LetterSpace]4\[LetterSpace]k2*Bn[t]), reaction\[LetterSpace]40 -> cell*reaction\[LetterSpace]40\[LetterSpace]k1*Bn[t], reaction\[LetterSpace]41 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]41\[LetterSpace]V, Bcp[t], reaction\[LetterSpace]41\[LetterSpace]Km], reaction\[LetterSpace]42 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]42\[LetterSpace]V, Bnp[t], reaction\[LetterSpace]42\[LetterSpace]Km], reaction\[LetterSpace]43 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]43\[LetterSpace]V, Ccp[t], reaction\[LetterSpace]43\[LetterSpace]Km], reaction\[LetterSpace]44 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]44\[LetterSpace]V, Pcp[t], reaction\[LetterSpace]44\[LetterSpace]Km], reaction\[LetterSpace]45 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]45\[LetterSpace]V, PCnp[t], reaction\[LetterSpace]45\[LetterSpace]Km], reaction\[LetterSpace]46 -> cell*reaction\[LetterSpace]46\[LetterSpace]k1*PCn[t], reaction\[LetterSpace]47 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]47\[LetterSpace]V, PCcp[t], reaction\[LetterSpace]47\[LetterSpace]Km], reaction\[LetterSpace]48 -> cell*function\[LetterSpace]3[reaction\[LetterSpace]48\[LetterSpace]Vs, Bn[t], reaction\[LetterSpace]48\[LetterSpace]n, reaction\[LetterSpace]48\[LetterSpace]K], reaction\[LetterSpace]49 -> cell*reaction\[LetterSpace]49\[LetterSpace]k1*Mr[t], reaction\[LetterSpace]5 -> cell*function\[LetterSpace]1[reaction\[LetterSpace]5\[LetterSpace]k, Mc[t]], reaction\[LetterSpace]50 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]50\[LetterSpace]V, Mr[t], reaction\[LetterSpace]50\[LetterSpace]Km], reaction\[LetterSpace]51 -> cell*function\[LetterSpace]1[reaction\[LetterSpace]51\[LetterSpace]k, Mr[t]], reaction\[LetterSpace]52 -> cell*(reaction\[LetterSpace]52\[LetterSpace]k1*Rc[t] - reaction\[LetterSpace]52\[LetterSpace]k2*Rn[t]), reaction\[LetterSpace]53 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]53\[LetterSpace]V, Rc[t], reaction\[LetterSpace]53\[LetterSpace]Km], reaction\[LetterSpace]54 -> cell*reaction\[LetterSpace]54\[LetterSpace]k1*Rc[t], reaction\[LetterSpace]55 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]55\[LetterSpace]V, Rn[t], reaction\[LetterSpace]55\[LetterSpace]Km], reaction\[LetterSpace]56 -> cell*reaction\[LetterSpace]56\[LetterSpace]k1*Rn[t], reaction\[LetterSpace]6 -> cell*reaction\[LetterSpace]6\[LetterSpace]k1*Mc[t], reaction\[LetterSpace]7 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]7\[LetterSpace]V, Cc[t], reaction\[LetterSpace]7\[LetterSpace]Km], reaction\[LetterSpace]8 -> cell*function\[LetterSpace]2[reaction\[LetterSpace]8\[LetterSpace]V, Ccp[t], reaction\[LetterSpace]8\[LetterSpace]Km], reaction\[LetterSpace]9 -> cell*function\[LetterSpace]3[parameter\[LetterSpace]0000097, Bn[t], reaction\[LetterSpace]9\[LetterSpace]n, reaction\[LetterSpace]9\[LetterSpace]K] }; parameters = { reaction\[LetterSpace]0\[LetterSpace]vsb -> 1.8, reaction\[LetterSpace]0\[LetterSpace]K -> 1.0, reaction\[LetterSpace]0\[LetterSpace]m -> 2.0, reaction\[LetterSpace]1\[LetterSpace]k -> 0.32, reaction\[LetterSpace]2\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]3\[LetterSpace]V -> 1.4, reaction\[LetterSpace]3\[LetterSpace]Km -> 1.006, reaction\[LetterSpace]4\[LetterSpace]k1 -> 0.8, reaction\[LetterSpace]4\[LetterSpace]k2 -> 0.4, reaction\[LetterSpace]5\[LetterSpace]k -> 3.2, reaction\[LetterSpace]6\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]7\[LetterSpace]V -> 1.2, reaction\[LetterSpace]7\[LetterSpace]Km -> 1.006, reaction\[LetterSpace]8\[LetterSpace]V -> 1.4, reaction\[LetterSpace]8\[LetterSpace]Km -> 0.3, reaction\[LetterSpace]9\[LetterSpace]n -> 2.0, reaction\[LetterSpace]9\[LetterSpace]K -> 0.6, reaction\[LetterSpace]10\[LetterSpace]k -> 1.2, reaction\[LetterSpace]11\[LetterSpace]V -> 3.4, reaction\[LetterSpace]11\[LetterSpace]Km -> 0.3, reaction\[LetterSpace]12\[LetterSpace]V -> 9.6, reaction\[LetterSpace]12\[LetterSpace]Km -> 1.006, reaction\[LetterSpace]13\[LetterSpace]k1 -> 0.8, reaction\[LetterSpace]13\[LetterSpace]k2 -> 0.4, reaction\[LetterSpace]14\[LetterSpace]V -> 2.4, reaction\[LetterSpace]14\[LetterSpace]Km -> 1.006, reaction\[LetterSpace]15\[LetterSpace]V -> 1.4, reaction\[LetterSpace]15\[LetterSpace]Km -> 0.3, reaction\[LetterSpace]16\[LetterSpace]k1 -> 0.8, reaction\[LetterSpace]16\[LetterSpace]k2 -> 0.4, reaction\[LetterSpace]17\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]18\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]19\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]20\[LetterSpace]Vs -> 2.2, reaction\[LetterSpace]20\[LetterSpace]n -> 2.0, reaction\[LetterSpace]20\[LetterSpace]K -> 0.6, reaction\[LetterSpace]21\[LetterSpace]V -> 2.4, reaction\[LetterSpace]21\[LetterSpace]Km -> 1.006, reaction\[LetterSpace]22\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]23\[LetterSpace]k1 -> 1.0, reaction\[LetterSpace]23\[LetterSpace]k2 -> 0.2, reaction\[LetterSpace]24\[LetterSpace]V -> 1.3, reaction\[LetterSpace]24\[LetterSpace]Km -> 0.4, reaction\[LetterSpace]25\[LetterSpace]V -> 2.0, reaction\[LetterSpace]25\[LetterSpace]Km -> 0.4, reaction\[LetterSpace]26\[LetterSpace]V -> 2.2, reaction\[LetterSpace]26\[LetterSpace]Km -> 0.3, reaction\[LetterSpace]27\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]28\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]29\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]40\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]41\[LetterSpace]V -> 0.2, reaction\[LetterSpace]41\[LetterSpace]Km -> 0.1, reaction\[LetterSpace]42\[LetterSpace]V -> 0.4, reaction\[LetterSpace]42\[LetterSpace]Km -> 0.1, reaction\[LetterSpace]43\[LetterSpace]V -> 0.2, reaction\[LetterSpace]43\[LetterSpace]Km -> 0.1, reaction\[LetterSpace]44\[LetterSpace]V -> 0.6, reaction\[LetterSpace]44\[LetterSpace]Km -> 0.1, reaction\[LetterSpace]45\[LetterSpace]V -> 0.2, reaction\[LetterSpace]45\[LetterSpace]Km -> 0.1, reaction\[LetterSpace]46\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]47\[LetterSpace]V -> 0.2, reaction\[LetterSpace]47\[LetterSpace]Km -> 0.1, reaction\[LetterSpace]48\[LetterSpace]Vs -> 1.6, reaction\[LetterSpace]48\[LetterSpace]n -> 2.0, reaction\[LetterSpace]48\[LetterSpace]K -> 0.6, reaction\[LetterSpace]49\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]50\[LetterSpace]V -> 1.6, reaction\[LetterSpace]50\[LetterSpace]Km -> 0.4, reaction\[LetterSpace]51\[LetterSpace]k -> 1.7, reaction\[LetterSpace]30\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]31\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]32\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]33\[LetterSpace]V -> 1.4, reaction\[LetterSpace]33\[LetterSpace]Km -> 0.3, reaction\[LetterSpace]34\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]35\[LetterSpace]V -> 3.0, reaction\[LetterSpace]35\[LetterSpace]Km -> 0.3, reaction\[LetterSpace]36\[LetterSpace]V -> 1.4, reaction\[LetterSpace]36\[LetterSpace]Km -> 1.006, reaction\[LetterSpace]37\[LetterSpace]V -> 3.0, reaction\[LetterSpace]37\[LetterSpace]Km -> 0.3, reaction\[LetterSpace]38\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]39\[LetterSpace]V -> 1.6, reaction\[LetterSpace]39\[LetterSpace]Km -> 0.3, reaction\[LetterSpace]52\[LetterSpace]k1 -> 0.8, reaction\[LetterSpace]52\[LetterSpace]k2 -> 0.4, reaction\[LetterSpace]53\[LetterSpace]V -> 4.4, reaction\[LetterSpace]53\[LetterSpace]Km -> 0.3, reaction\[LetterSpace]54\[LetterSpace]k1 -> 0.02, reaction\[LetterSpace]55\[LetterSpace]V -> 0.8, reaction\[LetterSpace]55\[LetterSpace]Km -> 0.3, reaction\[LetterSpace]56\[LetterSpace]k1 -> 0.02, cell -> 1.0 }; assignments = { function\[LetterSpace]2[V_,substrate_,Km_] -> (substrate*V)/(Km + substrate), function\[LetterSpace]3[Vs_,B_,n_,K_] -> (B^n*Vs)/(B^n + K^n), function\[LetterSpace]5[length_,tt_] -> Ceiling[Sin[0.001 + (Pi*tt)/length]/2], function\[LetterSpace]0[vsb_,K_,m_,Bn_] -> (K^m*vsb)/(Bn^m + K^m), function\[LetterSpace]1[k_,mRNA_] -> k*mRNA, parameter\[LetterSpace]0000097 -> 2.4 + 0.6000000000000001*parameter\[LetterSpace]0000096, parameter\[LetterSpace]0000096 -> function\[LetterSpace]5[12, t] }; events = { }; speciesAnnotations = { }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { Bc'[t] == 1.0*reaction\[LetterSpace]1 +1.0*reaction\[LetterSpace]41 -1.0*reaction\[LetterSpace]3 -1.0*reaction\[LetterSpace]4 -1.0*reaction\[LetterSpace]34, Bcp'[t] == 1.0*reaction\[LetterSpace]3 -1.0*reaction\[LetterSpace]18 -1.0*reaction\[LetterSpace]35 -1.0*reaction\[LetterSpace]41, Bn'[t] == 1.0*reaction\[LetterSpace]4 +1.0*reaction\[LetterSpace]42 -1.0*reaction\[LetterSpace]23 -1.0*reaction\[LetterSpace]36 -1.0*reaction\[LetterSpace]40, Bnp'[t] == 1.0*reaction\[LetterSpace]36 -1.0*reaction\[LetterSpace]19 -1.0*reaction\[LetterSpace]37 -1.0*reaction\[LetterSpace]42, Cc'[t] == 1.0*reaction\[LetterSpace]5 +1.0*reaction\[LetterSpace]43 -1.0*reaction\[LetterSpace]7 -1.0*reaction\[LetterSpace]13 -1.0*reaction\[LetterSpace]28, Ccp'[t] == 1.0*reaction\[LetterSpace]7 -1.0*reaction\[LetterSpace]8 -1.0*reaction\[LetterSpace]30 -1.0*reaction\[LetterSpace]43, In'[t] == 1.0*reaction\[LetterSpace]23 -1.0*reaction\[LetterSpace]38 -1.0*reaction\[LetterSpace]39, Mb'[t] == 1.0*reaction\[LetterSpace]0 -1.0*reaction\[LetterSpace]2 -1.0*reaction\[LetterSpace]24, Mc'[t] == 1.0*reaction\[LetterSpace]20 -1.0*reaction\[LetterSpace]6 -1.0*reaction\[LetterSpace]25, Mp'[t] == 1.0*reaction\[LetterSpace]9 -1.0*reaction\[LetterSpace]22 -1.0*reaction\[LetterSpace]26, Mr'[t] == 1.0*reaction\[LetterSpace]48 -1.0*reaction\[LetterSpace]49 -1.0*reaction\[LetterSpace]50, PCc'[t] == 1.0*reaction\[LetterSpace]13 +1.0*reaction\[LetterSpace]47 -1.0*reaction\[LetterSpace]14 -1.0*reaction\[LetterSpace]16 -1.0*reaction\[LetterSpace]32, PCcp'[t] == 1.0*reaction\[LetterSpace]14 -1.0*reaction\[LetterSpace]15 -1.0*reaction\[LetterSpace]31 -1.0*reaction\[LetterSpace]47, PCn'[t] == 1.0*reaction\[LetterSpace]16 +1.0*reaction\[LetterSpace]45 -1.0*reaction\[LetterSpace]21 -1.0*reaction\[LetterSpace]23 -1.0*reaction\[LetterSpace]46, PCnp'[t] == 1.0*reaction\[LetterSpace]21 -1.0*reaction\[LetterSpace]17 -1.0*reaction\[LetterSpace]33 -1.0*reaction\[LetterSpace]45, Pc'[t] == 1.0*reaction\[LetterSpace]10 +1.0*reaction\[LetterSpace]44 -1.0*reaction\[LetterSpace]12 -1.0*reaction\[LetterSpace]13 -1.0*reaction\[LetterSpace]27, Pcp'[t] == 1.0*reaction\[LetterSpace]12 -1.0*reaction\[LetterSpace]11 -1.0*reaction\[LetterSpace]29 -1.0*reaction\[LetterSpace]44, Rc'[t] == 1.0*reaction\[LetterSpace]51 -1.0*reaction\[LetterSpace]52 -1.0*reaction\[LetterSpace]53 -1.0*reaction\[LetterSpace]54, Rn'[t] == 1.0*reaction\[LetterSpace]52 -1.0*reaction\[LetterSpace]55 -1.0*reaction\[LetterSpace]56 }; 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]}]