  In this work we try to use the Monte Carlo algorithm, metropolis, to study
the behavior of 2D magnetic systems; honeycomb, hexagonal and square lattices.
In this study we use Ising Model of magnetism, with considering only nearest
neighbors in the Hamiltonian and with the Ferromagnetic correlation constant.
The code was written in Payton. In this calculation, by comparing the results
obtained for all lattices, we will show that the paramagnetic to ferromagnetic
phase transition depends on the number of neighbor and it is a direct
dependence.
