Construct this figure:

> x:=cos(t):
y:=sin(t):
x1:=3/4*x+1/4*cos(3*t):y1:=3/4*y-1/4*sin(3*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=yellow,axes=none,scaling=constrained);

maple