The highly accurate block-grid method for solving Laplace’s boundary value problems on polygons is developed for nonanalytic boundary conditions of the first kind. The quadrature approximation of the integral representations ...