The introduced method in this study consists of reducing a system of infinite boundary integro-differential equations (IBI-DE) into a system of algebraic equations, by expanding the unknown functions, as a series in terms of Laguerre polynomials with unknown coefficients. Properties of these polynomials and operational matrix of integration are first presented. Finally, two examples illustrate the simplicity and the effectiveness of the proposed method have been presented.
