An efficient numerical method based on Fibonacci polynomials to solve fractional differential equations

The Fibonacci sequence is significant because of the so-called golden ratio, which describes predictable patterns for everything. Fibonacci polynomials are related to Fibonacci numbers, and in this paper we extend their applicability by using them to solve fractional differential equations (FDEs) and systems of fractional differential equations (SFDEs). With the help of the Riemann–Liouville fractional integral operator for the fractional-order hybrid function of block-pulse functions and the Fibonacci polynomials defined in this paper, the solution of the considered FDE and SFDE is reduced to a system of algebraic equations, which can be solved by Newton’s iterative method. The fractional order is obtained by transforming x into xα, with α>0. Compared to other models, our method in some situations is better by twelve orders of magnitude. There are situations when we get the exact solution. The presented method proves to be simple and effective in solving nonlinear problems with given initial values.