An integral equation method for the numerical solution of a Dirichlet problem for second-order elliptic equations with variable coefficients

We develop a numerical approximation involving boundary integral techniques for the solution of the Dirichlet problem for second-order elliptic equations with variable coefficients. Using the concept of a parametrix, the problem is reduced to a boundary-domain integral equation to be solved for two unknown densities. Via a change of variables based on shrinkage of the boundary curve of the solution domain a parameterised system of boundary-domain integrals is obtained. It is shown how to write the singularities in this system in an explicit form such that boundary integral techniques can be applied for analysis and discretisation. An effective discretisation involving the Nystrom method is given, together with numerical experiments showing that the proposed approach can be turned into a practical working method.