IdentityMat(R: RING, N: INT): MAT
/**/ Id := IdentityMat(QQ,3); Id; matrix(QQ, [[1, 0, 0], [0, 1, 0], [0, 0, 1]]) /**/ type(Id[1,1]); RINGELEM /**/ RingOf(Id[1,1]); QQ