Construct the circles with center at
(cos(t),sin(t)) and and passing through
(2cos(t)/3+cos(2t)/3,2sin(t)/3+sin(2t)/3) with
t ranging over [0,2p]. restart;

> 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:=k*Pi/50;

> w:=evalf(m):

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

>