Purpose. The aim of this paper is numerical modeling of the slope stability of limestone quarry in dry and wet cases. Methodology. The study of slope stability was performed by the limit equilibrium method (LEM) through the SLIDE software, and numerical methods through PLAXIS and FLAC. Findings. In this study, the presence of water in the marl layer has a heavy influence on the safety factor (SF), which allows for instability of the slope. Originality. The originality of this work is based on two different approaches to determine the safety factor in the dry and wet cases through three methods, the limit equilibrium, finite element and finite difference methods (FEM, FDM). Practical value. This study has shown that numerical modeling of the slope stability using three methods, LEM, FDM and FEM can be expected to produce similar results. The different methods indicate that the finite difference method in our case is the most efficient and can be used as a practical check for the other methods.