The main goal of this paper is to solve a class of Darboux problems by converting them into the two-dimensional nonlinear Volterra integral equation of the second kind. The scheme approximates the solution of these integral equations using the discrete Galerkin method together with local radial basis functions, which use a small set of data instead of all points in the solution domain. We also employ the Gauss–Legendre integration rule on the influence domains of shape functions to compute the local integrals appearing in the method. Since the scheme is constructed on a set of scattered points and does not require any background meshes, it is meshless. The error bound and the convergence rate of the presented method are provided. Some illustrative examples are included to show the validity and efficiency of the new technique. Furthermore, the results obtained demonstrate that this method uses much less computer memory than the method established using global radial basis functions.
