Fuchsia and master integrals for splitting functions from differential equations in QCD

We report on the recent progress in reducing differential equations for Feynman master integrals to canonical form with the help of a method proposed by Roman Lee. For the first time, we present Fuchsia --- our open-source implementation of the Lee algorithm written in Python using mathematical routines of a free computer algebra system SageMath. We demonstrate Fuchsia by reducing differential equations for NLO contributions to splitting functions in QCD, which contain both loops and legs integrals.