Efficient numerical method for TDS-spectrum simulation is developed. TDS-modeling is required only integration of a ODE-system of low order. Desorption and diffusion activation energies interplay leads to several TDS-peaks. The effect of diffusion and desorption parameters on TDS-spectrum is illustrated.