In this paper, the properties of Legendre wavelets are introduced and are utilized to reduce linear Fredholm integro-differential equations system to a set of algebraic equations, rather than making the matrix sparse. In the end, for showing efficiency of the method, we use some numerical examples and compare with other methods.