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

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

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

6 -6 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

%% syntaxe : x_min x_max y_min y_max z_min z_max quadrilleXYZ
0 5 0 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-y-1)* 0.5 * Cos (x+y+1) + 2
end
} def

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

/pas 0.25 def
1 pas 4 0 pas 3 {f} surfaceparam3d