The main aim of this work is to present a numerical method based on Bernoulli wavelet basis functions for solving a class of fractional order integro-differential equations with weakly singular kernel. To this aim, we introduce the Riemann–Liouville integral of the basis functions and use them to reduce the main problem to a system of algebraic equations. Finally, an example is solved by the proposed scheme and the results show the efficiency and high accuracy of our new technique.
