 of the algorithm), as
          follows:
          = 'A': use all poles with non-negative imaginary parts
                 and all midpoints;
          = 'P': use part of the poles and midpoints.

</PRE>
<B>Input/Output Parameters</B>
<PRE>
  N       (input) INTEGER
          The order of the system.  N &gt;= 0.

  M       (input) INTEGER
          The column size of the matrix B.  M &gt;= 0.

  P       (input) INTEGER
          The row size of the matrix C.  P &gt;= 0.

  RANKE   (input) INTEGER
         