In this study, a polygonal approach is suggested to generalize the notion of the con…dence region of the univariate probability density function for the bivariate probability density function. The equal density approach is used to demonstrate that con…dence regions can be polygonal shapes. The bisection method is the preferred method in …nding the equal density value that reveals the desired con…dence coe¢ cient. Con…dence regions estimate not only bivariate unimodal probability functions but also bivariate multimodal probability functions. An approach is enhanced to estimate these con…dence regions for probability density functions which are de…ned as rectangular, polygonal and in…nite expanse areas. In order to show the applicable of the proposed method, four di¤erent examples are analyzed. The results show that the con…dence region is found no matter how complex the distribution function. In addition, the proposed method gives more e¢ cient results for multimodal probability density functions.