proc g3grid data=reflib.nums out=numdef;
   grid y*x=z / axis1=-5 to 5 by .5
                axis2=-5 to 5 by .5;
run;