  A new subtraction procedure for removal both ultraviolet and infrared
divergences in Feynman integrals is proposed. This method is developed for
computation of QED corrections to the electron anomalous magnetic moment. The
procedure is formulated in the form of a forest formula with linear operators
that are applied to Feynman amplitudes of UV-divergent subgraphs. The
contribution of each Feynman graph that contains propagators of electrons and
photons is represented as a finite Feynman-parametric integral. Application of
the developed method to the calculation of 2-loop and 3-loop contributions is
described.
