This paper investigates a numerical method for solving two-dimensional nonlinear Fredholm integral equations of the second kind on non-rectangular domains. The scheme utilizes the shape functions of the moving least squares (MLS) approximation constructed on scattered points as a basis in the discrete collocation method. The MLS methodology is an effective technique for approximating unknown functions which involves a locally weighted least square polynomial fitting. The proposed method is meshless, since it does not need any background mesh or cell structures and so it is independent of the geometry of the domain. The scheme reduces the solution of two-dimensional nonlinear integral equations to the solution of nonlinear systems of algebraic equations. The error analysis of the proposed method is provided. The efficiency and accuracy of the new technique are illustrated by several numerical examples.
