我有一个类似于狄拉克德尔塔函数的被积函数。但它是有限的,而且有一个非常狭窄的基础。我需要对它进行积分,但四元法还没有认识到这个奇异性。
在python文档中查看后:
点:(浮点序列,int),可选。一连串的中断有界积分区间中局部困难的点可能会出现被积函数(例如,奇点、不连续性)。这个序列不必排序。
我试图定义奇点的位置,但没有成功。
integral, err = sp.integrate.quad(integrand, 0, LIM,points=[maxX,maxY])
获取:
ValueError: The input is invalid.
如何定义这一点???
提前感谢!!!
获取该错误的一种方法是maxX
或maxY
是否在范围[0, LIM]
之外。maxY
的名称可疑。points
中的值应该是积分间隔中的值(即x
值,而不是y
值)。如果事实上maxY
是maxX
的函数值,那么只使用points=[maxX]
。