A vector iterative technique is developed for the extraction of eigenpairs related to the solution of finite element problems. The algorithm consists of using inverse iteration and conjugate gradient methods so as to obtain the solution vector associated to the smallest eigenvalue. Eigensolutions are sequentially calculated by replacing the coefficient matrix in the problem equilibrium equation using a deflation technique. The extensive usage of this technique, introduces multiple eigenvalue in the coefficient matrix, requiring a procedure to combine both methods. Also, a study is performed to find the appropriate starting vector to be used with methods. The algorithm has been implemented and the results of some example solutions are given that yield insight into its predictive capabilities.