{ diff(x1, t) = ((-1)*k2*x1^2/(x1^2 + k6^2) + 1*(1 - k17)*(k7 + (1 - x3)*k1*x2^2/(x2^2 + k9^2)*x1^2/(x1^2 + k13^2))*x5 + 1*k12*(k10 + k11^8/(x5^8 + k11^8)) + (-1)*k15*x1^2/(x1^2 + k14^2) + (-1)*k19*x1^4/(x1^4 + k20^4) + 1*k18*x6/(x6 + 1/100))/k21,
  diff(x2, t) = (1*k8*x1 + (-1)*k3*x2)/k21,
  diff(x3, t) = (1*k4*x1^4/(x1^4 + k16^4)*(1 - x3) + (-1)*k5*x3)/k21,
  diff(x4, t) = 0/k22,
  diff(x5, t) = (1*k2*x1^2/(x1^2 + k6^2) + (-1)*(1 - k17)*(k7 + (1 - x3)*k1*x2^2/(x2^2 + k9^2)*x1^2/(x1^2 + k13^2))*x5 + (-1)*k17*(k7 + (1 - x3)*k1*x2^2/(x2^2 + k9^2)*x1^2/(x1^2 + k13^2))*x5)/k23,
  diff(x6, t) = (1*k19*x1^4/(x1^4 + k20^4) + (-1)*k18*x6/(x6 + 1/100) + 1*k17*(k7 + (1 - x3)*k1*x2^2/(x2^2 + k9^2)*x1^2/(x1^2 + k13^2))*x5)/k24 }