The axis-symmetric modeling technique is based on expanding vector fields in cylindrical harmonics and computing the response on a two-dimensional cross-section separately for each azimuthal harmonic, significantly reducing computational costs. However, it has limitations when dealing with dipoles placed away from the symmetry axis due to challenges in the expansion of angular modes. To address this, we propose a reformulated axis-symmetric model based on the Fourier expansion of the delta function distribution concerning the azimuthal variable. This model is validated using standard Mie theory for off-axis dipoles and applied to study multiple-emitter-plasmon interactions. The emission properties of a non-cooperative ensemble near a plasmonic nanoparticle are observed to scale with the number of emitters considered, N. Notably, a Dicke effect-like superradiance (N2-dependence) is observed when a spatially disordered ensemble of dipoles oscillates collectively inside a plasmonic dimer gap. This kind of high-level cooperative quantum phenomenon is of high interest in fields such as quantum optics and light-harvesting.