    On exit, the leading P-by-NR part of this array contains
          the state/output matrix Cr of the reduced order system.

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

  D       (input/output) DOUBLE PRECISION array, dimension (LDD,M)
          On entry, the leading P-by-M part of this array must
          contain the input/output matrix of the original system.
          On exit, the leading P-by-M part of this array contains
          the input/output matrix Dr o