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