A method is outlined for calculating scattered fields from inhomogeneous penetrable objects using a coupled finite element-integral equation solution. The finite element equation can efficiently model fields in penetrable and inhomogeneous regions, while the integral equation exactly models fields on the finite element mesh boundary and in the exterior region. By decoupling the interior finite element and exterior integral equation meshes, considerable flexibility is found in both the number of field expansion points as well as their density. Only the nonmetal portions of the object need be modeled using a finite element expansion; exterior perfect conducting surfaces are modeled using an integral equation with a single unknown field since E(tan) is identically zero on these surfaces. Numerical convergence, accuracy, and stability at interior resonant frequencies am studied in detail.