> with(plots): with(plottools):
> C:=cone([0,0,5],1,-5,color=green):#display(C);
> d:=3:a:=sqrt(Pi)/2:q:=[0+d,0,5]:
> p1:=[a+d,-a,0]:p2:=[a+d,a,0]:p3:=[-a+d,a,0]:p4:=[-a+d,-a,0]:
> P:=polygon([q,p1,p2]),polygon([q,p2,p3]),polygon([q,p3,p4]),polygon([q,p4,p1]):#display(C,P);
> Disc:=plot([sqrt(1-x^2),-sqrt(1-x^2)],x=-1..1,color=blue,filled=true):#Disc;
> Cir:=circle([0,0],1,color=blue,thickness=3):
> Sq:=polygon([[a,-a],[a,a],[-a,a],[-a,-a]],color=blue,style=line,thickness=3):
> S:=polygon([[a,-a],[a,a],[-a,a],[-a,-a]],color=blue):#display(S);
> frames:=array(0..16):
> for i from 0 to 16 do
> h:=i/16:scale:=1-h:DT:=transform((x,y)->[x*scale,d+y*scale,h*5]):
> ST:=transform((x,y)->[d+x*scale,d+y*scale,h*5]):
> CT:=transform((x,y)->[x*scale,y*scale,h*5]):
> SqT:=transform((x,y)->[d+x*scale,y*scale,h*5]):
> frames[i]:=display(C,P,DT(Disc),ST(S),CT(Cir),SqT(Sq),polygon([[-2,-2,h*5],[-2,5,h*5],[5,5,h*5],[5,-2,h*5]],style=line,thickness=2,color=blue)): od:
> display(seq(frames[i],i=0..16),insequence=true,scaling=constrained,orientation=[60,75]);
>