{ diff(x1, t) = (1*k10*k1*x1*(1 - x1/k2) + (-1)*k10*k3*x1*x2 + (-1)*k10*k8*x1)/k10,
  diff(x2, t) = (1*k10*k9*x2*x3 + (-1)*k10*k4*x2 + (-1)*k10*k5*x1*x2)/k10,
  diff(x3, t) = (1*k10*k6*x3*(1 - x3/k7) + (-1)*k10*k9*x2*x3)/k10 }