In this talk, we explore the renormalization of the Ising model on hierarchical lattices, which leads to the recurrence relation of the corresponding partition functions. We will see that the singularity points of free energies lie on the Julia set of the renormalization transformation, which correspond to the points of phase transitions.