We present a family of unconditionally stable algorithms, based on the Suzuki product-formula approach, that solve the time-dependent Maxwell equations in systems with spatially varying permittivity and permeability. Salient features of these algorithms are illustrated by computing the density of states and by simulating the propagation of light in a two-dimensional photonic material.