This study introduces a polyharmonic framework for analyzing the electromagnetic (EM) field generated by an oscillating point charge near a dispersive bulk of size comparable to the wavelength under study. We critically evaluate traditional approaches such as Liénard-Wiechert, Landau, and Raimond, and propose a Fourier representation of sources that simplifies numerical implementation and enhances analytical clarity. Our method effectively addresses the limitations of conventional models and is applicable to both relativistic and non-relativistic scenarios. It includes the oscillating point dipole fields, providing a comprehensive understanding of the EM field behavior. The Finite Element Method (FEM) is employed for numerical analysis, demonstrating the method’s adaptability to complex geometries. While offering significant insights, this study acknowledges certain limitations and outlines directions for future research.