  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.
