  A method to calculate the hole spectral function in the discrete part of the
spectrum is suggested within the natural orbital representation of the one-body
density matrix of $A$-nucleon system using its relationship with the overlap
functions of the eigenstates in the $(A-1)$-nucleon system.
