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