Exponentially convergent semi-discretization is proposed for multilayered waveguides.
Midpoint integration with linear finite elements is key to exponential convergence.
Finite layers are discretized using a mesh stretched into complex plane.
Half-space is modeled using an optimally graded mesh.
Implementation requires minimal modification to existing FE-based codes.