A continuation method is developed for solving the inverse medium scattering problem of time harmonic Maxwell equations in R 3 . By using multifrequency scattering data, our reconstruction algorithm first employs the Born approximation for an initial guess and proceeds via recursive linearization on the wavenumber k. At each linearization step, one forward and one adjoint state of the Maxwell equations are solved. Numerical examples are also presented and discussed.