A matrix method for Fractional Sturm-Liouville problems on bounded domain