dyn_sys:=proc(  A::matrix,
                n::integer
              )
   local x,p,j,k,h;

   x:=array(1..4,[seq(array(0..n),i=1..4)]);
   x[1][0]:=vector([1,1]);
   x[2][0]:=vector([1,-1]);
   x[3][0]:=vector([-1,1]);
   x[4][0]:=vector([-1,-1]);

   for k to n do
      for j to 4 do
         x[j][k]:=linalg[multiply](A,x[j][k-1])
      od
   od;

   h:=(max(seq(seq(x[j][k][2],j=1..4),k=0..n))
       -min(seq(seq(x[j][k][2],j=1..4),k=0..n)))*0.03;

   p:=array(1..4,[seq(array(1..n),i=1..4)]);

   for j to 4 do
      for k to n do
         p[j][k]:=plottools[arrow](x[j][k-1],matadd(x[j][k],x[j][k-1],
                  1.0,-1),.003,h,.1)
      od
   od;

   plots[display]({seq(seq(p[j][k],j=1..4),k=1..n)});

end;