  We show how to use the multiple histogram method to combine canonical
ensemble Monte Carlo simulations made at different temperatures and densities.
The method can be applied to study systems of particles with arbitrary
interaction potential and to compute the thermodynamic properties over a range
of temperatures and densities. The calculation of the Helmholtz free energy
relative to some thermodynamic reference state enables us to study phase
coexistence properties. We test the method on the Lennard-Jones fluids for
which many results are available.
