我想在Octave中使用MCMC算法以最大精度计算以下表达式：“1 / e”。在阅读了一些教程后，我找到了计算π的公式，但我不明白它是如何工作的。
octave:2> S=1e7; a=rand(S,2); 4*mean(sum(a.*a,2)<1)
ans = 3.1418
有人可以解释并帮助我提示，如何使用这样的东西来计算'e'的值？
提前致谢。
答案 0 :(得分：2)
这是the dartboard method用于估算pi的应用。基本上你是在创建一个Sx2矩阵（把它想象成S个（x，y）坐标），所有这些矩阵的值都在0到1之间，所以几何上在1x1平方内。然后，您将对x和y值进行平方并将它们相加以获得距离原点的点的距离平方。 <1将所有这些距离转换为0或1，具体取决于该点是否位于以原点为中心的半径1的四分之一圆内。这个二进制数组的平均值是＆＃34; darts＆＃34;的比率。在总抛出的四分之一圈内击中，这是其面积的近似值。乘以4，您可以估算半径为1的整圆，其精确面积等于pi。
进行谷歌搜索会以类似的方式显示这个（希望）有用的出版物来计算e：Monte Carlo estimations of e