> t:=k*2*Pi/100;
> with(linalg):
> d:=norm([cos(t),sin(t)]-[1,0],2);
> plot([[x,y,s=0..2*Pi]$k=1..100],color=red, scaling=constrained,axes=none);
> x:=cos(t)+d*cos(s);y:=sin(t)+d*sin(s);
> t:=k*2*Pi/100;
> with(linalg):
> d:=sin(t);
> plot([[x,y,s=0..2*Pi]$k=1..100],color=red, scaling=constrained,axes=none);
> x:=3*cos(t)+d*cos(s);y:=3*sin(t)+d*sin(s);
> t:=k*2*Pi/100;
> with(linalg):
> d:=norm([2*cos(t)+cos(-2*t),2*sin(t)+sin(-2*t)]-[3*cos(t),3*sin(t)],2);
> plot([[x,y,s=0..2*Pi]$k=1..100],color=red, scaling=constrained,axes=none);
> x:=4*cos(t)+d*cos(s);y:=4*sin(t)+d*sin(s);
> t:=k*2*Pi/100;
> with(linalg):
> d:=norm([3*cos(t)+cos(-3*t),3*sin(t)+sin(-3*t)]-[4*cos(t),4*sin(t)],2);
> plot([[x,y,s=0..2*Pi]$k=1..100], color=red,scaling=constrained,axes=none);
> plot([[(1-k/25)*cos(t),k/25*sin(t),t=0..2*Pi]$k=0..25],
color=red,scaling=constrained,axes=none);
> x:=cos(t);
> y:=sin(t);
> m:=[[x,0],[0,y]];
> t:=n*Pi/50;
> plot([m$n=1..100],color=green,
> scaling=constrained,axes=none);
> x:=cos(t);y:=sin(t);
> m:=[[x,y],[cos(2*t),sin(2*t)]];
> t:=n*Pi/200;
> plot([m$n=1..400],color=yellow,scaling=constrained,axes=none);
> restart;
> x:=cos(t);y:=sin(t);
> m:=[[x,y],[cos(3*t),sin(3*t)]];
> t:=n*Pi/150;
> plot([m$n=1..300],color=gold,scaling=constrained,axes=none);
> x:=cos(t);y:=sin(t);
> m:=[x+t*sin(s),y+t*cos(s),s=-Pi..Pi];
> t:=n*Pi/20;
> w:=evalf(m);
> plot([w$n=1..100],color=wheat,scaling=constrained,axes=none);
