The temporal behavior of electric fields in arbitrary double-negative planar slabs is systematically investigated in this paper, from both analytical and numerical perspectives. Concerning infinite slabs, a set of exact expressions for an exponential current excitation is derived through an efficient complex analysis, and an integrated study of surface polariton frequencies is performed. Subsequently, the significant case of a source with a random spatial profile is explored in order to obtain rigorous relations for the field and transient phenomena damping time with respect to problem parameters. On the other hand, a robust finite-difference time-domain methodology is introduced for the comprehensive examination of finite slabs, whose numerical simulations dictate the adoption of a resonatorlike discipline. This inevitable, yet very instructive, convention is physically justified by the almost perfect surface mode reflections at the edges of the slab. In this manner, the proposed formulation reveals a prominent increase in the excited polariton amplitude, relative to the corresponding infinite arrangements, which leads to larger transient times.