>> n=7


>> A=diag(ones(1,n))        **staci hociaku maticu n*n, ide len o to aby bola vytvorena

A =

     1     0     0     0     0     0     0
     0     1     0     0     0     0     0
     0     0     1     0     0     0     0
     0     0     0     1     0     0     0
     0     0     0     0     1     0     0
     0     0     0     0     0     1     0
     0     0     0     0     0     0     1


>> for ii=0:(n-1) A(ii+1,(ii+1):end)=1:(n-ii); A((ii+1):end,ii+1)=1:(n-ii); end
>> A

A =

     1     2     3     4     5     6     7
     2     1     2     3     4     5     6
     3     2     1     2     3     4     5
     4     3     2     1     2     3     4
     5     4     3     2     1     2     3
     6     5     4     3     2     1     2
     7     6     5     4     3     2     1