To evaluate ESR spectrum at finite temperatures for specified spatial configurations of spins is very important issue to study quantum spin systems. Although a direct numerical estimation of the Kubo formula provides exact data, the application is limited to small size of the system because of the restriction of the computer capacity. The method of the Fourier transform of the autocorrelation function improved the restriction. As an extension of the method, we propose a new method for numerical calculation of the ESR spectrum from the time evolution of the magnetization by making use of the Wiener-Khinchin theorem.