{ diff(x1, t) = (1*k79*k3^k4/(k3^k4 + x16^k4)*(k78*(k1*x13 + k2) + k5*x9^k6/(k7^k6 + x9^k6)) + (-1)*k79*k8*x1/(k9 + x1))/k79,
  diff(x2, t) = (1*k79*k10*x1 + (-1)*k79*(k11*x2 - k12*x3) + (-1)*k79*k13*x2/(k14 + x2))/k79,
  diff(x3, t) = (1*k79*(k11*x2 - k12*x3) + (-1)*k79*k15*x3/(k16 + x3))/k79,
  diff(x4, t) = (1*k79*k17*x12^k18/(k19^k18 + x12^k18)*k20^k21/(k20^k21 + x3^k21) + (-1)*k79*k22*x4/(k23 + x4))/k79,
  diff(x5, t) = (1*k24*k79*x4 + (-1)*k79*((-k26)*x6 + k25*x5) + (-1)*k79*(1 - k78)*k27*x5/(k29 + x5) + (-1)*k28*k79*x5/(k29 + x5))/k79,
  diff(x6, t) = (1*k79*((-k26)*x6 + k25*x5) + (-1)*k79*(1 - k78)*k30*x6/(k32 + x6) + (-1)*k31*k79*x6/(k32 + x6))/k79,
  diff(x7, t) = (1*k79*k33*x6^k34/(k35^k34 + x6^k34) + (-1)*k79*k36*x7/(k37 + x7))/k79,
  diff(x8, t) = ((-1)*k79*(k39*x8 - k40*x9) + (-1)*k79*k41*x8/(k42 + x8) + 1*k79*k38*x7)/k79,
  diff(x9, t) = (1*k79*(k39*x8 - k40*x9) + (-1)*k79*k43*x9/(k44 + x9))/k79,
  diff(x10, t) = (1*k79*(k78*k45*x13 + (k78*k46 + k47)*k48^k50/(k48^k50 + x6^k50))*k49^k51/(k49^k51 + x3^k51) + (-1)*k79*k52*x10/(k53 + x10))/k79,
  diff(x11, t) = ((-1)*k79*(k55*x11 - k56*x12) + 1*k79*k54*x10 + (-1)*k79*k57*x11/(k58 + x11))/k79,
  diff(x12, t) = (1*k79*(k55*x11 - k56*x12) + (-1)*k79*k59*x12/(k60 + x12))/k79,
  diff(x13, t) = (1*(1 - k78)*k61*k79 + (-1)*k64*k78*x13*k79 + (-1)*k79*k63*x13/(k62 + x13))/k79,
  diff(x14, t) = (1*k78*k65*x13*k79 + 1*k79*k67*x3^k66/(k68^k66 + x3^k66) + (-1)*k79*k69*x14/(k70 + x14))/k79,
  diff(x15, t) = (1*k79*k71*x14 + (-1)*k79*(k72*x15 - k73*x16) + (-1)*k79*k74*x15/(k75 + x15))/k79,
  diff(x16, t) = (1*k79*(k72*x15 - k73*x16) + (-1)*k79*k76*x16/(k77 + x16))/k79 }