Scipy nquad 给出了错误的输出

问题描述

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

从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))

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

解决方法

nquad 期望被积函数的 first 参数是对应于最内积分的变量。如果您将 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

python scipy scipy