  The scattering of electromagnetic pulses is described using a non-singular
boundary integral method to solve directly for the field components in the
frequency domain, and Fourier transform is then used to obtain the complete
space-time behavior. This approach is stable for wavelengths both small and
large relative to characteristic length scales. Amplitudes and phases of field
values can be obtained accurately on or near material boundaries. Local field
enhancement effects due to multiple scattering of interest to applications in
microphotonics are demonstrated.
