A method of determining the eigenvalues and eigenvectors of structure vibrations is proposed which is based on the reciprocity of the orthogonality theorem and the eigenvectors with respect to the stiffness and inertia matrices. The method is applied to the generalized problem of eigenvalues and eigenvectors (arbitrary inertia matrix) and makes it possible to calculate all the eigenmodes. It contains simple operations of elementary Gaussian transformation types. The cumulation of errors from one stage to the next during the iterative process is eliminated.