Numerical solution of the glow curve differential equations

The three basic simultaneous differential equations describing the phenomena of thermoluminescence and thermally stimulated conductivity are solved numerically for any given set of trapping parameters and without any additional assumptions, thus enabling us to simulate the TL and TSC phenomena. The method employs a change of variable by which the high sensitivity of the solution to small numerical errors is overcome. Use is then made of an improved Runge-Kutta method for the solution of the new set of differential equations. Previous methods for extracting information from glow curves are tested using these calculated curves as data. The often used assumption | ṅ_c | ≪ | ṅ |, | ṁ | is explored and found to be correct for most of the temperature range of interest.