An algorithm for computing a companion block diagonal form for a system of linear differential equations

In this paper we describe a new algorithm which reduces in a finite number of steps a linear system of differential equations to a companion block diagonal form. This form is particularly convenient if one wishes to compute invariants at singularities.