. mat a = I(3) . mat list a symmetric a[3,3] c1 c2 c3 r1 1 r2 0 1 r3 0 0 1