Montecarlo Simulation of Ising Model

A simple simulation for 1,2 and 3 dimensional Ising model.

Montecarlo simulations show the phase transition from a disordered state to an ordered one in 2D and 3D while in the 1D case no phase transition is allowed for temperature different from zero.
The parameters T ( temperature in kb units) , J ( exchange constant ) and h ( external field ) can be adjusted freely.
Here J > 0 gives ferromagnetic order while J < 0 gives antiferromagnetic order.