Variational and Numerical Approximations for Higher Order Fractional Sturm-Liouville Problems