Ising Model — Metropolis Algorithm

2D ferromagnetic Ising model on a square lattice. At the critical temperature T_c = 2/ln(1+√2) ≈ 2.269 (Onsager 1944), the system undergoes a second-order phase transition.

Magnetization: — | Energy: —