2D Ising model, Metropolis updating.
$$ S=K\sum_{\langle i,j\rangle} s_i s_j $$ The critical coupling \(K_c \approx 0.44 \)