It's not ``integrate.quad``, which still gives the right result and no warnings. The problem is that numpy functions are giving warnings for integrals that are converging fine, because it's not happy when you put in large numbers (something like integrating ``x**-2`` from zero to inf ).

Figuring these out is quite time-consuming, but I'll go figure out what happens exactly for the first warning now.