Numerical Solution of Algebraic Riccati Equations