I'm looking for interesting applications of companion matrices. I can also use the Frobenius Normal Form.
I already covered the Cayley-Hamilton Theorem and the application to linearly recursive sequences and high-order scalar linear differential equations.
Companion matrices can also be used to find the roots of a polynomial equation: the eigenvalues of the Companion Matrix are the roots of the polynomial equation. It is a very robust way of finding the roots of an equation, but not the most efficient (computationally).