  An algorithm for the systematic analytical approximation of multi-scale
Feynman integrals is presented. The algorithm produces algebraic expressions as
functions of the kinematical parameters and mass scales appearing in the
Feynman integrals, allowing for fast numerical evaluation. The results are
valid in all kinematical regions, both above and below thresholds, up to in
principle arbitrary orders in the dimensional regulator. The scope of the
algorithm is demonstrated by presenting results for selected two-loop
three-point and four-point integrals with an internal mass scale that appear in
the two-loop amplitudes for Higgs+jet production.
