(* 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 = { s1[t], s10[t], s11[t], s12[t], s13[t], s14[t], s15[t], s16[t], s17[t], s18[t], s19[t], s2[t], s20[t], s21[t], s22[t], s23[t], s24[t], s25[t], s26[t], s27[t], s28[t], s29[t], s3[t], s30[t], s31[t], s32[t], s33[t], s34[t], s35[t], s36[t], s37[t], s38[t], s39[t], s4[t], s40[t], s41[t], s42[t], s43[t], s44[t], s45[t], s46[t], s47[t], s48[t], s49[t], s5[t], s50[t], s51[t], s52[t], s6[t], s7[t], s8[t], s9[t] }; initialValues = { s1[0] == 0.5, s10[0] == 1.0, s11[0] == 1.0, s12[0] == 1.0, s13[0] == 1.0, s14[0] == 1.0, s15[0] == 1.0, s16[0] == 1.0, s17[0] == 1.0, s18[0] == 1.0, s19[0] == 1.0, s2[0] == 0.5, s20[0] == 1.0, s21[0] == 1.0, s22[0] == 1.0, s23[0] == 1.0, s24[0] == 1.0, s25[0] == 1.0, s26[0] == 1.0, s27[0] == 1.0, s28[0] == 1.2, s29[0] == 1.2, s3[0] == 0.8, s30[0] == 1.5, s31[0] == 1.5, s32[0] == 1.2, s33[0] == 1.2, s34[0] == 1.2, s35[0] == 1.2, s36[0] == 1.5, s37[0] == 1.5, s38[0] == 1.2, s39[0] == 1.2, s4[0] == 0.8, s40[0] == 1.2, s41[0] == 1.2, s42[0] == 1.2, s43[0] == 1.2, s44[0] == 1.2, s45[0] == 1.2, s46[0] == 1.2, s47[0] == 1.2, s48[0] == 1.5, s49[0] == 1.8, s5[0] == 0.8, s50[0] == 2.0, s51[0] == 2.2, s52[0] == 2.5, s6[0] == 0.8, s7[0] == 1.0, s8[0] == 1.0, s9[0] == 1.0 }; rates = { re1, re10, re11, re12, re13, re14, re15, re16, re17, re18, re19, re2, re20, re21, re22, re23, re24, re25, re26, re27, re28, re29, re3, re30, re31, re32, re33, re34, re35, re36, re37, re38, re39, re4, re40, re41, re42, re43, re44, re45, re46, re47, re48, re49, re5, re50, re51, re52, re53, re54, re55, re56, re57, re58, re59, re6, re60, re61, re62, re63, re64, re65, re66, re67, re68, re69, re7, re70, re71, re72, re73, re74, re75, re76, re77, re78, re79, re8, re80, re81, re82, re83, re84, re85, re86, re87, re88, re9 }; rateEquations = { re1 -> kass\[LetterSpace]re1*s1[t] - kdiss\[LetterSpace]re1*s3[t], re10 -> kass\[LetterSpace]re10*s5[t] - kdiss\[LetterSpace]re10*s7[t], re11 -> kass\[LetterSpace]re11*s6[t] - kdiss\[LetterSpace]re11*s7[t], re12 -> kass\[LetterSpace]re12*s8[t] - kdiss\[LetterSpace]re12*s9[t], re13 -> -(kdiss\[LetterSpace]re13*s10[t]) + kass\[LetterSpace]re13*s8[t], re14 -> -(kdiss\[LetterSpace]re14*s11[t]) + kass\[LetterSpace]re14*s8[t], re15 -> -(kdiss\[LetterSpace]re15*s12[t]) + kass\[LetterSpace]re15*s8[t], re16 -> kass\[LetterSpace]re16*s6[t] - kdiss\[LetterSpace]re16*s9[t], re17 -> -(kdiss\[LetterSpace]re17*s13[t]) + kass\[LetterSpace]re17*s8[t], re18 -> -(kdiss\[LetterSpace]re18*s13[t]) + kass\[LetterSpace]re18*s5[t], re19 -> -(kdiss\[LetterSpace]re19*s11[t]) + kass\[LetterSpace]re19*s5[t], re2 -> kass\[LetterSpace]re2*s1[t] - kdiss\[LetterSpace]re2*s4[t], re20 -> kass\[LetterSpace]re20*s14[t] - kdiss\[LetterSpace]re20*s15[t], re21 -> -(kdiss\[LetterSpace]re21*s14[t]) + kass\[LetterSpace]re21*s8[t], re22 -> kass\[LetterSpace]re22*s15[t] - kdiss\[LetterSpace]re22*s16[t], re23 -> kass\[LetterSpace]re23*s15[t] - kdiss\[LetterSpace]re23*s17[t], re24 -> kass\[LetterSpace]re24*s15[t] - kdiss\[LetterSpace]re24*s18[t], re25 -> kass\[LetterSpace]re25*s15[t] - kdiss\[LetterSpace]re25*s19[t], re26 -> -(kdiss\[LetterSpace]re26*s16[t]) + kass\[LetterSpace]re26*s9[t], re27 -> -(kdiss\[LetterSpace]re27*s17[t]) + kass\[LetterSpace]re27*s9[t], re28 -> -(kdiss\[LetterSpace]re28*s18[t]) + kass\[LetterSpace]re28*s9[t], re29 -> kass\[LetterSpace]re29*s11[t] - kdiss\[LetterSpace]re29*s19[t], re3 -> kass\[LetterSpace]re3*s1[t] - kdiss\[LetterSpace]re3*s5[t], re30 -> kass\[LetterSpace]re30*s20[t] - kdiss\[LetterSpace]re30*s21[t], re31 -> kass\[LetterSpace]re31*s15[t] - kdiss\[LetterSpace]re31*s20[t], re32 -> kass\[LetterSpace]re32*s21[t] - kdiss\[LetterSpace]re32*s22[t], re33 -> kass\[LetterSpace]re33*s21[t] - kdiss\[LetterSpace]re33*s23[t], re34 -> kass\[LetterSpace]re34*s21[t] - kdiss\[LetterSpace]re34*s24[t], re35 -> kass\[LetterSpace]re35*s21[t] - kdiss\[LetterSpace]re35*s25[t], re36 -> kass\[LetterSpace]re36*s21[t] - kdiss\[LetterSpace]re36*s26[t], re37 -> kass\[LetterSpace]re37*s21[t] - kdiss\[LetterSpace]re37*s27[t], re38 -> kass\[LetterSpace]re38*s16[t] - kdiss\[LetterSpace]re38*s22[t], re39 -> kass\[LetterSpace]re39*s17[t] - kdiss\[LetterSpace]re39*s22[t], re4 -> kass\[LetterSpace]re4*s2[t] - kdiss\[LetterSpace]re4*s4[t], re40 -> kass\[LetterSpace]re40*s17[t] - kdiss\[LetterSpace]re40*s23[t], re41 -> kass\[LetterSpace]re41*s18[t] - kdiss\[LetterSpace]re41*s23[t], re42 -> kass\[LetterSpace]re42*s17[t] - kdiss\[LetterSpace]re42*s25[t], re43 -> kass\[LetterSpace]re43*s16[t] - kdiss\[LetterSpace]re43*s24[t], re44 -> kass\[LetterSpace]re44*s18[t] - kdiss\[LetterSpace]re44*s25[t], re45 -> kass\[LetterSpace]re45*s28[t] - kdiss\[LetterSpace]re45*s29[t], re46 -> kass\[LetterSpace]re46*s30[t] - kdiss\[LetterSpace]re46*s31[t], re47 -> kass\[LetterSpace]re47*s32[t] - kdiss\[LetterSpace]re47*s33[t], re48 -> kass\[LetterSpace]re48*s34[t] - kdiss\[LetterSpace]re48*s35[t], re49 -> kass\[LetterSpace]re49*s36[t] - kdiss\[LetterSpace]re49*s37[t], re5 -> kass\[LetterSpace]re5*s2[t] - kdiss\[LetterSpace]re5*s5[t], re50 -> kass\[LetterSpace]re50*s38[t] - kdiss\[LetterSpace]re50*s39[t], re51 -> kass\[LetterSpace]re51*s40[t] - kdiss\[LetterSpace]re51*s41[t], re52 -> kass\[LetterSpace]re52*s42[t] - kdiss\[LetterSpace]re52*s43[t], re53 -> kass\[LetterSpace]re53*s44[t] - kdiss\[LetterSpace]re53*s45[t], re54 -> kass\[LetterSpace]re54*s46[t] - kdiss\[LetterSpace]re54*s47[t], re55 -> kass\[LetterSpace]re55*s22[t] - kdiss\[LetterSpace]re55*s28[t], re56 -> kass\[LetterSpace]re56*s24[t] - kdiss\[LetterSpace]re56*s28[t], re57 -> kass\[LetterSpace]re57*s23[t] - kdiss\[LetterSpace]re57*s32[t], re58 -> kass\[LetterSpace]re58*s24[t] - kdiss\[LetterSpace]re58*s32[t], re59 -> kass\[LetterSpace]re59*s24[t] - kdiss\[LetterSpace]re59*s38[t], re6 -> kass\[LetterSpace]re6*s2[t] - kdiss\[LetterSpace]re6*s6[t], re60 -> kass\[LetterSpace]re60*s23[t] - kdiss\[LetterSpace]re60*s42[t], re61 -> kass\[LetterSpace]re61*s25[t] - kdiss\[LetterSpace]re61*s42[t], re62 -> kass\[LetterSpace]re62*s25[t] - kdiss\[LetterSpace]re62*s46[t], re63 -> kass\[LetterSpace]re63*s25[t] - kdiss\[LetterSpace]re63*s32[t], re64 -> kass\[LetterSpace]re64*s23[t] - kdiss\[LetterSpace]re64*s46[t], re65 -> kass\[LetterSpace]re65*s23[t] - kdiss\[LetterSpace]re65*s44[t], re66 -> kass\[LetterSpace]re66*s25[t] - kdiss\[LetterSpace]re66*s44[t], re67 -> kass\[LetterSpace]re67*s23[t] - kdiss\[LetterSpace]re67*s40[t], re68 -> kass\[LetterSpace]re68*s25[t] - kdiss\[LetterSpace]re68*s36[t], re69 -> kass\[LetterSpace]re69*s21[t] - kdiss\[LetterSpace]re69*s30[t], re7 -> kass\[LetterSpace]re7*s7[t] - kdiss\[LetterSpace]re7*s8[t], re70 -> kass\[LetterSpace]re70*s21[t] - kdiss\[LetterSpace]re70*s48[t], re71 -> kass\[LetterSpace]re71*s21[t] - kdiss\[LetterSpace]re71*s49[t], re72 -> kass\[LetterSpace]re72*s21[t] - kdiss\[LetterSpace]re72*s51[t], re73 -> kass\[LetterSpace]re73*s21[t] - kdiss\[LetterSpace]re73*s50[t], re74 -> kass\[LetterSpace]re74*s24[t] - kdiss\[LetterSpace]re74*s34[t], re75 -> kass\[LetterSpace]re75*s45[t] - kdiss\[LetterSpace]re75*s52[t], re76 -> kass\[LetterSpace]re76*s31[t] - kdiss\[LetterSpace]re76*s52[t], re77 -> kass\[LetterSpace]re77*s49[t] - kdiss\[LetterSpace]re77*s52[t], re78 -> kass\[LetterSpace]re78*s50[t] - kdiss\[LetterSpace]re78*s52[t], re79 -> kass\[LetterSpace]re79*s51[t] - kdiss\[LetterSpace]re79*s52[t], re8 -> kass\[LetterSpace]re8*s3[t] - kdiss\[LetterSpace]re8*s7[t], re80 -> kass\[LetterSpace]re80*s48[t] - kdiss\[LetterSpace]re80*s52[t], re81 -> kass\[LetterSpace]re81*s37[t] - kdiss\[LetterSpace]re81*s52[t], re82 -> kass\[LetterSpace]re82*s41[t] - kdiss\[LetterSpace]re82*s52[t], re83 -> kass\[LetterSpace]re83*s35[t] - kdiss\[LetterSpace]re83*s52[t], re84 -> kass\[LetterSpace]re84*s29[t] - kdiss\[LetterSpace]re84*s52[t], re85 -> kass\[LetterSpace]re85*s47[t] - kdiss\[LetterSpace]re85*s52[t], re86 -> kass\[LetterSpace]re86*s43[t] - kdiss\[LetterSpace]re86*s52[t], re87 -> kass\[LetterSpace]re87*s39[t] - kdiss\[LetterSpace]re87*s52[t], re88 -> kass\[LetterSpace]re88*s33[t] - kdiss\[LetterSpace]re88*s52[t], re9 -> kass\[LetterSpace]re9*s4[t] - kdiss\[LetterSpace]re9*s7[t] }; parameters = { kass\[LetterSpace]re1 -> 1.0, kass\[LetterSpace]re10 -> 1.0, kass\[LetterSpace]re11 -> 1.0, kass\[LetterSpace]re12 -> 1.0, kass\[LetterSpace]re13 -> 1.0, kass\[LetterSpace]re14 -> 1.0, kass\[LetterSpace]re15 -> 1.0, kass\[LetterSpace]re16 -> 1.0, kass\[LetterSpace]re17 -> 1.0, kass\[LetterSpace]re18 -> 1.0, kass\[LetterSpace]re19 -> 1.0, kass\[LetterSpace]re2 -> 1.0, kass\[LetterSpace]re20 -> 1.0, kass\[LetterSpace]re21 -> 1.0, kass\[LetterSpace]re22 -> 1.0, kass\[LetterSpace]re23 -> 1.0, kass\[LetterSpace]re24 -> 1.0, kass\[LetterSpace]re25 -> 1.0, kass\[LetterSpace]re26 -> 1.0, kass\[LetterSpace]re27 -> 1.0, kass\[LetterSpace]re28 -> 1.0, kass\[LetterSpace]re29 -> 1.0, kass\[LetterSpace]re3 -> 1.0, kass\[LetterSpace]re30 -> 1.0, kass\[LetterSpace]re31 -> 1.0, kass\[LetterSpace]re32 -> 1.0, kass\[LetterSpace]re33 -> 1.0, kass\[LetterSpace]re34 -> 1.0, kass\[LetterSpace]re35 -> 1.0, kass\[LetterSpace]re36 -> 1.0, kass\[LetterSpace]re37 -> 1.0, kass\[LetterSpace]re38 -> 1.0, kass\[LetterSpace]re39 -> 1.0, kass\[LetterSpace]re4 -> 1.0, kass\[LetterSpace]re40 -> 1.0, kass\[LetterSpace]re41 -> 1.0, kass\[LetterSpace]re42 -> 1.0, kass\[LetterSpace]re43 -> 1.0, kass\[LetterSpace]re44 -> 1.0, kass\[LetterSpace]re45 -> 1.0, kass\[LetterSpace]re46 -> 1.0, kass\[LetterSpace]re47 -> 1.0, kass\[LetterSpace]re48 -> 1.0, kass\[LetterSpace]re49 -> 1.0, kass\[LetterSpace]re5 -> 1.0, kass\[LetterSpace]re50 -> 1.0, kass\[LetterSpace]re51 -> 1.0, kass\[LetterSpace]re52 -> 1.0, kass\[LetterSpace]re53 -> 1.0, kass\[LetterSpace]re54 -> 1.0, kass\[LetterSpace]re55 -> 1.0, kass\[LetterSpace]re56 -> 1.0, kass\[LetterSpace]re57 -> 1.0, kass\[LetterSpace]re58 -> 1.0, kass\[LetterSpace]re59 -> 1.0, kass\[LetterSpace]re6 -> 1.0, kass\[LetterSpace]re60 -> 1.0, kass\[LetterSpace]re61 -> 1.0, kass\[LetterSpace]re62 -> 1.0, kass\[LetterSpace]re63 -> 1.0, kass\[LetterSpace]re64 -> 1.0, kass\[LetterSpace]re65 -> 1.0, kass\[LetterSpace]re66 -> 1.0, kass\[LetterSpace]re67 -> 1.0, kass\[LetterSpace]re68 -> 1.0, kass\[LetterSpace]re69 -> 1.0, kass\[LetterSpace]re7 -> 1.0, kass\[LetterSpace]re70 -> 1.0, kass\[LetterSpace]re71 -> 1.0, kass\[LetterSpace]re72 -> 1.0, kass\[LetterSpace]re73 -> 1.0, kass\[LetterSpace]re74 -> 1.0, kass\[LetterSpace]re75 -> 1.0, kass\[LetterSpace]re76 -> 1.0, kass\[LetterSpace]re77 -> 1.0, kass\[LetterSpace]re78 -> 1.0, kass\[LetterSpace]re79 -> 1.0, kass\[LetterSpace]re8 -> 1.0, kass\[LetterSpace]re80 -> 1.0, kass\[LetterSpace]re81 -> 1.0, kass\[LetterSpace]re82 -> 1.0, kass\[LetterSpace]re83 -> 1.0, kass\[LetterSpace]re84 -> 1.0, kass\[LetterSpace]re85 -> 1.0, kass\[LetterSpace]re86 -> 1.0, kass\[LetterSpace]re87 -> 1.0, kass\[LetterSpace]re88 -> 1.0, kass\[LetterSpace]re9 -> 1.0, kdiss\[LetterSpace]re1 -> 1.0, kdiss\[LetterSpace]re10 -> 1.0, kdiss\[LetterSpace]re11 -> 1.0, kdiss\[LetterSpace]re12 -> 1.0, kdiss\[LetterSpace]re13 -> 1.0, kdiss\[LetterSpace]re14 -> 1.0, kdiss\[LetterSpace]re15 -> 1.0, kdiss\[LetterSpace]re16 -> 1.0, kdiss\[LetterSpace]re17 -> 1.0, kdiss\[LetterSpace]re18 -> 1.0, kdiss\[LetterSpace]re19 -> 1.0, kdiss\[LetterSpace]re2 -> 1.0, kdiss\[LetterSpace]re20 -> 1.0, kdiss\[LetterSpace]re21 -> 1.0, kdiss\[LetterSpace]re22 -> 1.0, kdiss\[LetterSpace]re23 -> 1.0, kdiss\[LetterSpace]re24 -> 1.0, kdiss\[LetterSpace]re25 -> 1.0, kdiss\[LetterSpace]re26 -> 1.0, kdiss\[LetterSpace]re27 -> 1.0, kdiss\[LetterSpace]re28 -> 1.0, kdiss\[LetterSpace]re29 -> 1.0, kdiss\[LetterSpace]re3 -> 1.0, kdiss\[LetterSpace]re30 -> 1.0, kdiss\[LetterSpace]re31 -> 1.0, kdiss\[LetterSpace]re32 -> 1.0, kdiss\[LetterSpace]re33 -> 1.0, kdiss\[LetterSpace]re34 -> 1.0, kdiss\[LetterSpace]re35 -> 1.0, kdiss\[LetterSpace]re36 -> 1.0, kdiss\[LetterSpace]re37 -> 1.0, kdiss\[LetterSpace]re38 -> 1.0, kdiss\[LetterSpace]re39 -> 1.0, kdiss\[LetterSpace]re4 -> 1.0, kdiss\[LetterSpace]re40 -> 1.0, kdiss\[LetterSpace]re41 -> 1.0, kdiss\[LetterSpace]re42 -> 1.0, kdiss\[LetterSpace]re43 -> 1.0, kdiss\[LetterSpace]re44 -> 1.0, kdiss\[LetterSpace]re45 -> 1.0, kdiss\[LetterSpace]re46 -> 1.0, kdiss\[LetterSpace]re47 -> 1.0, kdiss\[LetterSpace]re48 -> 1.0, kdiss\[LetterSpace]re49 -> 1.0, kdiss\[LetterSpace]re5 -> 1.0, kdiss\[LetterSpace]re50 -> 1.0, kdiss\[LetterSpace]re51 -> 1.0, kdiss\[LetterSpace]re52 -> 1.0, kdiss\[LetterSpace]re53 -> 1.0, kdiss\[LetterSpace]re54 -> 1.0, kdiss\[LetterSpace]re55 -> 1.0, kdiss\[LetterSpace]re56 -> 1.0, kdiss\[LetterSpace]re57 -> 1.0, kdiss\[LetterSpace]re58 -> 1.0, kdiss\[LetterSpace]re59 -> 1.0, kdiss\[LetterSpace]re6 -> 1.0, kdiss\[LetterSpace]re60 -> 1.0, kdiss\[LetterSpace]re61 -> 1.0, kdiss\[LetterSpace]re62 -> 1.0, kdiss\[LetterSpace]re63 -> 1.0, kdiss\[LetterSpace]re64 -> 1.0, kdiss\[LetterSpace]re65 -> 1.0, kdiss\[LetterSpace]re66 -> 1.0, kdiss\[LetterSpace]re67 -> 1.0, kdiss\[LetterSpace]re68 -> 1.0, kdiss\[LetterSpace]re69 -> 1.0, kdiss\[LetterSpace]re7 -> 1.0, kdiss\[LetterSpace]re70 -> 1.0, kdiss\[LetterSpace]re71 -> 1.0, kdiss\[LetterSpace]re72 -> 1.0, kdiss\[LetterSpace]re73 -> 1.0, kdiss\[LetterSpace]re74 -> 1.0, kdiss\[LetterSpace]re75 -> 1.0, kdiss\[LetterSpace]re76 -> 1.0, kdiss\[LetterSpace]re77 -> 1.0, kdiss\[LetterSpace]re78 -> 1.0, kdiss\[LetterSpace]re79 -> 1.0, kdiss\[LetterSpace]re8 -> 1.0, kdiss\[LetterSpace]re80 -> 1.0, kdiss\[LetterSpace]re81 -> 1.0, kdiss\[LetterSpace]re82 -> 1.0, kdiss\[LetterSpace]re83 -> 1.0, kdiss\[LetterSpace]re84 -> 1.0, kdiss\[LetterSpace]re85 -> 1.0, kdiss\[LetterSpace]re86 -> 1.0, kdiss\[LetterSpace]re87 -> 1.0, kdiss\[LetterSpace]re88 -> 1.0, kdiss\[LetterSpace]re9 -> 1.0, c1 -> 1.0, c2 -> 1.0, default -> 1.0 }; assignments = { }; events = { }; speciesAnnotations = { }; reactionAnnotations = { }; units = { {"time" -> "", "metabolite" -> "substance", "extent" -> "substance"} }; (* Time evolution *) odes = { s1'[t] == -1.0*re1 -1.0*re2 -1.0*re3, s10'[t] == 1.0*re13 , s11'[t] == 1.0*re14 +1.0*re19 -1.0*re29, s12'[t] == 1.0*re15 , s13'[t] == 1.0*re17 +1.0*re18 , s14'[t] == 1.0*re21 -1.0*re20, s15'[t] == 1.0*re20 -1.0*re22 -1.0*re23 -1.0*re24 -1.0*re25 -1.0*re31, s16'[t] == 1.0*re22 +1.0*re26 -1.0*re38 -1.0*re43, s17'[t] == 1.0*re23 +1.0*re27 -1.0*re39 -1.0*re40 -1.0*re42, s18'[t] == 1.0*re24 +1.0*re28 -1.0*re41 -1.0*re44, s19'[t] == 1.0*re25 +1.0*re29 , s2'[t] == -1.0*re4 -1.0*re5 -1.0*re6, s20'[t] == 1.0*re31 -1.0*re30, s21'[t] == 1.0*re30 -1.0*re32 -1.0*re33 -1.0*re34 -1.0*re35 -1.0*re36 -1.0*re37 -1.0*re69 -1.0*re70 -1.0*re71 -1.0*re72 -1.0*re73, s22'[t] == 1.0*re32 +1.0*re38 +1.0*re39 -1.0*re55, s23'[t] == 1.0*re33 +1.0*re40 +1.0*re41 -1.0*re57 -1.0*re60 -1.0*re64 -1.0*re65 -1.0*re67, s24'[t] == 1.0*re34 +1.0*re43 -1.0*re56 -1.0*re58 -1.0*re59 -1.0*re74, s25'[t] == 1.0*re35 +1.0*re42 +1.0*re44 -1.0*re61 -1.0*re62 -1.0*re63 -1.0*re66 -1.0*re68, s26'[t] == 1.0*re36 , s27'[t] == 1.0*re37 , s28'[t] == 1.0*re55 +1.0*re56 -1.0*re45, s29'[t] == 1.0*re45 -1.0*re84, s3'[t] == 1.0*re1 -1.0*re8, s30'[t] == 1.0*re69 -1.0*re46, s31'[t] == 1.0*re46 -1.0*re76, s32'[t] == 1.0*re57 +1.0*re58 +1.0*re63 -1.0*re47, s33'[t] == 1.0*re47 -1.0*re88, s34'[t] == 1.0*re74 -1.0*re48, s35'[t] == 1.0*re48 -1.0*re83, s36'[t] == 1.0*re68 -1.0*re49, s37'[t] == 1.0*re49 -1.0*re81, s38'[t] == 1.0*re59 -1.0*re50, s39'[t] == 1.0*re50 -1.0*re87, s4'[t] == 1.0*re2 +1.0*re4 -1.0*re9, s40'[t] == 1.0*re67 -1.0*re51, s41'[t] == 1.0*re51 -1.0*re82, s42'[t] == 1.0*re60 +1.0*re61 -1.0*re52, s43'[t] == 1.0*re52 -1.0*re86, s44'[t] == 1.0*re65 +1.0*re66 -1.0*re53, s45'[t] == 1.0*re53 -1.0*re75, s46'[t] == 1.0*re62 +1.0*re64 -1.0*re54, s47'[t] == 1.0*re54 -1.0*re85, s48'[t] == 1.0*re70 -1.0*re80, s49'[t] == 1.0*re71 -1.0*re77, s5'[t] == 1.0*re3 +1.0*re5 -1.0*re10 -1.0*re18 -1.0*re19, s50'[t] == 1.0*re73 -1.0*re78, s51'[t] == 1.0*re72 -1.0*re79, s52'[t] == 1.0*re75 +1.0*re76 +1.0*re77 +1.0*re78 +1.0*re79 +1.0*re80 +1.0*re81 +1.0*re82 +1.0*re83 +1.0*re84 +1.0*re85 +1.0*re86 +1.0*re87 +1.0*re88 , s6'[t] == 1.0*re6 -1.0*re11 -1.0*re16, s7'[t] == 1.0*re8 +1.0*re9 +1.0*re10 +1.0*re11 -1.0*re7, s8'[t] == 1.0*re7 -1.0*re12 -1.0*re13 -1.0*re14 -1.0*re15 -1.0*re17 -1.0*re21, s9'[t] == 1.0*re12 +1.0*re16 -1.0*re26 -1.0*re27 -1.0*re28 }; 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]}]