Abstract
This paper presents a new pseudospectral technique to solve three dimensional integral equations (3D-IEs). The shifted Legendre Gauss-Lobatto collocation method is investigated to approximate the 3D-IEs. The main idea in the novel algorithm is to reduce the 3D-IEs to systems of algebraic equations. The applicability of the present method is examined by several numerical examples. By means of these numerical examples, we ensure that the present technique is a simple and very accurate numerical scheme for solving 3D-IEs.