Process points with PROC G3GRID. The SPLINE option specifies the bivariate spline method for the data set interpolation.

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