We address the problem of existence, uniqueness and approximation of the solution of a large class of fractional integro-differential equations with boundary conditions. Based on fixed point techniques with Boyd-Wong’s type conditions, the existence and uniqueness of the solution is studied. This, together with the use of certain types of biorthogonal systems, allows us to propose a method to approximate the solution that is tested with several examples.