In the present work, a numerical technique based on the hybrid Bernoulli
functions is proposed for solving a general class of fractional differential equations. The
Riemann-Liouville integral of the basis functions is introduced and used along with the
suitable collocation points to transform the main problem into a system of nonlinear
algebraic equations. Moreover, an example is included to confirm the efficiency and
accuracy of the method.