  successful exit;
          &lt; 0:  if INFO = -i, the i-th argument had an illegal
                value.

</PRE>
<A name="Numerical Aspects"><B><FONT SIZE="+1">Numerical Aspects</FONT></B></A>
<PRE>
  The algorithm requires 0(N*K*( P + Q )) floating point operations.

</PRE>
<A name="Method"><B><FONT SIZE="+1">Method</FONT></B></A>
<PRE>
  The Householder transformations and modified hyperbolic rotations
  computed by SLICOT Library routine MB02CU are applied to the
  corresponding parts of the matrices 