... for the construction of fractional Sturm–Liouville operators. The proposed algorithm is developed using operational integration matrices of generalized shifted Chebyshev polynomials. We introduce a new bound on the coefficients of the shifted. Chebyshev polynomials subsequently employed to establish an upper bound...