Meshfree techniques based on infinitely smooth radial kernels have the great potential to2 provide spectrally accurate function approximations with irregular domain in high dimen-3 sions. The maximum accuracy can mostly be found when the RBF shape parameter is small,4 i.e., when the radial kernel is relatively smooth. However, as the shape parameter goes to5 zero, the standard RBF interpolant matrix will be very ill-conditioned. The ill-conditioning6 can be alleviated using alternate bases. One of these alternative bases is the Hilbert–Schmidt7 SVD basis. The Hilbert–Schmidt SVD approach suggests a stable mechanism for replacing8 a set of near-flat kernels with scattered centres to a well-conditioned base for exactly the9 same space. In this work, the Gaussian Hilbert–Schmidt SVD basis functions method is pre-10 sented to numerically solve the linear two-dimensional Fredholm integral equations of the11 second kind. The method estimates the solution by the discrete collocation method based12 on Gaussian Hilbert–Schmidt SVD basis functions constructed on a set of scattered points.13 The emerged integrals in the scheme are approximately computed by the Gauss–Legendre14 quadrature rule. This approach reduces the problem under study to a linear system of alge-15 braic equations which can be solved easily via applying an appropriate numerical technique.16 Also, the convergence of the proposed approach is established. Finally, numerical results are17 compared with standard RBF method to indicate the accuracy and efficiency of the suggested18 approach
