A quadrature method for Cauchy singular integral equations with index \(-1\) (Q2902203)

A new quadrature method for first kind Cauchy singular integral equations with index \(-1\) on the interval (\(-1,1\)) is proposed. The method uses a Lagrange type polynomial interpolation operator based on the zeros of orthonormal polynomials with respect to the Jacobi weights. The integral equation is approached in a pair of Zygmund-type weighted spaces and the existence and uniqueness of the solution is proved. The polynomial approximation of the solution is constructed by solving a determined and well-conditioned linear system such that its condition number tends to the condition number associated to the linear bounded operator that describes the integral equation in operatorial form. The convergence of the method is proved by providing an error estimate. The theoretical results are tested and confirmed on five numerical examples illustrating the accuracy and the evolution of the condition numbers.