proc sgplot data=iris;
title 'Histogram of petal length with density curves';
histogram p_length;
density p_length;
density p_length / type=kernel;
  keylegend / location=inside position=topright;
run;
quit;
title;