A numerical method of evaluating electromagnetic fields in a generalized anisotropic medium

A transition matrix method is commonly used to deal with the problems of either plane-wave scattering from or the Green's function of a generalized anisotropic medium. This method, although rigorous analytically, introduces numerical breakdown, when the layers are electrically thick and the waves are evanescent. A variable transformation method is developed to deal with the exponentially-growing terms associated with exponential-matrix method. The proposed scheme is suitable for the numerical analysis of generalized anisotropic layers including ferrites, magneto-plasmas, chiral layers, and bianisotropic layers.>