Computing the fundamental solutions for equations of electrodynamics

An analytical method for computing a generalized solution of a linear symmetric hyperbolic system is suggested. This method is based on the 3D Fourier transformation, matrix computations, a special approximation (regularization) of the Dirac delta function and 3D inverse Fourier transform of generalized functions. The implementation of this method has been done in MATLAB for computing the fundamental solutions of differential equations of electrodynamics in bi-anisotropic materials. Computational experiments have confirmed the robustness of the method.