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