Solution of fractional-order differential equations based on the operational matrices of new fractional Bernstein functions

An algorithm for approximating solutions to fractional differential equations (FDEs) in a modified new Bernstein polynomial basis is introduced. Writing x→x^α(0<α<1) in the operational matrices of Bernstein polynomials, the fractional Bernstein polynomials are obtained and then transformed into matrix form. Furthermore, using Caputo fractional derivative, the matrix form of the fractional derivative is constructed for the fractional Bernstein matrices. We convert each term of the problem to the matrix form by means of fractional Bernstein matrices. A basic matrix equation which corresponds to a system of fractional equations is utilized, and a new system of nonlinear algebraic equations is obtained. The method is given with some priori error estimate. By using the residual correction procedure, the absolute error can be estimated. Illustrative examples are included to demonstrate the validity and applicability of the presented technique.