mplectic factorizations.
      BIT, 43 (4), pp. 775-790, 2003.

</PRE>
<A name="Numerical Aspects"><B><FONT SIZE="+1">Numerical Aspects</FONT></B></A>
<PRE>
  The algorithm requires 16*( M - K )*N + ( 26*K - 4 )*K*N floating
  point operations if STRAB = 'Z' and additional ( 12*K + 2 )*K*N
  floating point operations if STRAB = 'N'.

</PRE>

<A name="Comments"><B><FONT SIZE="+1">Further Comments</FONT></B></A>
<PRE>
  None
</PRE>

<A name="Example"><B><FONT SIZE="+1">Example</FONT></B></A>
<P>
<B>Program Te