The paper deals with two essential and related closely processes involved in the finite element slope stability analysis in two-dimensional problems, i.e. computation of the factors of safety (FOS) and location of the critical slide surfaces. A so-called phi-v inequality, sin 0 >= 1 - 2v is proved for any elasto-plastic material satisfying Mohr-Coulomb's yield criterion. In order to obtain an FOS of high precision with less calculation and a proper distribution of plastic zones in the critical equilibrium state, it is stated that the Poisson's ratio v should be adjusted according to the principle that the phi-v inequality always holds as reducing the strength parameters c and phi. While locating the critical slide surface represented by the critical slide line (CSL) under the plane strain condition, an initial value problem of a system of ordinary differential equations defining the CSL is formulated. A robust numerical solution for the initial value problem based on the predictor-corrector method is given in conjunction with the necessary and sufficient condition ensuring the convergence of solution. A simple example, the kinematic solution of which is available, and a challenging example from a hydraulic project in construction are analysed to demonstrate the effectiveness of the proposed procedures. Copyright (c) 2005 John Wiley & Sons, Ltd.