I need to find the final state (i.e. the state at $t\to+\infty$) of the following ODE system:
$$\begin{eqnarray\*}\frac{dA}{dt}&=&-aAX\\
\frac{dB}{dt}&=&-bBX\\
\frac{dX}{dt}&=&-(aA+bB)X\end{eqnarray\*}$$
where $A$, $B$, and $X$ are non-negative variables (with known initial values at $t=0$: $A(0)\ge 0$, $B(0)\ge 0$, $X(0)>0$), $a$ and $b$ are positive constants.


The ODE system appears to have no analytical solution (at least, both me and the automatic analytical ODE solvers have failed to find the solution).
But is it possible to find a formula for $A(+\infty)$ and $B(+\infty)$ without having a full formula for $A(t)$ and $B(t)$?


P.S. A "physical" interpretation: the ODE describes the process of irreversible competitive binding. $A$ and $B$ are the amounts of the competing free substances, $X$ is the amount of free "receptors", $a$ and $b$ are the reaction constants ($\approx$reaction speeds). So the game is about what runs out first - the substances or the free receptors. The required $A(+\infty)$ and $B(+\infty)$ give the final ratio of bound substances.


The main pitfall here is that both situations are possible: there can be more substances than receptors ($A(0)+B(0)>X(0)$) or vice versa ($A(0)+B(0)<X(0)$).



Is it possible to find the final state of the bilinear ODE system?
EDIT: It's obvious there can be no equilibria with all $A,B,X > 0$, and $A,B,X$ are decreasing. The only possibilities are a limit where $X=0$ or a limit where $A=B=0$. As noted, $A+B-X$ is invariant, so the first occurs, with $A(\infty) + B(\infty) = A(0) + B(0) - X(0)$, if $A(0) + B(0) - X(0) \ge 0$ while the latter occurs, with $X(\infty) = -A(0) - B(0) + X(0)$, if $A(0) + B(0) - X(0) \le 0$.


For further analysis of the first case, consider $A$ and $B$ as functions of $X$. The system becomes 


$$ \eqalign{\dfrac{dA}{dX} &= \dfrac{aA}{aA + bB}\cr
 \dfrac{dB}{dX} &= \dfrac{bB}{aA + bB}\cr}$$


This has implicit solutions
$$ \eqalign{c\_1 A^{b/a} &+ A - X = c\_2 \cr
 B &= c\_1 A^{b/a}\cr} $$
We can determine $c\_1$ and $c\_2$ from the initial conditions, and then take
$X=0$ to get the limiting values.


