> m:=[[-3,sin(t)],[cos(t),sin(t)],[cos(t)-cos(2*t),sin(t)-sin(2*t)]]:
> t:=k*Pi/50-Pi/2:
> plot([m$k=0..50],scaling=constrained,axes=none,color=red);