Abstract A computational procedure for the analysis of thermoluminescence (TL) glow peaks based on a simple model for TL, assuming transitions via the conduction band has been described. The method takes into consideration the full shape of the glow peak and utilizes the curve fitting technique and a general order kinetic equation. A consistent value for trap depth (E) corresponding to the observed peak temperature ( T m) is first calculated for chosen values of other parameters and is used to generate a theoretical glow curve which is then compared with the experimental one. By consideration of the various restraints the parameters are altered to obtain a good fit. A set of parameters which produce a glow curve differing only slightly from the experimental one in the least-square sense is taken as the final set to describe the TL process. TL parameters obtained in the case of MgO by application of this method are reported.