In nonlinear finite element stability analysis of structures, the foremost necessary procedure is the computation to precisely locate a singular equilibrium point, at which the instability occurs. The present study describes global and local procedures for the computation of stability points including bifurcation points and limit points. The starting point, at which the procedure will be initiated, may be close to or arbitrarily far away from the target point. It may also be an equilibrium point or non-equilibrium point. Apart from the usual equilibrium path, bypass and homotopy path are proposed as the global path to the stability point. A local iterative method is necessary, when it is inspected that the computed path point is sufficiently close to the stability point.