  We propose a new Monte Carlo algorithm for the free energy calculation based
on configuration space sampling. We implement this algorithm for Ising model.
Comparison with the exact free energy shows an excellent agreement. We analyse
the properties of this algorithm and compares it with Wang-Landau algorithm
which samples in energy space. This method is applicable to classical
statistical models. The possibility of extending it to quantum systems is
discussed.
