A new algorithm is introduced to determine continuous deep energy level distributions from broad DLTS spectra. It is based on a series expansion of the energy level distribution and uses the DLTS spectrum and its temperature derivatives. A digitized transient is not necessary. The mathematical derivation of the method is given, and it is applied to simulated DLTS spectra in order to test accuracy and energy resolution. Special attention is paid to methods of noise suppression. An iterative correction of the obtained energy level distribution allows further improvement in accuracy.