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.