OPTIONS PS=66 LS=72 NODATE NONUMBER;

/*
FILENAME GRAPH 'c:\papers\khuri\cont1.eps';
GOPTIONS HSIZE=3.5IN VSIZE=4IN 
FTEXT=SWISSB HTEXT=.3CM
GSFNAME=GRAPH GSFLEN=132 GSFMODE=replace 
DEVICE=psepsf;
*/
/*
FILENAME GRAPH 'C:\papers\khuri\cont2.eps';
GOPTIONS HSIZE=4IN VSIZE=4IN 
FTEXT=SWISSB HTEXT=.3CM
GSFNAME=GRAPH GSFLEN=132 GSFMODE=replace 
DEVICE=psepsf;
*/

* GOPTIONS FTEXT=SWISSB HTEXT=.25CM;

*****************************;
*** MAKING A CONTOUR PLOT ***;
*****************************;

DATA cont;
   DO x1 = -1 TO 1 BY .1; * .05;
   DO x2 = -1 TO 1 BY .1; * .05;
      spv = 5 + 4.5*(x1**4 + x2**4 - x1**2 - x2**2 + .5*(x1**2)*(x2**2));
      OUTPUT;
   END; END;

PROC GCONTOUR DATA=cont;
     PLOT x2*x1=spv   / LEVELS= 3.3 4.3 4.8 5.3 5.8 6.8 
          LLEVELS = 33 3 1 2 34 4   
          CLEVELS=BLACK BLACK BLACK BLACK BLACK BLACK ;

PROC G3D DATA=cont;
     PLOT x2*x1=spv / GRID ROTATE=60 
                      ZMIN=3.2 ZTICKNUM=5
                      XTICKNUM=5 YTICKNUM=5; 
RUN;