Numerical determination for stability of Hessenberg matrix

In this paper, the reduction with a scaling strategy from Hessenberg matrix to Frobenius-like form is given. An estimation of backward error propagation and numerical tests demonstrate that the reduction to Frobenius-like form and the operation on quasi-Routh array are numerically stable and they can be addressed for numerical determination of matrix stability.