A method is presented to solve the real algebraic Riccati equation -

, where

and

. The solution for the corresponding eigenvalue problem

, where

is a Hamiltonian matrix, is computed by an algorithm similar to the QR algorithm. Special symplectic matrices are used for the transformation of

such that the Hamiltonian form is preserved during the computations.