Abstract
A higher-order Nystrom scheme is developed for discretizing the scalar potential integral equation (SPIE) for analyzing electromagnetic field/wave interactions on penetrable scatterers. The unknown scalar potential and its normal derivative are expanded using higher-order Lagrange interpolation functions. Inserting this expansion into the SPIE and point-testing the resulting equation yields a matrix system that is solved for the unknown expansion coefficients. Numerical results demonstrate that this matrix system is well-conditioned regardless of the discretization order and mesh density, and that its solution has higher-order convergence.