(* 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 = { \[LetterSpace]101[t], \[LetterSpace]105[t], \[LetterSpace]11[t], \[LetterSpace]129[t], \[LetterSpace]13[t], \[LetterSpace]147[t], \[LetterSpace]15[t], \[LetterSpace]153[t], \[LetterSpace]174[t], \[LetterSpace]181[t], \[LetterSpace]19[t], \[LetterSpace]198[t], \[LetterSpace]21[t], \[LetterSpace]25[t], \[LetterSpace]27[t], \[LetterSpace]5[t], \[LetterSpace]75[t], \[LetterSpace]79[t], \[LetterSpace]84[t], \[LetterSpace]9[t], \[LetterSpace]96[t], \[LetterSpace]99[t], species\[LetterSpace]1[t], species\[LetterSpace]10[t], species\[LetterSpace]11[t], species\[LetterSpace]12[t], species\[LetterSpace]13[t], species\[LetterSpace]14[t], species\[LetterSpace]15[t], species\[LetterSpace]16[t], species\[LetterSpace]17[t], species\[LetterSpace]18[t], species\[LetterSpace]19[t], species\[LetterSpace]2[t], species\[LetterSpace]20[t], species\[LetterSpace]21[t], species\[LetterSpace]22[t], species\[LetterSpace]23[t], species\[LetterSpace]24[t], species\[LetterSpace]25[t], species\[LetterSpace]26[t], species\[LetterSpace]27[t], species\[LetterSpace]28[t], species\[LetterSpace]29[t], species\[LetterSpace]3[t], species\[LetterSpace]30[t], species\[LetterSpace]31[t], species\[LetterSpace]4[t], species\[LetterSpace]5[t], species\[LetterSpace]6[t], species\[LetterSpace]7[t], species\[LetterSpace]8[t], species\[LetterSpace]9[t] }; initialValues = { \[LetterSpace]101[0] == 0.0, \[LetterSpace]105[0] == 0.0, \[LetterSpace]11[0] == 0.0, \[LetterSpace]129[0] == 0.0, \[LetterSpace]13[0] == 0.0, \[LetterSpace]147[0] == 0.0, \[LetterSpace]15[0] == 0.0, \[LetterSpace]153[0] == 0.0, \[LetterSpace]174[0] == 0.0, \[LetterSpace]181[0] == 0.0, \[LetterSpace]19[0] == 0.0, \[LetterSpace]198[0] == 0.0, \[LetterSpace]21[0] == 0.0, \[LetterSpace]25[0] == 0.0, \[LetterSpace]27[0] == 0.0, \[LetterSpace]5[0] == 0.0, \[LetterSpace]75[0] == 0.0, \[LetterSpace]79[0] == 0.0, \[LetterSpace]84[0] == 0.0, \[LetterSpace]9[0] == 0.0, \[LetterSpace]96[0] == 0.0, \[LetterSpace]99[0] == 0.0, species\[LetterSpace]1[0] == 0.0, species\[LetterSpace]10[0] == 0.0, species\[LetterSpace]11[0] == 0.0, species\[LetterSpace]12[0] == 0.0, species\[LetterSpace]13[0] == 0.0, species\[LetterSpace]14[0] == 0.0, species\[LetterSpace]15[0] == 0.0, species\[LetterSpace]16[0] == 0.0, species\[LetterSpace]17[0] == 0.0, species\[LetterSpace]18[0] == 0.0, species\[LetterSpace]19[0] == 0.0, species\[LetterSpace]2[0] == 0.0, species\[LetterSpace]20[0] == 0.0, species\[LetterSpace]21[0] == 0.0, species\[LetterSpace]22[0] == 0.0, species\[LetterSpace]23[0] == 0.0, species\[LetterSpace]24[0] == 0.0, species\[LetterSpace]25[0] == 0.0, species\[LetterSpace]26[0] == 0.0, species\[LetterSpace]27[0] == 0.0, species\[LetterSpace]28[0] == 0.0999999951844375, species\[LetterSpace]29[0] == 0.0, species\[LetterSpace]3[0] == 0.0, species\[LetterSpace]30[0] == 0.0, species\[LetterSpace]31[0] == 0.0, species\[LetterSpace]4[0] == 0.0, species\[LetterSpace]5[0] == 0.0, species\[LetterSpace]6[0] == 0.0, species\[LetterSpace]7[0] == 0.0, species\[LetterSpace]8[0] == 0.0, species\[LetterSpace]9[0] == 0.0 }; rates = { \[LetterSpace]103, \[LetterSpace]125, \[LetterSpace]132, \[LetterSpace]139, \[LetterSpace]149, \[LetterSpace]156, \[LetterSpace]164, \[LetterSpace]172, \[LetterSpace]178, \[LetterSpace]192, \[LetterSpace]194, \[LetterSpace]29, \[LetterSpace]31, \[LetterSpace]33, \[LetterSpace]35, \[LetterSpace]37, \[LetterSpace]39, \[LetterSpace]41, \[LetterSpace]43, \[LetterSpace]77, \[LetterSpace]86, 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]57, reaction\[LetterSpace]58, reaction\[LetterSpace]59, reaction\[LetterSpace]6, reaction\[LetterSpace]60, reaction\[LetterSpace]61, reaction\[LetterSpace]62, reaction\[LetterSpace]63, reaction\[LetterSpace]64, reaction\[LetterSpace]65, reaction\[LetterSpace]66, reaction\[LetterSpace]67, reaction\[LetterSpace]68, reaction\[LetterSpace]69, reaction\[LetterSpace]7, reaction\[LetterSpace]70, reaction\[LetterSpace]8, reaction\[LetterSpace]9 }; rateEquations = { \[LetterSpace]103 -> \[LetterSpace]1*(-(\[LetterSpace]103\[LetterSpace]k2*\[LetterSpace]105[t]) + \[LetterSpace]103\[LetterSpace]k1*\[LetterSpace]101[t]*\[LetterSpace]99[t]), \[LetterSpace]125 -> \[LetterSpace]1*function\[LetterSpace]3[\[LetterSpace]125\[LetterSpace]k, \[LetterSpace]96[t], \[LetterSpace]99[t], \[LetterSpace]125\[LetterSpace]km], \[LetterSpace]132 -> \[LetterSpace]1*function\[LetterSpace]3[\[LetterSpace]132\[LetterSpace]k, \[LetterSpace]96[t], \[LetterSpace]105[t], \[LetterSpace]132\[LetterSpace]km], \[LetterSpace]139 -> \[LetterSpace]139\[LetterSpace]k1*\[LetterSpace]129[t], \[LetterSpace]149 -> \[LetterSpace]1*\[LetterSpace]149\[LetterSpace]k1*\[LetterSpace]129[t]^2*\[LetterSpace]147[t], \[LetterSpace]156 -> \[LetterSpace]156\[LetterSpace]k1*\[LetterSpace]147[t] - \[LetterSpace]156\[LetterSpace]k2*\[LetterSpace]5[t], \[LetterSpace]164 -> \[LetterSpace]164\[LetterSpace]k1*\[LetterSpace]153[t], \[LetterSpace]172 -> -(\[LetterSpace]172\[LetterSpace]k2*\[LetterSpace]174[t]) + \[LetterSpace]172\[LetterSpace]k1*\[LetterSpace]25[t], \[LetterSpace]178 -> -(\[LetterSpace]178\[LetterSpace]k2*\[LetterSpace]181[t]) + \[LetterSpace]178\[LetterSpace]k1*\[LetterSpace]27[t], \[LetterSpace]192 -> \[LetterSpace]1*\[LetterSpace]192\[LetterSpace]k1*\[LetterSpace]181[t]*\[LetterSpace]96[t], \[LetterSpace]194 -> \[LetterSpace]1*(-(\[LetterSpace]194\[LetterSpace]k2*\[LetterSpace]198[t]) + \[LetterSpace]194\[LetterSpace]k1*\[LetterSpace]174[t]*\[LetterSpace]96[t]), \[LetterSpace]29 -> \[LetterSpace]29\[LetterSpace]k1*\[LetterSpace]3*species\[LetterSpace]27[t], \[LetterSpace]31 -> \[LetterSpace]3*function\[LetterSpace]1[\[LetterSpace]9[t], \[LetterSpace]31\[LetterSpace]Km, \[LetterSpace]31\[LetterSpace]V], \[LetterSpace]33 -> \[LetterSpace]3*(-(\[LetterSpace]33\[LetterSpace]k2*\[LetterSpace]13[t]) + \[LetterSpace]33\[LetterSpace]k1*\[LetterSpace]11[t]^2*\[LetterSpace]9[t]), \[LetterSpace]35 -> \[LetterSpace]35\[LetterSpace]k1*species\[LetterSpace]13[t]*\[LetterSpace]15[t]^3, \[LetterSpace]37 -> \[LetterSpace]3*\[LetterSpace]37\[LetterSpace]k1*\[LetterSpace]19[t]^2*\[LetterSpace]5[t], \[LetterSpace]39 -> \[LetterSpace]3*function\[LetterSpace]1[\[LetterSpace]19[t], \[LetterSpace]39\[LetterSpace]Km, \[LetterSpace]39\[LetterSpace]V], \[LetterSpace]41 -> \[LetterSpace]3*\[LetterSpace]41\[LetterSpace]k1*\[LetterSpace]15[t]*\[LetterSpace]21[t], \[LetterSpace]43 -> \[LetterSpace]3*(\[LetterSpace]43\[LetterSpace]k1*\[LetterSpace]15[t]*\[LetterSpace]25[t] - \[LetterSpace]43\[LetterSpace]k2*\[LetterSpace]27[t]), \[LetterSpace]77 -> \[LetterSpace]1*function\[LetterSpace]2[parameter\[LetterSpace]1, \[LetterSpace]77\[LetterSpace]k, \[LetterSpace]75[t], \[LetterSpace]75[t]], \[LetterSpace]86 -> \[LetterSpace]1*(\[LetterSpace]86\[LetterSpace]k1*\[LetterSpace]79[t]*\[LetterSpace]84[t]^2 - \[LetterSpace]86\[LetterSpace]k2*\[LetterSpace]96[t]), reaction\[LetterSpace]1 -> \[LetterSpace]1*function\[LetterSpace]4[reaction\[LetterSpace]1\[LetterSpace]v], reaction\[LetterSpace]10 -> reaction\[LetterSpace]10\[LetterSpace]k1*\[LetterSpace]1*\[LetterSpace]84[t], reaction\[LetterSpace]11 -> reaction\[LetterSpace]11\[LetterSpace]k1*\[LetterSpace]1*\[LetterSpace]75[t], reaction\[LetterSpace]12 -> reaction\[LetterSpace]12\[LetterSpace]k1*\[LetterSpace]1*\[LetterSpace]79[t], reaction\[LetterSpace]13 -> \[LetterSpace]1*function\[LetterSpace]4[reaction\[LetterSpace]13\[LetterSpace]v], reaction\[LetterSpace]14 -> \[LetterSpace]1*function\[LetterSpace]4[reaction\[LetterSpace]14\[LetterSpace]v], reaction\[LetterSpace]15 -> reaction\[LetterSpace]15\[LetterSpace]k1*\[LetterSpace]1*\[LetterSpace]147[t], reaction\[LetterSpace]16 -> reaction\[LetterSpace]16\[LetterSpace]k1*\[LetterSpace]1*\[LetterSpace]101[t], reaction\[LetterSpace]17 -> \[LetterSpace]1*function\[LetterSpace]1[\[LetterSpace]129[t], reaction\[LetterSpace]17\[LetterSpace]Km, reaction\[LetterSpace]17\[LetterSpace]V], reaction\[LetterSpace]18 -> reaction\[LetterSpace]18\[LetterSpace]k1*\[LetterSpace]1*species\[LetterSpace]2[t], reaction\[LetterSpace]19 -> -(reaction\[LetterSpace]19\[LetterSpace]k2*species\[LetterSpace]1[t]) + reaction\[LetterSpace]19\[LetterSpace]k1*\[LetterSpace]15[t], reaction\[LetterSpace]2 -> \[LetterSpace]1*function\[LetterSpace]5[reaction\[LetterSpace]2\[LetterSpace]k, reaction\[LetterSpace]2\[LetterSpace]k1, species\[LetterSpace]30[t]], reaction\[LetterSpace]20 -> reaction\[LetterSpace]20\[LetterSpace]k1*\[LetterSpace]1*\[LetterSpace]99[t], reaction\[LetterSpace]21 -> reaction\[LetterSpace]21\[LetterSpace]k1*species\[LetterSpace]2[t] - reaction\[LetterSpace]21\[LetterSpace]k2*\[LetterSpace]11[t], reaction\[LetterSpace]22 -> \[LetterSpace]1*function\[LetterSpace]4[reaction\[LetterSpace]22\[LetterSpace]v], reaction\[LetterSpace]23 -> reaction\[LetterSpace]23\[LetterSpace]k1*\[LetterSpace]1*species\[LetterSpace]3[t], reaction\[LetterSpace]24 -> \[LetterSpace]1*function\[LetterSpace]1[species\[LetterSpace]4[t], reaction\[LetterSpace]24\[LetterSpace]Km, reaction\[LetterSpace]24\[LetterSpace]V], reaction\[LetterSpace]25 -> \[LetterSpace]3*function\[LetterSpace]1[species\[LetterSpace]18[t], reaction\[LetterSpace]25\[LetterSpace]Km, reaction\[LetterSpace]25\[LetterSpace]V], reaction\[LetterSpace]26 -> \[LetterSpace]3*function\[LetterSpace]1[species\[LetterSpace]20[t], reaction\[LetterSpace]26\[LetterSpace]Km, reaction\[LetterSpace]26\[LetterSpace]V], reaction\[LetterSpace]27 -> \[LetterSpace]1*function\[LetterSpace]3[reaction\[LetterSpace]27\[LetterSpace]k, \[LetterSpace]96[t], species\[LetterSpace]3[t], reaction\[LetterSpace]27\[LetterSpace]km], reaction\[LetterSpace]28 -> \[LetterSpace]1*function\[LetterSpace]3[reaction\[LetterSpace]28\[LetterSpace]k, \[LetterSpace]96[t], species\[LetterSpace]5[t], reaction\[LetterSpace]28\[LetterSpace]km], reaction\[LetterSpace]29 -> reaction\[LetterSpace]29\[LetterSpace]k1*\[LetterSpace]3*species\[LetterSpace]17[t], reaction\[LetterSpace]3 -> \[LetterSpace]1*function\[LetterSpace]6[reaction\[LetterSpace]3\[LetterSpace]k, reaction\[LetterSpace]3\[LetterSpace]k1, species\[LetterSpace]30[t], species\[LetterSpace]16[t], species\[LetterSpace]23[t]], reaction\[LetterSpace]30 -> -(reaction\[LetterSpace]30\[LetterSpace]k2*species\[LetterSpace]17[t]) + reaction\[LetterSpace]30\[LetterSpace]k1*species\[LetterSpace]3[t], reaction\[LetterSpace]31 -> \[LetterSpace]1*(-(reaction\[LetterSpace]31\[LetterSpace]k2*species\[LetterSpace]5[t]) + reaction\[LetterSpace]31\[LetterSpace]k1*species\[LetterSpace]3[t]*\[LetterSpace]101[t]), reaction\[LetterSpace]32 -> reaction\[LetterSpace]32\[LetterSpace]k1*\[LetterSpace]3*species\[LetterSpace]19[t], reaction\[LetterSpace]33 -> reaction\[LetterSpace]33\[LetterSpace]k1*species\[LetterSpace]4[t], reaction\[LetterSpace]34 -> reaction\[LetterSpace]34\[LetterSpace]k1*\[LetterSpace]1*species\[LetterSpace]4[t]^2*\[LetterSpace]147[t], reaction\[LetterSpace]35 -> reaction\[LetterSpace]35\[LetterSpace]k1*\[LetterSpace]3*species\[LetterSpace]20[t]^2*\[LetterSpace]5[t], reaction\[LetterSpace]36 -> \[LetterSpace]3*(-(reaction\[LetterSpace]36\[LetterSpace]k2*species\[LetterSpace]21[t]) + reaction\[LetterSpace]36\[LetterSpace]k1*species\[LetterSpace]18[t]*\[LetterSpace]11[t]^2), reaction\[LetterSpace]37 -> reaction\[LetterSpace]37\[LetterSpace]k1*species\[LetterSpace]6[t], reaction\[LetterSpace]38 -> reaction\[LetterSpace]38\[LetterSpace]k1*species\[LetterSpace]14[t], reaction\[LetterSpace]39 -> \[LetterSpace]1*function\[LetterSpace]7[reaction\[LetterSpace]39\[LetterSpace]k1, \[LetterSpace]174[t], species\[LetterSpace]15[t]], reaction\[LetterSpace]4 -> \[LetterSpace]1*function\[LetterSpace]5[reaction\[LetterSpace]4\[LetterSpace]k, reaction\[LetterSpace]4\[LetterSpace]k1, species\[LetterSpace]30[t]], reaction\[LetterSpace]40 -> reaction\[LetterSpace]40\[LetterSpace]k1*\[LetterSpace]1*species\[LetterSpace]7[t], reaction\[LetterSpace]41 -> -(reaction\[LetterSpace]41\[LetterSpace]k2*species\[LetterSpace]22[t]) + reaction\[LetterSpace]41\[LetterSpace]k1*species\[LetterSpace]7[t], reaction\[LetterSpace]42 -> \[LetterSpace]1*function\[LetterSpace]5[reaction\[LetterSpace]42\[LetterSpace]k, reaction\[LetterSpace]42\[LetterSpace]k1, species\[LetterSpace]30[t]], reaction\[LetterSpace]43 -> \[LetterSpace]1*function\[LetterSpace]7[reaction\[LetterSpace]43\[LetterSpace]k1, species\[LetterSpace]1[t], \[LetterSpace]174[t]], reaction\[LetterSpace]44 -> \[LetterSpace]3*(-(reaction\[LetterSpace]44\[LetterSpace]k2*species\[LetterSpace]23[t]) + reaction\[LetterSpace]44\[LetterSpace]k1*species\[LetterSpace]22[t]*\[LetterSpace]5[t]), reaction\[LetterSpace]45 -> \[LetterSpace]1*(reaction\[LetterSpace]45\[LetterSpace]k1*species\[LetterSpace]3[t]*species\[LetterSpace]7[t] - reaction\[LetterSpace]45\[LetterSpace]k2*species\[LetterSpace]8[t]), reaction\[LetterSpace]46 -> \[LetterSpace]3*(-(reaction\[LetterSpace]46\[LetterSpace]k2*species\[LetterSpace]25[t]) + reaction\[LetterSpace]46\[LetterSpace]k1*species\[LetterSpace]24[t]*\[LetterSpace]25[t]), reaction\[LetterSpace]47 -> \[LetterSpace]1*function\[LetterSpace]5[reaction\[LetterSpace]47\[LetterSpace]k, reaction\[LetterSpace]47\[LetterSpace]k1, species\[LetterSpace]30[t]], reaction\[LetterSpace]48 -> \[LetterSpace]3*(reaction\[LetterSpace]48\[LetterSpace]k1*species\[LetterSpace]18[t]*species\[LetterSpace]22[t]^2 - reaction\[LetterSpace]48\[LetterSpace]k2*species\[LetterSpace]26[t]), reaction\[LetterSpace]49 -> reaction\[LetterSpace]49\[LetterSpace]k1*\[LetterSpace]1*species\[LetterSpace]10[t]*\[LetterSpace]96[t], reaction\[LetterSpace]5 -> \[LetterSpace]1*function\[LetterSpace]5[reaction\[LetterSpace]5\[LetterSpace]k, reaction\[LetterSpace]5\[LetterSpace]k1, species\[LetterSpace]30[t]], reaction\[LetterSpace]50 -> -(reaction\[LetterSpace]50\[LetterSpace]k2*species\[LetterSpace]10[t]) + reaction\[LetterSpace]50\[LetterSpace]k1*species\[LetterSpace]25[t], reaction\[LetterSpace]51 -> \[LetterSpace]1*function\[LetterSpace]7[reaction\[LetterSpace]51\[LetterSpace]k1, species\[LetterSpace]9[t], \[LetterSpace]174[t]], reaction\[LetterSpace]52 -> reaction\[LetterSpace]52\[LetterSpace]k1*species\[LetterSpace]24[t] - reaction\[LetterSpace]52\[LetterSpace]k2*species\[LetterSpace]9[t], reaction\[LetterSpace]53 -> \[LetterSpace]3*(reaction\[LetterSpace]53\[LetterSpace]k1*species\[LetterSpace]18[t]*species\[LetterSpace]28[t] - reaction\[LetterSpace]53\[LetterSpace]k2*species\[LetterSpace]30[t]), reaction\[LetterSpace]54 -> \[LetterSpace]3*(-(reaction\[LetterSpace]54\[LetterSpace]k2*species\[LetterSpace]30[t]) + reaction\[LetterSpace]54\[LetterSpace]k1*species\[LetterSpace]28[t]*\[LetterSpace]9[t]), reaction\[LetterSpace]55 -> \[LetterSpace]3*(reaction\[LetterSpace]55\[LetterSpace]k1*species\[LetterSpace]21[t]*species\[LetterSpace]28[t] - reaction\[LetterSpace]55\[LetterSpace]k2*species\[LetterSpace]29[t]), reaction\[LetterSpace]56 -> \[LetterSpace]3*(-(reaction\[LetterSpace]56\[LetterSpace]k2*species\[LetterSpace]29[t]) + reaction\[LetterSpace]56\[LetterSpace]k1*species\[LetterSpace]28[t]*\[LetterSpace]13[t]), reaction\[LetterSpace]57 -> \[LetterSpace]3*(reaction\[LetterSpace]57\[LetterSpace]k1*species\[LetterSpace]26[t]*species\[LetterSpace]28[t] - reaction\[LetterSpace]57\[LetterSpace]k2*species\[LetterSpace]29[t]), reaction\[LetterSpace]58 -> -(reaction\[LetterSpace]58\[LetterSpace]k2*species\[LetterSpace]11[t]) + reaction\[LetterSpace]58\[LetterSpace]k1*species\[LetterSpace]22[t]^2*\[LetterSpace]9[t], reaction\[LetterSpace]59 -> reaction\[LetterSpace]59\[LetterSpace]k1*species\[LetterSpace]11[t]*species\[LetterSpace]28[t] - reaction\[LetterSpace]59\[LetterSpace]k2*species\[LetterSpace]29[t], reaction\[LetterSpace]6 -> \[LetterSpace]1*function\[LetterSpace]4[reaction\[LetterSpace]6\[LetterSpace]v], reaction\[LetterSpace]60 -> \[LetterSpace]3*(reaction\[LetterSpace]60\[LetterSpace]k1*species\[LetterSpace]23[t]*species\[LetterSpace]28[t] - reaction\[LetterSpace]60\[LetterSpace]k2*species\[LetterSpace]29[t]), reaction\[LetterSpace]61 -> \[LetterSpace]1*(-(reaction\[LetterSpace]61\[LetterSpace]k2*species\[LetterSpace]12[t]) + reaction\[LetterSpace]61\[LetterSpace]k1*species\[LetterSpace]7[t]*\[LetterSpace]99[t]), reaction\[LetterSpace]62 -> -(reaction\[LetterSpace]62\[LetterSpace]k2*species\[LetterSpace]13[t]) + reaction\[LetterSpace]62\[LetterSpace]k1*\[LetterSpace]11[t]^3*\[LetterSpace]19[t]^3, reaction\[LetterSpace]63 -> -(reaction\[LetterSpace]63\[LetterSpace]k2*species\[LetterSpace]14[t]) + reaction\[LetterSpace]63\[LetterSpace]k1*species\[LetterSpace]20[t]^3*\[LetterSpace]11[t]^3, reaction\[LetterSpace]64 -> reaction\[LetterSpace]64\[LetterSpace]k1*\[LetterSpace]3*\[LetterSpace]9[t], reaction\[LetterSpace]65 -> reaction\[LetterSpace]65\[LetterSpace]k1*\[LetterSpace]3*species\[LetterSpace]18[t], reaction\[LetterSpace]66 -> \[LetterSpace]1*function\[LetterSpace]4[reaction\[LetterSpace]66\[LetterSpace]v], reaction\[LetterSpace]67 -> reaction\[LetterSpace]67\[LetterSpace]k1*\[LetterSpace]1*species\[LetterSpace]15[t], reaction\[LetterSpace]68 -> reaction\[LetterSpace]68\[LetterSpace]k1*species\[LetterSpace]13[t]*species\[LetterSpace]31[t]^3, reaction\[LetterSpace]69 -> reaction\[LetterSpace]69\[LetterSpace]k1*species\[LetterSpace]14[t]*species\[LetterSpace]31[t]^3, reaction\[LetterSpace]7 -> \[LetterSpace]3*(-(reaction\[LetterSpace]7\[LetterSpace]k2*species\[LetterSpace]16[t]) + reaction\[LetterSpace]7\[LetterSpace]k1*\[LetterSpace]11[t]*\[LetterSpace]5[t]), reaction\[LetterSpace]70 -> reaction\[LetterSpace]70\[LetterSpace]k1*species\[LetterSpace]15[t] - reaction\[LetterSpace]70\[LetterSpace]k2*species\[LetterSpace]31[t], reaction\[LetterSpace]8 -> -(reaction\[LetterSpace]8\[LetterSpace]k2*\[LetterSpace]21[t]) + reaction\[LetterSpace]8\[LetterSpace]k1*\[LetterSpace]99[t], reaction\[LetterSpace]9 -> reaction\[LetterSpace]9\[LetterSpace]k1*\[LetterSpace]1*\[LetterSpace]96[t] }; parameters = { parameter\[LetterSpace]1 -> 0.01, reaction\[LetterSpace]31\[LetterSpace]k1 -> 1.0, reaction\[LetterSpace]31\[LetterSpace]k2 -> 0.1, reaction\[LetterSpace]32\[LetterSpace]k1 -> 0.0492, reaction\[LetterSpace]33\[LetterSpace]k1 -> 16.6, reaction\[LetterSpace]34\[LetterSpace]k1 -> 1000.0, reaction\[LetterSpace]35\[LetterSpace]k1 -> 255.068, reaction\[LetterSpace]36\[LetterSpace]k1 -> 1.6, reaction\[LetterSpace]36\[LetterSpace]k2 -> 1.6, reaction\[LetterSpace]37\[LetterSpace]k1 -> 0.16, reaction\[LetterSpace]38\[LetterSpace]k1 -> 0.232, reaction\[LetterSpace]39\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]40\[LetterSpace]k1 -> 0.232, reaction\[LetterSpace]41\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]41\[LetterSpace]k2 -> 0.2, reaction\[LetterSpace]42\[LetterSpace]k -> 2*^-05, reaction\[LetterSpace]42\[LetterSpace]k1 -> 0.00055, reaction\[LetterSpace]43\[LetterSpace]k1 -> 0.05, reaction\[LetterSpace]44\[LetterSpace]k1 -> 1.0, reaction\[LetterSpace]44\[LetterSpace]k2 -> 0.05288, reaction\[LetterSpace]45\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]45\[LetterSpace]k2 -> 0.1, reaction\[LetterSpace]46\[LetterSpace]k1 -> 2.9, reaction\[LetterSpace]46\[LetterSpace]k2 -> 0.2, reaction\[LetterSpace]47\[LetterSpace]k -> 0.0001, reaction\[LetterSpace]47\[LetterSpace]k1 -> 0.0022, reaction\[LetterSpace]48\[LetterSpace]k1 -> 1.6, reaction\[LetterSpace]48\[LetterSpace]k2 -> 1.6, reaction\[LetterSpace]49\[LetterSpace]k1 -> 1900.0, reaction\[LetterSpace]50\[LetterSpace]k1 -> 1.0, reaction\[LetterSpace]50\[LetterSpace]k2 -> 0.01, reaction\[LetterSpace]51\[LetterSpace]k1 -> 0.5, reaction\[LetterSpace]52\[LetterSpace]k1 -> 0.05, reaction\[LetterSpace]52\[LetterSpace]k2 -> 3.0, reaction\[LetterSpace]53\[LetterSpace]k1 -> 0.463, reaction\[LetterSpace]53\[LetterSpace]k2 -> 0.102, reaction\[LetterSpace]54\[LetterSpace]k1 -> 0.463, reaction\[LetterSpace]54\[LetterSpace]k2 -> 0.102, reaction\[LetterSpace]55\[LetterSpace]k1 -> 0.2, reaction\[LetterSpace]55\[LetterSpace]k2 -> 0.2, reaction\[LetterSpace]56\[LetterSpace]k1 -> 0.2, reaction\[LetterSpace]56\[LetterSpace]k2 -> 0.2, reaction\[LetterSpace]57\[LetterSpace]k1 -> 0.2, reaction\[LetterSpace]57\[LetterSpace]k2 -> 0.2, reaction\[LetterSpace]58\[LetterSpace]k1 -> 1.6, reaction\[LetterSpace]58\[LetterSpace]k2 -> 1.6, reaction\[LetterSpace]59\[LetterSpace]k1 -> 0.2, reaction\[LetterSpace]59\[LetterSpace]k2 -> 0.2, reaction\[LetterSpace]60\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]60\[LetterSpace]k2 -> 0.1, reaction\[LetterSpace]61\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]61\[LetterSpace]k2 -> 0.1, reaction\[LetterSpace]62\[LetterSpace]k1 -> 1.0, reaction\[LetterSpace]62\[LetterSpace]k2 -> 1.0, reaction\[LetterSpace]63\[LetterSpace]k1 -> 1.0, reaction\[LetterSpace]63\[LetterSpace]k2 -> 1.0, reaction\[LetterSpace]64\[LetterSpace]k1 -> 0.005, reaction\[LetterSpace]65\[LetterSpace]k1 -> 0.005, reaction\[LetterSpace]66\[LetterSpace]v -> 2*^-05, reaction\[LetterSpace]67\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]68\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]69\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]70\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]70\[LetterSpace]k2 -> 0.1, \[LetterSpace]29\[LetterSpace]k1 -> 0.0492, \[LetterSpace]31\[LetterSpace]Km -> 40.0, \[LetterSpace]31\[LetterSpace]V -> 2.34, \[LetterSpace]33\[LetterSpace]k1 -> 1.6, \[LetterSpace]33\[LetterSpace]k2 -> 1.6, \[LetterSpace]35\[LetterSpace]k1 -> 0.232, \[LetterSpace]37\[LetterSpace]k1 -> 255.068, \[LetterSpace]39\[LetterSpace]Km -> 0.53, \[LetterSpace]39\[LetterSpace]V -> 3.51, \[LetterSpace]41\[LetterSpace]k1 -> 0.2, \[LetterSpace]43\[LetterSpace]k1 -> 2.9, \[LetterSpace]43\[LetterSpace]k2 -> 0.2, \[LetterSpace]77\[LetterSpace]k -> 9.45, \[LetterSpace]86\[LetterSpace]k1 -> 0.03333, \[LetterSpace]86\[LetterSpace]k2 -> 0.03333, \[LetterSpace]103\[LetterSpace]k1 -> 1.0, \[LetterSpace]103\[LetterSpace]k2 -> 0.1, \[LetterSpace]125\[LetterSpace]k -> 1000.0, \[LetterSpace]125\[LetterSpace]km -> 0.0318, \[LetterSpace]132\[LetterSpace]k -> 3.51, \[LetterSpace]132\[LetterSpace]km -> 0.53, \[LetterSpace]139\[LetterSpace]k1 -> 16.6, \[LetterSpace]149\[LetterSpace]k1 -> 1000.0, \[LetterSpace]156\[LetterSpace]k1 -> 0.156, \[LetterSpace]156\[LetterSpace]k2 -> 0.156, \[LetterSpace]164\[LetterSpace]k1 -> 0.16, \[LetterSpace]172\[LetterSpace]k1 -> 0.1, \[LetterSpace]172\[LetterSpace]k2 -> 0.01, \[LetterSpace]178\[LetterSpace]k1 -> 1.0, \[LetterSpace]178\[LetterSpace]k2 -> 0.01, \[LetterSpace]192\[LetterSpace]k1 -> 1900.0, \[LetterSpace]194\[LetterSpace]k1 -> 8.69, \[LetterSpace]194\[LetterSpace]k2 -> 0.01, reaction\[LetterSpace]1\[LetterSpace]v -> 0.00146, reaction\[LetterSpace]2\[LetterSpace]k -> 0.0001, reaction\[LetterSpace]2\[LetterSpace]k1 -> 0.031, reaction\[LetterSpace]3\[LetterSpace]k -> 0.0001, reaction\[LetterSpace]3\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]4\[LetterSpace]k -> 0.000228, reaction\[LetterSpace]4\[LetterSpace]k1 -> 0.0285, reaction\[LetterSpace]5\[LetterSpace]k -> 2*^-05, reaction\[LetterSpace]5\[LetterSpace]k1 -> 0.00055, reaction\[LetterSpace]6\[LetterSpace]v -> 0.0125, reaction\[LetterSpace]7\[LetterSpace]k1 -> 1.0, reaction\[LetterSpace]7\[LetterSpace]k2 -> 0.05288, reaction\[LetterSpace]8\[LetterSpace]k1 -> 0.156, reaction\[LetterSpace]8\[LetterSpace]k2 -> 0.336, reaction\[LetterSpace]9\[LetterSpace]k1 -> 0.027778, reaction\[LetterSpace]10\[LetterSpace]k1 -> 0.027778, reaction\[LetterSpace]11\[LetterSpace]k1 -> 0.027778, reaction\[LetterSpace]12\[LetterSpace]k1 -> 0.03333, reaction\[LetterSpace]13\[LetterSpace]v -> 0.0156, reaction\[LetterSpace]14\[LetterSpace]v -> 0.01183, reaction\[LetterSpace]15\[LetterSpace]k1 -> 0.1266, reaction\[LetterSpace]16\[LetterSpace]k1 -> 0.065, reaction\[LetterSpace]17\[LetterSpace]Km -> 3.51, reaction\[LetterSpace]17\[LetterSpace]V -> 0.53, reaction\[LetterSpace]18\[LetterSpace]k1 -> 0.232, reaction\[LetterSpace]19\[LetterSpace]k1 -> 0.2333, reaction\[LetterSpace]19\[LetterSpace]k2 -> 1.8056, reaction\[LetterSpace]20\[LetterSpace]k1 -> 0.2, reaction\[LetterSpace]21\[LetterSpace]k1 -> 0.1, reaction\[LetterSpace]21\[LetterSpace]k2 -> 0.2, reaction\[LetterSpace]22\[LetterSpace]v -> 0.04528, reaction\[LetterSpace]23\[LetterSpace]k1 -> 0.2, reaction\[LetterSpace]24\[LetterSpace]Km -> 3.51, reaction\[LetterSpace]24\[LetterSpace]V -> 0.53, reaction\[LetterSpace]25\[LetterSpace]Km -> 40.0, reaction\[LetterSpace]25\[LetterSpace]V -> 2.34, reaction\[LetterSpace]26\[LetterSpace]Km -> 0.53, reaction\[LetterSpace]26\[LetterSpace]V -> 3.51, reaction\[LetterSpace]27\[LetterSpace]k -> 1000.0, reaction\[LetterSpace]27\[LetterSpace]km -> 0.0318, reaction\[LetterSpace]28\[LetterSpace]k -> 3.51, reaction\[LetterSpace]28\[LetterSpace]km -> 0.53, reaction\[LetterSpace]29\[LetterSpace]k1 -> 0.2, reaction\[LetterSpace]30\[LetterSpace]k1 -> 0.156, reaction\[LetterSpace]30\[LetterSpace]k2 -> 0.336, \[LetterSpace]1 -> 1.0, \[LetterSpace]3 -> 1.0 }; assignments = { function\[LetterSpace]2[TGFbeta_,k_,s_,s2_] -> k*s*s2*TGFbeta, function\[LetterSpace]6[k_,k1_,m_,m2_,m3_] -> (k + k1*m)/(1 + m2 + m3), function\[LetterSpace]5[k_,k1_,m_] -> k + k1*m, function\[LetterSpace]7[k1_,substrate_,modifier_] -> k1*(1 + modifier)*substrate, function\[LetterSpace]3[k_,rec_,s_,km_] -> (k*rec*s)/(km + s), function\[LetterSpace]1[substrate_,Km_,V_] -> (substrate*V)/(Km + substrate), function\[LetterSpace]4[v_] -> v, parameter\[LetterSpace]2 -> 2*species\[LetterSpace]16[t] + species\[LetterSpace]2[t] + \[LetterSpace]11[t], parameter\[LetterSpace]4 -> species\[LetterSpace]1[t] + \[LetterSpace]15[t] + \[LetterSpace]181[t] + \[LetterSpace]27[t], parameter\[LetterSpace]5 -> species\[LetterSpace]19[t] + \[LetterSpace]105[t] + \[LetterSpace]129[t] + \[LetterSpace]13[t] + \[LetterSpace]153[t] + \[LetterSpace]19[t] + \[LetterSpace]21[t] + \[LetterSpace]9[t] + \[LetterSpace]99[t], parameter\[LetterSpace]6 -> species\[LetterSpace]16[t] + species\[LetterSpace]18[t] + species\[LetterSpace]19[t] + species\[LetterSpace]21[t] + species\[LetterSpace]23[t] + species\[LetterSpace]26[t] + species\[LetterSpace]27[t] + species\[LetterSpace]6[t] + \[LetterSpace]13[t] + \[LetterSpace]147[t] + \[LetterSpace]153[t] + \[LetterSpace]5[t] + \[LetterSpace]9[t], parameter\[LetterSpace]9 -> species\[LetterSpace]6[t] + \[LetterSpace]147[t] + \[LetterSpace]153[t], parameter\[LetterSpace]10 -> species\[LetterSpace]11[t] + species\[LetterSpace]16[t] + 2*species\[LetterSpace]18[t] + species\[LetterSpace]19[t] + species\[LetterSpace]21[t] + species\[LetterSpace]23[t] + 2*species\[LetterSpace]26[t] + species\[LetterSpace]27[t] + \[LetterSpace]13[t] + \[LetterSpace]5[t] + \[LetterSpace]9[t], parameter\[LetterSpace]8 -> species\[LetterSpace]17[t] + species\[LetterSpace]18[t] + species\[LetterSpace]20[t] + species\[LetterSpace]21[t] + species\[LetterSpace]26[t] + species\[LetterSpace]27[t] + species\[LetterSpace]3[t] + species\[LetterSpace]4[t] + species\[LetterSpace]5[t] + species\[LetterSpace]6[t] + species\[LetterSpace]8[t], parameter\[LetterSpace]7 -> species\[LetterSpace]18[t] + species\[LetterSpace]20[t] + species\[LetterSpace]21[t] + species\[LetterSpace]26[t] + species\[LetterSpace]4[t] + species\[LetterSpace]6[t], parameter\[LetterSpace]3 -> \[LetterSpace]174[t] + \[LetterSpace]181[t] + \[LetterSpace]198[t] + \[LetterSpace]25[t] + \[LetterSpace]27[t] }; events = { }; speciesAnnotations = { }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "", "extent" -> ""} }; (* Time evolution *) odes = { \[LetterSpace]101'[t] == 1.0*\[LetterSpace]132 +1.0*reaction\[LetterSpace]2 +1.0*reaction\[LetterSpace]28 -1.0*\[LetterSpace]103 -1.0*reaction\[LetterSpace]16 -1.0*reaction\[LetterSpace]31, \[LetterSpace]105'[t] == 1.0*\[LetterSpace]103 -1.0*\[LetterSpace]132, \[LetterSpace]11'[t] == 1.0*reaction\[LetterSpace]21 -2.0*\[LetterSpace]33 -1.0*reaction\[LetterSpace]7 -2.0*reaction\[LetterSpace]36 -3.0*reaction\[LetterSpace]62 -3.0*reaction\[LetterSpace]63, \[LetterSpace]129'[t] == 1.0*\[LetterSpace]125 +1.0*\[LetterSpace]132 -1.0*\[LetterSpace]139 -2.0*\[LetterSpace]149 -1.0*reaction\[LetterSpace]17, \[LetterSpace]13'[t] == 1.0*\[LetterSpace]33 -1.0*reaction\[LetterSpace]56, \[LetterSpace]147'[t] == 1.0*reaction\[LetterSpace]14 -1.0*\[LetterSpace]149 -1.0*\[LetterSpace]156 -1.0*reaction\[LetterSpace]15 -1.0*reaction\[LetterSpace]34, \[LetterSpace]15'[t] == -3.0*\[LetterSpace]35 -1.0*\[LetterSpace]41 -1.0*\[LetterSpace]43 -1.0*reaction\[LetterSpace]19, \[LetterSpace]153'[t] == 1.0*\[LetterSpace]149 -1.0*\[LetterSpace]164, \[LetterSpace]174'[t] == 1.0*\[LetterSpace]172 +1.0*reaction\[LetterSpace]3 -1.0*\[LetterSpace]194 -1.0*reaction\[LetterSpace]39, \[LetterSpace]181'[t] == 1.0*\[LetterSpace]178 -1.0*\[LetterSpace]192, \[LetterSpace]19'[t] == 3.0*\[LetterSpace]35 +1.0*\[LetterSpace]139 +3.0*reaction\[LetterSpace]68 -2.0*\[LetterSpace]37 -1.0*\[LetterSpace]39 -3.0*reaction\[LetterSpace]62, \[LetterSpace]198'[t] == 1.0*\[LetterSpace]194 , \[LetterSpace]21'[t] == 1.0*\[LetterSpace]39 +1.0*reaction\[LetterSpace]8 +2.0*reaction\[LetterSpace]32 -1.0*\[LetterSpace]41, \[LetterSpace]25'[t] == -1.0*\[LetterSpace]43 -1.0*\[LetterSpace]172 -1.0*reaction\[LetterSpace]46, \[LetterSpace]27'[t] == 1.0*\[LetterSpace]43 -1.0*\[LetterSpace]178, \[LetterSpace]5'[t] == 1.0*\[LetterSpace]29 +1.0*\[LetterSpace]156 +1.0*reaction\[LetterSpace]32 -1.0*\[LetterSpace]37 -1.0*reaction\[LetterSpace]7 -1.0*reaction\[LetterSpace]35 -1.0*reaction\[LetterSpace]44, \[LetterSpace]75'[t] == 1.0*reaction\[LetterSpace]1 +2.0*reaction\[LetterSpace]12 -2.0*\[LetterSpace]77 -1.0*reaction\[LetterSpace]11, \[LetterSpace]79'[t] == 1.0*\[LetterSpace]77 -1.0*\[LetterSpace]86 -1.0*reaction\[LetterSpace]12, \[LetterSpace]84'[t] == 1.0*reaction\[LetterSpace]6 -2.0*\[LetterSpace]86 -1.0*reaction\[LetterSpace]10, \[LetterSpace]9'[t] == 1.0*\[LetterSpace]37 +1.0*\[LetterSpace]164 -1.0*\[LetterSpace]31 -1.0*\[LetterSpace]33 -1.0*reaction\[LetterSpace]54 -1.0*reaction\[LetterSpace]58 -1.0*reaction\[LetterSpace]64, \[LetterSpace]96'[t] == 1.0*\[LetterSpace]86 -1.0*\[LetterSpace]192 -1.0*\[LetterSpace]194 -1.0*reaction\[LetterSpace]9 -1.0*reaction\[LetterSpace]49, \[LetterSpace]99'[t] == 1.0*reaction\[LetterSpace]13 +1.0*reaction\[LetterSpace]17 -1.0*\[LetterSpace]103 -1.0*\[LetterSpace]125 -1.0*reaction\[LetterSpace]8 -1.0*reaction\[LetterSpace]20 -1.0*reaction\[LetterSpace]61, species\[LetterSpace]1'[t] == 1.0*reaction\[LetterSpace]4 +1.0*reaction\[LetterSpace]19 -1.0*reaction\[LetterSpace]43, species\[LetterSpace]10'[t] == 1.0*reaction\[LetterSpace]50 -1.0*reaction\[LetterSpace]49, species\[LetterSpace]11'[t] == 1.0*reaction\[LetterSpace]58 -1.0*reaction\[LetterSpace]59, species\[LetterSpace]12'[t] == 1.0*reaction\[LetterSpace]61 , species\[LetterSpace]13'[t] == 1.0*reaction\[LetterSpace]62 -1.0*\[LetterSpace]35 -1.0*reaction\[LetterSpace]68, species\[LetterSpace]14'[t] == 1.0*reaction\[LetterSpace]63 -1.0*reaction\[LetterSpace]38 -1.0*reaction\[LetterSpace]69, species\[LetterSpace]15'[t] == 1.0*reaction\[LetterSpace]66 -1.0*reaction\[LetterSpace]67 -1.0*reaction\[LetterSpace]70, species\[LetterSpace]16'[t] == 1.0*reaction\[LetterSpace]7 , species\[LetterSpace]17'[t] == 2.0*\[LetterSpace]29 +1.0*reaction\[LetterSpace]26 +1.0*reaction\[LetterSpace]30 -1.0*reaction\[LetterSpace]29, species\[LetterSpace]18'[t] == 1.0*reaction\[LetterSpace]35 +1.0*reaction\[LetterSpace]37 -1.0*reaction\[LetterSpace]25 -1.0*reaction\[LetterSpace]36 -1.0*reaction\[LetterSpace]48 -1.0*reaction\[LetterSpace]53 -1.0*reaction\[LetterSpace]65, species\[LetterSpace]19'[t] == 1.0*\[LetterSpace]31 -1.0*reaction\[LetterSpace]32, species\[LetterSpace]2'[t] == 1.0*reaction\[LetterSpace]5 -1.0*reaction\[LetterSpace]18 -1.0*reaction\[LetterSpace]21, species\[LetterSpace]20'[t] == 1.0*reaction\[LetterSpace]33 +1.0*reaction\[LetterSpace]38 +1.0*reaction\[LetterSpace]69 -1.0*reaction\[LetterSpace]26 -2.0*reaction\[LetterSpace]35 -3.0*reaction\[LetterSpace]63, species\[LetterSpace]21'[t] == 1.0*reaction\[LetterSpace]36 -1.0*reaction\[LetterSpace]55, species\[LetterSpace]22'[t] == 1.0*reaction\[LetterSpace]41 -1.0*reaction\[LetterSpace]44 -2.0*reaction\[LetterSpace]48 -2.0*reaction\[LetterSpace]58, species\[LetterSpace]23'[t] == 1.0*reaction\[LetterSpace]44 -1.0*reaction\[LetterSpace]60, species\[LetterSpace]24'[t] == -1.0*reaction\[LetterSpace]46 -1.0*reaction\[LetterSpace]52, species\[LetterSpace]25'[t] == 1.0*reaction\[LetterSpace]46 -1.0*reaction\[LetterSpace]50, species\[LetterSpace]26'[t] == 1.0*reaction\[LetterSpace]48 -1.0*reaction\[LetterSpace]57, species\[LetterSpace]27'[t] == 1.0*reaction\[LetterSpace]25 -1.0*\[LetterSpace]29, species\[LetterSpace]28'[t] == -1.0*reaction\[LetterSpace]53 -1.0*reaction\[LetterSpace]54 -1.0*reaction\[LetterSpace]55 -1.0*reaction\[LetterSpace]56 -1.0*reaction\[LetterSpace]57 -1.0*reaction\[LetterSpace]59 -1.0*reaction\[LetterSpace]60, species\[LetterSpace]29'[t] == 1.0*reaction\[LetterSpace]55 +1.0*reaction\[LetterSpace]56 +1.0*reaction\[LetterSpace]57 +1.0*reaction\[LetterSpace]59 +1.0*reaction\[LetterSpace]60 , species\[LetterSpace]3'[t] == 1.0*reaction\[LetterSpace]22 +1.0*reaction\[LetterSpace]24 -1.0*reaction\[LetterSpace]23 -1.0*reaction\[LetterSpace]27 -1.0*reaction\[LetterSpace]30 -1.0*reaction\[LetterSpace]31 -1.0*reaction\[LetterSpace]45, species\[LetterSpace]30'[t] == 1.0*reaction\[LetterSpace]53 +1.0*reaction\[LetterSpace]54 , species\[LetterSpace]31'[t] == 1.0*reaction\[LetterSpace]70 -3.0*reaction\[LetterSpace]68 -3.0*reaction\[LetterSpace]69, species\[LetterSpace]4'[t] == 1.0*reaction\[LetterSpace]27 +1.0*reaction\[LetterSpace]28 -1.0*reaction\[LetterSpace]24 -1.0*reaction\[LetterSpace]33 -2.0*reaction\[LetterSpace]34, species\[LetterSpace]5'[t] == 1.0*reaction\[LetterSpace]31 -1.0*reaction\[LetterSpace]28, species\[LetterSpace]6'[t] == 1.0*reaction\[LetterSpace]34 -1.0*reaction\[LetterSpace]37, species\[LetterSpace]7'[t] == 1.0*reaction\[LetterSpace]42 -1.0*reaction\[LetterSpace]40 -1.0*reaction\[LetterSpace]41 -1.0*reaction\[LetterSpace]45 -1.0*reaction\[LetterSpace]61, species\[LetterSpace]8'[t] == 1.0*reaction\[LetterSpace]45 , species\[LetterSpace]9'[t] == 1.0*reaction\[LetterSpace]47 +1.0*reaction\[LetterSpace]52 -1.0*reaction\[LetterSpace]51 }; 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]}]