The aim of this article is to study in a rigorous way at the computational level the critical
properties of 2D Ising Model. We have used the Metropolis algorithm for a square lattice with
periodic boundary conditions in a ferromagnet without external field. We will also address the
effect of an external magnetic field, and lattices with triangular and hexagonal symmetry.
2015/2016