Abstract-This paper demonstrates the one-dimensional computational results of the propagation of Gaussian electromagnetic pulse through dielectric slabs of finite thickness with variation in permittivity. The numerical approach used is the characteristic-based method solving the time-domain Maxwell curl equations involved with nonuniform permittivity. In the numerical model, all dielectric slabs are assumed to be isotropic, lossless, and linear. The permittivity of dielectric slab may increase or decrease linearly or sinusoidally. The numerical permittivity is finely discretized such that the variation between two adjacent grids is so small that the non-uniform permittivity is assumed to be piecewise continuous and consequently can be modeled as an individual block. The numerical results of various electric fields, both in the time-and frequency-domain, are presented and compared based on the dielectric slab of constant permittivity for close investigating the effects of the non-uniform permittivity distribution on the electromagnetic fields. It is also shown that under certain arrangement of Gaussian electromagnetic pulse and dielectric slab thickness the pattern of field propagation, reflection and transmission, can be reproduced in different time scales and frequency ranges.