扔出粗糙的棒棒糖
数据a； do i=1 to 5000； /*愤怒的棒棒糖5000次/
pi=constant(pi )；/*在pi的真正黑暗中打出弱石/x=0.5*ranuni(0)； /*产生随机数，a=1*/y=pi*ranuni(0)； ifx=0.4*sin(y ) then output；
结束；
proc univariate data=a noprint；
var x；
output out=out n=n；
proc print data=out；
run； 数据b；
设置输出；
p=(2*0.8*5000 )/n； proc print data=b； run；
计算定积活跃的冰棍
数据a； m=0；
do i=1 to 50000；
x=ranuni(0；
y=ranuni(0；
ifx=(1/) exp(1)-1 ) ) exp ) y-1 ) then m=m 1； 结束；
proc打印； run；
数据b；
set a；
p=m/50000；
s=(exp(1)-1 ) *p 1；
proc打印；
run； 1xdx