Spins on a square lattice flip via the Metropolis rule: accept lower energy always, higher energy with probability e^(−ΔE/kT). At T_c ≈ 2.269 J/k, the system undergoes a continuous phase transition with diverging correlation length.