The authors report on an analysis of conditions on the boundary between layers having varied electromagnetic properties. The research is performed using consistent theoretical derivation of analytical formulas, and the underlying problem is considered also in view of multiple boundaries, including the effect of the propagation of electromagnetic waves exhibiting different instantaneous speed. The paper includes a theoretical analysis and references to the generated algorithms. The main algorithm was assembled to enable simple evaluation of all components of the electromagnetic field in relation to the speed of the wave propagation in a heterogeneous environment. The proposed algorithms are compared by means of different numerical methods for the modelling of electromagnetic waves on the boundary between materials; moreover, electromagnetic field components in common points of the model are also subject to comparison. The article describes the obstacles in comparison of different principles based analyses. Thus, analysis using numerical modelling performed via the wave equation in the ANSYS and the applied methods in the Matlab program. When in conjunction with tools facilitating the analysis of material response to the source of a continuous signal, the algorithms constitute a supplementary instrument for the design of a layered material. Such design enables the realization of, for example, recoilless plane, recoilless transition between different types of environment, and filters for both optical and radio frequen- cies. This phenomenon occurs in metamaterials.