Solution of Integral Equations by a Chebyshev Expansion Method
  A new spectral type method for solving the one dimensional quantum-mechanical Lippmann-Schwinger integral equation in configuration space is described. The radial interval is divided into partitions, not necessarily of equal length. Two independent local solutions of the integral equation are obtained in each interval via Clenshaw-Curtis quadrature in terms of Chebyshev Polynomials. The local solutions are then combined into a global solution by solving a matrix equation for the coefficients. This matrix is sparse and the equation is easily soluble. The method shows excellent numerical stability, as is demonstrated by several numerical examples. 
G. H. Rawitscher, I. Koltracht and R. A. Gonzales (U. Conn)