A trick that is standard in my little world is this: the matrix
$$
M =
\left(
\begin{array}{rr}
-2 & -4 \\
4 & 6
\end{array}
\right)
$$
has trace $4$ and determinant $4.$ The characteristic roots satisfy $\lambda^2 - 4 \lambda + 4 = 0.$ The Cayley-Hamilton Theorem (if this is not familiar, see the ADDENDUM) says that
$$ a_{n+2} = 4 a_{n+1} - 4 a_n, $$
$$ b_{n+2} = 4 b_{n+1} - 4 b_n. $$
It is easy enough to confirm these with direct calculations.
Because of the repeated characteristic value $2,$ we get $a_n = A 2^n + B n 2^n,$ with $b_n = C 2^n + D n 2^n.$
Calculating the first few of each to solve for the coefficients, we get
$$ a_n = 2^n - 2n 2^n, \; \; \; \; b_n = 2n 2^n. $$
ADDENDUM:
Not everyone has seen Cayley-Hamilton. I did say it could be confirmed by straightforward calculation:
Suppose we have the system
$$ \color{blue}{ a_{n+1} = \alpha a_n + \beta b_n,}$$
$$ \color{blue}{ b_{n+1} = \gamma a_n + \delta b_n.} $$
We will find $a_{n+2}$ in two slightly different ways.
$$ a_{n+2} = \alpha a_{n +1} + \beta b_{n +1} = \alpha(\alpha a_n + \beta b_n) + \beta ( \gamma a_n + \delta b_n) = (\alpha^2 + \beta \gamma) a_n +(\alpha \beta + \beta \delta) b_n $$
Let me go straight to this, define
$$ \Psi = (\alpha + \delta) a_{n+1} - (\alpha \delta - \beta \gamma) a_n, $$
$$ \Psi = (\alpha + \delta)( \alpha a_n + \beta b_n) - (\alpha \delta - \beta \gamma) a_n, $$
$$ \Psi = (\alpha^2 + \alpha \delta) a_n + (\alpha \beta + \beta \delta)b_n - (\alpha \delta - \beta \gamma) a_n, $$
$$ \Psi = (\alpha^2 + \beta \gamma) a_n + (\alpha \beta + \beta \delta)b_n. $$
From
$$ a_{n+2} = (\alpha^2 + \beta \gamma) a_n +(\alpha \beta + \beta \delta) b_n $$
we find
$$ a_{n+2} = \Psi, $$
or
$$ \color{blue}{ a_{n+2} = (\alpha + \delta) a_{n+1} - (\alpha \delta - \beta \gamma) a_n.} $$
An analogous calculation works for $b_{n+2}= (\alpha + \delta) b_{n+1} - (\alpha \delta - \beta \gamma) b_n .$