> restart;

> x:=cos(t);

> y:=sin(t);

> r:=((x-1)^2+y^2)^(1/2):

> m:=[x+r*sin(s), y+r*cos(s), s=-Pi..Pi]:

> t:=2*n*Pi/100;

> w:=evalf(m):

> plot([w$n=1..100],color=red,axes=none,scaling=constrained);

>