The present paper studies the numerical computation of the extreme eigenvalues of a n x n real symmetric matrix A by the means of the Newton?s approximate method for the characteristic polynomial P_A (\lambda) An iterative algorithm is also presented,involving the computation of a trace of an appropriate matrix, instead of using the evaluation of P_A (\lambda) and its derivative. Numerical examples solved with this algorithm are to be found within as well.