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