{ diff(x1, t) = (-1)*k78*(x1*x2*k1 - x22*k2)/k78,
  diff(x2, t) = (-1)*k78*(x1*x2*k1 - x22*k2)/k78,
  diff(x3, t) = (-1)*k78*(x23*x3*k4 - x24*k5)/k78,
  diff(x4, t) = ((-1)*k78*(x23*x4*k6 - x25*k7) + (-1)*k78*(x6*x4*k17 - x26*k18))/k78,
  diff(x5, t) = (1*k78*x25*k8 + 1*k78*x26*k70 + (-1)*k78*(x5*x7*k9 - x27*k10) + (-1)*k78*(x8*x5*k11 - x28*k12) + 1*k78*x28*k13 + (-1)*k78*(x5*x10*k25 - x31*k26) + 1*k78*x31*k27)/k78,
  diff(x6, t) = ((-1)*k78*(x6*x4*k17 - x26*k18) + 1*k78*x26*k70 + 1*k78*x32*k16)/k78,
  diff(x7, t) = (-1)*k78*(x5*x7*k9 - x27*k10)/k78,
  diff(x8, t) = ((-1)*k78*(x8*x5*k11 - x28*k12) + (-1)*k78*(x8*x30*k61 - x29*k62))/k78,
  diff(x9, t) = (1*k78*x28*k13 + 1*k78*x29*k63 + (-1)*k78*(x9*x12*k14 - x32*k15) + 1*k78*x32*k16 + (-1)*k78*(x9*x13*k19 - x33*k20) + (-1)*k78*(x15*x9*k22 - x34*k23) + 1*k78*x34*k24)/k78,
  diff(x10, t) = (-1)*k78*(x5*x10*k25 - x31*k26)/k78,
  diff(x11, t) = (1*k78*x31*k27 + (-1)*k78*(x11*x18*k28 - x37*k29) + (-1)*k78*(x11*x19*k30 - x38*k31) + 1*k78*x38*k32)/k78,
  diff(x12, t) = (-1)*k78*(x9*x12*k14 - x32*k15)/k78,
  diff(x13, t) = ((-1)*k78*(x9*x13*k19 - x33*k20) + 1*k78*x33*k21 + (-1)*k78*(x13*x30*k66 - x35*k67) + (-1)*k78*(x13*x17*k68 - x36*k69))/k78,
  diff(x14, t) = 1*k78*x33*k21/k78,
  diff(x15, t) = (-1)*k78*(x15*x9*k22 - x34*k23)/k78,
  diff(x16, t) = 1*k78*x34*k24/k78,
  diff(x17, t) = ((-1)*k78*(x13*x17*k68 - x36*k69) + 1*k78*(x48*k64 - x17*k65))/k78,
  diff(x18, t) = (-1)*k78*(x11*x18*k28 - x37*k29)/k78,
  diff(x19, t) = (-1)*k78*(x11*x19*k30 - x38*k31)/k78,
  diff(x20, t) = (1*k78*x38*k32 + (-1)*k78*(x20*k33 - x57*k34))/k78,
  diff(x21, t) = ((-1)*k79*(x57*x21*k35/k71 - x39*k36) + (-1)*k79*(x21*x58*k39/k71 - x41*k40) + (-1)*k79*(x21*x40*k43/k71 - x42*k44))/k79,
  diff(x22, t) = (1*k78*(x1*x2*k1 - x22*k2) + (-1)*k78*x22*k3)/k78,
  diff(x23, t) = (1*k78*x22*k3 + (-1)*k78*(x23*x3*k4 - x24*k5) + (-1)*k78*(x23*x4*k6 - x25*k7) + 1*k78*x25*k8)/k78,
  diff(x24, t) = 1*k78*(x23*x3*k4 - x24*k5)/k78,
  diff(x25, t) = (1*k78*(x23*x4*k6 - x25*k7) + (-1)*k78*x25*k8)/k78,
  diff(x26, t) = (1*k78*(x6*x4*k17 - x26*k18) + (-1)*k78*x26*k70)/k78,
  diff(x27, t) = 1*k78*(x5*x7*k9 - x27*k10)/k78,
  diff(x28, t) = (1*k78*(x8*x5*k11 - x28*k12) + (-1)*k78*x28*k13)/k78,
  diff(x29, t) = (1*k78*(x8*x30*k61 - x29*k62) + (-1)*k78*x29*k63)/k78,
  diff(x30, t) = ((-1)*k78*(x8*x30*k61 - x29*k62) + 1*k78*x29*k63 + (-1)*k78*(x13*x30*k66 - x35*k67) + 1*k78*(x55*x56*k59 - x30*k60))/k78,
  diff(x31, t) = (1*k78*(x5*x10*k25 - x31*k26) + (-1)*k78*x31*k27)/k78,
  diff(x32, t) = (1*k78*(x9*x12*k14 - x32*k15) + (-1)*k78*x32*k16)/k78,
  diff(x33, t) = (1*k78*(x9*x13*k19 - x33*k20) + (-1)*k78*x33*k21)/k78,
  diff(x34, t) = (1*k78*(x15*x9*k22 - x34*k23) + (-1)*k78*x34*k24)/k78,
  diff(x35, t) = 1*k78*(x13*x30*k66 - x35*k67)/k78,
  diff(x36, t) = 1*k78*(x13*x17*k68 - x36*k69)/k78,
  diff(x37, t) = 1*k78*(x11*x18*k28 - x37*k29)/k78,
  diff(x38, t) = (1*k78*(x11*x19*k30 - x38*k31) + (-1)*k78*x38*k32)/k78,
  diff(x39, t) = 1*k79*(x57*x21*k35/k71 - x39*k36)/k79,
  diff(x40, t) = (1*k79*(x58*x58*k41/k71 - x40*k42) + (-1)*k79*(x21*x40*k43/k71 - x42*k44) + (-1)*k79*(x40*x43*k45/k71 - x44*k46))/k79,
  diff(x41, t) = 1*k79*(x21*x58*k39/k71 - x41*k40)/k79,
  diff(x42, t) = 1*k79*(x21*x40*k43/k71 - x42*k44)/k79,
  diff(x43, t) = (-1)*k79*(x40*x43*k45/k71 - x44*k46)/k79,
  diff(x44, t) = (1*k79*(x40*x43*k45/k71 - x44*k46) + (-1)*k79*x44*k47)/k79,
  diff(x45, t) = (1*k79*x44*k47 + (-1)*k79*(x45*x46*k51/k71 - x47*k52) + 1*k79*x47*k53 + (-1)*k79*(x45*x49*k48/k71 - x50*k49) + 1*k79*x50*k50)/k79,
  diff(x46, t) = (-1)*k79*(x45*x46*k51/k71 - x47*k52)/k79,
  diff(x47, t) = (1*k79*(x45*x46*k51/k71 - x47*k52) + (-1)*k79*x47*k53)/k79,
  diff(x48, t) = (1*k79*x47*k53 + (-1)*k78*(x48*k64 - x17*k65))/k79,
  diff(x49, t) = (-1)*k79*(x45*x49*k48/k71 - x50*k49)/k79,
  diff(x50, t) = (1*k79*(x45*x49*k48/k71 - x50*k49) + (-1)*k79*x50*k50)/k79,
  diff(x51, t) = (1*k79*x50*k50 + (-1)*k78*(x51*k54 - x52*k55))/k79,
  diff(x52, t) = (1*k78*(x51*k54 - x52*k55) + (-1)*k78*(x52*x53*k56 - x54*k57) + 1*k78*x54*k58)/k78,
  diff(x53, t) = (-1)*k78*(x52*x53*k56 - x54*k57)/k78,
  diff(x54, t) = (1*k78*(x52*x53*k56 - x54*k57) + (-1)*k78*x54*k58)/k78,
  diff(x55, t) = (1*k78*x54*k58 + (-1)*k78*(x55*x56*k59 - x30*k60))/k78,
  diff(x56, t) = (-1)*k78*(x55*x56*k59 - x30*k60)/k78,
  diff(x57, t) = (1*k78*(x20*k33 - x57*k34) + (-1)*k79*(x57*x21*k35/k71 - x39*k36) + (-1)*k79*(x57*x57*k37/k71 - x58*k38))/k79,
  diff(x58, t) = (1*k79*(x57*x57*k37/k71 - x58*k38) + (-1)*k79*(x58*x58*k41/k71 - x40*k42) + (-1)*k79*(x21*x58*k39/k71 - x41*k40))/k79 }