Octave e的倍频程近似

Octave e的倍频程近似,octave,mcmc,Octave,Mcmc,我想用MCMC算法在倍频程中以最大精度计算以下表达式:“1/e”。在阅读了一些教程之后,我找到了一个计算π的公式,但我不明白它是如何工作的 octave:2> S=1e7; a=rand(S,2); 4*mean(sum(a.*a,2)<1) ans = 3.1418 octave:2>S=1e7;a=兰特(S,2);4*mean(sum(a.*a,2)这是估算pi的一个应用程序。本质上,您创建的是一个Sx2矩阵(将其视为(x,y)坐标的S个数),所有值都在0和1之间,几何上在1

我想用MCMC算法在倍频程中以最大精度计算以下表达式:“1/e”。在阅读了一些教程之后,我找到了一个计算π的公式,但我不明白它是如何工作的

octave:2> S=1e7; a=rand(S,2); 4*mean(sum(a.*a,2)<1)
ans = 3.1418 
octave:2>S=1e7;a=兰特(S,2);4*mean(sum(a.*a,2)这是估算pi的一个应用程序。本质上,您创建的是一个Sx2矩阵(将其视为(x,y)坐标的S个数),所有值都在0和1之间,几何上在1x1平方之内。然后,您将x和y值平方,并将它们相加,以获得点到原点的距离平方。