Construct this figure:

>

> restart;

> x1:=cos(t);

> y1:=sin(t);

> x3:=cos(3*t);

> y3:=sin(3*t);

> m:=[[x1,y1],[x3,y3]];

> t:=k*Pi/50;

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

>