我有点困在一个函数上,我试图通过scipy,python进行数值积分。在
为了简单起见,我将函数定义为:integral f(x,y)= SUM[double integral(ax+by)dxdy]
a和b是常数,但对于每个积分的方程,它们是不同的。我已经将每个函数分别积分,然后将所有积分的结果求和,但是这需要花费大量的时间来计算,这对于我正在尝试实现的目标并不理想。在
是否有一种方法可以通过扩展和来同时集成整个函数:
^{pr2}$
然后将带有(a,b)元组列表的函数传递给scipy的dblquad函数?在
目前,我正在努力在文献中找到任何与此相关的东西。在
*编辑
我提供了一个示例代码,以更清楚地显示我想要实现的目标:import sys
import re
import math
from scipy.integrate import dblquad
def f((x,y),variables):
V=0
for v in variables:
a,b=v
V=V+ax+by
return (V)
def integral(x_max,y_max,variables):
return dblquad(f, 0, y_max, lambda x: 0, lambda x: x_max,args=variables)
def main():
variables=[(1,2),(3,4),(5,6)] #example variables. The length