input mp-solid

outcolor:=0.5[jaune,white];
incolor:=0.5[vert,white];

figureespace(-10u,-10u,10u,10u);
Initialisation(500,30,20,25);
a=20;b=36;h=8;
zz=-h/2;
ObjetNew1((sqrt(1+((zz+h)**2)/4),0,zz+h)    
  for m=a downto 0:
    for n=if m=a:b-1 else:b fi downto 1:%if k=a:b-1 else:b fi downto 0:
      ,(sqrt(1+((zz+h*m/a)**2)/4)*cosd(n*(360 div b)),sqrt(1+((zz+h*m/a)**2)/4)*sind(n*(360 div b)),zz+h*m/a)
    endfor
  endfor
  )(4,a*b,b+a*b,b-1+a*b,a*b-1
  for m=a downto 1:
    for n=if m=a:m*b-1 else: m*b fi downto (m-1)*b+2:
      ,4,n,b+n,(b-1)+n,n-1
    endfor
    ,4,m*b+1,(m+1)*b,m*b,(m-1)*b+1
  endfor
  );
Ferme1:=false;
AffichageObjet1;
finespace;
end