An efficient, finite temperature, recursion method (RM) is introduced for c
alculations of the mean-square relative displacements sigma(2) in multiple
scattering (MS) XAFS Debye-Waller factors. Instead of calculating total pro
jected densities of modes, the calculations are based on a double delta-fun
ction representation. Results for the Debye-Waller factors are found to be
in agreement with equation-of-motion (EM) method to within about 10% percen
t for all MS paths.