Square lattice Ising model with H = −J Σ s_i s_j. The critical temperature is T_c = 2J/ln(1+√2) ≈ 2.269J (Onsager 1944). Wolff cluster updates dramatically reduce critical slowing down near T_c.