A new efficient method named matrix method is presented based on the traditional calculation methods of stabilizing force. The method is derived from the thrust force method and designed to address the deficiency of traditional method. It cannot only get the best pile location with minimum stabilizing force when is applied to the design calculation of single row anti slide piles, but also be applied to the optimal calculation of multi row anti slide piles. The cases comparing results show the method is more reasonable. For major projects of landslide control, the combination of matrix method and finite element method can make the result safer and more accurate.