                ( 1 ),
     H  = I - tau *u *u',    u  = ( v )
      i          i  i  i      i   (  i)

  where v  is an (M-P+I-2)-vector.  The components of v  are stored
         i                                             i
  in the i-th row of A, beginning from the location i+1, and tau
                                                                i
  is stored in TAU(i).

</PRE>
<A name="Numerical Aspects"><B><FONT SIZE="+1">Numerical Aspects</FONT></B></A>
<PRE>
  The algorithm is backward stable.