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