%@AUTEUR: Maxime Chupin %@DATE: 31 mars 2008 \documentclass[a4paper]{article} \usepackage{pst-solides3d} \usepackage{xcolor} \usepackage[margin=2cm]{geometry} \definecolor{moblue}{rgb}{0.06,0.247,0.662} \pagestyle{empty} \begin{document} \multido{\iRotZ=1+5}{73}{% \begin{pspicture}(-6,-6)(6,6) \psset{SphericalCoor=true,viewpoint=12 4 19,Decran=40,linewidth=0.1pt} \psset{solidmemory} \psset{lightsrc=55 5 20} \defFunction[algebraic]{F}(t){t}{}{3/16*(2*t*t-16)} \defFunction[algebraic]{g}(t) {(2/3*((3-cos(1.668))*cos(t)+cos(1.668)*cos(t)*cos(3*t)+sin(t)*sin(3*t)))} {(2/3*((3-cos(1.668))*sin(t)+cos(1.668)*sin(t)*cos(3*t)-cos(t)*sin(3*t)))} {(2/3*sin(1.668)*(1-cos(3*t)))} \defFunction[algebraic]{g'}(t) {((2*(1.291152019431635*cos(t)*sin(3*t)+3.097050673143878*sin(t)*cos(3*t)-3.097050673143878*sin(t)))/3)} {((2*(1.291152019431635*sin(t)*sin(3*t)-3.097050673143878*cos(t)*cos(3*t)+3.097050673143878*cos(t)))/3)} {(1.990558883170573*sin(3*t))} \defFunction[algebraic]{g''}(t) {((2*(-10.58230403886327*sin(t)*sin(3*t)+6.970506731438783*cos(t)*cos(3*t)-3.097050673143878*cos(t)))/3)} {((2*(10.58230403886327*cos(t)*sin(3*t)+6.970506731438783*sin(t)*cos(3*t)-3.097050673143878*sin(t)))/3)} {(5.97167664951172*cos(3*t))} \codejps{ /RotRad \iRotZ\space 180 div 3.14159 mul def} \psSolid[object=courbe, range=0 RotRad, ngrid=200 16, fillcolor=orange, function=g, name=A1, action=none, r=0.03] \psSolid[object=cone, h=2.07,r=0.666667, % fillcolor=moblue!70, hue=0 1 0.5 1, name=C2, ngrid=1 30, action=none] \codejps{% /RotAxeCone \iRotZ\space 3 mul def /RotZC \iRotZ\space 90 add def /C20{ C2 {0 0 RotAxeCone rotateOpoint3d } solidtransform } def /C21{ C20 {-84.5 0 0 rotateOpoint3d } solidtransform } def /C22{ C21 {0 -2.06459 0.67 translatepoint3d} solidtransform } def /C23{ C22 {0 0 RotZC rotateOpoint3d} solidtransform } def } \psSolid[object=datfile,file=spherecone,name=CS,action=none] \psSolid[object=fusion,base=A1 C23 CS] \composeSolid \end{pspicture} \newpage} \end{document}