(eng): A numerical procedure is proposed to solve a matrix polynomial frequently encountered in discrete-time control and signal processing. The algorithm is based on simple rewriting of the original equation in terms of a reduced Sylvester matrix. In contrast to previously published methods, it does not make use of elementary polynomial operations. Moreover and most notably, it is numerically reliable.