7-8

> restart;

> x:=cos(t)/(2-(cos(t))^2):

> y:=sin(t)*cos(t)/(2-(cos(t))^2):

> a:=plot([x,y,t=-Pi..Pi],color=red,scaling=constrained):

> a:;

> x1:=diff(x,t):

> y1:=diff(y,t):

> x2:=diff(x1,t):

> y2:=diff(y1,t):

> f:=x1^2+y1^2:

> g:=x1*y2-y1*x2:

> r:=f^(3/2)/abs(g):

> xc:=x-f/g*y1:

> yc:=y+f/g*x1:

> with(plots):

Warning, the name changecoords has been redefined

> c:=animate([xc+r*cos(s),yc+r*sin(s),s=0..2*Pi],t=-Pi..Pi,color=green,scaling=constrained,frames=30):

> c:

> display(a,c):

> d:=animate([(1-s)*xc+s*x,(1-s)*yc+s*y,s=0..1],t=-Pi..Pi,color=brown,scaling=constrained,frames=30,title="press button here"):

> d:

> display(a,c,d,view=[-1..1,-2..2],axes=none);

>

>