The Rayleigh quotient iterative algorithm with computational complexity O(N ) was originally suggested for finding the eigenvalues of an N-by-N matrix with fast convergence rate. In this paper, a modified payleigh quotient iterative (MRQI) algorithm is proposed to track the desired roots, i.e. the roots close to the unit circle of spectral polynomials. The MRQI algorithm associated with a zero suppression technique, called the parallel Rayleigh Quotient iterative (PRQI) algotithm, is further proposed to assure that rooting processors converge to different desired principal roots of spectral polynomials. The PROI alglrithm only with computational complexity O(N) tracks the nonstationary roots of spectral polynomials even starting from random initializations. Simulations show the better tracking performance of the suggested algorithm compared to that of the gradient Newton algorithm[11].