An Algorithm for the Symbolic Reduction of Multi-loop Feynman Integrals via Generating Functions

We develop a generating-function formulation for the symbolic reduction of multi-loop Feynman integrals. In this framework, integration-by-parts identities are rewritten as differential equations for sector-wise generating functions, so the reduction problem can be studied in a non-commutative algebra of differential operators rather than only through relations among individual integrals. This viewpoint leads to an iterative algorithm that generates candidate equations, extracts symbolic reduction rules, updates the active rule set, and tests completeness on the lattice of integral indices. We illustrate the method with the sunset topology, planar and non-planar massless double-box topologies, representative subsectors, and a degenerate example in which the top sector contains no master integral. Together, these examples show how symbolic reduction rules, descendant equations, and completeness criteria can be organized within a single algebraic framework.