pectively, of the converged eigenvalues;
                this error is unlikely to appear;
          = N+1:  if T has almost reciprocal eigenvalues; perturbed
                values were used to solve Lyapunov equations, but
                the matrix T, if given (for FACT = 'F'), is
                unchanged.

</PRE>
<A name="Method"><B><FONT SIZE="+1">Method</FONT></B></A>
<PRE>
  The condition number of the Riccati equation is estimated as

  cond = ( norm(Theta)*norm(A) + norm(inv(Omega))*norm(Q) +
    