Rational Chebyshev bases and Galerkin method are used to obtain the approximate solution of a
system of high-order integro-differential equations on the interval [0,∞). This method is based on
replacement of the unknown functions by their truncated series of rational Chebyshev expansion.
Test examples are considered to show the high accuracy, simplicity, and efficiency of this method.
