Convergence of a numerical method for solving a hypersingular integral equation on a segment with the use of piecewise linear approximations on a nonuniform grid (Q2358661)

The following linear hypersingular integral equation is considered \[ \int\limits_{a}^{b}\frac{g(y)dy}{(x-y)^2}+\int\limits_{a}^{b}B(x,y)g(y)dy=f(x), \quad x\in(a,b). \tag{1} \] Here the right-hand side \(f(x)\) is a Hölder-continuous function on the segment \([a,b]\) and the kernel \(B(x,y)\) can have a polar singularity of the form \[ B(x,y)=\frac{B^*(x,y)}{|x-y|^{\alpha}}, \quad \alpha\in[0,1), \] where \(B^*(x,y)\) is a Hölder-continuous function of two variables on the set \([a,b]\times[a,b]\). The unknown function \(g(x)\) is sought in the class of functions satisfying the following conditions \[ g\in C[a,b]\cap C^1(a,b), \quad g'(x)=\frac{\gamma^*(x)}{\sqrt{(b-x)(x-a)}}, \quad g(a)=g(b)=0, \] where \(\gamma^*\) is a Hölder-continuous function on the segment \([a,b]\). The first integral in equation (1) is hypersingular and is treated in the sense of the Hadamard principal value. The considered equation is solved by a numerical scheme based on piecewise linear approximation of the unknown function. The uniform convergence of the obtained numerical solutions to the exact solution is proved as the grid pitch and the radius of the neighborhood in which the regularization is performed simultaneously tend to zero.