Computes the integral with Gauss-Legendre quadrature using m nodes in Python,computesquadrature,## module ga
文章由Byrx.net分享于2019-03-23 10:03:45
Computes the integral with Gauss-Legendre quadrature using m nodes in Python,computesquadrature,## module ga
## module gaussQuad''' I = gaussQuad(f,a,b,m). Computes the integral of f(x) from x = a to b with Gauss-Legendre quadrature using m nodes.'''from gaussNodes import *def gaussQuad(f,a,b,m): c1 = (b + a)/2.0 c2 = (b - a)/2.0 x,A = gaussNodes(m) sum = 0.0 for i in range(len(x)): sum = sum + A[i]*f(c1 + c2*x[i]) return c2*sum
评论关闭