On the reduction of matrix functions

Straightforward reduction of matrix functions without explicit use of the Lagrange-Sylvester (LS) theorem and of the characteristic roots is achieved by expansions in terms of Lucas polynomials. Use of the LS theorem within the framework of the Lucas polynomial method and, furthermore, of circulant matrices leads to substantial simplifications.