goptions reset=all; 
%let n=1000; %let r=20;
Axis1 label=(angle=90 f=centb h=1.2) order= -3 -2 -1 0 1 2 3;
Axis2 label=(angle=90 f=centb h=1.2);


data anno; 
  xsys='2'; ysys='2'; 
  y=-3; x = (1-.2*exp(-1*Y*Y)); function="move"; output; 
  function = "draw"; size=3;
do Y=-3 to 3 by .1; 
  x = (1-.2*exp(-1*Y*Y)); output; 
end;
 
data a; 
do group = 1 to &r; 
  t&n=0; Xt=0; Sn=0; output; t&n=1/&n; output; 
  do tt = 1 to &n; t=tt; et = normal(9018275); Sn +et;
    Xt = Sn/sqrt(&n); 
    t&n   = t/&n;   Label t&n   = "0  <=   t/&n    <=  1";
  output; t=tt+1; t&n=t/&n; if t le &n then output; end; end; 
  proc print; run; 
proc gplot annotate=anno; 
plot Xt*t&n=group/nolegend vaxis=axis1;  
label Xt = "X(t,&n)";  
symbol1 v=none i=join r=&r; 
title "&r realizations of X(t,&n)"; 
run;