Ferromagnetic Ising model. Tc = 2/ln(1+√2) ≈ 2.269. Watch the phase transition. Toggle Wolff cluster algorithm.