Generalization of Leverrier's algorithm to polynomial matrices of arbitrary degree

Leverrier's algorithm to find the determinant of ( sI_{r}-A ) is generalized to the polynomial matrices of arbitrary degree which are in the form I_{r}s^{n} + H_{n-1}s^{n-1} + ... + H_{1}s + H_{0} .