* SAS script to simulate two sets of
  random IQ scores (mu=100, sigma=15);

options linesize=70 nodate pageno=1;
data random_iq;
   title "Simulated IQ scores";
   do i = 1 to 5;
      x = rand("normal", 100, 15);
      y = rand("normal", 100, 15);
      output;
   end;
   drop i;

proc corr;
proc gplot;
   plot y * x;

run;
quit;