One of the main difficulties in studying quantum field theory, in the perturbative regime, is the calculation of D-dimensional Feynman integrals. In general, one introduces the so-called Feynman parameters and, associated with them, the cumbersome parametric integrals. Solving these integrals beyond the one-loop level can be a difficult task. The negative-dimensional integration method (NDIM) is a technique whereby such a problem is dramatically reduced. We present the calculation of two-loop integrals in three different cases: scalar ones with three different masses, massless with arbitrary tensor rank, with and N insertions of a two-loop diagram.