Abstract. We calculate the pressure and density of polarized non-relativistic systems of two-component fermions coupled via a contact interaction at finite temperature. For the unpolarized one-dimensional system with an attractive interaction, we perform a thirdorder lattice perturbation theory calculation and assess its convergence by comparing with hybrid Monte Carlo. In that regime, we also demonstrate agreement with real Langevin. For the repulsive unpolarized one-dimensional system, where there is a so-called complex phase problem, we present lattice perturbation theory as well as complex Langevin calculations. For our studies, we employ a Hubbard-Stratonovich transformation to decouple the interaction and automate the application of Wick's theorem for perturbative calculations, which generates the diagrammatic expansion at any order. We find excellent agreement between the results from our perturbative calculations and stochastic studies in the weakly interacting regime. In addition, we show predictions for the strong coupling regime as well as for the polarized one-dimensional system. Finally, we show a first estimate for the equation of state in three dimensions where we focus on the polarized unitary Fermi gas.