Generate the contour after default interpolation.

proc gcontour data=numdef;
   plot y*x=z / 
      haxis=axis1 
      vaxis=axis1;
   run;
quit;