[ Eq(Derivative(x1, t), ((-1)*k86*(k1*x2*x1 - k2*x3) + (-1)*k86*(k32*x20*x1 - k33*x25))/k86),
  Eq(Derivative(x2, t), (-1)*k86*(k1*x2*x1 - k2*x3)/k86),
  Eq(Derivative(x3, t), (1*k86*(k1*x2*x1 - k2*x3) + (-1)*k86*(k3*x3*x8 - k4*x9) + (-1)*k86*(k28*x3*x21 - k29*x23) + (-1)*k86*(k40*x27*x3 - k41*x28))/k86),
  Eq(Derivative(x4, t), ((-1)*k86*k5*x9*x4/(k13 + x4) + (-1)*k86*k17*x24*x4/(k18 + x4) + (-1)*k86*k23*x22*x4/(k24 + x4) + (-1)*k86*k36*x26*x4/(k37 + x4) + (-1)*k86*k55*x34*x4/(k56 + x4) + (-1)*k86*k61*x36*x4/(k62 + x4) + (-1)*k86*k69*x37*x4/(k70 + x4) + (-1)*k86*k77*x41*x4/(k78 + x4))/k86),
  Eq(Derivative(x5, t), (1*k86*k5*x9*x4/(k13 + x4) + 1*k86*k17*x24*x4/(k18 + x4) + 1*k86*k23*x22*x4/(k24 + x4) + 1*k86*k36*x26*x4/(k37 + x4) + 1*k86*k55*x34*x4/(k56 + x4) + 1*k86*k61*x36*x4/(k62 + x4) + 1*k86*k69*x37*x4/(k70 + x4) + 1*k86*k77*x41*x4/(k78 + x4))/k86),
  Eq(Derivative(x6, t), (1*k86*k5*x9*x4/(k13 + x4) + (-1)*k86*(k7*x6*x10 - k8*x12) + 1*k86*k17*x24*x4/(k18 + x4) + 1*k86*k23*x22*x4/(k24 + x4) + 1*k86*k36*x26*x4/(k37 + x4) + (-1)*k86*(k46*x27*x6 - k45*x31) + 1*k86*k47*x12*x27 + 1*k86*k55*x34*x4/(k56 + x4) + 1*k86*k61*x36*x4/(k62 + x4) + 1*k86*k69*x37*x4/(k70 + x4) + 1*k86*k77*x41*x4/(k78 + x4))/k86),
  Eq(Derivative(x7, t), ((-1)*k86*k6*x9*x7/(k14 + x7) + (-1)*k86*k30*x24*x7/(k31 + x7) + (-1)*k86*k25*x22*x7/(k26 + x7) + (-1)*k86*k38*x26*x7/(k39 + x7) + (-1)*k86*k57*x34*x7/(k58 + x7) + (-1)*k86*k63*x36*x7/(k64 + x7) + (-1)*k86*k71*x37*x7/(k72 + x7) + (-1)*k86*k79*x41*x7/(k80 + x7))/k86),
  Eq(Derivative(x8, t), ((-1)*k86*(k3*x3*x8 - k4*x9) + (-1)*k86*(k34*x25*x8 - k35*x26) + (-1)*k86*(k53*x23*x8 - k54*x34) + (-1)*k86*(k67*x24*x8 - k68*x37))/k86),
  Eq(Derivative(x9, t), 1*k86*(k3*x3*x8 - k4*x9)/k86),
  Eq(Derivative(x10, t), (1*k86*k6*x9*x7/(k14 + x7) + (-1)*k86*(k7*x6*x10 - k8*x12) + 1*k86*k30*x24*x7/(k31 + x7) + 1*k86*k25*x22*x7/(k26 + x7) + 1*k86*k38*x26*x7/(k39 + x7) + 1*k86*k47*x12*x27 + 1*k86*k57*x34*x7/(k58 + x7) + 1*k86*k63*x36*x7/(k64 + x7) + 1*k86*k71*x37*x7/(k72 + x7) + 1*k86*k79*x41*x7/(k80 + x7))/k86),
  Eq(Derivative(x11, t), (1*k86*k6*x9*x7/(k14 + x7) + 1*k86*k30*x24*x7/(k31 + x7) + 1*k86*k25*x22*x7/(k26 + x7) + 1*k86*k38*x26*x7/(k39 + x7) + 1*k86*k57*x34*x7/(k58 + x7) + 1*k86*k63*x36*x7/(k64 + x7) + 1*k86*k71*x37*x7/(k72 + x7) + 1*k86*k79*x41*x7/(k80 + x7))/k86),
  Eq(Derivative(x12, t), (1*k86*(k7*x6*x10 - k8*x12) + (-1)*k86*(k12*x12 - k27*x18) + (-1)*k86*k44*x12*x27 + (-1)*k86*k47*x12*x27 + (-1)*k86*(k49*x12*x27 - k50*x32) + (-1)*k86*k89*x12)/k86),
  Eq(Derivative(x13, t), (-1)*k86*k9*x12*x13/k86),
  Eq(Derivative(x14, t), 1*k86*k9*x12*x13/k86),
  Eq(Derivative(x15, t), (1*k86*k9*x12*x13 + (-1)*k86*(k10*x15 - k11*x16) + (-1)*k86*k88*x15)/k86),
  Eq(Derivative(x16, t), 1*k86*(k10*x15 - k11*x16)/k86),
  Eq(Derivative(x17, t), ((-1)*k86*(k16*x23*x17 - k15*x24) + (-1)*k86*(k21*x20*x17 - k22*x22) + (-1)*k86*(k59*x25*x17 - k60*x36) + (-1)*k86*(k65*x34*x17 - k66*x37) + (-1)*k86*(k75*x39*x17 - k76*x41))/k86),
  Eq(Derivative(x18, t), (1*k86*(k12*x12 - k27*x18) + (-1)*k86*(k51*x18*x27 - k52*x33))/k86),
  Eq(Derivative(x19, t), (-1)*k86*(k19*x19*x21 - k20*x20)/k86),
  Eq(Derivative(x20, t), (1*k86*(k19*x19*x21 - k20*x20) + (-1)*k86*(k21*x20*x17 - k22*x22) + (-1)*k86*(k32*x20*x1 - k33*x25))/k86),
  Eq(Derivative(x21, t), ((-1)*k86*(k28*x3*x21 - k29*x23) + (-1)*k86*(k19*x19*x21 - k20*x20))/k86),
  Eq(Derivative(x22, t), 1*k86*(k21*x20*x17 - k22*x22)/k86),
  Eq(Derivative(x23, t), (1*k86*(k28*x3*x21 - k29*x23) + (-1)*k86*(k16*x23*x17 - k15*x24) + (-1)*k86*(k53*x23*x8 - k54*x34) + (-1)*k86*(k90*x27*x23 - k91*x35))/k86),
  Eq(Derivative(x24, t), (1*k86*(k16*x23*x17 - k15*x24) + (-1)*k86*(k67*x24*x8 - k68*x37))/k86),
  Eq(Derivative(x25, t), (1*k86*(k32*x20*x1 - k33*x25) + (-1)*k86*(k34*x25*x8 - k35*x26) + (-1)*k86*(k42*x27*x25 - k43*x29) + (-1)*k86*(k59*x25*x17 - k60*x36))/k86),
  Eq(Derivative(x26, t), 1*k86*(k34*x25*x8 - k35*x26)/k86),
  Eq(Derivative(x27, t), ((-1)*k86*(k40*x27*x3 - k41*x28) + (-1)*k86*(k42*x27*x25 - k43*x29) + (-1)*k86*(k46*x27*x6 - k45*x31) + (-1)*k86*(k49*x12*x27 - k50*x32) + (-1)*k86*(k51*x18*x27 - k52*x33) + (-1)*k86*k87*x27 + (-1)*k86*(k90*x27*x23 - k91*x35))/k86),
  Eq(Derivative(x28, t), 1*k86*(k40*x27*x3 - k41*x28)/k86),
  Eq(Derivative(x29, t), 1*k86*(k42*x27*x25 - k43*x29)/k86),
  Eq(Derivative(x30, t), 1*k86*k44*x12*x27/k86),
  Eq(Derivative(x31, t), 1*k86*(k46*x27*x6 - k45*x31)/k86),
  Eq(Derivative(x32, t), 1*k86*(k49*x12*x27 - k50*x32)/k86),
  Eq(Derivative(x33, t), 1*k86*(k51*x18*x27 - k52*x33)/k86),
  Eq(Derivative(x34, t), (1*k86*(k53*x23*x8 - k54*x34) + (-1)*k86*(k65*x34*x17 - k66*x37))/k86),
  Eq(Derivative(x35, t), 1*k86*(k90*x27*x23 - k91*x35)/k86),
  Eq(Derivative(x36, t), 1*k86*(k59*x25*x17 - k60*x36)/k86),
  Eq(Derivative(x37, t), (1*k86*(k65*x34*x17 - k66*x37) + 1*k86*(k67*x24*x8 - k68*x37))/k86),
  Eq(Derivative(x38, t), 0/k86),
  Eq(Derivative(x39, t), (1*k86*(k73*x42*x40 - k74*x39) + (-1)*k86*(k75*x39*x17 - k76*x41))/k86),
  Eq(Derivative(x40, t), (-1)*k86*(k73*x42*x40 - k74*x39)/k86),
  Eq(Derivative(x41, t), 1*k86*(k75*x39*x17 - k76*x41)/k86),
  Eq(Derivative(x42, t), (-1)*k86*(k73*x42*x40 - k74*x39)/k86) ]