Using operational matrix of two-dimensional Bernstein polynomials for solving two-dimensional integral equations of fractional order

In this paper, we first preset two-dimensional Bernstein polynomials operational matrix (2DBPOM) for two-dimensional fractional integral equations (2DFIE). Then, we introduce a computational method based on 2DBPOM to reduce the 2DFIE to those of solving a system of algebraic equations. Also, we give the error analysis. Finally, this method is applied for solving three examples to illustrate the applicability and accuracy of presented method.