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