An efficient algorithm for solving higher-order fractional Sturm–Liouville eigenvalue problems

In this paper, we present a simple and efficient computational algorithm for solving eigenvalue problems of high fractional-order differential equations with variable coefficients. The method of solution is based on utilizing the series solution to convert the governing fractional differential equation into a linear system of algebraic equations. Then, the eigenvalues can be calculated by finding the roots of the corresponding characteristic polynomial. Notice that this class of eigenvalue problems is very promising to the solution of linear fractional partial differential equations (FPDE). The numerical results demonstrate reliability and efficiency of the proposed algorithm. Based on our simulations some theoretical conjectures are reported.