The numerical solution of Cauchy singular integral equations with additional fixed singularities

In this paper we propose a quadrature method for the numerical solution of Cauchy singular integral equations with additional fixed singularities. The unknown function is approximated by a weighted polynomial which is the solution of a finite dimensional equation obtained discretizing the involved integral operators by means of a Gauss-Jacobi quadrature rule. Stability and convergence results for the proposed procedure are proved. Moreover, we prove that the linear systems one has to solve, in order to determine the unknown coefficients of the approximate solutions, are well conditioned. The efficiency of the proposed method is shown through some numerical examples.