Numerical solutions of the randall-wilkins and one trap one recombination models for first order kinetic

Randall-Wilkins and One Trap One Recombination (otor) models have been proposed to explain thermoluminescence emission and it should be emphasized that each model has its own allowed charge carrier transitions, trapping parameters and differential equations set. The equations are generally first or higher order linear differential equations with constant coefficients and their numerical solutions are an initial value problem. From this point on, numerical solutions of the thermoluminescence equations have been effectively used. In this paper the models were solved, numerically by using Euler and Runge-Kutta methods on Mathematica 8.0. In this work, although the fastest result calculated by Explicit Euler method, the most accurate results were calculated Linearly Implicit Euler method.