Fast and accurate computation of dominant modes for optical fibers with arbitrary profile through Pade' approximants