The aim of this paper is the accurate and efficient analysis surface-plasmon-mode resonances a graphene microdisk stack in terahertz range. By means suitable generalized boundary conditions Fourier series expansion, problem formulated terms sets one-dimensional integral equations vector Hankel transform domain for harmonics surface current densities. In virtue Helmholtz decomposition, unknowns ...