“…From their set, we will then separate the LM fields. The wave equation has particular solutions in the form of waves [9] with the field E x = Ψ(y) exp(iβ κ z), where the function Ψ(y) is determined by the relations Ψ(y) = u (1) κ exp[iκ 1 (y + d)] + v (1) κ exp[−iκ 1 (y + d)], y < −d; u (2) κ exp[−iκ 2 (y − d)] + v (2) κ exp[iκ 2 (y − d)], y > d (1) for |y| > d. Here, β κ = k 2 n 2 1 − κ 2 1 are the propagation constants of modes and κ 1 and κ 2 are the transverse wave numbers in the substrate and the superstrate, respectively (κ 2 = k 2 (n 2 2 − n 2 1 ) + κ 2 1 ). Since κ 2 is expressed via κ 1 , in what follows we will assume that all quantities are functions of κ = κ 1 .…”