#### Abstract

Göhmann, Klümper and Seel derived the multiple integral formula of the density matrix of the XXZ Heisenberg chain at finite temperatures. We have applied the high temperature expansion (HTE) method to isotropic case of their formula in a finite magnetic field and obtained coefficients for several short-range correlation functions. For example, we have succeeded to obtain the coefficients of the HTE of the 3rd neighbor correlation function < σ jσ z j+3 > for zero magnetic field up to the order of 25. These results expand our previous results on the emptiness formation probability [Z.Tsuboi, M.Shiroishi, J. Phys.A: Math. Gen. 38(2005) L363] to more general correlation functions. MSC: 82B23; 45G15; 82B20; 82B80 PACS: 75.10.Jm, 02.30.Ik, 05.70.-a, 05.30.-d