A new numerical method for scattering from inhomogeneous bodies is presented. In particular, the 2D case of a TM-polarizated incident wave scattered by an infinite cylinder is considered. The scattered field is sought in two different domains. The first one is a bounded region inside the scattering body with an inhomogeneous permittivity epsilon(x, y). The second one is an unbounded homogeneous region outside the scatterer. An approximate solution for the scattered field inside the scatterer is sought by applying the QTSM technique. The method of discrete sources is used to approximate the scattered field in the unbounded region outside the scattering body. A comparison of the numerical solution with an analytic solution is performed. (C) 2002 Elsevier Science Ltd. All rights reserved.