Numerical modeling of multimode doped optical waveguides