  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.
