autocrop
%% PhS
%%
%% échelle et étendue de l'image
30 setxunit
-1 9 setxrange
-5 3 setyrange

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%% Positionnement de la Caméra
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

/ScreenDist 0.1 def
/ZoomFactor_x 100 def
/ZoomFactor_y 100 def

5 -8 4 SetCamPos
-1 1.1 -0.3 SetCamVec
0.03 0.03 1 SetCamUp
ComputeCamera

%% traces

1.5 setlinewidth
190 255 div 190 255 div 200 255 div setrgbcolor

-4%% syntaxe : x_min x_max y_min y_max z_min z_max quadrilleXYZ
0 5 -5 5 0 0 quadrilleXYZ

2 setlinewidth
/arrowscale {2 2} def

%% syntaxe : xyz_min xyz_max lu axesRVB
%% (étendue du pointille, longueur des vecteurs)
%% trace des vecteurs unitaires et des axes
-2 5 1 axesRVB

%% la fonction z = f (x, y)
/f { % x y
2 dict begin
   /y exch def
   /x exch def
#rpn# Cos (x^2 + y^2) * 4 * Exp (-x^2 - y^2)
end
} def

1.2 setlinewidth
/coul1 {120 255 div 70 255 div 9 255 div setrgbcolor} def
coul1

/pas 0.25 def
0 pas 4 -4 pas 4 {f} surfaceparam3d