From Sven Kreiss
This is the result of my 2D Ising Monte Carlo simulation. The code is written in Python in less than 70 lines and all videos shown here were simulated on a laptop in less than two minutes. The last image is a bit hard to see, but it is a plot of energy versus iteration number.
First with spin +1/2 and -1/2 below and above critical temperature.
This simulation can simulate arbitrary spins. Here, I chose the range -10 to 10.
As a check, one can use this algorithm to confirm the spin +1/2 and spin -1/2 case from above.