{ diff(x1, t) = (1*k39*k1*x1*(1 - x1/k2) + (-1)*k39*k3*x5/(x5 + k4)*x2*x1/(k7 + x1)*(k5 + k6*(1 - k5)/(x3 + k6)))/k39,
  diff(x2, t) = (1*k39*k8*x6*x1/(x6*x1 + k9)*(k10 + k11*(1 - k10)/(x3 + k11)) + (-1)*k39*k12*x2 + 1*k39*(piecewise(k29/k38, and(and(piecewise(t - 120*ceiling(t/120), xor(t < 0, 120 < 0), t - 120*floor(t/120)) >= 0, piecewise(t - 120*ceiling(t/120), xor(t < 0, 120 < 0), t - 120*floor(t/120)) <= k38), t < 360), 0) + piecewise(k37/k38, and(and(and(piecewise(t - 120*ceiling(t/120), xor(t < 0, 120 < 0), t - 120*floor(t/120)) >= 0, piecewise(t - 120*ceiling(t/120), xor(t < 0, 120 < 0), t - 120*floor(t/120)) <= k38), t > 1440), t < 1800), 0) + piecewise(k37/k38, and(and(and(piecewise(t - 120*ceiling(t/120), xor(t < 0, 120 < 0), t - 120*floor(t/120)) >= 0, piecewise(t - 120*ceiling(t/120), xor(t < 0, 120 < 0), t - 120*floor(t/120)) <= k38), t > 2880), t < 3240), 0) + piecewise(k37/k38, and(and(and(piecewise(t - 120*ceiling(t/120), xor(t < 0, 120 < 0), t - 120*floor(t/120)) >= 0, piecewise(t - 120*ceiling(t/120), xor(t < 0, 120 < 0), t - 120*floor(t/120)) <= k38), t > 4320), t < 4680), 0) + piecewise(k37/k38, and(and(and(piecewise(t - 120*ceiling(t/120), xor(t < 0, 120 < 0), t - 120*floor(t/120)) >= 0, piecewise(t - 120*ceiling(t/120), xor(t < 0, 120 < 0), t - 120*floor(t/120)) <= k38), t > 5760), t < 6120), 0)))/k39,
  diff(x3, t) = (1*k39*k13 + 1*k39*k5*x1 + (-1)*k39*k15*x3)/k39,
  diff(x4, t) = (1*k39*k16*x2 + (-1)*k39*k17*x4)/k39,
  diff(x5, t) = (1*k39*k18 + 1*k39*k40 + (-1)*k39*k21*x5)/k39,
  diff(x6, t) = (1*k39*k41 + (-1)*k39*k26*x6)/k39 }