利用python编写蒙特卡罗计算定积分

Question

计算定积分编写蒙特卡罗方法计算定积分S= ∫2 1 1/xdx的程序,并计算自然常数e的值

Answer 1

示例代码如下：
1.蒙特卡罗方法计算定积分代码

import random
def monte_carlo_integration(func, a, b, n):
    s = 0
    for i in range(n):
        x = random.uniform(a, b)
        s += func(x)
    return s / n * (b - a)
……

2.定义被积函数代码

……
def f(x):
    return 1 / x
……

3.计算定积分及计算自然常数e并输出的代码

……
S = monte_carlo_integration(f, 1, 2, 100000)
e = monte_carlo_integration(lambda x: 1 / x, 1, 100, 100000)
print("定积分S= ∫2 1 1/xdx ≈", S)
print("自然常数e ≈", e)

希望能帮到你，加油~~~