  A new approach is presented to evaluate multi-loop integrals, which appear in
the calculation of cross-sections in high-energy physics. It relies on a fully
numerical method and is applicable to a wide class of integrals with various
mass configurations. As an example, the computation of two-loop planar and
non-planar box diagrams is shown. The results are confirmed by comparisons with
other techniques, including the reduction method, and by a consistency check
using the dispersion relation.
