This paper proposes a numerical finite-element-based model aimed at optimizing narrow-vein stope stability. This model combines empirical and numerical methods to develop a sequence, which may determine an acceptable stope safety factor. A stope stability analysis was conducted through the Mathews stability graph method, which requires two factors: the hydraulic radius (HR) and stability number (N'). The Mathews stability graph method is used to assess the stability of an underground design. Variations in stope dimensions are estimated by changing the HR and Factor A within the N', which is determined through numerical methods. The results of the numerical simulation indicate that the HR increases with an increase in stope dimensions, while Factor A maintains an inverse relationship with the maximum stress induced on the excavation walls. This document demonstrates the potential of combining empirical and numerical methods in stope design optimization, especially when developed in small narrow vein mines.