We review the EMSolve code, a 3D finite element program for the time-dependent Maxwell's equations. At the core of the program is a class library of arbitrary-order H(curl) and H(div) basis functions, these basis functions have been developed for tetrahedral, hexahedral, and prism elements. We have experimented with a variety of degrees-of-freedom for these bases, and we have developed nearly optimal basis from the point of view of matrix conditioning. A current research issue is higher-order time integration, and we have preliminary evidence that symplectic Runge-Kutta methods are well suited for Maxwell's equations.