We derive an algorithm for automatic calculation of perturbative p-fun
ctions and anomalous dimensions in any local quantum field theory with
canonical kinetic terms. The infrared rearrangement is performed by i
ntroducing a common mass parameter in all the propagator denominators.
We provide a set of explicit formulae for all the necessary scalar in
tegrals up to three loops. (C) 1998 Published by Elsevier Science B.V.