Reducing differential equations for multiloop master integrals

We present an algorithm of the reduction of the differential equations for master integrals the Fuchsian form with the right-hand side matrix linearly depending on dimensional regularization parameter $\epsilon$. We consider linear transformations of the functions column which are rational in the variable and in $\epsilon$. Apart from some degenerate cases described below, the algorithm allows one to obtain the required transformation or to ascertain irreducibility to the form required. Degenerate cases are quite anticipated and likely to correspond to irreducible systems.