Computer codes for computational assistance in the study of asymptotic stability

AbstractThe Routh-Hurwitz, Liénard-Chipart and Stodola criteria are reviewed and then combined into an algorithm in order to represent efficiently Liapunov's sufficient conditions for asymptotic stability. Three computer subprograms, which were written in CDC Extended FORTRAN IV to encode this algorithm, are presented. To assist the user further, one of the routines computes the coefficients of a characteristic polynomial by the method of Leverrier. As a demonstration of the computational assistance provided by this software package, the nine-dimensional domain of asymptotic stability of a system of four o.d.e.s is estimated; these equations have been put forth by other authors to simulate the moose and wolf interactions on Isle Royale, Michigan