A fast numerical method for harmonic equation based on natural boundary integral

Based on the coupling of the natural boundary integral method and the finite elements method, we mainly investigate the numerical solution of Neumann problem of harmonic equation in an exterior elliptic. Using our trigonometric wavelets and Galerkin method, there obtained stiffness matrix is symmetrical and circulant, which lead us to a fast numerical method based on fast Fourier transform. Furthermore, we do not need to compute the entries of the stiffness matrix. On the other hand, we prove that the numerical solution possesses exponential convergence rate. Especially, examples state that our method still has good accuracy for small j when the solution u ₀(θ) is almost singular.