We present new method for the numerical reconstruction of the variable refractive index of multi-layered circular weakly guiding dielectric waveguides using the measurements of the propagation constants of their eigenwaves. Our numerical examples show stable reconstruction of the dielectric permittivity function ε for random noise level using these measurements.