Robust Numerical Electromagnetic Eigenfunction Expansion Algorithms