The integral representations of the solution around the vertices of the interior reentered angles (on the “singular” parts) are approximated by the composite midpoint rule when the boundary functions are from These approximations are connected with the 9-point approximation of Laplace's equation on each rectangular grid on the “nonsingular” part of the polygon by the fourth-order gluing operator. It is proved that the uniform error is of order where and is the mesh step. For the -order derivatives () of the difference between the approximate and the exact solutions, in each “ singular” part order is obtained; here is the distance from the current point to the vertex in question and is the value of the interior angle of the th vertex. Numerical results are given in the last section to support the theoretical results.
A fourth order block-hexagonal grid approximation for the solution of Laplace’s equation with singularities
