我想在三角形上对函数进行数值积分,类似于
import scipy.integrate as integrate
inside = lambda x: integrate.quad(lambda x,y: 1, 0, x, args=(x))[0]
outside = integrate.quad(inside, 0, 1)[0]
print(outside)
0.5
但是使用scipy.integrate.fixed_quad
函数(其积分阶数n
作为参数)。然而,当我写
inside = lambda x: integrate.fixed_quad(lambda x,y: 1, 0, x, args=(x), n=5)
print(inside(5))
Traceback (most recent call last): File "", line 1, in File "", line 1, in File "/Users/username/anaconda/lib/python3.5/site-packages/scipy/integrate/ quadrature.py", line 82, in fixed_quad return (b-a)/2.0 * np.sum(w*func(y, *args), axis=0), None
TypeError: () argument after * must be an iterable, not int
我不知道我做错了什么,因为我正在遵循 scipy.integrate.fixed_quad 上的文档.
最佳答案
问题在于您对 args
、args=(x)
的定义。它应该作为元组传递,因此您需要添加一个额外的逗号以使其成为元组:
inside = lambda x: integrate.fixed_quad(lambda x,y: 1, 0, x, args=(x,), n=5)
然后
inside(5)
产量
(5.0, None)
线路
integrate.quad(lambda x,y: 1, 0, x, args=(x))[0]
与quad
一样工作,检查args
是否是一个元组;如果没有,则进行转换(直接从源代码中获取):
if not isinstance(args, tuple):
args = (args,)
在 fixed_quad
中,情况并非如此,这就是您在一种情况而非两种情况下收到错误的原因。
关于python - scipy.integrate.fixed_quad 可以计算具有函数边界的积分吗?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/42270420/