Algorithms for approximating finite Hilbert transform with end-point singularities and its derivatives (Q645708)

Algorithms are proposed for the numerical evaluation of Cauchy principal value integrals \[ \int^1_{-1} \{\omega(t) f(t)/(t- x)\}\,dt,\quad -1< x< 1, \] with a weight function with Jacobi type singularities \(\omega(t)= (1-t)^\alpha (1+t)^\beta\), where \(\alpha=\pm 1/2\), for a given function \(f(t)\) and Hadamard finite-part integrals \[ \int^1_{-1} \{\omega(t) f(t)/(t- x)^2\}\,dt. \] The function \(f\) is interpolated by using a finite sum of Chebyshev polynomials. The scheme used in this paper gives uniform approximations, namely the errors are bounded independently of \(x\), and is quite efficient for smooth \(f\). Further the approximation of hyper-singular integrals \[ \int^1_{-1} \{\omega(t) f(t)/(t- x)^n\}\,dt,\quad n\geq 3 \] and uniform convergence are discussed. Numerical examples are also given to demonstrate the performance of the present scheme.