Leverrier's algorithm for singular systems

An algorithm is presented which allows the computation of the transfer function of a singular system from its state-space description, without inverting a polynomial matrix. This algorithm is an extension of the Leverrier algorithm for the more general case of singular systems and it reduces the computational cost.