In this paper we present numerical analysis of the phase transition of the area-interaction model, which is a standard model of Statistical Mechanics. The theoretical results are based on a recent paper by Dereudre \& Houdebert \cite{dereudre_houdebert_2018_SharpTransitionWR} which provides a complete phase diagram except on a bounded (implicit) domain. With our numerical analysis we give an approximative explicit description of this domain. %This region is related to the theoretical function $β\mapsto \widetilde{z}_c^a(β, 1)$ which is percolation threshold of the model, for given $β$. %We provide a numerical approximation of this function in order to fully understand the phase diagram. Furthermore our numerical results confirm the still unproven conjecture stating that non-uniqueness holds if and only if $z= β$ is large enough, with a value of the threshold obtained from the simulation of $β_c \simeq 1.726$.