is real. A more general and simpler algorithm is developed which converges to the matrix sign function of
whatever the spectrum of
may be, and thus answers the question. The numerical stability of the algorithms is discussed and some results are quoted.