[
    $rapport:=4, angle3:=0, courbe:=Nil,
    $angle1:=arcsin(rapport/(rapport+1)), $angle2:=arcsin(1/(rapport+1)),
    $H:=6, $R:=H*tan(angle1), $h:=7,$r:=h*tan(angle2),
  
    $C1:=cone(Origin,-H*vecK,R,45),
    $x:=angle1+angle2,
    $u:=[exp(i*angle3)*sin(x),-cos(x)],
    $C2:=cone(Origin, h*u,r,30),
    $A:=h*[sin(angle1),-cos(angle1)]/cos(angle2),
    for k from 1 to 72 do
	$angle1':=angle3*rapport,
	$u:=[exp(i*angle3)*sin(x),-cos(x)],
	A':=rot3d(A,[Origin,vecK],angle3),
	A'':=rot3d(A',[Origin,-u],angle1'),
	Insert(courbe,A''),Inc(angle3,pi/36)
    od,
    L:=lineTube(courbe,0.15,12),
    S:=Sphere(Origin,Norm(A),30,30),
    Build3D(
	[1,Rgb(0.74,1,0.73), C1],
	[1, aliceblue+0.6*i,S],
	[1,Rgb(0.74,0.73,1), C2],
	[1,red-i, L]
    ),
    Display3D()
]