In order to prove that solutions of linear equations with constant coefficients grow at most exponentially, you can use a multidimensional version of Gronwall's lemma.
Remember indeed that the solution
)
of a

-th order linear differential equation is a component of the solution
)
of a

-dimensional linear equation. Then an exponential bound on the euclidean norm of
)
gives an exponential bound on
)
. Here's how to obtain such a bound.
If
=AY(t)+B)
where

are

-matrices, we have
\|\leq\|A\|\|Y(t)\|+\|B\|)
for every

, so that
-Y(0)\|\leq \|A\| \int_0^t \|Y(s)\| ds + \|B\| t)
, hence
\|\leq \|Y(0)\|+\|B\|t+\|A\|\int_0^t \|Y(s)\| ds)
. And now you can apply the one-dimensional Gronwall lemma "in integral form" (quoting the
wikipedia) to the function
\|)
. The only online reference I could find for the final result is this french
exercise sheet (Exercice 1.4), but you don't need it to find the proof, just apply Gronwall's lemma like I wrote and do an integration by parts to compute the integral. You should get:
\|\leq \|Y(0)\|e^{\|A\|t}+\frac{\|B\|}{\|A\|}(e^{\|A\|t}-1))
.
This is what you need.