In this article, a numerical method is presented for solving non-linear multi-order fractional differential equations (FDEs) using orthonormal Boubaker polynomials as basis functions. The operational matrix is constructed for the product of basis functions. The fractional derivative is considered in the Caputo sense. Using this technique, non-linear multi-order FDEs converted into a system of non-linear algebraic equations, which can then be solved numerically. The proposed method is implemented using MATLAB for numerical analysis. Using numerical results, validity, and accuracy of the proposed method is given and compared the results with existing literature
Read full abstract