input mp-solid

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

%
figureespace(-100u,-100u,100u,100u);
Initialisation(500,75,20,10);
creux:=false;
%nb:=18;subh:=12;
nb:=20;subh:=16;
%----Le soufre
outcolor:=jaune;
Objetsphere1("R=3.56");%S1
TR:=(0,0,20.1);
Objetsphere2("R=3.56");%S2
%----L'oxygène
transformation:=true;
outcolor:=rouge;
TR:=(0,0,15);
angy:=-109.5;
angz:=120;
Objetsphere3("R=2.17");
angz:=0;
Objetsphere4("R=2.17");
angz:=-120;
Objetsphere5("R=2.17");
nb:=10;subh:=10;
%---La liaison simple O-S
outcolor:=rouge;
TR:=(0,0,7.5);
angz:=120;
Objetcylindre7("r=0.5","h=7.5");
outcolor:=jaune;
TR:=(0,0,0);
Objetcylindre8("r=0.5","h=7.5");
%Les liaisons double O-S
%1ere liaison
outcolor:=rouge;
TR:=(0.75,0,7.5);
angz:=0;
Objetcylindre9("r=0.5","h=7.5");
outcolor:=jaune;
TR:=(0.75,0,0);
Objetcylindre10("r=0.5","h=7.5");
outcolor:=rouge;
TR:=(-0.75,0,7.5);
Objetcylindre11("r=0.5","h=7.5");
outcolor:=jaune;
TR:=(-0.75,0,0);
Objetcylindre12("r=0.5","h=7.5");
%2eme liaison
outcolor:=rouge;
TR:=(0.75,0,7.5);
angz:=-120;
Objetcylindre13("r=0.5","h=7.5");
outcolor:=jaune;
TR:=(0.75,0,0);
Objetcylindre14("r=0.5","h=7.5");
outcolor:=rouge;
TR:=(-0.75,0,7.5);
Objetcylindre15("r=0.5","h=7.5");
outcolor:=jaune;
TR:=(-0.75,0,0);
Objetcylindre16("r=0.5","h=7.5");
transformation:=false;
%---La liaison S-S
TR:=(0,0,0);
angx:=0;angy:=0;angz:=0;
outcolor:=jaune;
Objetcylindre6("r=0.5","h=20.1");
%%%affichage
TR:=(0,0,0);
angx:=0;angy:=0;angz:=0;
nbobj:=16;
DessineFusion;
TraceAxesD(15,15,25);
finespace;
end