高斯-勒让德积分(Gauss-Legendre integral)的原理及Python实现

Astrologer__

已于 2022-12-22 13:41:03 修改

阅读量5.5k

点赞数 6

分类专栏： python 笔记文章标签： python

于 2021-04-26 22:52:11 首次发布

本文链接：https://blog.csdn.net/qq_19905867/article/details/116176450

版权

该博客介绍了如何利用Python实现高斯Legendre公式进行数值积分。通过导入`scipy.special.roots_legendre`计算节点和权重，然后应用这些节点和权重求解函数在[0, 1]区间上的积分。示例代码展示了计算π的近似值，结果为3.1415926535897927，证明了算法的准确性。文章还引用了多篇参考资料，包括书籍、在线计算器和数学世界网站。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

算法

#%%
from scipy.special import roots_legendre
import numpy as np

#%%
f = lambda x: 4./(1+x*x) #被积函数

N = 100 #取100个样本点
a, b = 0, 1 #积分上下限

#%%
x, wi = roots_legendre(N)
xi =</