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 the A-nucleon system using its relationship with the overlap functions of the eigenstates in the (A-1)-nucleon system.