代码闲暇时间，将代码过程重要的一些代码收藏起来，下面代码是关于Gauss-Legendre integration in Python的代码，希望对小伙伴也有用处。
''' I = gaussQuad2(f,xc,yc,m).
Gauss-Legendre integration of f(x,y) over a
quadrilateral using integration order m.
{xc},{yc} are the corner coordinates of the quadrilateral.
'''
from numpy import zeros,dot
def gaussQuad2(f,x,y,m):
def jac(x,y,s,t):
J = zeros((2,2))
def map(x,y,s,t):
N = zeros(4)
xCoord = dot(N,x)
yCoord = dot(N,y)
return xCoord,yCoord
s,A = gaussNodes(m)
sum = 0.0
for i in range(m):
for j in range(m):
xCoord,yCoord = map(x,y,s[i],s[j])
return sum免责声明：本站发布的内容（图片、视频和文字）以原创、转载和分享为主，文章观点不代表本网站立场，如果涉及侵权请联系站长邮箱：is@yisu.com进行举报，并提供相关证据，一经查实，将立刻删除涉嫌侵权内容。
