{ diff(x1, t) = (1*k3*k7/(k6 + x41) + (-1)*k3*k47*x1*(1 + piecewise(1, and(piecewise(t - 24*ceiling(t/24), xor(t < 0, 24 < 0), t - 24*floor(t/24)) >= 6, piecewise(t - 24*ceiling(t/24), xor(t < 0, 24 < 0), t - 24*floor(t/24)) <= 18), 0)/(1 + piecewise(1, and(piecewise(t - 24*ceiling(t/24), xor(t < 0, 24 < 0), t - 24*floor(t/24)) >= 6, piecewise(t - 24*ceiling(t/24), xor(t < 0, 24 < 0), t - 24*floor(t/24)) <= 18), 0)))/(k46 + x1))/k3,
  diff(x2, t) = (1*k3*k10*x1*(1 + k9*x27)/(k8 + x41) + (-1)*k3*k12*x2/(k11 + x2))/k3,
  diff(x3, t) = (1*k15*x2*k14*(1 + x15/(1 + x15))*(1 + k13*x27) + (-1)*k3*k17*x3/(k16 + x3))/k3,
  diff(x4, t) = 1*1/k18*(x3 - x4)/k4,
  diff(x5, t) = (1*k24*(x11 + k19)/(k20*(1 + (x7/k21)^k23))*(1 + k22*x34) + 1*k25*x33/x11 + (-1)*k4*k26*x5)/k4,
  diff(x6, t) = (1*k4*k27*x5^k28 + (-1)*k4*k29*x6 + (-1)*k30*x6 + 1*k31*x7)/k4,
  diff(x7, t) = (1*k30*x6 + (-1)*k31*x7 + (-1)*k32*x7*(1 + x16))/k5,
  diff(x8, t) = (1*k36*x7^k35/(k33^k35 + x7^k35)*(1 + k34*x24) + (-1)*k4*k37*x8)/k4,
  diff(x9, t) = (1*k4*k38*x8 + (-1)*(k39*x9 - k40*x10) + (-1)*k4*k43*x9)/k4,
  diff(x10, t) = (1*(k39*x9 - k40*x10) + (-1)*k5*k41*x10 + 1*k5*k42*x11 + (-1)*k5*k44*x10)/k5,
  diff(x11, t) = (1*k5*k41*x10 + (-1)*k5*k42*x11 + (-1)*k5*k45*x11 + (-1)*(k67*x11*x16 - k68*x17))/k5,
  diff(x12, t) = (1*k4*k49*(k50 - x12)/(k48 + k50 - x12) + 1*k4*k52*x13/(k51 + x13) + (-1)*k54*x20*x12/(k53 + x12) + (-1)*k4*k56*x12/(k55 + x12))/k4,
  diff(x13, t) = ((-1)*k4*k52*x13/(k51 + x13) + 1*k4*k58*x18*x14/(k57 + x14))/k4,
  diff(x14, t) = (1*k4*k56*x12/(k55 + x12) + (-1)*k4*k58*x18*x14/(k57 + x14))/k4,
  diff(x15, t) = (1*k4*k60*x12/(k59 + x12) + (-1)*k4*k61*x15)/k4,
  diff(x16, t) = (1*k4*k63*x12*(k64 - x16)/(k62 + k64 - x16) + (-1)*k4*k66*x16/(k65 + x16) + (-1)*(k67*x11*x16 - k68*x17))/k4,
  diff(x17, t) = (1*(k67*x11*x16 - k68*x17) + (-1)*k5*k69*x17)/k5,
  diff(x18, t) = (1*k70*x17 + (-1)*k4*k71*x18)/k4,
  diff(x19, t) = 1*k3*1/k72*(piecewise(1, and(piecewise(t - 24*ceiling(t/24), xor(t < 0, 24 < 0), t - 24*floor(t/24)) >= 6, piecewise(t - 24*ceiling(t/24), xor(t < 0, 24 < 0), t - 24*floor(t/24)) <= 18), 0) - x19)/k3,
  diff(x20, t) = 1*k3*1/k73*(x19 - x20)/k3,
  diff(x21, t) = (1*k4*k74*(1 + x25) + (-1)*k4*k75*x21)/k4,
  diff(x22, t) = ((-1)*k4*k76*x22 + (-1)*k77*x22 + 1*k78*x23 + 1*k4*k85*x21)/k4,
  diff(x23, t) = (1*k77*x22 + (-1)*k78*x23 + 1*k5*k79*x24 + (-1)*k80*x24*x16)/k5,
  diff(x24, t) = ((-1)*k5*k79*x24 + (-1)*k5*k81*x24 + 1*k86*x23*(1 + x16))/k5,
  diff(x25, t) = (1*k82*x24 + (-1)*k4*k83*x25)/k4,
  diff(x26, t) = ((-1)*k4*k84*x26 + 1*k134*x25*x33*1/(1 + (x7/k135)^k136))/k4,
  diff(x27, t) = 1*1/k87*(x34 - x27)/k3,
  diff(x28, t) = (1*k4*k92*(1 + k91*x4/(k89 + x4))*(k90 - x28)/(k88 + k90 - x28) + (-1)*k4*k99*x28/(k98 + x28) + (-1)*k4*k102*x28*x4 + 1*k106*x32)/k4,
  diff(x29, t) = (1*k4*k97*(1 + k96*x4/(k95 + x4))*(k93 - x29)/(k94 + k93 - x29) + (-1)*k4*k101*x29/(k100 + x29) + (-1)*k4*k103*x29*x4 + 1*k107*x33)/k4,
  diff(x30, t) = (1*k4*k103*x29*x4 + (-1)*k105*x30)/k4,
  diff(x31, t) = (1*k4*k102*x28*x4 + (-1)*k104*x31)/k4,
  diff(x32, t) = (1*k104*x31 + (-1)*k106*x32)/k5,
  diff(x33, t) = (1*k105*x30 + (-1)*k107*x33)/k5,
  diff(x34, t) = (1*k4*k114*x35 + (-1)*k4*k115*x34)/k4,
  diff(x35, t) = (1*k111*(1 - k110*x33/(k108 + x33))*(1 - k112*x10/(k109 + x10))*(1 + x38) + (-1)*k4*k113*x35)/k4,
  diff(x36, t) = (1*k118*(1 - k117*x32/(k116 + x32)) + (-1)*k4*k119*x36)/k4,
  diff(x37, t) = (1*k4*k120*x36 + (-1)*k4*k121*x37 + (-1)*k4*k122*x37*x34)/k4,
  diff(x38, t) = (1*k4*k122*x37*x34 + (-1)*k4*k133*x38)/k4,
  diff(x39, t) = ((-1)*k3*k124*x39/(k123 + x39) + 1*k3*k129*(1 + k128*x3/(k127 + x3))*(k125 - x39)/(k126 + k125 - x39) + (-1)*k3*k130*x39*x3 + 1*k3*k132*x41)/k3,
  diff(x40, t) = (1*k3*k130*x39*x3 + (-1)*k3*k131*x40)/k3,
  diff(x41, t) = (1*k3*k131*x40 + (-1)*k3*k132*x41)/k3 }