  An algorithm for the reduction of massive Feynman integrals with any number
of loops and external momenta to a minimal set of basic integrals is proposed.
The method is based on the new algorithm for evaluating tensor integrals,
representation of generalized recurrence relations for a given kind of
integrals as a linear system of PDEs and the reduction of this system to a
standard form using algorithms proposed by G. Reid. Basic integrals reveal as
parametric derivatives of the system in the standard form and the number of
basic integrals in the minimal set is determined by the dimension of the
solution space of the system of PDEs.
