Create the scatter plots. On the PLOT statement, the PBSPLINE option fits the spline curves to the scatter points.

proc sgscatter data=sashelp.iris(where=(species="Virginica"));
  title "Multi-Celled Spline Curve for Species Virginica";
  plot (sepallength sepalwidth)*(petallength petalwidth)
       / pbspline;
run;