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