Use scipy.integrate.quad
Use scipy.integrate.quad

When df is given use scipy.integrate.quad (which is a wrapper for the Fortran library QUADPACK) to compute the contour integrals that are not integers.

rparini/cxroots,rparini/cxroots