Iterative numerical computation of the electromagnetic fields inside weakly nonlinear infinite dielectric cylinders of arbitrary cross sections using the distorted-wave Born approximation