In this paper, we introduce a numerical method for the solution of two-dimensional Fredholm integral equations. The method is based on interpolation by Gaussian radial basis function based on Legendre–Gauss–Lobatto nodes and weights. Numerical examples are presented and results are compared with the analytical solution to demonstrate the validity and applicability of the method.
