python求定积分_python使用蒙特卡洛积分法（Monte-Carlo integration）求定积分

最新推荐文章于 2024-08-24 22:00:01 发布

weixin_39933484

最新推荐文章于 2024-08-24 22:00:01 发布

阅读量2.5k

点赞数 3

文章标签： python求定积分

本文链接：https://blog.csdn.net/weixin_39933484/article/details/113452541

版权

### 使用蒙特卡洛投点法计算定积分
import random
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt 
count=0
p=0
for i in range(10000):
    count+=1
    point_x=random.random()
    point_y=random.random()
    value=point_x/25+1/5
    if point_y<=value:
        p+=1
print('p is :',p/count)

p is : 0.2204

### 使用平均法求定积分
import numpy as np
point_list=np.linspace(0,1,10001)
value=0
for i in range(10000):
    x0=point_list[i]
    x1=point_list[i+1]
    v=0.5*(x1-x0)*(x0/25+0.2+x1/25+0.2)
    value+=v
print('p is :',value)

p is : 0.2200000000000001