     contain the system output matrix C.
          On exit, the leading P-by-N part of this array contains
          the intermediary output matrix C1+BETA*D1*K (see METHOD).

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

  D       (input/output) DOUBLE PRECISION array, dimension (LDD,M)
          On entry, if JOBD = 'D', the leading P-by-M part of this
          array must contain the system direct input/output
     