質問

I have been learning about QUADPACK and I used scipy.integrate.quad to calculate an integral from 0 to infinite. It gave a very good result, but now I want to know which integration method (QAGI, QAWF,etc.) the software has applied. Is there any way for printing the technique? Does the software apply some decision tree? Thanks in advance for your time. Regards.

役に立ちましたか?

解決

As scipy is open source, you can actually read the code for integrate.quad, which says:

For finite integration limits, the integration is performed using a Clenshaw-Curtis method which uses Chebyshev moments. ... If one of the integration limits is infinite, then a Fourier integral is computed (assuming w neq 0).

The call in _quad for an infinite bound is to _quadpack._qagie.

ライセンス: CC-BY-SA帰属
所属していません StackOverflow
scroll top