y=lab.triu(x)
returns the upper triangular part of x,
the other elements of y are set to 0.
lab.triu(x,k)
returns the elements on and above the k-th diagonal of x as non-zero. k=0 is the main diagonal, k>0 is above the main diagonal and k<0
is below the main diagonal.