On the numerical method for solving a hypersingular integral equation with the computation of the solution gradient (Q2434667)

The author considers a numerical scheme for solving the hypersingular integral equation \[ \int\int_\Sigma g(x)K_r(x^0,x)d\sigma_x=f(x^0),\quad x^0\in\Sigma, \] where \(\sigma\) is a closed smooth surface of the class \(C^4\), the boundary of the bounded domain \(\Omega_0\subset \mathbb R^3\). The kernel is given by the formula \[ K_r(x^0,x)=\frac \partial {\partial n(x^0)}.\frac\partial {\partial n(x)}F(x^0-x) \quad x,x^0 \in\mathbb R^3, \,x\neq x^0,\quad F(x)=\frac 1 {4\pi| x|}, \] where \(\frac \partial {\partial n(x)}\) is the derivation along the normal to the surface \(\Sigma\) at the point \(x\in \Sigma\). The numerical method is constructed by the use of quadratic functions of the type of the method of discrete singularities with a regularization for use of an irregular grid. The convergence of the method is proven.