Abstract: In order to analyze the electromagnetic scattered properties of a 2-D given body with the TEz wave incidence, a meshless time-domain algorithm subject to the incident TEz wave is developed. Based on the clouds of points, a weighted least square technique is introduced to approximate the spatial derivatives at the discrete points, and then an approximate Riemann solver is used for computing the numerical flux related to the governing equations. Physical quantities at the midpoint are determined by the quadratic function reconstruction instead of the traditional linear function reconstruction so as to increase the numerical accuracy. After that, an explicit four-stage Runge-Kutta scheme is employed in time-marching. By using the developed algorithm, scattered fields and bistatic radar cross sections (RCS) are calculated for the cases of a 2-D aerodynamic shape with an external store. Numerical results show that all the factors that include the direction of the incident wave and the size of the external store as well as the relative position between the external store and the main body can affect the total electromagnetic scattered properties.