In contrast to the classical boundary layer approach which is not uniquely solvable at the irregular wave numbers, the modified layer method always has a unique solution to the exterior boundary problem for the vector Helmholtz equation. We compare the stability of the two methods numerically.