%% %% Ein Beispiel der DANTE-Edition %% %% 7. Auflage %% %% Beispiel 35-00-4 auf Seite 818. %% %% Copyright (C) 2016 Herbert Voss %% %% It may be distributed and/or modified under the conditions %% of the LaTeX Project Public License, either version 1.3 %% of this license or (at your option) any later version. %% %% See http://www.latex-project.org/lppl.txt for details. %% %% %% ==== % Show page(s) 1 %% %% \documentclass[]{exaarticle} \pagestyle{empty} \setlength\textwidth{193.16928pt} \setlength\parindent{0pt} %StartShownPreambleCommands \usepackage{pst-solides3d} %StopShownPreambleCommands \begin{document} \psset{viewpoint=80 -130 20 rtp2xyz,Decran=50,lightsrc=viewpoint} \begin{pspicture}(-7,-8)(7,8) {\psset{object=courbe,r=0,linecolor=blue,resolution=360,range=-5 5} \multido{\rA=-4.0+0.5}{16}{% \defFunction[algebraic]{Fx}(x)% {x}% {\rA}% { 10*(x^3+x*(\rA)^4-x/5)*Euler^(-x^2-(\rA)^2)+Euler^(-((x-1.225)^2+(\rA)^2)) } \psSolid[function=Fx](0,5 \rA\space sub,0)} \multido{\rA=-4.0+0.5}{16}{% \defFunction[algebraic]{Fy}(y)% {\rA}% {y}% { 10*((\rA)^3+(\rA)*y^4-(\rA)/5)*Euler^(-(\rA)^2-y^2)+Euler^(-(((\rA)-1.225)^2+y^2)) } \psSolid[function=Fy](5 \rA\space sub ,0,0)}} \psSurface[ngrid=.15 .15,fillcolor=green!30,incolor=gray!30, linewidth=0.4pt,algebraic](-4,-4)(4,4)% { 10*(x^3+x*y^4-x/5)*Euler^(-x^2-y^2)+Euler^(-((x-1.225)^2+y^2)) } \axesIIID[linecolor=gray](0,0,0)(7,7,5) \end{pspicture} %the equation in postfix notation % x 3 exp x y 4 exp mul add x 5 div sub 10 mul % 2.729 x dup mul y dup mul add neg exp mul % 2.729 x 1.225 sub dup mul y dup mul add neg exp add \end{document}