There several different ways to handle
[itex]x_1'= a- ax_1- ax_2[/itex]
[itex]x_2'= b- bx_1- bx_2[/itex]
The more "elementary" method would be to differentiate the first equation again:
[itex]x_1''= -ax_1'- ax_2'= ax_1'- a(b- bx_1- bx_2)= (a+ ab)x_1+ abx_2- ab[itex]
but, from the first equation, multiplied by b, [itex]abx_2= ab- bx_1'- abx_1[itex]so <br />
[itex]x_1''= (a+ ab)x_1+ (ab- bx_1'- abx_1)- ab[/itex]<br />
[itex]x_1''+ bx_1'- ax_1= 0[/itex].<br />
That has characteristic equation [itex]r^2+ br- a= 0[/itex] with roots [itex]r= \frac{-b\pm\sqrt{b^2+ 4a}}{2}[/itex] so that [itex]x_1(t)= e^{-bt}(C_1e^{\sqrt{b^2+ 4a}t}+ C_2e^{-\sqrt{b^2+ 4a}t}[/itex] if [itex]b^2+ 4a>0[/itex], [itex]x_1(t)= e^{-bt}(C_1 cos(\sqrt{-b^2- 4a}t)+ C_2 sin(\sqrt{-b^2- 4a}t)[/itex] if [itex]b^2+ 4a< 0[/itex], or [itex]x_1(t)= e^{-bt}(C_1+ C_2t)[/itex] if [itex]b^2+ 4a= 0[/itex].<br />
<br />
In any case, [itex]x_2= (ab- bx_1'- abx_1)/ab[/itex].<br />
<br />
A more "sophisticated" method would be to write the differential equations as a matrix equation:<br />
[tex]\begin{pmatrix}x_1 \\ x_2\end{pmatrix}'= \begin{pmatrix}-a & -a \\ -b & -b\end{pmatrix}\begin{pmatrix}x_1 \\ x_2 \end{pmatrix}+ \begin{pmatrix}a \\ b\end{pmatrix}[/tex]<br />
<br />
The characteristic equation for that matrix is <br />
[tex]\left|\begin{array}{cc}-a- \lambda & -a \\ -b & -b- \lambda\end{array}\right|= (-a- \lambda)(-b- \lambda)- ab= \lambda^2+ (a+ b)\lambda= 0[/tex]<br />
which has roots [itex]\lambda= 0[/itex] and [itex]\lambda= -a- b[itex].<br />
We can find the eigenvectors corresponding to those eigenvectors and so diagonalize the matrix,[/itex][/itex][/itex][/itex][/itex][/itex]