In this paper, we focus on obtaining an approximate solution for multi-dimensional Fredholm integral equations of second kind. An expansion method is used for treatment multi-dimensional Fredholm integral equation of second kind. This method reduces multi-dimensional integral equation to a partial differential equation. After constructing boundary conditions, this partial differential equation reduces to algebraic equation that can be solved easily with any of the usual methods. Furthermore some theorems are proved for convergence analysis. Finally, for showing the efficiency of the method we use some numerical examples