  An algorithm for calculation of three-loop propagator diagrams in HQET, based
on integration-by-parts recurrence relations, is constructed and implemented as
a REDUCE package Grinder, and in Axiom.
