Generate the surface plot.

proc g3d data=iris;
note j=r f="Albany AMT/bo" "Species:  : c=green "Virginica   "
         j=r    c=red "Versicolor      "
         j=r    c=blue "Setosa      ";
 scatter PetalLength*PetalWidth=SepalLength/
      color=color
      shape=shape;
run;
quit;