Pregunta

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.

¿Fue útil?

Solución

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.

Licenciado bajo: CC-BY-SA con atribución
No afiliado a StackOverflow
scroll top