A relativistic harmonic confinement model for quarks and a similar current confinement model for gluons have been used to obtain an equation of state for quark-gluon plasma. Such models may be deduced from QCD under certain approximations, by considering small quantum fluctuations about a background field. At high temperatures a T7 dependence of pressure and energy density is obtained with relativistic harmonic mode of confinement.