用随机数法求 Ln38 的近似值

如题，

用随机数法求 Ln38 的近似值

我寻思着直接用：
data a ;
y=log(38);
put y;
run;
这多简单啊~~

可是老师出的题目要求用随机数法求

言外之意就是要用到ranuni(seed)

请高手赐教！

由于ln38=ln38-ln1
因此，ln38是1/x在（1，38）区间的积分
故可以用随机模拟的方法给出如下程序：
m=10000
x=runif(m,min=1,max=38)
theta=(38-1)*mean(1/x)
theta
the=log(38)
the
> > [1] 3.642904
> > [1] 3.637586
输出结果完全符合。

tsclw 发表于 2010-11-24 19:11
由于ln38=ln38-ln1
因此，ln38是1/x在（1，38）区间的积分
故可以用随机模拟的方法给出如下程序：
m=10000
x=runif(m,min=1,max=38)
theta=(38-1)*mean(1/x)
theta
the=log(38)
the
> > [1] 3.642904
> > [1] 3.637586
输出结果完全符合。

一语惊醒梦中人，，，虽然你的程序看不懂，，，哈哈  ，，

多多交流啊

哈哈，不好意思，我用R编的，转化成SAS的就行了