The solution to the problem of electromagnetic (EM) wave scattering from a set of smallsize impedance particles of arbitrary shape is derived by the asymptotic approach. Particles are located in a homogeneous domain with constant ɛ 0 and μ 0 . The solution is derived under the condition b → 0, where b is the characteristic size of the particle; further, the number M(b) of particles tends to infinity at a specific rate. The regularising procedure consists of the derivation of the explicit form of a solution that excludes the necessity to solve the respective integral equation for determination of the fields at the surface of particles and thus avoids integrating the Green function derivatives, which are in the kernel of this boundary integral equation. Practical application of this approach yields an ability to create media with the desired inhomogeneous distribution of effective refractive index n and magnetic permeability μ(x). Explicit analytical formulas are derived for these physical parameters and supported by computations.This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.