  A method is introduced to calculate the UV-divergent parts at one-loop level
in dimensional regularization. The method is based on the recursion, and the
basic integrals are just the scaleless integrals after the recursive reduction,
which involve no other momentum scales except the loop momentum itself. The
method can be easily implemented in any symbolic computer language, and an
implementation in Mathematica is ready to use.
