RunDec: a Mathematica package for running and decoupling of the strong   coupling and quark masses
  In this paper the formulae are collected which are needed for the computation of the strong coupling constant and quark masses at different energy scales and for different number of active flavours. All equations contain the state-of-the-art QCD corrections up to three- and sometimes even four-loop order. For the practical implementation {\tt Mathematica} is used and a package containing useful procedures is provided. 
K.G. Chetyrkin, J.H. Kuehn and M. Steinhauser