Regular, Singular and Hypersingular Integrals over Fractal Contours

: The paper is devoted to the approximate calculation of Riemann deﬁnite integrals, singular and hypersingular integrals over closed and open non-rectiﬁable curves and fractals. The conditions of existence for the Riemann deﬁnite integrals over non-rectiﬁable curves and fractals are provided. We give a deﬁnition of a singular integral over non-rectiﬁable curves and fractals which generalizes the known one. We deﬁne hypersingular integrals over non-rectiﬁable curves and fractals. We construct quadratures for the calculation of Riemann deﬁnite integrals, singular and hypersingular integrals over non-rectiﬁable curves and fractals and the corresponding error estimates for various classes of functions. Singular and hypersingular integrals are deﬁned up to an additive constant (or a combination of constants) that are subject to a convention that depends on the actual problem being solved. We illustrate our theoretical results with numerical examples for Riemann deﬁnite integrals, singular integrals and


Introduction 1.Review of Approximate Methods for Calculating Singular and Hypersingular Integrals
Starting the middle of the last century, the methods of singular integral equations (SIE) and then hypersingular integral equations (HIE) have been increasingly used in the study and modeling of various problems in physics, natural science and technology: in aerodynamics, electrodynamics, elasticity theory, nuclear and atomic physics, geophysics, and mathematical physics.
Analytical methods for solving singular and hypersingular integral equations are known only for very special cases (see [1] for singular integral equations and [2], Ref. [3] for hypersingular integral equations).Thus, numerical methods are broadly employed to solve singular and hypersingular integral equations.
The development of approximate methods for solving SIE started in the 1950s.The number of publications devoted to approximate methods for solving SIE, their generalizations and related Riemann and Hilbert boundary problems has not diminished since.The main approximate methods for solving SIE are presented in [4,5], where one can find an extensive bibliography.The main approximate methods for solving HIE can be found in the publications [6][7][8][9][10].
The major basic component of any approach to solving SIE and HIE is an efficient approximate method for evaluating the corresponding singular (SI) and hypersingular integrals (HI).It substantiates the need for development of numerical methods for evaluating SI and HI.
There are numerous publications devoted to numerical evaluation of SI and HI over smooth curves (see [4,8,[11][12][13] and the literature therein).The numerical approach to the evaluation of SI and HI over fractals is different.Even though there are few papers devoted to solutions of SIE and HIE on the prefractals of the Cantor perfect set and the Sierpinski "carpet" [14,15], and the research on existence and uniqueness of SIE solutions on fractals is being actively developed [16], the works devoted to generic numerical methods for calculation of SI and HI and solution of SIE and HIE over fractals are unknown to the authors.
The demand for approximate calculations of Riemann definite integrals, singular and hypersingular integrals over non-rectifiable curves and fractals is determined by numerous applications in physics and engineering that are based on integral equations over fractals.In particular, modeling of fractal systems in electrodynamics, ultrahigh frequency technology, antennas [17][18][19][20], in solid state physics and geophysics [21,22] is an active and growing field.Further progress in this field is hampered by the lack of numerical methods for such problems.This work should make up for this deficiency.We pay special attention to numerical calculations of singular and hypersingular integrals over non-rectifiable curves and fractals and their peculiarities.
This paper is devoted to approximate methods for the calculation of integralsincluding singular and hypersingular integrals-on non-rectifiable curves and fractals.We recall the major definitions and summarize some known facts concerning the integrals over fractals.Next, we shall discuss a few general cases of such integrals, introduce approximate quadrature formulas and give their error estimates.Then, we introduce the cases of singular and hypersingular integrals.Finally, we illustrate the results with a few simple numerical examples.

Definitions
Let L be a contour on the complex plane.Let A = [a, b] or A = L. Definition 1. Class of Hölder functions H α (M; A), (0 < α ≤ 1) consists of functions f (x) defined on A and satisfying at all points x and x of this set the inequality | f (x ) Definition 2. The class W r (M; A) consists of functions defined on A, continuous and having continuous derivatives up to (r − 1)-th order inclusive and piece-wise continuous r-th order derivative satisfying on this set the inequality | f (r) (x)| ≤ M, r ≥ 1. Definition 3. The class W r H α (M; A) consists of functions f (x) belonging to the class W r (M; A) and satisfying an additional condition f (r) 3, . . ., in the sense of Cauchy-Hadamard principal value is defined as the limit: where ξ(v) is a function satisfying the following conditions: (1) the limit exists; (2) ξ(v) has a continuous p − 1 degree derivative at a neighborhood of zero.
Throughout the paper, we shall use the Whitney extension of continuous functions.For the reader's convenience, we provide the corresponding statement.
Let F ⊂ R n be a closed set, H(α, F) is a space of functions defined on F and satisfying the Hölder condition with an exponent α.The following statement holds.Lemma 1 ([24]).The linear Whitney extension operator E 0 maps the space H(α, F) to the space H(α, R n ) continuously.The norm of the mapping is bounded by a constant that does not depend on the set F.

Definitions of Regular Integrals
Consider the integral where Γ is a non-rectifiable curve on a complex plane.In the literature, there are several definitions of integrals over non-rectifiable curves and fractals.Here are some of them.
We will start with an informal definition of the integral (1).Let g(x) be a one-to-one and continuous mapping of the interval [0, 1] onto the curve Γ.Then, the evaluation of the integral (1) reduces to the evaluation of the integral 1 0 f (g(x))dg(x). (2) First, we have to solve the issue of the existence of the integral (2).V. Kondurar [25] obtained the following result.

Statement 1. The Stieltjes contour integral
Kats [26] notes that the Kondurar class {H α , H β } is not geometrically invariant since conditions for existence of the integral depend on the choice of parametrization z = z(x).The generalization of the Statement 1 for functions of bounded Φ variation is also provided in [26].Definition 5. Let Φ(x) be a continuous increasing function for x ≥ 0, Φ(0) = 0, f (z) be the function on the curve Γ.The quantity is called the Φ variation of the function f .The supremum is taken with respect to all possible partitions τ = {z Let V Φ (Γ) stand for the class of functions of bounded Φ variation, i.e., all functions for which Φ variation is finite.Moreover CV Φ (Γ) is the subclass of V Φ (Γ) which consists of continuous functions.
A number of statements about the existence of the integral (1) has been provided in [26].Here, we recall the generalization of Statement 1.
The definition of integral over non-rectifiable curves on the complex plane called geometric approximation by Kats [26] is as follows.If a function is given in the neighborhood of Γ, for large enough n, the curve can be approximated by a polygon Γ n with n edges and we shall assume This approach has been investigated in [30][31][32].
Kats [33] presented the definition of the integral over an arc Γ Γ f (z)dz based on approximation of f (z) with algebraic polynomials {p n (z)} : where P n (z) is anti-derivative of the polynomial p n (z).
It has been shown that the limit exists if (1) The sequence p n (z) converges to f (z) in the metric of the Hölder space H α , 0 < α ≤ 1; (2) α > d(Γ) − 1, d(Γ)-is the upper metric dimension of the curve Γ; (3) The arc Γ does not twist into a spiral at the ends.
In [34,35], the integral Γ f (z)dz is defined as a generalized function with support Γ.
Let us assume that Γ is non-rectifiable or fractal contour in complex plane and D + is the domain bounded by Γ.
Recall [36] that the operator ∂z is defined by the formula By Formula (4), evaluation of the integral (3) over the closed curve Γ reduces to evaluation of the integral in the right-hand side of (4) over region D + bounded with the curve Γ.
The extension of function f (z), z ∈ Γ is carried out with the Whitney extension operator [24].When f ∈ H α , the Whitney operator extends each continuous function f (z) on Γ to a continuous function f (z, z) on the complex plane C. Moreover Remark 1.To evaluate integral of the function f ∈ C 1 (D + ) over a smooth curve Γ, Vekua [36] proposed the formula By analogy, with Formula (4), this formula can be used to define integrals on non-rectifiable curves and fractals where f -the Whitney extension [24], p. 205.This representation seems to be more convenient to construct quadrature formulas.

Definitions of Singular and Hypersingular Integrals
Consider the integral Γ f (z)dz, assuming that the function f has a pole at some point on the contour Γ.In case of such a pole we can represent f as a product of regular and singular functions f = f 0 v, here f 0 ∈ H α (Γ), f 0 (t) = 0, and Here, D + is the domain bounded by the contour Γ.
Now, we can give a definition for a singular integral.
Consider the singular integral This leads us to the following statement.

Definition 8 ([37]
).A singular integral S Γ f over a closed non-rectifiable curve is defined by , where f (z) is a Whitney extension for f .
Using Equation ( 5) and repeating the arguments by Mironova [37], we have the following definition of singular integrals.Definition 9. A singular integral S Γ f over a closed non-rectifiable curve is defined by , where f (z) is the Whitney extension for f .
Consider the hypersingular integral 1 πi Γ f (τ) (τ−t) p dτ, t ∈ Γ, p = 2, 3, . . . .In the case of a smooth closed curve Γ, a hypersingular integral on the complex plane C is defined by Using the last formula, the definition follows.Definition 10.A hypersingular integral over a non-rectifiable curve Γ is defined by where g(z)− is a Whitney extension for g(τ, t).
Using Formula (5), we have the following definition of hypersingular integrals.

Definition 11.
A hypersingular integral over a non-rectifiable curve Γ is defined by where g(z)− is a Whitney extension for g(τ, t).
Singular and hypersingular integrals with a singularity of the order higher than the one given above can be considered as distributions (generalized functions) (see [38]).Within this approach, the singular and hypersingular integrals can be defined up to an additive constant, which depends on the particular method of regularization.Definitions of the singular and hypersingular integrals ( 8)-( 11) also, in fact, rely upon some particular methods of fractal boundary regularization.For instance, for t ∈ Γ we assign the integral 1 πi Γ f (τ) τ−t dτ the value of f (t), as if Γ would have been a smooth contour.Here, we assume that f (t) is an analytic function in D + .There are, however, other methods of contour regularization that might be more appropriate when solving some particular physical or engineering problems.Here, we shall discuss one more approach to such alternative regularization.
Assume that Γ is a piece-wise smooth continuous contour and t is one of its "corner" points.Then, the integral 1 where α is the angle between left and right tangents to Γ at point t.If Γ is a fractal contour and we approximate it with some sequence of piece-wise smooth prefractals Γ n , and starting from some n > n 0 , t ∈ Γ n , we can define the angle between the right and the left tangents to the prefractal α n (t) on each of the prefractals.Now, we can distinguish the two following cases: (a) ∃ α ∞ (t) = lim n→∞ α n (t); (b) there is no such limit.In the first case, we can introduce the following definition of the singular integral: Definition 12.A singular integral S Γ f over a fractal contour Γ approximated by a sequence of prefractals Γ n such that t ∈ Γ n is defined by where f (z) is a Whitney extension for f , provided the limit exists, where α n (t) is the angle between the right and the left tangents to the prefractal Γ n at the point t.
This definition emphasizes that the singular integral evaluation depends on the choice of the prefractal sequence that we use to approximate the idealized fractal contour.It also exhibits the fact that there is no unique value that could be assigned as a value of singular integral over the fractal curve, and only in some exceptional cases such a unique value-the one given by Definition 12-can be found and substantiated under certain assumptions.

Remark 2.
A similar situation takes place when defining hypersingular integrals.

Approximate Calculation of the Stieltjes Integral
In this section, we construct quadrature formulas to calculate the Stieltjes integral over non-rectifiable curves and fractals.
The most natural way is to calculate the Stieltjes integral (1) using known quadrature formulas, which is applied for the integral (2).Using (2), however, requires an explicit form of the function z = g(x).Constructing it in general form is very complicated, if not impossible.Therefore, using (2) in approximate calculations is difficult.So, the more natural way is to employ the Definition 6-or an equivalent Equation ( 5)-for complex functions and the corresponding vector fields.
In this section, we present the method for calculation of integrals of continuous functions over non-rectifiable curves and fractals under the condition that the integral over non-rectifiable curve is understood in the sense of Definition 6.
We shall assume-without loss of generality-that the curve Γ is a closed contour.If the curve Γ does not form a closed contour, we can connect its start and end points with a smooth curve and perform two separate calculations for a closed partially fractal contour and its smooth part.
For example, let us assume that a non-rectifiable curve Γ with the start point B and the end point C lies in the plane of the complex variable z, in the region Ω = [−1, 1] 2 .Let f (z) be an analytic function defined on Ω. Connect the points B and C with polygon γ lying in D. Any piece-wise smooth function suitable for quadrature formula construction can be selected. Let Construct the cubature formula to calculate the integral (4 Here, is the total residual term (the error).
There are two contributions to the error of the method: an error of the rectangle cubature formula over ∆ kl squares that do not have common points with Γ and an error over the boundary Γ R nn ( f ) = q nn ( f ) + r nn ( f ) .
The contribution r n,n ( f ) to the residue corresponds to changing the domain of integration from l=0 * stands for a union of the tiles we use in the cubature rule.
Obviously, the first component q nn ( f ) does not exceed A/n α (the error of the rectangle cubature formula for double integral, and α is the Hölder class exponent).
To estimate the second term, it is necessary to estimate the number of squares ∆ kl , k, l = 0, 1, . . ., n − 1 covering the curve Γ.In doing so, we employ fractal measure (Minkowski's measure).Definition 13 ([39]).Let Ω be a compact set on the plane.Divide the plane onto squares with sides ε > 0. Let N ε (Ω) stand for the number of squares that have intersections with Ω.The quantity is called the upper metric dimension of the set Ω.
It follows from the definition that for the cubature formula in question Γ) .Here, ε = 1 n .The maximum error of the rectangular quadrature rule in the square ∆ kl does not exceed M/n 2+α .If the square ∆ kl intersects the region D+ and is not included in the number of the squares to be summed in (6), then its lack in the cubature Formula ( 6) introduces the error of the order M/n 2 .Thus, the quantity of the second component of the error of the cubature Formula (6) does not exceed CMn d(Γ) /n 2 = CM/n 2−d(Γ) , where d(Γ) is the metric dimension of the fractal Γ.
Therefore the error of ( 6) is ).So, we have come up with the following theorem.
), where d(Γ) is metric dimension of the set Γ, and C and M are independent positive constants.Now, let us look at the case of a smoother function f , where, for instance, we assume that f (z) ∈ W r+1 H α (M, D+ ).We will use the same notation.Let the region D+ lie inside a square ∂z will be approximated with interpolation polynomial P r,r ( ∂ f (z) ∂z , ∆ k,l ) constructed with respect to r + 1 equidistant points.
We will calculate the integral Γ f (z)dz by the quadrature formula Here, is determined the same way as in Equation ( 6).
Employing the estimates of the constructive theory of functions [40], we have ).
Thus, we have proved the following assertion. Here, l=0 * * means summation over the squares ∆ k,l for which the measure of intersection with ∆ i,j is zero and their centers (v k , v l ) belong to Ω.In this case, the error of the cubature is Now, we shall assume that the function f (z) is not defined in D + .In this case, to construct cubature formulas we employ the Whitney extension of the function f (z) onto region D + .Finally, we obtain the following cubature formulas. Let Then, the cubature formula reads Knowing that ).We now consider another approach to calculating integrals over non-rectifiable curves and fractals by calculating the contour integral directly.
Let us calculate the integral Γ f (z)dz.On contour Γ, we choose n equidistant points Let us estimate the error: Parameters h and n are related as nh d(Γ) = l, where d(Γ) is the fractal dimensionality of the contour Γ and l is its length in fractal measure.
Finally, we have Remark 3. As the nodes of the quadrature we can choose any points of the complex plane that deviate from Γ by no more than h, it does not affect the error estimate.

Approximate Calculation of Singular Integrals
In this section, we shall use the same notation.Let f (z) We shall construct cubature formulas to calculate singular integrals on the class of functions W 1 H α (M, D + ) based on the following formula Let the region D+ be within the square Construct the cubature formula

∑
l=0 * means summation over indexes (k, l) so that the nodes ūkl , k, l = 0, 1, . . ., n − 1 are within the region D+ , |R nn ( f )| is the error of the cubature Formula (11), and r nn ( f ) has the same meaning as in Equation (6).We call the squares ∆ kl marked if they are entirely in the region D+ , (∆ kl \∂∆ kl ) ∩ Γ = ∅ .Let us estimate the error of the cubature Formula (11), which consists of three components: first, the accuracy of calculation of the integrals defined in marked squares; second, in the cubature Formula (11), we do not count the squares, whose centers are outside the region D+ ; third, when constructing Formula (11), we disregard the fact that if the distance between the center of the square ∆ kl and the boundary of the region D+ is less than A/n, then ∆ kl ∩ D+ = ∆ kl .We estimate each of the components.Let t ∈ ∆ kl . Obviously, Here and below, C stands for positive constants that do not depend on n.
Let us estimate the second contribution to the error of cubature (11).This estimate essentially depends on the topology of the fractal Γ.
During the cubature construction, the region D + has been covered by squares with the side h = A/n.Let us enumerate the squares crossed by Γ counterclockwise starting from the one with the singularity point t.We shall call these squares marked.
Consider two limiting cases.First, we shall assume that the fractal is covered by a rather coarse grid or that its structure approaches the structure of a smooth curve.In this case, the number of tiles crossed by the contour in the vicinity of the singularity is limited to just the two nearest tiles.In the second case, the grid is dense enough to cover the finest variations in the fractal curve and those variations scale, so most of the tiles that neighbor the singularity are crossed by the contour.
In the first case, the distance from the squares K s in the vicinity of the singularity point t and the singularity can be estimated as where α > 0, K s is the s-th marked square and.Then, if the m-th marked square is missed in the summation, the lack of its contribution leads to the error The squares more distant from the singularity contribute O(n −2 ) to the error.So, the second component of the error can be estimated as where m is the number of marked squares.Let L be the Hausdorff measure of the fractal Γ, then mh d(Γ) = L, and we can conclude In the second case, as before, we start marking the squares from K 0 which contains the singularity point.K 0 is surrounded by eight squares ∆ kl , the next layer contains 16 squares, and, generally, the i-th layer from K 0 contains 8i tiles.Now, we shall assume that all the squares around K 0 up to the r-th layer are marked.As they do not contribute to the sum, they introduce the error So, we can conclude that each layer introduces the error O(1/n), and we have to estimate the number of such "boundary" layers r.For this purpose, let us compare the number of marked cells m = 4r(r + 1) with characteristic length L = m( 1 n ) d(Γ) of the fractal.Obviously, r ∼ m 1/2 = (Ln d(Γ) ) 1/2 .Thus, the second contribution to the error has the order of Combining the obtained estimates leads us to the following assertion.
The cubature Formula (11) has the error Here, d(Γ) is upper metric dimension of the set Γ, C 1 and C 2 are independent positive constants.Now, consider the case of a smoother function.Let f ∈ W r+1 H α (M, D + ).Let the region D+ be within the square Ω * = [−A, A] 2 .We will use the foregoing notation.We will approximate function ∂z in each square ∆ k,l , k, l = 0, 1, . . ., n − 1 by the interpolating polynomial P r,r ( ∂ f (z) ∂z , ∆ k,l ) constructed with respect to r + 1 equidistant nodes.Let us construct the cubature formula Repeating the foregoing arguments (given standard estimates of constructive function theory [40]) leads us to the following assertion.
The cubature Formula ( 13) has the error Here, d(Γ) is upper metric dimension of the set Γ, C 1 and C 2 are independent positive constants.

Approximate Calculation of Hypersingular Integrals
In this section, we study the calculations of hypersingular integrals Γ f (τ) (τ − t) p dτ , p = 2, 3, . . .over non-rectifiable curves and fractals.We use the same notation as in previous sections.
Let f (z) ∈ W p+r−1 H α (M, D + ).We shall construct cubature formulas for the calculation of hypersingular integrals on the classes of functions W p+r−1 H α (M, D + ) based on the following formula (see Definition 10) Let the region D+ be within the square Here, we use the same notation for the summation symbols as before.
Repeating the foregoing arguments (given standard estimates of constructive function theory [40]) leads us to the following assertion.
The cubature Formula (15) has the error Here, d(Γ) is upper metric dimension of the set Γ, and C 1 and C 2 are independent positive constants.

Numerical Illustrations
In this section, we give a few examples of numerical calculations for integrals over fractals.As test fractals, we choose the Koch snowflake and its analog based on scaled squares (see Figure 1a,b).We shall refer them to as fractals "a" and "b".The corresponding fractal dimensions for the fractals "a As integrands, we use five functions.For examples of regular functions, we use As singular function integrands, we employ Finally, for an example of hypersingular integral, we use the following integrand We perform two versions of calculations for the corresponding contour integrals: the first is based on Formula (7) and its analog for the singular case (11), and the second corresponds to a direct calculation of the integral over the corresponding prefractals with mid-point or trapezoidal rules.The results of calculations for the integral Γ f 1,2 (z)dz for regular functions are given in Table 1.As the function f 1 is linear in both arguments, the midpoint rule gives an exact result for both approaches.The numerical results for a given order of the prefractal match exactly.The results for function f 2 obtained by direct calculation and by Formula (7) also agree well.This demonstrates the correctness of the integral over fractals definitions that we have been using.
The case of singular functions is less trivial.Even an elementary direct calculation of an integral over a contour meets some unexpected complications.In the example above we chose the step of the quadrature as the length of the corresponding prefractal side.This choice does not guarantee convergence to a correct result in some contour geometries.Consider a sequence of contours composed of the prefractals for the fractal Figure 1 shifted so that the lowest side of each prefractal is centered at the singularity point z = 0 + 0i (Figure 2).The sides of the prefractal scale are 1/3 k , where k is the order of the prefractal.In order to calculate the singular integral over the contour, we employ the trapezoidal rule on each side of the prefractal.Even though the integration step goes to zero, any fixed order quadrature formula has a constant error contribution in the vicinity of the singular point.In order to compensate for this effect, we either have to use a sufficiently high order quadrature or further subdivide the sides of the prefractal.This effect has little to do with the fractal-like structure of the contour and can be also observed with piece-wise smooth contours of certain geometries.We attract the readers' attention to this peculiarity only as a warning that the experience in numerical calculations of principal value integrals over real ranges does not always translate directly to their complex contour integral counterparts.The results of calculations are summarized in Table 2.
Table 2. Contour integral calculation results for function f 3 (z) = 1/z over the contours shown in Figure 2. The first column indicates the order of the corresponding prefractal, and the second column contains the length of the corresponding prefractal side.The calculated real part of the integral does not exceed 10 −10 and is not indicated.Take note that even though the prefractal side length approaches zero, there is a non-vanishing error in the evaluation of the singular integral.This error goes down, however, with more accurate quadrature formulas being used.Here i stands for the imaginary unit.As we have mentioned above, the singularity contribution πi f (t) when Formula ( 11) is employed strongly depends on the approximation of the fractal and can vary between 0 and 2πi f (t).We have constructed the prefractal contours so that the singularity always strikes at the middle of one of the straight sides of the prefractal, which corresponds to the πi f (t), as if the fractal is approximated by a smooth contour.Here, we emphasize one more time that this choice is a result of a somewhat arbitrary convention that we make when performing singular integral calculations over a fractal.For instance, if we shift the prefractals in Figure 2 so that the singularity at z = 0 hits the vertex of the prefractals rather than a midpoint of the side, the limiting fractal will be exactly the same, as the side length goes to zero, but the value of the integral would change from πi to πi/2.

Order
The results of calculations for f 4 (z, z) are summarized in Table 3.As the function is not analytical, unlike the previous example, the contribution of the double integral in Equation ( 11) is not trivial.Direct calculations are preformed using an eight-point Gauss-Legendre rule, which guarantees eight significant digits in our case.The results of direct calculations of the integral agree well with the calculations performed on the base of Formula (11).Again, the values that we report here are based on the same convention as in the previous example.Another enlightening example of the delicacy of singular integrals over fractals can be given by calculating the integrals of functions f 3 (z) and f 4 (z, z) over fractal "b".Both approaches to the calculation of the integral give identical results provided the singularity contribution in Formula ( 11) is chosen correctly.But an infinitesimal-in the infinite prefractal order limits-shift of the contour to place the singularity at one of the nearest sides of the prefractals changes the result to πi.
The numerical results for γ b f 4 (z, z)dz are given in Table 4. Finally, we give an example of a hypersingular integral calculation.We calculate the integral using Definition 12 The calculations have been performed with parameters a 0 = 1, a 1 = 3 and a 2 = 1.The results of calculations are given in Table 5.If no regularization is applied, the direct calculation of the contour integral at the left-hand side of Equation ( 16) is not feasible as the real part of the integral rapidly diverges.The imaginary part of the directly evaluated integral, however, is finite and evaluates to the values close to πia 1 /3 consistent with the results evaluated from ( 16) and (17).16) over the fractal "b" prefractals (Koch snowflake).The first column indicates the order of the corresponding prefractal, and the second column shows the grid step for the evaluation of the double integral in (17).Here i stands for the imaginary unit.

Discussion and Conclusions
We have discussed definitions of regular, singular and hypersingular contour integrals over non-rectifiable curves and fractals.One of the main observation is in the difference between singular (hypersingular) integrals over piece-wise smooth curves and fractals.Let us emphasize this difference.
Consider an integral where Γ is a closed bounded curve.Suppose that Γ is a smooth curve and t ∈ Γ.The integral Γ f (τ) τ−t dτ is a singular integral, and some regularization is required.The standard regularization method is It is known [1], that for a smooth closed curve dτ is a definite integral.To calculate it, we use standard quadrature formulas.
In some form or other, this scheme can be applied to the construction of quadrature formulas for the calculation of singular integrals.Now, let Γ be a piece-wise smooth curve, and t * be a point where there is no tangent to the given curve.For t = t * , the regularization of the integral Γ where α is an angle between left and right tangents to the curve γ at the point t * .For non-rectifiable curves, the construction described above is not applicable.There is, at least, a countable set of points where the curve Γ has no tangent lines.Moreover, left and right tangent lines may have different angles between them at different points.There are also curves with no tangent line at any point.
To calculate singular integral over a non-rectifiable curve or fractal, we implement the regularization similar to the regularization for a piece-wise smooth curve Γ f (τ) The first integral is calculated by Stokes's formula and the Whitney extension where f (z) is the Whitney extension for f .So, The second integral f (t) Γ dτ τ−t calculation depends on t.If t is a finite decimal or binary and, for a large enough N, it is included in the list of the vertices of the nth prefractals (n ≥ N), then we assume f (t) Γ dτ τ−t = iα N f (t), here α N is an angle between the left and the right tangent lines at the point t of the N-th order prefractal.If the condition is not fulfilled, any value between 0 and 2π can be ascribed to the integral Γ dτ τ−t .Singular and hypersingular integrals are particular cases of generalized functions.In the theory of generalized functions [38,41], there is a known statement that all the regularizations differ by a constant.Within this approach, it is legitimate to defineaccording to Mironova [37]-a singular integral as where f (z) is the Whitney extension for f .In this work, we define a singular integral by using Formula (18).This approach is substantiated by the representation of the singularity point t.If t is presented as an infinite fraction, then, when solving a particular problem, its value has to be approximated by a finite decimal fraction.Therefore, by choosing an approximate representation of t and the prefractal sequence, which approximates the fractal curve, the researcher ascribes the value to the singular integral dτ τ−t according to the problem being solved.(It seems that further generalizations of this construction are also possible.For instance, singular and hypersingular integrals over fractals could be treated as stochastic objects with distributions depending on the fractal curve.This approach, however, is subject to future research).
The stability of quadrature and cubature formulas for one-and multi-dimensional singular integrals has been studied in the monograph [11].Upper bounds of the errors for a number of cubature formulas have been obtained assuming an ε-perturbation of the integrands.Besides the upper bounds, for some cubatures, the expected values for the errors have also been given.These results can be easily transferred to the cubatures discussed in the article Similar arguments are applicable to quadrature formulas for hypersingular integrals.
In this paper, we have constructed quadrature and cubature formulas for the calculation of Riemann, singular and hypersingular integrals over non-rectifiable curves and fractals.Some quadrature and cubature formulas have been constructed based on various definitions of integrals over non-rectifiable curves and fractals.We obtained error estimates on classes of functions having derivatives of the first order satisfied the Hölder condition with α, 0 < α ≤ 1.
The obtained results show that having derivatives greater than the first order does not affect the accuracy of cubature formulas with rectangular grids.To increase the accuracy of cubature formulas, it is necessary to construct cubature formulas with several grids, which account for the boundary layer.A similar problem arises when we construct cubature formulas to calculate integrals over non-rectifiable curves and fractals based on the Whitney extension.This is caused by the feature of Whitney's extension: the extension u(z) of function f (z) defined on the boundary Γ of region D has derivatives of the first order in D\Γ and |grad u(z)| ≤ K(dist(z, Γ)) α−1 , K = const.
The authors intend to construct cubature formulas with variable grids accounting for the boundary layer in future works.

Theorem 2 . 1 ∑ k=0 n− 1 ∑
Let f (z) ∈ W r+1 H α ( D+ ).The cubature Formula(7) has the error |R nn ( f )| ≤ C(n −r−α + Mn −2+d(Γ) ), where d(Γ) is upper metric dimension of the set Γ, C and M are independent positive constants.From the estimate of |R n,n ( f )| given in Theorem 2 we see that the contributions to the error of the quadrature coming from the squares that enter the sum n−l=0 * and that coming from omitting the squares crossed by the region boundary are not equal.It would be just natural to modify the quadrature(7) so that the contributions equalize.For this purpose, we introduce a finer supplementary grid ∆ k,l which allows us to compensate for the slow decay of the second term and which contributes only in the near-boundary region.Cover the region Ω * by squares ∆

1 ∑
l=0 * means summation over indexes (k, l) so that the nodes ūkl , k, l = 0, 1, . . ., n − 1 lie in the region D+ , |R nn ( f )| is the error of the cubature formula when calculating the singular integral on the function f (z).

Figure 1 .
Figure 1.Fractals used in numerical examples: (a) an analogue of the Koch snowflake based on scaled squares; (b) the Koch snowflake.

Figure 2 .
Figure 2. A sequence of prefractals used in singular integral example.

Table 1 .
Convergence table for integrals Γ f i (z)dz of regular functions over fractal "a".Here i stands for the imaginary unit.

Table 3 .
Convergence table for integrals of a singular function f 4 (z, z) over the contours shown in Figure2.The first column indicates the order of the corresponding prefractal, and the second column contains the length of the corresponding prefractal side.

Table 4 .
(11)ergence table for integrals of a singular function f 4 (z, z) over the fractal "b" prefractals (Koch snowflake).The first column indicates the order of the corresponding prefractal, and the second column shows the grid step for the evaluation of the double integral(11).

Table 5 .
Convergence table for integrals of a hypersingular function (