The paper presents a new approach for the analysis of slope stability that is based on the numerical solution of a differential equation, which describes the thrust force distribution within the potential sliding mass. It is based on the evaluation of the thrust force value at the endpoint of the slip line. A coupled approximation of the slip and thrust lines is applied. The model is based on subdivision of the sliding mass into slices that are normal to the slip line and the equilibrium differential equation is obtained as the slice width approaches zero. Opposed to common iterative limit equilibrium procedures the present method is straightforward and gives an estimate of slope stability at the value of the safety factor prescribed in advance by standard requirements. Considering the location of the thrust line within the soil mass above the trial slip line eliminates the possible development of a tensile thrust force in the stable and critical states of the slope. The location of the upper boundary point of the thrust line is determined by the equilibrium of the upper triangular slice. The method can be applied to any smooth shape of a slip line, i.e., to a slip line without break points. An approximation of the slip and thrust lines by quadratic parabolas is used in the numerical examples for a series of slopes.