We propose here a computational approach that is very efficient in performing the whole diagonalization of a configuration interaction (CI) matrix. The efficiency of the approach has been established by calculating all the eigenvectors of a semi-empirical full CI matrix representing the pi structure of the 1,2-benzopentacene. This CI matrix, with a dimension of the order of a million, was diagonalized in a few hours with a standard PC.