Calculation of Infrared-Divergent Feynman Diagrams with Zero Mass Threshold

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.