A 2D Ising model of ferromagnetism using the Metropolis algorthim