
  In both cases, an orthogonal rotation U1' is constructed such that
  U1'*A1*U1 has the appropriate form. Let UT be an extension of U1
  to an N-by-N orthogonal matrix, using identity submatrices. Then A
  is replaced by UT'*A*UT and the contents of array U is U * UT.

</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 Text</B>
<PRE>
  None
</PRE>
<B>Program Data</B>
<PRE>
  None
</PR