A variational method to approximate the field of weakly guiding optical fibers by Laguerre-Gauss/Bessel expansion