%let name=random_scatter;
filename odsout '.';

/*
Using example from "SAS 9.3 Functions and CALL Routines: Reference"
https://books.google.com/books?id=GwTewpZva_YC&pg=PA16&lpg=PA16&dq=sas+%22Nonindependent+Random+Normal+Variables%22&source=bl&ots=4wnY8UQASD&sig=ge-vDf9xKWWxmE8sBQB1_FeJg5w&hl=en&sa=X&ved=0ahUKEwiGt-_yxPnQAhWDVywKHabHD4cQ6AEIJzAC#v=onepage&q=sas%20%22Nonindependent%20Random%20Normal%20Variables%22&f=false
*/


data normal;
seed1 = 11111;
seed2 = 22222;
seed3 = 33333;
do i = 1 to 10000;
 call rannor(seed1, x1);
 call rannor(seed2, x2);
 call rannor(seed3, x3);
 output;
 end;
run;


goptions device=png;
goptions xpixels=600 ypixels=600;
goptions noborder;

ODS LISTING CLOSE;
ODS HTML path=odsout body="&name..htm" 
 (title="SAS Pseudo-Random Numbers (SAS/Graph gplot)") 
 options(pagebreak='no') style=sasweb;

goptions htitle=15pt htext=11pt ftitle="albany amt/bold" ftext="albany amt" ctext=gray33;

symbol1 font='albany amt/unicode' value='25cf'x height=.70
 interpol=none color=A3174bf55 repeat=3;

axis1 color=graybb value=(color=gray55) label=none minor=none offset=(0,0);

title1 ls=1.5 "Nonindependent Random Normal Variables";
proc gplot data=normal;
plot x1*x3 /* x3*x2 x1*x2 */ / 
 haxis=axis1 vaxis=axis1 
 autohref chref=graydd
 autovref cvref=graydd
 des='' name="&name"; 
run;

quit;
ODS HTML CLOSE;
ODS LISTING;