Ferromagnetic Ising model on a square lattice. Metropolis algorithm samples the Boltzmann distribution. Critical temperature T_c = 2/ln(1+√2) ≈ 2.269 (Onsager 1944).