The total energy of the Heisenberg ferromagnetic system is calculated by the many-body Green's function method under random phase approximation when the correlation function 〈S+S-〉 is considered. A general expression of the total energy universally applicable to any spin quantum number S is obtained. Numerical results are computed. The energy is lower than the case when the correlation function is not included.