> with(plots): with(plottools):
> Disc:=plot([sqrt(1-x^2),-sqrt(1-x^2)],x=-1..1,color=blue,filled=true):#Disc;
> Disc2:=plot([sqrt(1-x^2),-sqrt(1-x^2)],x=-1..1,color=white,filled=true):#Disc2;
> frames:=array(0..16):
> for i from 0 to 16 do
> ss:=sqrt(1-(i/16)^2):cs:=i/16:
> sT:=transform((x,y)->[ss*x,ss*y]):cT:=transform((x,y)->[2+cs*x,cs*y]):
> frames[i]:=display(cT(Disc),sT(Disc2),Disc,view=[-1..3,-1..1]):od:
> display(seq(frames[16-i],i=0..16),seq(frames[i],i=1..16),insequence=true,scaling=constrained);
> display(frames[1],scaling=constrained);
>