A generalized donsker-kač formula to compute the fundamental modes in complex loaded waveguides