Computer Simulation
The equilibrium properties are obtained from the partition function
Z=? e-?E(i) where E(i) is the total energy of configuration i, ?=1/kT and the sum is over all 2N configurations
the sum over all configurations is like a 2N dimensional integral
the Metropolis algorithm is used to simulate the problem numerically and is a cornerstone of computational physics