> with(LinearAlgebra):
> M := Matrix([[1,0,0],[0,1,1],[1,1,1]]);

                               [1    0    0]
                               [           ]
                          M := [0    1    1]
                               [           ]
                               [1    1    1]

> M^2;

                            [1    0    0]
                            [           ]
                            [1    2    2]
                            [           ]
                            [2    2    2]

> M^3;

                            [1    0    0]
                            [           ]
                            [3    4    4]
                            [           ]
                            [4    4    4]

> #?identity
> I3 := Matrix(3,3,shape=identity);

                               [1    0    0]
                               [           ]
                         I3 := [0    1    0]
                               [           ]
                               [0    0    1]

> R:=I3 + M + M^2;

                               [3    0    0]
                               [           ]
                          R := [1    4    3]
                               [           ]
                               [3    3    4]

> map(z->signum(0,z,0),R);

                            [1    0    0]
                            [           ]
                            [1    1    1]
                            [           ]
                            [1    1    1]

>