SUMMARYThe boundary integral representation of second-order derivatives of the primary function involves secondorder (hypersingular) and third-order (supersingular) derivatives of the Green's function. By defining these highly singular integrals as a difference of boundary limits, interior minus exterior, the limiting values are shown to exist. With a Galerkin formulation, coincident and edge-adjacent supersingular integrals are separately divergent, but the sum is finite, while the individual hypersingular integrals are finite. Moreover, the cancellation of the supersingular divergent terms only requires a continuous interpolation of the surface potential, and there is no continuity requirement on the surface flux. The algorithm is efficient, the non-singular integrals vanish and the singular integrals are computed entirely analytically, and accurate values are obtained for smooth surfaces. However, it is shown that a