“…where ξ (s) is a given function, λ, μ are constants, and x(s) is an unknown function. The CSIEs have been solved via various numerical techniques such as using orthogonal Legendre polynomial [6], Lagrangian interpolation with Gauss-Jacobi mechanical quadrature [8], spline method [12,13], Galerkin technique [14], collocation method [15][16][17], application of Jacobi polynomials [18], using Chebyshev polynomials of the second kind [19], quadrature formula [20][21][22], reproducing kernel Hilbert space method [23,24], and other schemes [25][26][27]. Recently, several types of operational matrix methods with truncated series have been proposed for solving the integral and integro-differential equations (see [16,28]).…”