The equations governing the traffic of charge carriers during the filling, by ionizing radiation, of traps and luminescence centers in an insulator were numerically solved. A set of 4 simultaneous differential equations were used governing the time-dependent functions of concns. of e and holes in the conduction and valence bands and in traps and centers. The results are more general and accurate than those reported previously, since no assumptions concerning the proximity to equil. have to be made. The level of trapped charge n was calcd. at time t+T, and this is taken to be a good representation of the trapped charge d. Results were obtained for very high and very low dose rates of radiation. Exptl. findings of the dose dependence of thermoluminescence (TL) are susceptible to anal. by this approach. Under an appropriate choice of parameters, superlinearity emerges, although the results are not identical to those of the previous approxns. In addn., an important result from the anal. is the possible dependence of TL output on the dose rate for a const. total dose. Recent exptl. results of such a dependence on TL in quartz are in general accord with the numerical results.