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