plot([[n*cos(t),n*sin(t),t=0..2*Pi]$n=1..20],color=blue,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);

plot([[x,y,s=0..2*Pi]$k=1..100]

,color=plum,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);

plot([[x,y,s=0..2*Pi]$k=1..100],

color=gray,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);

plot([[x,y,s=0..2*Pi]$k=1..100],color=pink,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);

plot([[x,y,s=0..2*Pi]$k=1..100],

color=pink,scaling=constrained,axes=none);

restart;

plot([[(1-k/25)*cos(t),k/25*sin(t),t=0..2*Pi]

$k=0..25],color=blue,

scaling=constrained,axes=none);

　

　

restart;

x:=cos(t);

y:=sin(t);

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

t:=n*Pi/50;

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

restart;

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=blue,scaling=constrained,axes=none);

　

restart;

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

t:=n*Pi/150;

plot([m$n=1..300],color=pink,

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

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