A new time-domain numerical method for full-wave electromagnetic scattering and propagation is presented. The propagator method can be described as a numerical integral operator technique for moving the electromagnetic field through numerical space at successive increments in time. The numerical equations are obtained by discretizing electric and magnetic field propagator integrals. A primary feature of the propagator method is that all six electromagnetic-field components are calculated at each numerical grid point and all components are in time synchronization. Constant spatial and time increments are maintained throughout an inhomogeneous numerical space by employing an extrapolation procedure that ensures causality. A simple and effective first-order absorbing boundary condition, described as the null boundary condition, is introduced. Examples, provided in one-, two-, and three-dimensions, include plane wave reflection from and transmission through a planar boundary, and scattering and radar cross sections for canonical dielectric objects.