This paper is concerned with a numerical procedure for fractional Volterra integro-differential equations with weakly singular kernels. The fractional derivative is in the Caputo sense. In this study, Bernoulli polynomial of first kind is used and its matrix form is given. Then, the matrix form based on the collocation points is constructed for each term of the problem. Hence, the proposed scheme simplifies the problem to a system of algebraic equations. Error analysis is also investigated. Numerical examples are announced to demonstrate the validity of the method