Abstract.A high-order, leapfrog integrator is presented and analyzed for Maxwell's time domain equations. The integrator is shown to be quite robust in terms of its efficiency and memory demands. Moreover, the integrator retains the original simplicity of the traditional second-order, leapfrog integrator. Rigorous Fourier analyses are presented to quantify the dispersion, dissipation, and stability properties of the integrator. To limit the scope of the presentation, the 2x2, 2x4, 4x2, and 4x4 schemes are examined and compared. Here the first digit denotes temporal accuracy; the second digit denotes spatial accuracy. Numerical demonstrations, using the three-dimensional rectangular waveguide as the object under test, are provided to further substantiate the theoretical analysis.