The progress in two-dimensional materials has led to rapid experimental developments quantum plasmonics, where light is manipulated using plasmons. Although numerical methods can be used quantitatively describe plasmons spatially inhomogeneous systems, they are limited relatively small setups. Here, we present a novel semi-analytical method media within the framework of Random Phase Approximati...