Recently a new approximate analytical method based on simple physical reasoning for the calculation of surface plasmon frequencies in the complexes of metallic nanoparticles (MNP) (particularly spheres) has been proposed . The method called Eliminated Quadrupole Moment Approximation (EQMA) allows adequate description of existing experimental data concerning surface plasmons in the system of two coupled spheres . We present the results of numerical simulation based on Discrete Dipole Approximation (DDA) () performed in order to identify the limits of applicability of EQMA. It is shown that EQMA works very well for interparticle distances down to 0.1 of the diameter of the sphere. It is demonstrated also that obtaining of reliable numerical results with use of DDSCAT (simulation tool for DDA, ) requires presentation of the spheres as an ensemble of more than 500000 point dipoles, whereas in EQMA each of the spheres is substituted by only one dipole. The obvious advantage of EQMA is that it provides the resonance frequency for the given values of the parameters of the problem practically instantly while the DDSCAT requires several hours to obtain the same result.© (2010) COPYRIGHT SPIE--The International Society for Optical Engineering. Downloading of the abstract is permitted for personal use only.