如何定义integ.quad方法的不连续点序列

时间:2014-10-25 10:03:54

标签: python scipy numerical-integration quad

我有一个类似于狄拉克三角洲函数的被积函数。但它是有限的,基础非常狭窄。 我需要整合它,但quad方法还没有认识到这种奇点。

查看python文档后:

  

积分 :(浮点数序列,整数),可选。一系列的休息   有界积分区间内的局部困难点   可能发生被积函数(例如,奇点,不连​​续性)。该   序列不必排序。

我试图定义奇点所在的点,但没有成功。

integral, err = sp.integrate.quad(integrand, 0, LIM,points=[maxX,maxY])

获得:

ValueError: The input is invalid.

如何定义这一点???

提前致谢!!!

1 个答案:

答案 0 :(得分:0)

获取该错误的一种方法是maxXmaxY超出范围[0, LIM]maxY的名称是可疑的。 points中的值应为积分间隔中的值(即x值,而不是y值。如果实际上maxYmaxX处的函数值,那么只使用points=[maxX]