The dose dependence of the concentration of point defects in alkali-halides as well as other crystals, as exhibited by the dependence of the thermoluminescence (TL), optical absorption and ESR on the dose of non-ionizing UV excitation is studied using numerical simulation. The relevant set of coupled rate equations are first written and plausible sets of trapping parameters are chosen. Instead of using simplifying assumptions previously used for reaching conclusions concerning this dose behavior, exact numerical solutions have now been reached. Depending on the parameters chosen, different dose dependencies are seen. In some cases, linear dose dependence is reached in a broad range. Sublinear dose dependence, close to a D 1/2 dependence when D is the dose of excitation can be reached when retrapping is stronger than trapping in other traps stabilizing the defects. When strong competition between stabilizing traps takes place, an initial linear range is observed followed by strong superlinearity and an approach to saturation. All these behaviors have been observed experimentally in TL measurements as well as ESR and optical absorption in different materials. Similarities and dissimilarities to linear and non-linear dose dependencies obtained experimentally and by simulations when ionizing irradiation is used for excitation are discussed.