Two-dimensional shifted Legendre polynomials operational matrix method for solving the two-dimensional integral equations of fractional order

This work approximates the unknown functions based on the two-dimensional shifted Legendre polynomials operational matrix method (2D-SLPOM) for the numerical solution of two-dimensional fractional integral equations. The present method reduces these equations to a system of algebraic equations and then this system will be solved numerically by Newton's method. Moreover, an estimation of the error bound for this algorithm will be shown by preparing some theorems. Some examples are presented to demonstrate the validity and applicability of the proposed method with respect to the two-dimensional block pulse functions method (2D-BPFs) and two-dimensional Bernstein polynomials operational matrix method (2D-BPOM).