The reconstruction of cylindrically layered media is investigated in this article. The inverse problem is modeled using a source-type integral equation with a series of cylindrical waves as incidences, and a conventional Born iterative procedure is modified for solving the integral equation. In the modified iterative procedure, a conventional single-point approximation for the calculation of the field inside media is replaced by a multi-points approximation to improve the numerical stability of its solution. Numerical simulations for different permittivity distributions are demon- strated in terms of artificial scattering data with the procedure. The result shows that the procedure enjoys both accuracy and stability in the numerical computation.