Fuchsia: a tool for reducing differential equations for Feynman master integrals to epsilon form

We present $text{Fuchsia}$ $-$ an implementation of the Lee algorithm, which for a given system of ordinary differential equations with rational coefficients $partial_x,mathbf{f}(x,epsilon) = mathbb{A}(x,epsilon),mathbf{f}(x,epsilon)$ finds a basis transformation $mathbb{T}(x,epsilon)$, i.e., $mathbf{f}(x,epsilon) = mathbb{T}(x,epsilon),mathbf{g}(x,epsilon)$, such that the system turns into the epsilon form: $partial_x, mathbf{g}(x,epsilon) = epsilon,mathbb{S}(x),mathbf{g}(x,epsilon)$, where $mathbb{S}(x)$ is a Fuchsian matrix. A system of this form can be trivially solved in terms of polylogarithms as a Laurent series in the dimensional regulator $epsilon$. That makes the construction of the transformation $mathbb{T}(x,epsilon)$ crucial for obtaining solutions of the initial equations. In principle, $text{Fuchsia}$ can deal with any regular systems, however its primary task is to reduce differential equations for Feynman master integrals. It ensures that solutions contain only regular singularities due to the properties of Feynman integrals.