  Two-loop vertex Feynman diagrams with infrared and collinear divergences are
investigated by two independent methods. On the one hand, a method of
calculating Feynman diagrams from their small momentum expansion extended to
diagrams with zero mass thresholds is applied. On the other hand, a numerical
method based on a two-fold integral representation is used. The application of
the latter method is possible by using lightcone coordinates in the parallel
space. The numerical data obtained with the two methods are in impressive
agreement.
