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