%% 
%%  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}