> x := (e,t)-> (1+e*(cos(3*t)+cos(5*t)/3))*cos(t)+e*(3*sin(3*t)+5/3*sin(5*t))*sin(t):

> y := (e,t) -> (1+e*(cos(3*t)+cos(5*t)/3))*sin(t)-e*(3*sin(3*t)+5/3*sin(5*t))*cos(t):

> c := (e,t) -> x(e,t)/sqrt((x(e,t))^2+(y(e,t))^2): s:= (e,t) -> y(e,t)/sqrt((x(e,t))^2+(y(e,t))^2):

> with(plots):plotoptions("transparent=true"):

> animate({[(1+(.05)*(cos(3*(u*Pi/40))+cos(5*(u*Pi/40))/3))*cos((u*Pi)/40)-s*sin((u*Pi)/40),(1+(.05)*(cos(3*(u*Pi)/40)+cos(5*(u*Pi)/40)/3))*sin((u*Pi)/40)+s*cos((u*Pi)/40),s=-2..2],[x(.05,ss),y(.05,ss),ss=0..2*Pi]},u=0..79);

> animate({[s(.05,(u*Pi)/40)*x(.05,t)-c(.05,(u*Pi)/40)*y(.05,t),c(.05,(u*Pi)/40)*x(.05,t)+s(.05,(u*Pi)/40)*y(.05,t)-((c(.05,(u*Pi)/40)*x(.05,(u*Pi)/40)+s(.05,(u*Pi)/40)*y(.05,(u*Pi)/40))+(c(.05,(u*Pi)/40)*x(.05,((u+40)*Pi)/40)+s(.05,(u*Pi)/40)*y(.05,((u+40)*Pi)/40)))/2,t=0..2*Pi],[v,1,v=-1.3..1.3],[l,-1,l=-1.3..1.3]}, u=0..79,view=[-1.3..1.3,-1.3..1.3],axes=NONE,color=red,thickness=2);

>

>