In this paper an expansion method, based on twodimensional delta functions (2D-DFs), is developed to nd numerical solutions of two-dimensional linear Fredholm integral equations. The main characteristic behind this method is that this method reduce such problems to a system of algebraic equations. Since this approach does not need integration, all calculations can be easily implemented. Finally, we estimate the error of the method, and present two numerical examples to demonstrate the accuracy of the method.
