[step:Apply an elementary Gronwall estimate to bound every trajectory]Define $A := |x_0|+aT$. For the fixed index $k$, define $y_k: [t_0,t_1] \to [0,\infty)$ by $y_k(t) := |x_k(t)|$. The previous step gives
\begin{align*}
y_k(t) \leq A+b\int_{t_0}^{t} y_k(s)\,d\mathcal{L}^1(s)
\end{align*}
for every $t \in [t_0,t_1]$.
We use the following elementary integral estimate. If $y: [t_0,t_1]\to [0,\infty)$ is integrable and satisfies
\begin{align*}
y(t) \leq A+b\int_{t_0}^{t} y(s)\,d\mathcal{L}^1(s)
\end{align*}
for every $t \in [t_0,t_1]$, then
\begin{align*}
y(t) \leq Ae^{b(t-t_0)}
\end{align*}
for every $t \in [t_0,t_1]$.
Indeed, define $z: [t_0,t_1]\to [0,\infty)$ by
\begin{align*}
z(t):=A+b\int_{t_0}^{t} y(s)\,d\mathcal{L}^1(s).
\end{align*}
Then $z$ is absolutely continuous, $y(t)\leq z(t)$ for every $t$, and
\begin{align*}
z'(t)=by(t)\leq bz(t)
\end{align*}
for $\mathcal{L}^1$-almost every $t \in [t_0,t_1]$. Define $w: [t_0,t_1]\to \mathbb{R}$ by $w(t):=e^{-b(t-t_0)}z(t)$. Since $z$ is absolutely continuous and the exponential factor is smooth, $w$ is absolutely continuous, and for $\mathcal{L}^1$-almost every $t$,
\begin{align*}
w'(t)=e^{-b(t-t_0)}(z'(t)-bz(t))\leq 0.
\end{align*}
Therefore $w(t)\leq w(t_0)=A$ for every $t \in [t_0,t_1]$, so $z(t)\leq Ae^{b(t-t_0)}$. Since $y(t)\leq z(t)$, the estimate follows.
Applying this estimate to $y_k$ gives
\begin{align*}
|x_k(t)|=y_k(t)\leq Ae^{b(t-t_0)}\leq Ae^{bT}
\end{align*}
for every $t \in [t_0,t_1]$. Thus, with
\begin{align*}
R:=(|x_0|+aT)e^{bT},
\end{align*}
we have $|x_k(t)|\leq R$ for every $k \in \mathbb{N}$ and every $t \in [t_0,t_1]$.[/step]