A computationally simple and efficient method is proposed to determine the stability of linear discrete-time systems. The method directly utilises the system (state-transition) matrix and avoids the evaluation of the characteristic polynomial. Another interesting feature of this approach is that it can be extended to extract the eigenvalues of matrices, in a recursive fashion.