Abstract
A new matrix which is the product of a chain of factored matrices is constructed in this paper to find the matrix Routh array and the matrix continued fraction inversion for multivariate systems if the partial matrix coefficients are known. A linear transformation matrix ia established from tho matrix Routh array developed to transform the second Cauer matrix form in state-space coordinates to tho phase variable matrix form.