We compute the free energy density for pure non-Abelian gauge theory at high temperature and zero chemical potential. The three-loop result to o (~~) is where T is the temperature, C is the Riemann zeta function, g~ g(,ii)~i'2, ji is the MS renormalization scale, g(p) is the corresponding coupling constant, and d~ and CA are the dimension and Casimir number of the adjoint representation. We examine the sensitivity of this result to the choice of renormalization scale ,ii. We also give a result for the free energy of scalar d4 theory, correcting a result previously given in the literature. PACS number(s): ll.lO.Wx, 12.38.B~ PETER ARNOLD AND CHENGXING ZHAI 50 -