Computation of a large group of interior eigenvalues at the middle spectrum is an important problem for quantum many-body systems. where the level statistics provides characteristic signatures of quantum chaos. We propose an exact numerical method. dual applications of Chebyshev polynomials (DACP). https://www.thebanyantreers.shop/product-category/jug/