intalg.spad line 50 [edit on github]
algebraic Hermite redution.
HermiteIntegrate(f, ') returns [g, h] such that f = g' + h and h has a only simple finite normal poles.
HermiteIntegrate(f, ', d0) returns [g, h] such that f = g' + d0*g + h and h has a only simple finite normal poles. Note: d0 must be integral.