I'm using scipy.integrate.dblquad for a 2-dimensional integral. I'd like to know what underlying method it is using to calculate this. I can't see this in the docs. AFAICT QUADPACK is just for one dimensional integrals so it must be using something else. Thanks, John.