如何解决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
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。