  The partition function by means of the static path approximation (SPA) plus
the random-phase approximation (RPA) treatment can be written as a contour
integral form without solving the RPA equations for a separable interaction.
This method is an efficient way to evaluate numerically the partition function
for realistic calculations. As an illustration, we adopt the pairing model at
finite temperature.
