Abstract. One of the common methods to solve algebraic Riccati equations is using matrix sign function to calculate stable invariant subspaces. In this paper, we show that it is possible to compute the matrix sign function using the quadratically convergent Kovarik’s method. Since the explicit inversion of a matrix, in every iteration, must be computed, this will introduce unstability in numerical results.