蒙特卡洛方法应用

e的定义如下：

 

求解它有两种方法：
一是公式法，

取一个很大的数代入。

python代码如下:

x = 1024*1024
e = pow(1+1/x,x)
print("{}".format(e))

 

二是蒙特卡洛方法 (题目本身不难，重点的是这种思维模式)

 

import random as rd

scatterNum = 1024*1024 #撒点的总数
count = 0

for i in range(scatterNum):
    x = rd.uniform(1,2)
    y = rd.uniform(0,1)
    if x*y <= 1:
        count+=1
e = pow(2,scatterNum/count)
print("自然常数e：{}".format(e))