f the reduced order system.
          If NR = 0 and the given system is stable, then D contains
          the steady state gain of the system.

  LDD     INTEGER
          The leading dimension of array D.  LDD &gt;= MAX(1,P).

  RCOND   (output) DOUBLE PRECISION
          The reciprocal condition number of the matrix A22-g*I
          (see METHOD).

</PRE>
<B>Workspace</B>
<PRE>
  IWORK   INTEGER array, dimension 2*(N-NR)

  DWORK   DOUBLE PRECISION array, dimension 4*(N-NR)

</PRE>
<B>Error Indicator</B>
