A new method based on the Clenshaw-Curtis quadrature for the numerical solution of the integro-differential Schrodinger equation is investigated. The method shows that it converges quickly and the truncation errors decrease faster than any power of the inverse number of the Chebyshev support points. Discretization technique is presented in detail. Accompanying C^+^+ code for the Fredholm type method is available upon request.