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
输出结果完全符合。