> restart:

> x:=cos(t):

> y:=sin(t):

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

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

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

> t:=n*Pi/50:

> w:=evalf(m):

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

>