We present a functional-integral method that is suitable for obtaining nonperturbative solutions to the partition function of fermion systems. The density matrix is cast as an expansion about the classical path contributions, and the quantum corrections can be calculated systematically by using retarded propagators in the space holomorphic functions. It is shown that the semiclassical expansion amounts to a high-temperature expansion, where the first nonvanishing quantum correction is of O(${\mathrm{\ensuremath{\beta}}}^{2}$). Each term of the semiclassical expansion may contain the coupling parameter in all orders of perturbation theory. The method can also be extended to boson systems.