Using homotopy analysis method to find the eigenvalues of higher order fractional Sturm–Liouville problems

We utilize the homotopy analysis method to find eigenvalues of fractional Sturm–Liouville problems. Inasmuch as very few papers have been devoted to estimating eigenvalues of these kind of problems, this work enjoys a particular significance in many different branches of science. The convergence of the homotopy analysis method is also considered on the high order fractional Sturm–Liouville problem. The numerical results acknowledge the ability of the proposed method. Eigenvalues are computed within a couple of minutes CPU time at core i3, 2.7 GHz PC.