The development of an efficient numerical algorithm (computer code) capable of determining the unknown material composition ALMD shape of an arbitrary target from the measured electromagnetic waves in the far field region will enhance the capability of the defense radar system to defeat known evasive schemes this can be achieved by improving the proven iterative numerical algorithm, generalized pulsespectrum technique (GPST), by taking advantages of parallelism, vectorization, adaptive grids, domain decomposition, etc. The phase I effort is to develop an efficient computer code for calculating the time-harmonic electromagnetic waves scattered from a two-dimensional scatterer of arbitrary material composition and shape. This computer code is an important subroutine of GPST for solving the inverse scattering problems, and it can be 30 times faster than the existing finite difference time-zone method (with the same accuracy) by the estimation of the total floating point arithmetic operation count needed for the calculation.