The continuation power flow method combined with the Jacobi-Davidson method is presented to trace the critical eigenvalues for power system small signal stability analysis. The continuation power flow based on a predictor- corrector technique is applied to evaluate a continuum of steady state power flow solutions as system parameters change; meanwhile, the critical eigenvalues are found by the Jacobi-Davidson method, and thereby the trajectories of the critical eigenvalues, Hopf bifurcation and saddle node bifurcation points can also be found by the proposed method. The numerical simulations are studied in the IEEE 30-bus test system.