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

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

示例代码如下：
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)

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

不知道你这个问题是否已经解决, 如果还没有解决的话:

• 你可以看下这个问题的回答https://ask.csdn.net/questions/7577966

如果你已经解决了该问题, 非常希望你能够分享一下解决方案, 写成博客, 将相关链接放在评论区, 以帮助更多的人 ^-^