the amplitudes $c^t_j$ and $c^r_j$ as $E_{j}(z)=c_j^t e^{i\omega z/c}+c_j^r e^{-i\omega z/c}$ and $cB_{j}(z)=c_j^t e^{i\omega z/c}-c_j^r e^{-i\omega z/c}$. Therefore, electromagnetic fields at adjacent layers are related by $$\begin{aligned}
 \begin{pmatrix}
  E_{j+1} \cr B_{j+1}
 \end{pmatrix}
 = 
\begin{pmatrix}
  \cos(\frac{\omega d}{c}) & i c \sin(\frac{\omega d}{c}) \cr
  \frac{i}{c} \sin(\frac{\omega d}{c})  & \cos(\frac{\omega d}{c})
 \end{pmatrix}
 \begin{pmatrix}
  E'_{j} \cr B'_{j}
 \end{pmatrix}.
 \label{eq:To}\end{aligned}$$ We combine Eq. (\[eq:bc\]) and Eq. (\[eq:To\]) to define the transfer matrix $$\begin{aligned}
 T_j
 = 
\begin{pmatrix}
  \cos(\frac{\omega d}{c}) & i c \sin(\frac{\omega d}{c}) \cr
  \frac{i}{c} \sin(\frac{\omega d}{c})  & \cos(\frac{\omega d}{c})
 \end{pmatrix}
 \begin{pmatrix}
  1 & 0 \cr
  -\frac{\sigma^N_j}{\epsilon_0 c^2} & 1 
 \end{pmatrix}
\label{eq:tmatrix}\end{aligned}$$ that satisfies $$\begin{aligned}
 \begin{pmatrix}
  E_{j+1} \cr B_{j+1}
 \end{pmatrix}
 =T_j
 \begin{pmatrix}
  E_j \cr B_j
 \end{pmatrix}.
 \label{eq:tmatr}\end{aligned}$$ The matrix is obtained from Maxwell’s equations and is the product of two matrices. The first matrix expresses the propagation of light in the interlayer space and the second matrix represents the boundary condition at a graphene layer. Our approaches are based on the approximation that the electronic current is sufficiently localizing at a graphene layer. But in fact, because electronic wave function of the $\pi$-orbital is slightly spreading in a space between graphene layers, light propagation in the space is likely to be subjected to the spread 