A modernized version of Gibbs' use of the grand canonical ensemble

In the last chapter of his Introduction to statistical mechanics Gibbs introduces the idea of the grand canonical ensemble. He had previously determined the properties of an assembly or a phase containing a given number of systems by averaging the properties of the assembly over an ensemble of examples canonically distributed in phase, keeping the number of systems in the assembly fixed. This means of course constructing what we now call the partition function for the assembly by summing or integrating e−E/kT over the whole of the accessible phase space.