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.roneverhart.com/33-14-50R24-HD878-R-T-Special/