(* 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 = { 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], x29[t], x3[t], x30[t], x31[t], x32[t], x33[t], x34[t], x35[t], x36[t], x37[t], x38[t], x39[t], x4[t], x40[t], x41[t], x42[t], x43[t], x44[t], x45[t], x46[t], x47[t], x48[t], x49[t], x5[t], x50[t], x51[t], x52[t], x53[t], x54[t], x55[t], x56[t], x57[t], x58[t], x59[t], x6[t], x60[t], x61[t], x62[t], x63[t], x64[t], x65[t], x66[t], x67[t], x68[t], x69[t], x7[t], x70[t], x71[t], x72[t], x73[t], x74[t], x75[t], x76[t], x77[t], x78[t], x79[t], x8[t], x80[t], x81[t], x82[t], x83[t], x84[t], x85[t], x86[t], x87[t], x88[t], x89[t], x9[t], x90[t], x91[t], x92[t], x93[t], x94[t] }; initialValues = { x10[0] == 0.0, x11[0] == 0.0, x12[0] == 81000.0, x13[0] == 0.0, x14[0] == 12000.0, x15[0] == 0.0, x16[0] == 0.0, x17[0] == 0.0, x18[0] == 0.0, x19[0] == 0.0, x2[0] == 50000.0, x20[0] == 0.0, x21[0] == 0.0, x22[0] == 11000.0, x23[0] == 0.0, x24[0] == 26300.0, x25[0] == 0.0, x26[0] == 72000.0, x27[0] == 0.0, x28[0] == 0.0, x29[0] == 0.0, x3[0] == 0.0, x30[0] == 40000.0, x31[0] == 101000.0, x32[0] == 0.0, x33[0] == 0.0, x34[0] == 0.0, x35[0] == 0.0, x36[0] == 0.0, x37[0] == 0.0, x38[0] == 0.0, x39[0] == 0.0, x4[0] == 0.0, x40[0] == 0.0, x41[0] == 40000.0, x42[0] == 0.0, x43[0] == 0.0, x44[0] == 40000.0, x45[0] == 0.0, x46[0] == 0.0, x47[0] == 22000000.0, x48[0] == 0.0, x49[0] == 0.0, x5[0] == 0.0, x50[0] == 0.0, x51[0] == 0.0, x52[0] == 0.0, x53[0] == 40000.0, x54[0] == 0.0, x55[0] == 21000000.0, x56[0] == 0.0, x57[0] == 0.0, x58[0] == 0.0, x59[0] == 0.0, x6[0] == 0.0, x60[0] == 10000000.0, x61[0] == 0.0, x62[0] == 0.0, x63[0] == 0.0, x64[0] == 0.0, x65[0] == 0.0, x66[0] == 0.0, x67[0] == 0.0, x68[0] == 0.0, x69[0] == 0.0, x7[0] == 0.0, x70[0] == 0.0, x71[0] == 0.0, x72[0] == 0.0, x73[0] == 0.0, x74[0] == 0.0, x75[0] == 0.0, x76[0] == 0.0, x77[0] == 0.0, x78[0] == 0.0, x79[0] == 0.0, x8[0] == 0.0, x80[0] == 0.0, x81[0] == 0.0, x82[0] == 0.0, x83[0] == 0.0, x84[0] == 0.0, x85[0] == 0.0, x86[0] == 0.0, x87[0] == 0.0, x88[0] == 0.0, x89[0] == 0.0, x9[0] == 0.0, x90[0] == 0.0, x91[0] == 0.0, x92[0] == 0.0, x93[0] == 0.0, x94[0] == 0.0 }; rates = { v1, v10, v100, v101, v102, v103, v104, v105, v106, v107, v108, v109, v11, v110, v111, v112, v113, v114, v115, v116, v117, v118, v119, v12, v120, v121, v122, v123, v124, v125, v13, v14, v15, v16, v17, v18, v19, v2, v20, v21, v22, v23, v24, v25, v26, v27, v28, v29, v3, v30, v31, v32, v33, v34, v35, v36, v37, v38, v39, v4, v40, v41, v42, v43, v44, v45, v46, v47, v48, v49, v5, v50, v51, v52, v53, v54, v55, v56, v57, v58, v59, v6, v60, v61, v62, v63, v64, v65, v66, v67, v68, v69, v7, v70, v71, v72, v73, v74, v75, v76, v77, v78, v79, v8, v80, v81, v82, v83, v84, v85, v86, v87, v88, v89, v9, v90, v91, v92, v93, v94, v95, v96, v97, v98, v99 }; rateEquations = { v1 -> k1*x1*x2[t] - kr1*x3[t], v10 -> -(kr10*x10[t]) + k10*x16[t]*x6[t], v100 -> k58*x60[t]*x81[t] - kr58*x85[t], v101 -> k59*x85[t], v102 -> k6*x15[t] - kr6*x17[t], v103 -> k6*x32[t] - kr6*x63[t], v104 -> k6*x33[t] - kr6*x64[t], v105 -> -(kr6*x19[t]) + k6*x25[t], v106 -> k4*x12[t]*x25[t] - kr4*x88[t], v107 -> k5*x88[t], v108 -> -(kr6*x20[t]) + k6*x27[t], v109 -> k4*x12[t]*x27[t] - kr4*x89[t], v11 -> k2*x10[t]^2 - kr2*x11[t], v110 -> k5*x89[t], v111 -> -(kr6*x21[t]) + k6*x29[t], v112 -> k4*x12[t]*x29[t] - kr4*x90[t], v113 -> k5*x90[t], v114 -> k6*x34[t] - kr6*x65[t], v115 -> k4*x12[t]*x34[t] - kr4*x91[t], v116 -> k5*x91[t], v117 -> k6*x35[t] - kr6*x66[t], v118 -> k4*x12[t]*x35[t] - kr4*x92[t], v119 -> k5*x92[t], v12 -> k3*x11[t] - kr3*x8[t], v120 -> k6*x36[t] - kr6*x67[t], v121 -> k4*x12[t]*x36[t] - kr4*x93[t], v122 -> k5*x93[t], v123 -> k6*x37[t] - kr6*x68[t], v124 -> k4*x12[t]*x37[t] - kr4*x94[t], v125 -> k5*x94[t], v13 -> k13, v14 -> -(kr14*x17[t]) + k14*x14[t]*x8[t], v15 -> k15*x9[t], v16 -> k16*x15[t]*x22[t] - kr16*x23[t], v17 -> k17*x23[t]*x24[t] - kr17*x25[t], v18 -> k18*x25[t]*x26[t] - kr18*x27[t], v19 -> k19*x27[t] - kr19*x25[t]*x28[t], v2 -> k2*x3[t]^2 - kr2*x4[t], v20 -> -(kr20*x29[t]) + k20*x25[t]*x43[t], v21 -> -(kr21*x25[t]*x26[t]) + k21*x29[t], v22 -> k22*x15[t]*x31[t] - kr22*x32[t], v23 -> k23*x32[t] - kr23*x33[t], v24 -> k24*x22[t]*x33[t] - kr24*x34[t], v25 -> k25*x24[t]*x34[t] - kr25*x35[t], v26 -> k18*x26[t]*x35[t] - kr18*x36[t], v27 -> -(kr19*x28[t]*x35[t]) + k19*x36[t], v28 -> k28*x28[t]*x41[t] - kr28*x42[t], v29 -> k29*x42[t] - kr29*x43[t]*x45[t], v3 -> k3*x4[t] - kr3*x5[t], v30 -> -(kr20*x37[t]) + k20*x35[t]*x43[t], v31 -> -(kr21*x26[t]*x35[t]) + k21*x37[t], v32 -> k32*x35[t] - kr32*x15[t]*x38[t], v33 -> k33*x38[t] - kr33*x30[t]*x40[t], v34 -> k34*x25[t] - kr34*x15[t]*x30[t], v35 -> -(kr35*x22[t]*x24[t]) + k35*x30[t], v36 -> (Vm36*x40[t])/(Km36 + x40[t]), v37 -> k37*x33[t] - kr37*x15[t]*x40[t], v38 -> -(kr24*x39[t]) + k24*x22[t]*x40[t], v39 -> k37*x34[t] - kr37*x15[t]*x39[t], v4 -> k4*x12[t]*x23[t] - kr4*x7[t], v40 -> -(kr40*x38[t]) + k40*x24[t]*x39[t], v41 -> k41*x30[t]*x33[t] - kr41*x35[t], v42 -> k42*x44[t]*x45[t] - kr42*x46[t], v43 -> k43*x46[t], v44 -> k44*x45[t]*x47[t] - kr44*x48[t], v45 -> k45*x48[t], v46 -> k44*x45[t]*x49[t] - kr44*x50[t], v47 -> k47*x50[t], v48 -> -(kr48*x52[t]) + k48*x51[t]*x53[t], v49 -> k49*x52[t], v5 -> k5*x7[t], v50 -> k50*x49[t]*x53[t] - kr50*x54[t], v51 -> k49*x54[t], v52 -> k52*x51[t]*x55[t] - kr52*x56[t], v53 -> k53*x56[t], v54 -> k52*x51[t]*x57[t] - kr52*x58[t], v55 -> k55*x58[t], v56 -> k56*x59[t]*x60[t] - kr56*x61[t], v57 -> k57*x61[t], v58 -> k58*x57[t]*x60[t] - kr58*x62[t], v59 -> k59*x62[t], v6 -> k6*x2[t] - kr6*x6[t], v60 -> k60*x6[t], v61 -> k61*x16[t], v62 -> k60*x8[t], v63 -> -(kr16*x18[t]) + k16*x17[t]*x22[t], v64 -> -(kr17*x19[t]) + k17*x18[t]*x24[t], v65 -> -(kr18*x20[t]) + k18*x19[t]*x26[t], v66 -> k19*x20[t] - kr19*x19[t]*x69[t], v67 -> -(kr20*x21[t]) + k20*x19[t]*x71[t], v68 -> k21*x21[t] - kr21*x19[t]*x26[t], v69 -> k22*x17[t]*x31[t] - kr22*x63[t], v7 -> k7*x5[t], v70 -> k23*x63[t] - kr23*x64[t], v71 -> k24*x22[t]*x64[t] - kr24*x65[t], v72 -> k25*x24[t]*x65[t] - kr25*x66[t], v73 -> k18*x26[t]*x66[t] - kr18*x67[t], v74 -> k19*x67[t] - kr19*x66[t]*x69[t], v75 -> k28*x41[t]*x69[t] - kr28*x70[t], v76 -> k29*x70[t] - kr29*x71[t]*x72[t], v77 -> -(kr20*x68[t]) + k20*x66[t]*x71[t], v78 -> -(kr21*x26[t]*x66[t]) + k21*x68[t], v79 -> -(kr32*x17[t]*x38[t]) + k32*x66[t], v8 -> -(kr8*x15[t]) + k8*x14[t]*x5[t], v80 -> k34*x19[t] - kr34*x17[t]*x30[t], v81 -> -(kr37*x17[t]*x40[t]) + k37*x64[t], v82 -> -(kr37*x17[t]*x39[t]) + k37*x65[t], v83 -> k41*x30[t]*x64[t] - kr41*x66[t], v84 -> k42*x44[t]*x72[t] - kr42*x73[t], v85 -> k43*x73[t], v86 -> k44*x47[t]*x72[t] - kr44*x74[t], v87 -> k45*x74[t], v88 -> k44*x72[t]*x75[t] - kr44*x76[t], v89 -> k47*x76[t], v9 -> k7*x23[t], v90 -> k48*x53[t]*x77[t] - kr48*x78[t], v91 -> k49*x78[t], v92 -> k50*x53[t]*x75[t] - kr50*x79[t], v93 -> k49*x79[t], v94 -> k52*x55[t]*x77[t] - kr52*x80[t], v95 -> k53*x80[t], v96 -> k52*x77[t]*x81[t] - kr52*x82[t], v97 -> k55*x82[t], v98 -> k56*x60[t]*x83[t] - kr56*x84[t], v99 -> k57*x84[t] }; parameters = { Km36 -> 771977800000000.0, RT -> 50000.0, Vm36 -> 615.0325, k1 -> 0.002372521, k10 -> 3803.728, k11 -> 0.000480156, k12 -> 31.71871, k13 -> 0.4545611, k14 -> 6.370566*^-07, k15 -> 46468.78, k16 -> 0.0004021305, k17 -> 0.0003099213, k18 -> 0.004463938, k19 -> 349.772, k2 -> 0.000480156, k20 -> 5.17656*^-05, k21 -> 0.4722901, k22 -> 0.0001445554, k23 -> 420.3359, k24 -> 0.007178843, k25 -> 0.0006871213, k28 -> 9.826084*^-06, k29 -> 931.1092, k3 -> 31.71871, k32 -> 14.19908, k33 -> 10.96212, k34 -> 0.2467995, k35 -> 1.836058, k37 -> 29.34687, k4 -> 3.047285*^-05, k40 -> 7.409959*^-05, k41 -> 0.001522817, k42 -> 0.009688174, k43 -> 51.60945, k44 -> 0.001406622, k45 -> 6340.081, k47 -> 1632.425, k48 -> 0.0006874119, k49 -> 10.73099, k50 -> 0.0005464454, k52 -> 0.003826571, k53 -> 62181.84, k55 -> 1120.398, k56 -> 0.004700229, k57 -> 19.75184, k58 -> 0.0001714511, k59 -> 6.409354, k6 -> 0.0004123214, k60 -> 0.08693199, k61 -> 0.006499143, k7 -> 0.003011324, k8 -> 0.0005174108, kr1 -> 0.1146248, kr10 -> 171.6947, kr11 -> 0.5100538, kr12 -> 2.220991, kr14 -> 196.6479, kr16 -> 0.4509308, kr17 -> 2.524092, kr18 -> 11.1361, kr19 -> 5.84737*^-06, kr2 -> 0.5100538, kr20 -> 12.816, kr21 -> 1.714441*^-05, kr22 -> 0.6220457, kr23 -> 17.39321, kr24 -> 563.2135, kr25 -> 1.218132, kr28 -> 0.9683624, kr29 -> 0.0001096614, kr3 -> 2.220991, kr32 -> 5.54527*^-05, kr33 -> 1.788597*^-05, kr34 -> 0.0001283286, kr35 -> 0.0003866434, kr37 -> 5.477036*^-06, kr4 -> 0.1230832, kr40 -> 2.748877, kr41 -> 44.60169, kr42 -> 1.870396, kr44 -> 0.5985189, kr48 -> 1489.015, kr50 -> 9.954943, kr52 -> 19.85279, kr56 -> 1.229629, kr58 -> 0.1138168, kr6 -> 0.294324, kr8 -> 0.9058936, x1 -> 4962.0, c1 -> 1.0, c2 -> 1.0, c3 -> 4.3*^-06 }; assignments = { k5 -> Piecewise[{{1.55, RT/(1 + kr1/(k1*x1)) < 3100}, {0.2, RT/(1 + kr1/(k1*x1)) > 100000}}, 1.55 - (0.000013500000000000001*RT)/(1 + kr1/(k1*x1))], C -> RT/(1 + kr1/(k1*x1)), EGF\[LetterSpace]EGFR\[LetterSpace]act -> x11[t] + x15[t] + x17[t] + x18[t] + x19[t] + x20[t] + x21[t] + x23[t] + x25[t] + x27[t] + x29[t] + x32[t] + x33[t] + x34[t] + x35[t] + x36[t] + x37[t] + x5[t] + x63[t] + x64[t] + x65[t] + x66[t] + x67[t] + x68[t] + x7[t] + x8[t] + x88[t] + x89[t] + x90[t] + x91[t] + x92[t] + x93[t] + x94[t], SHC\[LetterSpace]P\[LetterSpace]t -> x33[t] + x34[t] + x35[t] + x36[t] + x37[t] + x38[t] + x39[t] + x40[t] + x64[t] + x65[t] + x66[t] + x67[t] + x68[t] + x91[t] + x92[t] + x93[t] + x94[t], Raf\[LetterSpace]act -> x45[t] + x46[t] + x48[t] + x50[t] + x72[t] + x73[t] + x74[t] + x76[t], ERK\[LetterSpace]PP -> x59[t] + x83[t], MEK\[LetterSpace]PP -> x51[t] + x77[t], Ras\[LetterSpace]GTP -> x28[t] + x42[t] + x69[t] + x70[t] }; events = { }; speciesAnnotations = { }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { x10'[t] == 1.0*v10 -2.0*v11, x11'[t] == 1.0*v11 -1.0*v12, x12'[t] == 1.0*v15 -1.0*v4 -1.0*v106 -1.0*v109 -1.0*v112 -1.0*v115 -1.0*v118 -1.0*v121 -1.0*v124, x13'[t] == 1.0*v61 , x14'[t] == -1.0*v8 -1.0*v14, x15'[t] == 1.0*v8 +1.0*v32 +1.0*v34 +1.0*v37 +1.0*v39 -1.0*v16 -1.0*v22 -1.0*v102, x16'[t] == -1.0*v10 -1.0*v61, x17'[t] == 1.0*v14 +1.0*v79 +1.0*v80 +1.0*v81 +1.0*v82 +1.0*v102 -1.0*v63 -1.0*v69, x18'[t] == 1.0*v5 +1.0*v9 +1.0*v63 -1.0*v64, x19'[t] == 1.0*v64 +1.0*v66 +1.0*v68 +1.0*v105 +1.0*v107 -1.0*v65 -1.0*v67 -1.0*v80, x2'[t] == 1.0*v13 -1.0*v1 -1.0*v6, x20'[t] == 1.0*v65 +1.0*v108 +1.0*v110 -1.0*v66, x21'[t] == 1.0*v67 +1.0*v111 +1.0*v113 -1.0*v68, x22'[t] == 1.0*v35 -1.0*v16 -1.0*v24 -1.0*v38 -1.0*v63 -1.0*v71, x23'[t] == 1.0*v16 -1.0*v4 -1.0*v9 -1.0*v17, x24'[t] == 1.0*v35 -1.0*v17 -1.0*v25 -1.0*v40 -1.0*v64 -1.0*v72, x25'[t] == 1.0*v17 +1.0*v19 +1.0*v21 -1.0*v18 -1.0*v20 -1.0*v34 -1.0*v105 -1.0*v106, x26'[t] == 1.0*v21 +1.0*v31 +1.0*v68 +1.0*v78 -1.0*v18 -1.0*v26 -1.0*v65 -1.0*v73, x27'[t] == 1.0*v18 -1.0*v19 -1.0*v108 -1.0*v109, x28'[t] == 1.0*v19 +1.0*v27 -1.0*v28, x29'[t] == 1.0*v20 -1.0*v21 -1.0*v111 -1.0*v112, x3'[t] == 1.0*v1 -2.0*v2, x30'[t] == 1.0*v33 +1.0*v34 +1.0*v80 -1.0*v35 -1.0*v41 -1.0*v83, x31'[t] == 1.0*v36 -1.0*v22 -1.0*v69, x32'[t] == 1.0*v22 -1.0*v23 -1.0*v103, x33'[t] == 1.0*v23 -1.0*v24 -1.0*v37 -1.0*v41 -1.0*v104, x34'[t] == 1.0*v24 -1.0*v25 -1.0*v39 -1.0*v114 -1.0*v115, x35'[t] == 1.0*v25 +1.0*v27 +1.0*v31 +1.0*v41 -1.0*v26 -1.0*v30 -1.0*v32 -1.0*v117 -1.0*v118, x36'[t] == 1.0*v26 -1.0*v27 -1.0*v120 -1.0*v121, x37'[t] == 1.0*v30 -1.0*v31 -1.0*v123 -1.0*v124, x38'[t] == 1.0*v32 +1.0*v40 +1.0*v79 -1.0*v33, x39'[t] == 1.0*v38 +1.0*v39 +1.0*v82 -1.0*v40, x4'[t] == 1.0*v2 -1.0*v3, x40'[t] == 1.0*v33 +1.0*v37 +1.0*v81 -1.0*v36 -1.0*v38, x41'[t] == 1.0*v43 +1.0*v85 -1.0*v28 -1.0*v75, x42'[t] == 1.0*v28 -1.0*v29, x43'[t] == 1.0*v29 -1.0*v20 -1.0*v30, x44'[t] == 1.0*v43 +1.0*v85 -1.0*v42 -1.0*v84, x45'[t] == 1.0*v29 +1.0*v45 +1.0*v47 -1.0*v42 -1.0*v44 -1.0*v46, x46'[t] == 1.0*v42 -1.0*v43, x47'[t] == 1.0*v51 +1.0*v93 -1.0*v44 -1.0*v86, x48'[t] == 1.0*v44 -1.0*v45, x49'[t] == 1.0*v45 +1.0*v49 -1.0*v46 -1.0*v50, x5'[t] == 1.0*v3 -1.0*v7 -1.0*v8, x50'[t] == 1.0*v46 -1.0*v47, x51'[t] == 1.0*v47 +1.0*v53 +1.0*v55 -1.0*v48 -1.0*v52 -1.0*v54, x52'[t] == 1.0*v48 -1.0*v49, x53'[t] == 1.0*v49 +1.0*v51 +1.0*v91 +1.0*v93 -1.0*v48 -1.0*v50 -1.0*v90 -1.0*v92, x54'[t] == 1.0*v50 -1.0*v51, x55'[t] == 1.0*v59 +1.0*v101 -1.0*v52 -1.0*v94, x56'[t] == 1.0*v52 -1.0*v53, x57'[t] == 1.0*v53 +1.0*v57 -1.0*v54 -1.0*v58, x58'[t] == 1.0*v54 -1.0*v55, x59'[t] == 1.0*v55 -1.0*v56, x6'[t] == 1.0*v6 -1.0*v10 -1.0*v60, x60'[t] == 1.0*v57 +1.0*v59 +1.0*v99 +1.0*v101 -1.0*v56 -1.0*v58 -1.0*v98 -1.0*v100, x61'[t] == 1.0*v56 -1.0*v57, x62'[t] == 1.0*v58 -1.0*v59, x63'[t] == 1.0*v69 +1.0*v103 -1.0*v70, x64'[t] == 1.0*v70 +1.0*v104 -1.0*v71 -1.0*v81 -1.0*v83, x65'[t] == 1.0*v71 +1.0*v114 +1.0*v116 -1.0*v72 -1.0*v82, x66'[t] == 1.0*v72 +1.0*v74 +1.0*v78 +1.0*v83 +1.0*v117 +1.0*v119 -1.0*v73 -1.0*v77 -1.0*v79, x67'[t] == 1.0*v73 +1.0*v120 +1.0*v122 -1.0*v74, x68'[t] == 1.0*v77 +1.0*v123 +1.0*v125 -1.0*v78, x69'[t] == 1.0*v66 +1.0*v74 -1.0*v75, x7'[t] == 1.0*v4 -1.0*v5, x70'[t] == 1.0*v75 -1.0*v76, x71'[t] == 1.0*v76 -1.0*v67 -1.0*v77, x72'[t] == 1.0*v76 +1.0*v87 +1.0*v89 -1.0*v84 -1.0*v86 -1.0*v88, x73'[t] == 1.0*v84 -1.0*v85, x74'[t] == 1.0*v86 -1.0*v87, x75'[t] == 1.0*v87 +1.0*v91 -1.0*v88 -1.0*v92, x76'[t] == 1.0*v88 -1.0*v89, x77'[t] == 1.0*v89 +1.0*v95 +1.0*v97 -1.0*v90 -1.0*v94 -1.0*v96, x78'[t] == 1.0*v90 -1.0*v91, x79'[t] == 1.0*v92 -1.0*v93, x8'[t] == 1.0*v7 +1.0*v12 -1.0*v14 -1.0*v62, x80'[t] == 1.0*v94 -1.0*v95, x81'[t] == 1.0*v95 +1.0*v99 -1.0*v96 -1.0*v100, x82'[t] == 1.0*v96 -1.0*v97, x83'[t] == 1.0*v97 -1.0*v98, x84'[t] == 1.0*v98 -1.0*v99, x85'[t] == 1.0*v100 -1.0*v101, x86'[t] == 1.0*v60 , x87'[t] == 1.0*v62 , x88'[t] == 1.0*v106 -1.0*v107, x89'[t] == 1.0*v109 -1.0*v110, x9'[t] == 1.0*v5 +1.0*v107 +1.0*v110 +1.0*v113 +1.0*v116 +1.0*v119 +1.0*v122 +1.0*v125 -1.0*v15, x90'[t] == 1.0*v112 -1.0*v113, x91'[t] == 1.0*v115 -1.0*v116, x92'[t] == 1.0*v118 -1.0*v119, x93'[t] == 1.0*v121 -1.0*v122, x94'[t] == 1.0*v124 -1.0*v125 }; 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]}]