function B = triuu( A )

     B = triu( A, 1 ) + eye( size( A ) );

return
