  Trager's Hermite reduction solves the integration problem for algebraic
functions via integral bases. A generalization of this algorithm to D-finite
functions has so far been limited to the Fuchsian case. In the present paper,
we remove this restriction and propose a reduction algorithm based on integral
bases that is applicable to arbitrary D-finite functions.
