Using a leading-order semiclassical approximation, we calculate the third-and fourth-order virial coefficients of nonrelativistic spin-1/2 fermions in a harmonic trapping potential in arbitrary spatial dimensions, and as functions of temperature, trapping frequency and coupling strength. Our simple, analytic results for the interaction-induced changes ∆b3 and ∆b4 agree qualitatively, and in some regimes quantitatively, with previous numerical calculations for the unitary limit of threedimensional Fermi gases. arXiv:1908.00070v1 [cond-mat.quant-gas]