KEYWORDSNon-polynomial spline; Finite di erence method; Group explicit iterative method; Elliptic partial di erential equation.Abstract. In this paper, we propose a new high-accuracy method based on nonpolynomial spline for the numerical solution to two-dimensional elliptic partial di erential equations. Using a non-polynomial spline approximation in x-direction and central di erence in y-direction, we obtain a new nine-point compact nite-di erence formulation. A four-point Group Explicit (GE) iterative scheme with an acceleration tool is then applied to the obtained system. The formulation procedure is presented in detail. The e ciency of the proposed method is then illustrated by some test problems. The numerical results are found to be in good agreement with the exact solutions.