In this paper, three numerical methods, namely the Legendre wavelet method, Taylor series expansion method and the improved block-Pulse functions method are implemented to solve linear systems of Fredholm integral equations. To test the validity of these methods, two numerical examples with known exact solution are presented. Numerical results show that the convergence and accuracy of these methods are in good agreement with analytical solution. However, according to comparison of these methods, we conclude that the Legendre wavelet method provides more accurate results.