One has a ODE system of the form
$\frac{d}{dt}\begin{bmatrix}x\\y\\z\end{bmatrix}=\begin{bmatrix}a&b&c\\d&e&f\\g&h&i\end{bmatrix}\begin{bmatrix}x\\y\\z\end{bmatrix}$
Then, for example, to solve this by forward difference Euler method,
how to derive numerical stability condition,
which relates parameters $a$ to $i$ with sample rate i.e. time step $k$ ?