Abstract
In this work, we present a straightforward and simple method of moments (MOM) solution procedure, with minimum mathematical manipulations, to solve the coupled integral equations for multiple, homogeneous and inhomogeneous, dielectric bodies of arbitrary shape directly in the time domain. The standard surface and volume integral equation formulations are used for homogeneous and inhomogeneous bodies, respectively. The numerical solution procedure does not involve a time-marching process as is usually adopted for time domain problems and seems to be one of the primary reasons for the late-time instabilities as a result of error accumulation. The present solution method is stable for a very long time as evidenced by several representative numerical examples presented for validation.