In this work we construct numerical schemes, using radial basis functions (RBFs), to solve Sturm--Liouville eigenvalue problems in which the so-called fractional derivatives are involved. We take into account the fact that the solutions to such problems might be nonsmooth and, consequently, we augment the space generated by RBF base functions with some fractional polynomials to overcome this problem. We show through some examples the effectiveness of our method.

CEMAT - Center for Computational and Stochastic Mathematics