{ diff(x1, t) = ((-1)*k86*(k1*x2*x1 - k2*x3) + (-1)*k86*(k32*x20*x1 - k33*x25))/k86,
  diff(x2, t) = (-1)*k86*(k1*x2*x1 - k2*x3)/k86,
  diff(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,
  diff(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,
  diff(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,
  diff(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,
  diff(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,
  diff(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,
  diff(x9, t) = 1*k86*(k3*x3*x8 - k4*x9)/k86,
  diff(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,
  diff(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,
  diff(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,
  diff(x13, t) = (-1)*k86*k9*x12*x13/k86,
  diff(x14, t) = 1*k86*k9*x12*x13/k86,
  diff(x15, t) = (1*k86*k9*x12*x13 + (-1)*k86*(k10*x15 - k11*x16) + (-1)*k86*k88*x15)/k86,
  diff(x16, t) = 1*k86*(k10*x15 - k11*x16)/k86,
  diff(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,
  diff(x18, t) = (1*k86*(k12*x12 - k27*x18) + (-1)*k86*(k51*x18*x27 - k52*x33))/k86,
  diff(x19, t) = (-1)*k86*(k19*x19*x21 - k20*x20)/k86,
  diff(x20, t) = (1*k86*(k19*x19*x21 - k20*x20) + (-1)*k86*(k21*x20*x17 - k22*x22) + (-1)*k86*(k32*x20*x1 - k33*x25))/k86,
  diff(x21, t) = ((-1)*k86*(k28*x3*x21 - k29*x23) + (-1)*k86*(k19*x19*x21 - k20*x20))/k86,
  diff(x22, t) = 1*k86*(k21*x20*x17 - k22*x22)/k86,
  diff(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,
  diff(x24, t) = (1*k86*(k16*x23*x17 - k15*x24) + (-1)*k86*(k67*x24*x8 - k68*x37))/k86,
  diff(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,
  diff(x26, t) = 1*k86*(k34*x25*x8 - k35*x26)/k86,
  diff(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,
  diff(x28, t) = 1*k86*(k40*x27*x3 - k41*x28)/k86,
  diff(x29, t) = 1*k86*(k42*x27*x25 - k43*x29)/k86,
  diff(x30, t) = 1*k86*k44*x12*x27/k86,
  diff(x31, t) = 1*k86*(k46*x27*x6 - k45*x31)/k86,
  diff(x32, t) = 1*k86*(k49*x12*x27 - k50*x32)/k86,
  diff(x33, t) = 1*k86*(k51*x18*x27 - k52*x33)/k86,
  diff(x34, t) = (1*k86*(k53*x23*x8 - k54*x34) + (-1)*k86*(k65*x34*x17 - k66*x37))/k86,
  diff(x35, t) = 1*k86*(k90*x27*x23 - k91*x35)/k86,
  diff(x36, t) = 1*k86*(k59*x25*x17 - k60*x36)/k86,
  diff(x37, t) = (1*k86*(k65*x34*x17 - k66*x37) + 1*k86*(k67*x24*x8 - k68*x37))/k86,
  diff(x38, t) = 0/k86,
  diff(x39, t) = (1*k86*(k73*x42*x40 - k74*x39) + (-1)*k86*(k75*x39*x17 - k76*x41))/k86,
  diff(x40, t) = (-1)*k86*(k73*x42*x40 - k74*x39)/k86,
  diff(x41, t) = 1*k86*(k75*x39*x17 - k76*x41)/k86,
  diff(x42, t) = (-1)*k86*(k73*x42*x40 - k74*x39)/k86 }