 reflected field is $E_{1-}$ and the field that is transmitted through the entire structure is $E^{\prime}_{(N+1)+}$. To compute the transmission spectra for the system, first we must construct the transfer matrix, $M$, from the individual interface matrices, $I_j$, and propagation matrices, $P_j$, where the index, $j$, specifies the layer in question [@Orfanidis], $$I_j=
\frac{1}{\tau_j}
\begin{pmatrix}
1            &  r_j \\[6pt]
r_j    &   1        \\
\end{pmatrix}$$

$$P_j=
\begin{pmatrix}
e^{2 \pi i n_j d_j \xi}    & 0                                         \\[6pt]
0                                    &   e^{-2 \pi i n_j d_j \xi}        \\
\end{pmatrix}$$

where $r_j$ and $\tau_j$ are the reflection and transmission coeffiecents, respectively, $$r_j = \frac{\mu_{j+1} n_{j} - \mu_{j} n_{j+1} }{\mu_{j+1} n_{j} + \mu_{j} n_{j+1} }$$ $$\tau_j = \frac{2 \mu_{j+1} n_{j}  }{\mu_{j+1} n_{j} + \mu_{j} n_{j+1} }$$ In scaled variables, the phase argument, $i k_j l_j$ becomes $2 \pi i n_j d_j \xi$. The frequency, $f$, becomes $\xi=f\Lambda/c_0$, where $c_0$ is the speed of light in vacuum. The incident and scattered field are related by,

$$\begin{pmatrix}
E_{1+}\\[6pt]
E_{1-}\\
\end{pmatrix}
=
\begin{pmatrix}
M_{11}  & M_{12} \\[6pt]
M_{21}  &   M_{22}           \\
\end{pmatrix}
\begin{pmatrix}
E^{\prime}_{(N+1)+}\\[6pt]
0\\
\end{pmatrix}$$

where, $$M =
reflected field is $ E_{1-}$ and the field that is transmitted through the integral social organization is $ E^{\prime}_{(N+1)+}$. To calculate the transmission spectra for the system, foremost we must construct the transfer matrix, $ M$, from the individual interface matrices, $ I_j$, and generation matrix, $ P_j$, where the index, $ j$, stipulate the layer in question [ @Orfanidis ], $ $ I_j= 
 \frac{1}{\tau_j } 
 \begin{pmatrix } 
 1             &   r_j \\[6pt ] 
 r_j     &    1         \\ 
 \end{pmatrix}$$ 

 $ $ P_j= 
 \begin{pmatrix } 
 e^{2 \pi i n_j d_j \xi }     & 0                                          \\[6pt ] 
 0                                     &    e^{-2 \pi i n_j d_j \xi }         \\ 
 \end{pmatrix}$$ 

 where $ r_j$ and $ \tau_j$ are the observation and transmission coeffiecents, respectively, $ $ r_j = \frac{\mu_{j+1 } n_{j } - \mu_{j } n_{j+1 } } { \mu_{j+1 } n_{j } + \mu_{j } n_{j+1 } } $ $ $ $ \tau_j = \frac{2 \mu_{j+1 } n_{j }   } { \mu_{j+1 } n_{j } + \mu_{j } n_{j+1 } } $ $ In scaled variable, the phase controversy, $ i k_j l_j$ becomes $ 2 \pi i n_j d_j \xi$. The frequency, $ f$, becomes $ \xi = f\Lambda / c_0 $, where $ c_0 $ is the speed of light in vacuum. The incident and confused field are related by, 

 $ $ \begin{pmatrix } 
 E_{1+}\\[6pt ] 
 E_{1-}\\ 
 \end{pmatrix } 
 = 
 \begin{pmatrix } 
 M_{11 }   & M_{12 } \\[6pt ] 
 M_{21 }   &    M_{22 }            \\ 
 \end{pmatrix } 
 \begin{pmatrix } 
 E^{\prime}_{(N+1)+}\\[6pt ] 
 0\\ 
 \end{pmatrix}$$ 

 where, $ $ M =
 revlected field is $E_{1-}$ and tme field that is transmmtted tgrough tfe entire structure is $E^{\primx}_{(N+1)+}$. Ti comklte the transmission rpectra flr the ststen, first we must construct bhe txaisfer matrix, $M$, nrom the ingividual interxaze matrices, $I_j$, and propagation matricqs, $P_j$, wnege the index, $j$, spebisies nht layer in question [@Orfanidis], $$I_j=
\fdac{1}{\tau_j}
\uegin{pmatrix}
1            &  r_j \\[6pt]
r_j    &   1        \\
\end{pmatrix}$$

$$P_j=
\begln{pmatrix}
e^{2 \kj i b_j d_j \xi}    & 0                                         \\[6pt]
0                                    &   e^{-2 \'i i n_t d_j \xi}        \\
\end{pmdtrix}$$

whrre $r_j$ and $\tau_m$ are thw reflection and tranvmission coeffiecegts, respewtnvely, $$r_j = \frac{\mu_{j+1} n_{j} - \my_{j} n_{j+1} }{\mu_{j+1} n_{j} + \nu_{j} n_{j+1} }$$ $$\vau_n = \fraf{2 \mn_{j+1} n_{j}  }{\mu_{j+1} h_{j} + \mu_{j} n_{j+1} }$$ In scaled variablex, eye phase argujent, $i k_t l_j$ becomes $2 \pi i n_j d_j \xi$. The frequenby, $f$, becomes $\xi=f\Lambda/c_0$, wheee $c_0$ is the speed of pight in dacuum. The incident and scattered field are relatad by,

$$\ueein{kmqtrix}
D_{1+}\\[6pt]
F_{1-}\\
\end{pmatrix}
=
\begin{pmatrix}
M_{11}  & M_{12} \\[6pt]
M_{21}  &   M_{22}           \\
\ekd{pmatrix}
\begin{pmctrix}
E^{\prime}_{(N+1)+}\\[6pt]
0\\
\enc{playtix}$$

where, $$M =
reflected field is $E_{1-}$ and the field transmitted the entire is $E^{\prime}_{(N+1)+}$. To the first we must the transfer matrix, from the individual interface matrices, $I_j$, propagation matrices, $P_j$, where the index, $j$, specifies the layer in question [@Orfanidis], \frac{1}{\tau_j} \begin{pmatrix} 1 & r_j \\[6pt] r_j & 1 \\ \end{pmatrix}$$ $$P_j= \begin{pmatrix} \pi n_j \xi} 0 \\[6pt] 0 & e^{-2 \pi i n_j d_j \xi} \\ \end{pmatrix}$$ where $r_j$ and $\tau_j$ the reflection and transmission coeffiecents, respectively, $$r_j = n_{j} - \mu_{j} n_{j+1} n_{j} + \mu_{j} n_{j+1} }$$ = \mu_{j+1} n_{j} n_{j} \mu_{j} }$$ In scaled the phase argument, $i k_j l_j$ becomes $2 \pi i n_j d_j \xi$. The frequency, $f$, becomes where $c_0$ speed of in The and scattered field by, $$\begin{pmatrix} E_{1+}\\[6pt] E_{1-}\\ \end{pmatrix} = M_{12} \\[6pt] M_{21} & M_{22} \\ \end{pmatrix} \begin{pmatrix} 0\\ \end{pmatrix}$$ $$M =
 reflected field is $E_{1-}$ and the fiEld that is tRansmIttEd tHrOugh The eNtire structure IS $E^{\prIme}_{(N+1)+}$. To compute the transmIssioN sPEctrA FoR the sYstem, fiRSt WE MusT cOnStrUcT ThE tranSfeR matrix, $m$, from the inDivIdUal interface MAtRices, $I_j$, and ProPagation matrIceS, $P_j$, wheRe The INdex, $j$, SpeCifieS the laYEr in quEstion [@OrfAnIDis], $$I_j=
\fRAc{1}{\tau_j}
\bEGIn{PmatRix}
1            &  r_j \\[6pt]
r_j    &   1        \\
\end{pmatrIX}$$

$$P_J=
\Begin{pmatrix}
e^{2 \pI i n_j d_j \Xi}    & 0                                         \\[6PT]
0                                    &   e^{-2 \PI I n_j D_j \xI}        \\
\end{pmatriX}$$

wHere $r_J$ And $\tau_j$ ARe THE RefLEction and tranSmission coeFFieCents, rEsPecTIvely, $$r_J = \frac{\Mu_{J+1} N_{j} - \mU_{j} n_{j+1} }{\mu_{j+1} n_{j} + \mu_{J} n_{j+1} }$$ $$\tAu_j = \frac{2 \mu_{J+1} n_{j}  }{\mu_{j+1} N_{J} + \mu_{j} n_{j+1} }$$ IN Scaled vAriablEs, tHe pHase ARgUmEnt, $I k_J L_j$ bECoMes $2 \PI i n_J d_j \xi$. The FrEqUency, $F$, becOMES $\Xi=f\LAmbDa/c_0$, wHere $c_0$ Is the speed of lIghT in vACuuM. The iNcideNt anD sCatteRed fieLd are ReLated by,

$$\begin{pmaTrix}
e_{1+}\\[6pt]
E_{1-}\\
\end{pmAtrIx}
=
\BegIn{PmatrIX}
M_{11}  & M_{12} \\[6pt]
M_{21}  &   m_{22}           \\
\enD{pmAtrix}
\beGin{pmatRIx}
E^{\PrIME}_{(n+1)+}\\[6pT]
0\\
\end{pmatrix}$$

where, $$M =
 reflected field is $E_{1- }$ and the  fiel d t hat  i s tr ansm itted throught he e ntire structure is $E^ {\pri me } _{(N + 1) +}$.To comp u te t hetr an smi ss i on  spec tra  for th e system,fir st  we must con s tr uct the tr ans fer matrix,$M$ , from  t hei ndivi dua l int erface matric es, $I_j$ ,a nd pro p agation m at rice s, $P_j$, where t h ei ndex, $j$, spe cifies  t h el a yer  in  question[@ Orfan i dis], $ $ I_ j = \fr a c{1}{\tau_j}\begin{pmat r ix} 
1        &  r _j \\ [6 p t]r_j    &  1       \\
\ end{pm a trix}$$ 
$$P_j= 
\begi n{p mat rix} e^ {2  \p ii  n_ j  d _j\ xi}     & 0                            \\[6pt ]
0                            &  e^ {-2 \pi i n_j d _j \ xi}       \ \\en d{ pmatr i x}$$
whe re$r_j$ a nd $\ta u _j$  a r e th e reflection and t ra n s mi ssion co effiec e nt s, respecti ve ly,  $$r _ j  = \f rac{ \ mu _{j+1} n _{j} - \m u_ {j} n_{ j+ 1} }{\ mu _{j +1}  n_{j }  + \ mu_{j}  n_{j+1}  }$$$ $\tau_j = \fra c {2 \mu_{j+1}n _{ j }  } {\mu _{j +1} n_{j} +  \mu _ {j}n_{j + 1}  }$ $  In s caled  v a ri a bles, the phase arg um ent, $ i k_j  l_j$ becomes  $2 \pi in _ j  d_j \xi $. T h ef requency, $f$,  beco mes $\xi=f \ Lambda/c _0$,where $c _0$ is th e speed of  li ght  in  va c u um . The inciden t andsc attered  fi eld are  re lat edby, 
$$\begin{ pmatrix} 
E _{ 1+ }\ \[6 pt]
E _ {1-}\\
\ en d{p ma tri x}
=\ begin{ pmatr ix}M_ {1 1 } & M_{12 }  \ \ [ 6pt] 
M _{ 21} & M_{22 }        \\\end{pmat rix } 
\be gi n{ pmatrix }
E^{\prime}_ {( N+1)+}\\[6 pt ]
0 \\
\en d { pmatrix} $$

where, $$M =
 reflected_field is_$E_{1-}$ and the field_that is_transmitted_through the_entire_structure is $E^{\prime}_{(N+1)+}$._To compute the_transmission spectra for the_system, first we_must_construct the transfer matrix, $M$, from the individual interface matrices, $I_j$, and propagation matrices,_$P_j$,_where the_index,_$j$,_specifies the layer in question_[@Orfanidis], $$I_j=
\frac{1}{\tau_j}
\begin{pmatrix}
1   _ _      & _r_j_\\[6pt]
r_j  _ &   1    _   \\
\end{pmatrix}$$

$$P_j=
\begin{pmatrix}
e^{2 \pi i_n_j d_j \xi}___ & 0 _     _     _      _     _    _ _    _  \\[6pt]
0_ _  _     _  _     __  ___ _  __ __&   e^{-2 \pi i_n_j_d_j \xi}     _  \\
\end{pmatrix}$$

where $r_j$_and_$\tau_j$ are the reflection_and transmission coeffiecents, respectively, $$r_j_= \frac{\mu_{j+1} n_{j} - \mu_{j} n_{j+1}_}{\mu_{j+1} n_{j}_+ \mu_{j}_n_{j+1} }$$ $$\tau_j = \frac{2 \mu_{j+1} n_{j}  }{\mu_{j+1} n_{j} +_\mu_{j} n_{j+1} }$$ In scaled variables,_the phase argument, $i_k_j l_j$_becomes_$2 \pi i_n_j_d_j \xi$._The frequency, $f$, becomes $\xi=f\Lambda/c_0$, where $c_0$_is the_speed of light in vacuum. The_incident and scattered field_are_related by,

$$\begin{pmatrix}
E_{1+}\\[6pt]
E_{1-}\\
\end{pmatrix}
=
\begin{pmatrix}
M_{11}  & M_{12} \\[6pt]
M_{21}_ &   M_{22} _    ___  \\
\end{pmatrix}
\begin{pmatrix}
E^{\prime}_{(N+1)+}\\[6pt]
0\\
\end{pmatrix}$$

where, $$M =