Process points with PROC G3GRID. SPLINE specifies the bivariate spline method for the data set interpolation. 
proc g3grid data=reflib.nums out=numspl;
   grid y*x=z / spline
                axis1=-5 to 5 by .5
                axis2=-5 to 5 by .5;
run;