> restart;

> m:=[[3*cos(t)+cos(3*t),3*sin(t)-sin(3*t)],[3*cos(t)-cos(3*t),3*sin(t)+sin(3*t)]];

> t:=n*Pi/50;

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

>