For any d-dimensional self-interacting fermionic model, all coefficients in the high-temperature expansion of its grand canonical partition function can be put in terms of multivariable Grassmann integrals. A new approach to calculate such coefficients, based on direct exploitation of the Grassmannian nature of fermionic operators, is presented. We apply the method to the soluble Hatsugai-Kohmoto model, reobtaining well-known results. (C) 1999 American Institute of Physics. [S0022-2488(99)00310-2].