Abstract. In this research, a new three-dimensional limit equilibrium method was developed. In the proposed method, the slip surface was assumed to be spherical and slices were assumed along the slip surface radius converging to a point. Moreover, force and moment equilibrium equations were used. In order to calculate Factor of Safety (FS) with the proposed method, a code was developed for two-and three-dimensional states. The model used the iteration (trial and error) method to simultaneously satisfy force and moment equilibrium equations. In order to ensure convergence of the numerical model developed in the three-dimensional state, the under-relaxation method was used in the trial and error operations. In examples, the value of FS obtained by the proposed method was 2.5 to 12 percent higher than those obtained by the other methods. In order to validate the proposed method and assess the functionality of the three-dimensional numerical model, a few examples were solved under di erent circumstances; the results were properly compliant with the results of other methods.