In this work, the thermal properties, such as total energy, pressure and free energy of an interacting electron gas system have been investigated. The ThomasFermi equation was numerically solved for a harmonically confined, two-dimensional gas system at finite temperatures. Numerical results were compared with analytical ones provided by the literature for a non-interacting case at finite temperatures and for a interacting case at T = 0 K. Excellent agreement was achieved for both cases. The results indicate that interactions affect the density of particles and hence, the density-dependent thermal functions significantly. © 2012 World Scientific Publishing Company.