High‐order methods for elliptic equations with variable coefficients

In this article, we give a simple method for developing finite difference schemes on a uniform square gird. We consider a general, two-dimensional, second-order, partial differential equation with variable coefficients. In the case of a nine-point scheme, we obtain the known results of Young and Dauwalder in a fairly elegant fashion. We show how this can be extended to obtain fourth-order schemes on thirteen points. We derive two such schemes which are attractive because they can be adapted quite easily bnto obtain formulas for gird points near the boundary. In addition to this, these formulas only require nine evaluations for the typical forcing function. Numerical examples are given to demonstrate the performance of one of the fourth-order schemes.