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.