We present a formalism based on the method of moment to solve the volume integral equation using tetrahedral (3-D) and triangular (2-D) elements. We introduce a regularization scheme to handle the strong singularity of the Green's tensor. This regularization scheme is extended to neighboring elements, which dramatically improves the accuracy and the convergence of the technique. Scattering by high-permittivity scatterers, like semiconductors, can be accurately computed. Furthermore, plasmon-polariton resonances in dispersive materials can also be reproduced.