Process points with PROC G3GRID. SMOOTH= specifies the smoothing parameter to use during spline interpolation.
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;