{ diff(x1, t) = (-1)*k14*(x1*x55 - k15*x3)/k115,
  diff(x2, t) = 1*k93*(x28*k12 - k94*x2)/k115,
  diff(x3, t) = (1*k14*(x1*x55 - k15*x3) + (-1)*k16*(x3*x4 - k17*x5) + (-1)*k18*x3)/k115,
  diff(x4, t) = ((-1)*k16*(x3*x4 - k17*x5) + (-1)*k87*(x48*x4 - k88*x12) + (-1)*k16*(x27*x4 - k17*x5))/k115,
  diff(x5, t) = (1*k16*(x3*x4 - k17*x5) + (-1)*k20*(x5 - k21*x6) + 1*k22*x6/(k23 + x6) + 1*k16*(x27*x4 - k17*x5) + (-1)*k19*x5)/k115,
  diff(x6, t) = (1*k20*(x5 - k21*x6) + (-1)*k22*x6/(k23 + x6) + (-1)*k24*(x6*x7 - k25*x8) + 1*k30*(x11 - k31*x6*x49) + (-1)*k36*(x6*x28 - k37*x29) + 1*(k40*x42 - k41*x6*x31) + (-1)*k86*x6 + 1*k20*(x17 - k21*x6))/k115,
  diff(x7, t) = ((-1)*k24*(x6*x7 - k25*x8) + 1*k34*x54/(k35 + x54))/k115,
  diff(x8, t) = (1*k24*(x6*x7 - k25*x8) + (-1)*(k26*x8 - k27*x9))/k115,
  diff(x9, t) = (1*(k26*x8 - k27*x9) + (-1)*k28*(x9*x10 - k29*x11))/k115,
  diff(x10, t) = ((-1)*k28*(x9*x10 - k29*x11) + 1*k32*(x49 - k33*x54*x10))/k115,
  diff(x11, t) = (1*k28*(x9*x10 - k29*x11) + (-1)*k30*(x11 - k31*x6*x49))/k115,
  diff(x12, t) = (1*k87*(x48*x4 - k88*x12) + (-1)*(k89*x12 - k90*x50))/k115,
  diff(x13, t) = 1*k91*(x34*k11 - k92*x13)/k115,
  diff(x14, t) = ((-1)*k43*x14*x49/(k44 + x14) + 1*k45*x53/(k46 + x53))/k115,
  diff(x15, t) = (1*k65*x44 + (-1)*k66*x15)/k115,
  diff(x16, t) = ((-1)*k47*x16*x53/(k48 + x16) + 1*k49*x52*(x46 + k51)/(x52 + k50))/k115,
  diff(x17, t) = (1*k19*x5 + (-1)*k20*(x17 - k21*x6))/k115,
  diff(x18, t) = ((-1)*k52*x18*x52/(k53 + x18) + 1*k57*x22)/k115,
  diff(x19, t) = (1*k52*x18*x52/(k53 + x18) + (-1)*k54*x19*x20 + (-1)*k52*x19*x52/(k53 + x19) + 1*k58*x21)/k115,
  diff(x20, t) = ((-1)*k54*x19*x20 + 1*k57*x22 + (-1)*k54*(x20*x45 - k55*x24) + 1*k58*x21 + (-1)*k81*x40*x20 + 1*k84*x23 + (-1)*k81*(x46*x20 - k82*x47) + 1*k85*x33)/k115,
  diff(x21, t) = (1*k54*x19*x20 + (-1)*k56*x21 + 1*k56*x24 + (-1)*k58*x21)/k115,
  diff(x22, t) = (1*k56*x21 + (-1)*k57*x22)/k115,
  diff(x23, t) = (1*k83*x33 + (-1)*k84*x23)/k115,
  diff(x24, t) = (1*k54*(x20*x45 - k55*x24) + (-1)*k56*x24)/k115,
  diff(x25, t) = ((-1)*k59*x25*x45/(k60 + x25) + 1*k61*x26/(k62 + x26))/k115,
  diff(x26, t) = (1*k59*x25*x45/(k60 + x25) + (-1)*k61*x26/(k62 + x26) + (-1)*k59*x26*x45/(k60 + x26) + 1*k61*x51/(k62 + x51))/k115,
  diff(x27, t) = (1*k18*x3 + (-1)*k16*(x27*x4 - k17*x5))/k115,
  diff(x28, t) = ((-1)*k36*(x6*x28 - k37*x29) + 1*k42*x31 + (-1)*k93*(x28*k12 - k94*x2))/k115,
  diff(x29, t) = (1*k36*(x6*x28 - k37*x29) + (-1)*k38*(x29 - k39*x42))/k115,
  diff(x30, t) = (1*k77*(x35*x41 - k78*x30) + (-1)*k79*x30/(k80 + x30) + 1*k84*x23)/k115,
  diff(x31, t) = (1*(k40*x42 - k41*x6*x31) + (-1)*k42*x31 + (-1)*k63*(x32*x31 - k64*x44) + 1*k66*x15)/k115,
  diff(x32, t) = ((-1)*k63*(x32*x31 - k64*x44) + 1*k70*x37)/k115,
  diff(x33, t) = (1*k81*x40*x20 + (-1)*k83*x33 + 1*k83*x47 + (-1)*k85*x33)/k115,
  diff(x34, t) = ((-1)*k67*(x35*x34 - k68*x36) + 1*k70*x37 + (-1)*k71*x34/(k72 + x34) + (-1)*k73*(x34*x38 - k74*x39) + 2*k76*x43 + (-1)*k91*(x34*k11 - k92*x13))/k115,
  diff(x35, t) = ((-1)*k67*(x35*x34 - k68*x36) + (-1)*k77*(x35*x41 - k78*x30) + 1*k66*x15)/k115,
  diff(x36, t) = (1*k67*(x35*x34 - k68*x36) + (-1)*k69*x36)/k115,
  diff(x37, t) = (1*k69*x36 + (-1)*k70*x37)/k115,
  diff(x38, t) = (1*k71*x34/(k72 + x34) + (-1)*k73*(x34*x38 - k74*x39))/k115,
  diff(x39, t) = (1*k73*(x34*x38 - k74*x39) + (-1)*k75*x39)/k115,
  diff(x40, t) = (1*k79*x30/(k80 + x30) + (-1)*k81*x40*x20 + (-1)*k79*x40/(k80 + x40) + 1*k85*x33)/k115,
  diff(x41, t) = (-1)*k77*(x35*x41 - k78*x30)/k115,
  diff(x42, t) = (1*k38*(x29 - k39*x42) + (-1)*(k40*x42 - k41*x6*x31))/k115,
  diff(x43, t) = (1*k75*x39 + (-1)*k76*x43)/k115,
  diff(x44, t) = (1*k63*(x32*x31 - k64*x44) + (-1)*k65*x44)/k115,
  diff(x45, t) = (1*k52*x19*x52/(k53 + x19) + (-1)*k54*(x20*x45 - k55*x24))/k115,
  diff(x46, t) = (1*k79*x40/(k80 + x40) + (-1)*k81*(x46*x20 - k82*x47))/k115,
  diff(x47, t) = (1*k81*(x46*x20 - k82*x47) + (-1)*k83*x47)/k115,
  diff(x48, t) = (-1)*k87*(x48*x4 - k88*x12)/k115,
  diff(x49, t) = (1*k30*(x11 - k31*x6*x49) + (-1)*k32*(x49 - k33*x54*x10))/k115,
  diff(x50, t) = 1*(k89*x12 - k90*x50)/k115,
  diff(x51, t) = (1*k59*x26*x45/(k60 + x26) + (-1)*k61*x51/(k62 + x51))/k115,
  diff(x52, t) = (1*k47*x16*x53/(k48 + x16) + (-1)*k49*x52*(x46 + k51)/(x52 + k50))/k115,
  diff(x53, t) = (1*k43*x14*x49/(k44 + x14) + (-1)*k45*x53/(k46 + x53))/k115,
  diff(x54, t) = (1*k32*(x49 - k33*x54*x10) + (-1)*k34*x54/(k35 + x54))/k115,
  diff(x55, t) = (-1)*k14*(x1*x55 - k15*x3)/k115 }