An algorithm for approximately diagonalizing a matrix

An algorithm is presented for selecting a sequence of similarity transformations which successively reduce the norm of the difference between the transformed matrix and a desired form. The algorithm is useful in reducing the cost of simulating the transient response of linear, time-invariant, dynamical systems on the digital computer and can be used to produce canonical forms.