The main objective of this project is to study phase transitions in Ising models on fractal lattices with different Hausdorff dimensions (d_H ≈ 1.892 and  d_H ≈ 1.595). The Ising model is presented for the cases of one-dimensional and two-dimensional lattices, and the computational tools used for simulation are explained. Subsequently, a numerical analysis is presented for the fractal case, utilizing Monte Carlo techniques and the Metropolis algorithm, to obtain a critical temperature value and a correlation length function. Analytical methods are also developed to provide a comparison with the numerically obtained results.