> restart;

> m:=[[1,0],[cos(2*Pi/3),sin(2*Pi/3)],[cos(4*Pi/3),sin(4*Pi/3)],[1,0]];

>

> w:=[]:

> t:=0.1:

> for k to 20 do w:=[op(w),m]:sm:=[m[2],m[3],m[1],m[2]]:m:=(1-t)*m+t*sm:od:

> plot(w,axes=none,scaling=constrained,color=blue);

>