Purpose - To investigate the numerical aspects of the electromagnetic scattering of a plane wave by a set of buried cylinders Approach - The cylindrical wave approach is employed. The analytical model is implemented in a Fortran code. The numerical aspects of the technique are presented, with particular emphasis on the numerical evaluation of the spectral integrals involved in the procedure. Findings - The implemented tool allows a fast computation of the electromagnetic field scattered by an arbitrary disposition of circular cylinders below an interface. Comparisons with the Finite Element Method are proposed, showing very good agreement. Originality and value - The computational advantages of the proposed technique are explained. The method can be used in a wide class of applications, for example in Ground Penetrating Radar tasks, microscopy, biomedical applications.