COMPUTATION OF COEFFICIENT FUNCTIONS FOR AN EXPANSION OF PARAPOSITRONIUM DECAY WIDTH USING REDUCE

A kind of asymptotic expansion for calculation of the (para)positronium decay rate is formulated. The method is completely relativistic and gauge invariant. The expansion goes over the inverse square of relative momentum $k$ of two photons in Euclidean region and the actual width is determined at $k^2=-m^2$, $m$ being the electron mass.