{ diff(x1, t) = (1*k33*k1*x1*(1 - x1/k2) + 1*k33*k10*x1/k2*x2/k9*(k2 - x1) + (-1)*k33*k3*x6/(x6 + k6)*(k5 + k6*(1 - k5)/(x4 + k6))*x3*x1/(k7 + x1))/k33,
  diff(x2, t) = ((-1)*k33*k10*x1/k2*x2/k9*(k2 - x1) + 1*k33*k8*x2*(1 - x2/k9) + (-1)*k33*k3*1/10*x6/(x6 + k6)*(k5 + k6*(1 - k5)/(x4 + k6))*x3*x2/(k7 + x2))/k33,
  diff(x3, t) = ((-1)*k33*k16*x3 + 1*k33*k32)/k33,
  diff(x4, t) = (1*k33*k17 + 1*k33*k18*x1 + 1*k33*k18*x2 + (-1)*k33*k20*x4)/k33,
  diff(x5, t) = (1*k33*k30*x3 + (-1)*k33*k31*x5)/k33,
  diff(x6, t) = (1*k33*k21 + (-1)*k33*k24*x6 + 1*k33*k22*x5/(x5 + k23))/k33,
  diff(x7, t) = (1*k33*k25*x5/(x5 + k26)*(k28*(1 - k27)/(x4 + k28) + k27) + (-1)*k33*k29*x7)/k33 }