A Numerical method for solving a class of fractional Sturm-Liouville eigenvalue problems