Generate the surface plot.

proc g3d data=sashelp.iris;
  scatter PetalLength*PetalWidth=SepalLength;
run;
quit;