In this paper, we consider an approximation to the Maxwell's eigenvalue problem based on a very weak formulation of two div-curl systems with complementary boundary conditions. We formulate each of these div-curl systems as a general variational problem with different test and trial spaces as done in [BP] , i.e., the solution space is (L2(D))3 and components in the test spaces are in subspaces of H1(D), the Sobolev space of order one on the computational domain D. A finite element least-squares approximation to these variational problems is used as a basis for the approximation. Using the structure of the continuous eigenvalue problem, a discrete approximation to the eigenvalues is set up involving only the approximation to either of the div-curl systems. We give some theorems which guarantee the convergence of the eigenvalues to those of the continuous problem without the occurrence of spurious values. Finally, some results of numerical experiments are given.