The light propagation in a structure with regularly placed metal nanoparticles is considered. The modeling is based on the Maxwell equations for the electromagnetic waves. We propose a new numerical technique for solving such problems. The well-known finite-difference time-domain (FDTD) simulation technique deals with the solution of the time-dependent wave equations. To the contrary, the new approach allows to find the stationary amplitudes of the electromagnetic oscillations directly from the linear algebraic equation system. An advantage of this method is its numerical stability, while finding an efficient method of the solution requires a further investigation. The resulting system of algebraic equations can be solved, e. g., by the Gaussian exclusion method. An essential limitation of the latter method is that the required computational resources increase very rapidly with the number of grid points. Finding a more efficient method could lead to a breakthrough in solving such physical problems, since the new numerical technique does not require any integration over time and, thus, could be really fast. As one of the possibilities, we have found that these algebraic equations can be solved by the steepest descent method, minimizing the squared l_2-norm of the equation system. The developed numerical technique is applicable to data-intense problems related to reflection, scattering and absorption of electromagnetic waves in various structures, including solar cells and hybrid liquid crystal cells as relevant examples considered in some detail.