Doubling the size of the matrices, we can solve second-order problems. Here is the constant-coefficient problem of three coupled harmonic oscillators.