import random
def f(x):
# 被积函数
return x**2
def Monte_Carlo_Methond(a, b):
# (a,b)是积分区间
sum = 0
y1 = f(a)
y2 = f(b)
# y1,y2是函数的上下界
for i in range(1000000):
x = random.uniform(a, b)
y = random.uniform(0, y2)
# (x,y)是在矩形内随机生成的点
if y <= f(x):
sum += 1
return (sum / 1000000) * (abs(b - a)) * (abs(max(y1, y2)))
print(Monte_Carlo_Methond(0, 1))
[Pytnon]用蒙塔卡罗面积法对x^2求定积分
最新推荐文章于 2026-08-17 11:10:18 发布
本文介绍了如何使用Python实现蒙特卡洛方法来估算函数∫₀¹ x² dx,通过100万次随机采样计算概率区域,进而得到积分近似值。

488

被折叠的 条评论
为什么被折叠?



