  In this paper, we present a finite-element-extended boundary condition
(FE-EBC) method for an efficient calculation of the electromagnetic wave
scattering from inhomogeneous magneto-dielectric objects. To this end, we apply
the hierarchical Legendre polynomial basis functions on large curved
inhomogeneous hexahedral elements and propose an efficient numerical algorithm
for a fast computation of the finite-element matrix entries in the presence of
the material and Jacobian inhomogeneities. Also, we present a multiple-harmonic
expansion on the surface boundary to represent the surface magnetic field
accurately. The accuracy, efficiency and convergence of the method are studied
through some numerical examples.
