Metropolis Monte Carlo simulation on a square lattice. H = −J Σ sisj − h Σ si, with J = 1 and periodic boundary conditions. Critical temperature Tc ≈ 2.269.