> restart:

> x1:=cos(t):

> y1:=sin(t):

> x2:=cos(3*t):

> y2:=sin(3*t):

> m:=[[x1,y1],[x2,y2]]:

> t:=n*Pi/50:

> plot([m$n=1..100],color=blue,scaling=constrained,axes=none);

>