Algorithm for differential equations for Feynman integrals in general dimensions

We present an algorithm for determining the minimal order differential equations associated to a given Feynman integral in dimensional or analytic regularisation. The algorithm is an extension of the Griffiths-Dwork pole reduction adapted to the case of twisted differential forms. In dimensional regularisation, we demonstrate the applicability of this algorithm by explicitly providing the inhomogeneous differential equations for the multiloop two-point sunset integrals: up to 20 loops for the equal mass case, the generic mass case at two- and three-loop orders. Additionally, we derive the differential operators for various infrared-divergent two-loop graphs. In the analytic regularisation case, we apply our algorithm for deriving a system of partial differential equations for regulated Witten diagrams, which arise in the evaluation of cosmological correlators of conformally coupled $\phi^4$ theory in four-dimensional de Sitter space.