  We have provided a concise introduction to the Ising model as one of the most
important models in statistical mechanics and in studying the phenomenon of
phase transition. The required theoretical background and derivation of the
Hamiltonian of the model have also been presented. We finally have discussed
the computational method and details to numerically solve the two- and
three-dimensional Ising problems using Monte Carlo simulations. The related
computer codes in both Python and Fortran, as well as a simulation trick to
visualize the spin lattice, have also been provided.
