{ diff(x1, t) = (1*k11*k12^k1/(k12^k1 + x10^k1) + (-1)*(k13/(k14 + x1) + k2)*x1*k8)/k8,
  diff(x2, t) = (1*k9*k10^k1/(k10^k1 + x10^k1) + (-1)*(k15/(k16 + x2) + k2)*x2*k8)/k8,
  diff(x3, t) = (1*k18*x1*k8 + (-1)*k2*x3*k8 + (-1)*k33*x3/(k34 + x3)*k8 + 1*k38*x4/(k37 + x4)*k8)/k8,
  diff(x4, t) = ((-1)*k2*x4*k8 + 1*k33*x3/(k34 + x3)*k8 + (-1)*k35*x4/(k36 + x4)*k8 + (-1)*k38*x4/(k37 + x4)*k8 + 1*k39*x5/(k40 + x5)*k8)/k8,
  diff(x5, t) = ((-1)*k2*x5*k8 + (-1)*(k19*x5*x8 - k20*x9)*k8 + 1*k35*x4/(k36 + x4)*k8 + (-1)*k39*x5/(k40 + x5)*k8 + (-1)*(1 + (k5 - 1)*ceiling(sin(t/k4*pi)*9/10))*k6*x5/(k41 + x5)*k8)/k8,
  diff(x6, t) = (1*k17*x2*k8 + (-1)*k2*x6*k8 + (-1)*k25*x6/(k26 + x6)*k8 + 1*k30*x7/(k29 + x7)*k8)/k8,
  diff(x7, t) = ((-1)*k2*x7*k8 + 1*k25*x6/(k26 + x6)*k8 + (-1)*k27*x7/(k28 + x7)*k8 + (-1)*k30*x7/(k29 + x7)*k8 + 1*k31*x8/(k32 + x8)*k8)/k8,
  diff(x8, t) = ((-1)*k2*x8*k8 + (-1)*(k19*x5*x8 - k20*x9)*k8 + 1*k27*x7/(k28 + x7)*k8 + (-1)*k31*x8/(k32 + x8)*k8 + (-1)*k42*x8/(k43 + x8)*k8)/k8,
  diff(x9, t) = (1*(k19*x5*x8 - k20*x9)*k8 + (-1)*k21*x9*k8 + (-1)*(k23*x9*k8 - k24*x10*k7))/k8,
  diff(x10, t) = ((-1)*k22*x10*k7 + 1*(k23*x9*k8 - k24*x10*k7))/k7 }