Thermally stimulated current analysis has been performed on silicon detectors irradiated with high fast neutron fluences. The complexity of the TSC spectra required the development of a numerical procedure to deconvolute peak components. The numerical code has been used to directly calculate trap parameters: a comparison between numerical and experimental results is shown. Twentyone peaks related to various lattice defects have been found for the highest fluence of irradiation of 10(14) n/cm(2).