python - Scipy nquad 给出错误的输出

Question

我正在尝试编写代码来解决

从1到2的积分从0到x的积分，xy^2 dydx = 31/15

我还使用 Wolfram|Alpha 验证了手动解决方案。但代码给出的输出为 3.1。

from scipy.integrate import nquad

def f(x, y):
    return x * y**2
def limits_y(x):
    return [0, x]
def limits_x():
    return [1, 2]

ans,err = nquad(f, [limits_y, limits_x])
print(round(ans,3))

知道为什么吗？我错过了什么吗？

Answer 1

nquad期望被积函数的第一个参数是与最内层积分对应的变量。如果您将 f 的定义更改为

def f(y, x):
    return x * y**2

然后nquad返回预期结果:

In [73]: def f(y, x):
    ...:     return x*y**2 
    ...:

In [74]: ans, err = nquad(f, [limits_y, limits_x])

In [75]: ans
Out[75]: 2.066666666666667