A new operational matrix for solving two-dimensional nonlinear integral equations of fractional order

In this paper, first, we derive the operational matrix of two-dimensional orthogonal triangular functions (2D-TFs) for two-dimensional fractional integrals. Then, we apply this operational matrix and properties of Two-dimensional orthogonal triangular functions to reduce two-dimensional fractional integral equations to a system of algebraic equations. Finally, in order to show the validity and efficiency, we present some numerical examples.