  We propose an algorithm for determining the irreducible polynomials over
finite fields, based on the use of the companion matrix of polynomials and the
generalized Jordan normal form of square matrices.
