python - 如何使用 Python 的 scipy.integrate.quad 评估多元函数的单积分？

Question

我正尝试使用 scipy.integrate.quad 将一个函数集成到 Python 中。这个特定的函数有两个参数。我只想整合一个论点。示例如下所示。

from scipy import integrate as integrate
def f(x,a):  #a is a parameter, x is the variable I want to integrate over
    return a*x

result = integrate.quad(f,0,1)

这个例子不起作用(你可能很清楚)，因为当我尝试它时 Python 提醒我:

TypeError: f() takes exactly 2 arguments (1 given)

我想知道当给定的函数通常是一个多变量函数时，如何使用 integrate.quad() 对单个变量进行积分，额外的变量提供参数函数。

Answer 1

在 the scipy documentation 中找到了答案.

您可以执行以下操作:

from scipy import integrate as integrate
def f(x,a):  #a is a parameter, x is the variable I want to integrate over
    return a*x

result = integrate.quad(f,0,1,args=(1,))

quad 方法中的 args=(1,) 参数将使 a=1 用于积分计算。

这也可以传递给具有两个以上变量的函数:

from scipy import integrate as integrate
def f(x,a,b,c):  #a is a parameter, x is the variable I want to integrate over
    return a*x + b + c

result = integrate.quad(f,0,1,args=(1,2,3))

这将使 a=1, b=2, c=3 用于积分评估。

对于要以这种方式集成的函数，重要的是要记住将要集成的变量置于函数的第一个参数上。