Abstract: The stability analysis of slope and foundation is one of the important contents in geotechnical engineering. The limit equilibrium method(LEM) is a common type of geotechnical engineering stability analysis method. However, it requires numerous assumptions and simplifications in the shape of the sliding surface and the forces acting on the slices. The discontinuity layout optimization(DLO)method is a soil stability analysis method that has gradually been developed in recent years. Based on the method of linear optimization of force and deformation, DLO can achieve results with arbitrary shapes of sliding surfaces, which can overcome some of the disadvantages of the limit equilibrium method. This paper uses the Mohr-Coulomb non-associative flow rule correction coefficient proposed by Davis to correct the soil strength parameters adopted by the DLO method and analyzes the influence of the soil strength parameters c and φ on the position and shape of the sliding surface. On this basis, the paper develops computational procedures for rock-soil stability analysis and uses a series of examples with the DLO method and the LEM method to carry out comparative analysis. The results show that the stability coefficient and sliding surface morphology obtained by the DLO method of Mohr-Coulomb are close to those from the LEM method. Compared with the LEM, the DLO method is not limited by the assumptions regarding sliding surface morphology and can obtain the optimal potential sliding surface position and the rupture pattern of the unstable strip block based on nonlinear programming, thereby providing the basis for slope reinforcement.