   2
    op(A)'*op(U)'*op(U)*op(A) - op(U)'*op(U) = -scale *op(F)'*op(F),

  where U and F are N-by-N upper triangular matrices, and

     F = R,                                  if M &gt;= N, or

     F = ( R ),    if LTRANS = .FALSE.,  or
         ( 0 )

     F = ( 0  Z ), if LTRANS = .TRUE.,       if M &lt; N.
         ( 0  R )

  The canonical equation is solved for U.

</PRE>
<A name="References"><B><FONT SIZE="+1">References</FONT></B></A>
<PRE>
  [1] Bartels, R.H. and Stewart, G.W.
      Solution of 