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