calculations. We propose and implement a new method to efficiently evaluate such integrals
in the physical region through the numerical integration of a suitable set of differential
equations, where the initial conditions are provided in the unphysical region via the sector
decomposition method. We present numerical results for a set of two-loop integrals, where
the non-planar ones complete the master integrals for gg→ γγ and\(q\overline {q}\)→ γγ …