  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.
