proc idmat lng . mat[][] .
   len mat[][] lng
   for i to lng
      len mat[i][] lng
      mat[i][i] = 1
   .
.
idmat 4 m[][]
print m[][]