This talk presents a verified numerical method for computing a subset of eigenvalues of real symmetric (Hermitian) matrices. The proposed approach efficiently guarantees the accuracy of selected eigenvalues for dense matrices. The theoretical foundation required for the verification is described, and numerical experiments are presented to demonstrate the effectiveness of the method.