bove. Then, the required triangular part of the result can be
  written as

     alpha*V + beta*tri( B )  + beta*(stri( B ))' +
              alpha*diag( V ) + beta*diag( tri( B ) ).

</PRE>
<A name="References"><B><FONT SIZE="+1">References</FONT></B></A>
<PRE>
  None.

</PRE>
<A name="Numerical Aspects"><B><FONT SIZE="+1">Numerical Aspects</FONT></B></A>
<PRE>
  The algorithm requires approximately

      2             2
     M x N + 1/2 x N x M

  operations.

</PRE>

<A name="Comments"><B><FONT SIZE="+1