Harmonic oscillator thermal density matrix: First-order differential equations for the position representation

A first-order partial differential equation is derived whose solution enables us to find straightforwardly the off-diagonal matrix elements in the position representation of the harmonic oscillator density operator. This approach constitutes an alternative to techniques that require advanced knowledge of mathematical and quantum mechanical results.