In this work, we present a comparative study of thermodynamic quantum equilibrium observables of spin-1 Heisenberg chain. As a theoretical approach, the modified spin-wave theory is chosen. From the numerical side, we consider finite-temperature Lanczos, kernel polynomial, and density matrix renormalization techniques. The results show general consistency of thermodynamic quantities, such as heat capacity and magnetic susceptibility, calculated within the approaches. For the modified spin-wave, the results show some inaccuracy especially at low temperature, failure in capturing the gapped features of spin-1 Heisenberg chain, while the numerical approaches illustrate a good achievement at $$T\rightarrow 0$$
T
→
0
within the proper parameters chosen.