您好,登錄后才能下訂單哦!
代碼閑暇時間,將代碼過程重要的一些代碼收藏起來,下面代碼是關于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進行舉報,并提供相關證據,一經查實,將立刻刪除涉嫌侵權內容。