“…Therefore, finding efficient numerical methods to approximate the solutions of these equations has become the main objective of many mathematicians. Some of these methods include Legendre wavelets [17], higher-order finite element method [18], generalized differential transform method [27], shifted Legendre polynomials [16,21,25], hybrid of block-pulse functions and shifted Legendre polynomials operational matrix method [31], Müntz-Legendre wavelets [32], fractional-order orthogonal Bernstein polynomials [38], delta functions operational matrix method [39], hybrid of block-pulse and parabolic functions [37], hat functions [35,40], two-dimensional orthonormal Bernstein polynomials [41][42][43], two-dimensional block-pulse operational matrix method [44], homotopy analysis method [47], Haar wavelet [4,49], orthonormal Bernoulli polynomials [52], shifted Jacobi polynomials [20,54,56], Bernstein polynomials [30,55], the second kind Chebyshev wavelets [51], etc. In this research study, some classes of two-dimensional nonlinear fractional integral equations of the second kind are considered in the following forms:…”