No label defined (Q4512227)

Let us consider the problem of calculating the electrostatic field in a region \(\Omega\doteq \mathbb{R}^3\setminus \overline S\) if the boundary value of the potential \(g(P), P\in S\) is given. Here \(S\subset \mathbb{R}^3\) is a smooth surface. This problem is reduced to the solution of a first kind integral equation with weak singularity in the kernel: \(\int_{S}K(P,M)\tau(M) dS_{M}=g(P)\), \(P\in S\). Let the surface \(S\) be given by the parametric equations \(x=x(\alpha,\beta), y=y(\alpha,\beta), z=z(\alpha,\beta)\), \((\alpha,\beta)\in [-1,1]\times [-1,1]\), and let the ranges for the parameters \(\alpha\) and \(\beta\) are divided into \(N_{\alpha}\) and \(N_{\beta}\) equal parts, respectively. As a result we obtain the following lattice \( \alpha_{k}\doteq H_{\alpha}k-1+([{N_{\alpha}-k\over N_{\alpha }}]-[{k\over N_{\alpha}}])\delta_{\alpha}\), \(k=0,\dots,N_{\alpha}\); \(\beta_{l}\doteq H_{\beta }l-1+([{N_{ \beta }-l\over N_{\beta }}]-[{l\over N_{\beta }}])\delta_{\beta }\), \(l=0,\dots,N_{\beta }\). Here \(H_{\alpha}\doteq 2/ N_{\alpha}\), \(H_{\beta}\doteq 2/ N_{\beta}\), \([\cdot]\) is the integer part of the expression. Using bicubic basis functions for the approximation of the density \(\tau(\alpha,\beta)\) the authors propose an algorithm for the approximate solution of the given integral equation.