p:=plot([[n*cos(t),n*sin(t),t=0..2*Pi]$n=1..20],color=red,axes=none):

restart;

x:=cos(t)+d*cos(s);y:=sin(t)+d*sin(s);

t:=k*2*Pi/100;

with(linalg):

d:=norm([cos(t),sin(t)]-[1,0],2);

p:=plot([[x,y,s=0..2*Pi]$k=1..100],color=red,scaling=constrained,axes=none):

restart;

x:=cos(t)+d*cos(s);y:=sin(t)+d*sin(s);

t:=k*2*Pi/100;

with(linalg):

d:=sin(t);

p:=plot([[x,y,s=0..2*Pi]$k=1..100],color=red,scaling=constrained,axes=none):

restart;

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);

p:=plot([[x,y,s=0..2*Pi]$k=1..100],color=red,scaling=constrained,axes=none):

restart;

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);

p:=plot([[x,y,s=0..2*Pi]$k=1..100],color=red,scaling=constrained,axes=none):

restart;

p:=plot([[(1-k/25)*cos(t),k/25*sin(t),t=0..2*Pi]$k=0..25],color=red,scaling=constrained,axes=none):

restart;

x:=cos(t);

y:=sin(t);

m:=[[x,0],[0,y]];

t:=n*Pi/50;

p:=plot([m$n=1..100],color=red,scaling=constrained,axes=none):

p;

restart;

x:=cos(t);

y:=sin(t);

m:=[[x,y],[cos(2*t),sin(2*t)]];

t:=n*Pi/200;

p:=plot([m$n=1..400],color=red,scaling=constrained,axes=none):

restart;

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

t:=n*Pi/150;

p:=plot([m$n=1..300],color=red,scaling=constrained,axes=none):

restart;

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):

p:=plot([m$n=1..100],color=red,scaling=constrained,axes=none):