估算定积分python代码

以下是使用蒙特卡洛方法估算定积分的Python代码示例：

import random

def estimate_integral(f, a, b, n):
    total = 0
    for _ in range(n):
        x = random.uniform(a, b)
        total += f(x)
    average = total / n
    integral = (b - a) * average
    return integral

# 示例函数：估算 y = x^2 在 0 到 1 上的定积分
def f(x):
    return x**2

a = 0  # 积分下限
b = 1  # 积分上限
n = 100000  # 投点数量

estimated_integral = estimate_integral(f, a, b, n)
print("Estimated integral:", estimated_integral)

这段代码使用蒙特卡洛方法估算了函数 y = x^2 在区间 [0, 1] 上的定积分。通过生成随机点并计算函数值的平均值，然后乘以区间长度，得到了定积分的估计值。