Product quasi-interpolation method for weakly singular integral equation eigenvalue problem (Q554580)

The numerical solution of the integral equation eigenvalue problem \[ \lambda\phi(x)=\int\limits_{-1}^{1}K(x,y)\phi(y)dy, \] for an eigenvalue \(\lambda\) and a corresponding eigenfunction \(\phi(x)\) is considered. The case when \(K(x,y)\) is a weakly singular kernel is discussed. \(K(x,y)\) has the form \[ K(x,y)=|x-y|^{-\alpha},\quad 0<\alpha<1, \quad \text{or} \quad K(x,y)=\log|x-y|. \] A discrete method of accuracy \(O(h^m)\) for this problem is constructed. The method is based first on improving the boundary behaviour of the exact solution \(\phi(x)\) with the help of a change of variables, and second on the product integration method based on a discrete quasi-interpolant of order \(m\). The error analysis associated to the suggested method is given, and several examples showing the numerical efficiency of such approach are presented.