Numerical curve fitting of glow peaks, using a computer, is shown to be more objective and accurate than the existing graphical curve fitting method in evaluating trap activation energies. Besides using the numerical curve fitting method for glow peaks of first-order kinetics, the method is extended to the case of second-order thermoluminescence or thermally stimulated current peaks. The method was checked both for numerically generated peaks and experimental thermoluminescence curves of NaCl and ZnS : Er3+ samples. The results are in good agreement with the given values of activation energies in the former case and values calculated by other means in the latter.