We derive the partition function of the one-body and two-body systems of classical noncommutative harmonic oscillator in two dimensions. Then, we employ the path integral approach to the quantum noncommutative harmonic oscillator and derive the partition function of the both systems at finite temperature.
