%macro ml(a,b,n);
data a;
do i=1 to &n;
u=uniform(0);
x=&a+&b*tan(cos(-1)*(u-0.5));
output;
end;%mend ml(a,b,n);
%ml(2,4,100);run;
这是SAS程序,a,b是该分布中的参数,n是产生随机数的个数,都由你自己指定.

2# harlon1976
谢谢！能不能用SAS画出柯西分布的概率图呢？