data one; array a a1-a3; array b b1-b3; array p p1-p3; array x x1-x3; array y y1-y3; do over a; a=exp(rannor(0)*.5); b=rannor(0); end; do j=1 to 10; t=rannor(0); do over p; y=0; p=1/(1+exp(-1.7*a*(t-b))); x=ranuni(0); if x le p then y = 1; else y=0; end; output; end; file 'c:\SASex\out1'; do over a; put a b; end; data two; set one; file 'C:\SASex\out2'; put '5'(y1-y3) (f1.0); run;