input mp-solid

vardef Transform(expr PT)=
  RotZ(RotX(RotY(PT+TR)))-TR
enddef;

transformation:=true;

%1'25" :(
figureespace(-10u,-10u,10u,10u);
Initialisation(1000,120,20,10);
%--les atomes
nb:=24;subh:=20;
outcolor:=noir;
Objetsphere1("R=7");%C
TR:=(0,0,16.1);
outcolor:=ciel;
Objetsphere2("R=3.2");%H1
angy:=109.5;
Objetsphere3("R=3.2");%H2
angz:=120;
Objetsphere4("R=3.2");%H2
angz:=-120;
Objetsphere5("R=3.2");%H2
%--les liaisons
TR:=(0,0,0);
angy:=0;
angz:=0;
nb:=12;subh:=6;
outcolor:=gris;
Objetcylindre6("r=0.25","h=16.1");
angy:=109.5;
Objetcylindre7("r=0.25","h=16.1");
angz:=120;
Objetcylindre8("r=0.25","h=16.1");
angz:=-120;
Objetcylindre9("r=0.25","h=16.1");
nbobj:=9;
DessineFusion;
finespace;
end