Abstract: A computation algorithm based on the Poincaré Mapping in combination with Pseudo-Arc Length Continuation Method is presented for calculating the unstable response with saddle-node bifurcation, and the singularity, which occurs using the general continuation method combined with Poincaré Mapping to follow the path, is also proved. A normalization equation can be introduced to avoid the singularity in the process of iteration, and a new iteration algorithm will be presented too. There will be two directions in which the path can be continued at each point, but only one can be used. The method of determining the direction will be presented in the paper. It can be concluded that this method is effective in analysis of nonlinear dynamic system with saddle-node bifurcations.