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

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