Generate the surface plot.
proc g3d data=reflib.hat;
   plot y*x=z;
run;
quit;