{ diff(x1, t) = (k1*(1 - x1/k5) - k2*x2)*x1,
  diff(x2, t) = ((-k3) + k4*x1)*x2 }