Abstract-This paper proposes a multilevel Green's function interpolation method (MLGFIM) to solve electromagnetic scattering from objects comprising both conductor and bi-isotropic objects using volume/surface integral equation (VSIE). Based on equivalence principle, the volume integral equation (VIE) in terms of volume electric and magnetic flux densities and surface integral equation (SIE) in terms of surface electric current density are first formulated for inhomogeneous bi-isotropic and conducting objects, respectively, and then are discretized using the method of moments (MoM). The MLGFIM is adopted to speed up the iterative solution of the resultant equation and reduces the memory requirement. Numerical examples are presented to show good accuracy and versatility of the proposed algorithm in dealing with a wide array of scattering problems.