We present a new high-order integral algorithm for the solution of scatteri
ng problems by heterogeneous bodies under TE radiation, Here, a scatterer i
s represented by a (continuously or discontinuously) varying refractive ind
ex n(x) within a two-dimensional (2-D) bounded region; solutions of the ass
ociated Helmholtz equation under given incident fields are then obtained by
high-order inversion of the Lippmann-Schwinger integral equation. The algo
rithm runs in O(N log(N)) operations, where N is the number of discretizati
on points. Our method provides highly accurate solutions in short computing
times, even for problems in which the scattering bodies contain complex ge
ometric singularities.