The analysis of Thermally Stimulated Depolarization Currents (TSDC) spectra in polymers requires the introduction of a distribution of relaxation times. Also, a decay of the induced polarization following a stretched exponential law can be used. Two different ways for introducing this effect in the TSDC equations are compared here together with the decomposition in elementary curves by a Simulated Annealing Monte Carlo procedure to find the distribution. It is found that the stretched exponential applied to computer generated curves corresponds in both cases to a nearly Debye process. Also the analysis of the alpha-peak of poly(epsilon-caprolactone) with the three approaches are compared.