In recent years, there is a considerable and growing interest in developing fast integral equation methods for solving Maxwell’s equations. Volume integral equations are a versatile technique to model inhomogeneous scattering objects. Numerical tests show excellent convergence properties of the technique for strongly inhomogeneous media with high dielectric contrast.