> restart;

> x:=cos(t);

> y:=sin(t);

> x1:=1/2*x+1/2*cos(t):y1:=1/2*y-1/2*sin(t):

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

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

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

> w:=evalf(m):

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

>