The program "sector_decomposition" can be used to compute numerically
the Laurent expansion of (divergent) multi-loop Feynman integrals.