An implementation of the Metropolis and Heat Bath algorithms for the 2D Ising model - View it on GitHub