**2. Auxiliary Results**

Thus, the objective is to obtain an integral representation of the density of a strictly stable law with a characteristic function Equation (1). Without losing generality we will further assume everywhere that *λ* = 1. A strictly stable law with a parameter *λ* = 1 is commonly called the standard strictly stable law. An abbreviated notation of the characteristic function is accepted for standard strictly stable laws *g*<sup>ˆ</sup>(*<sup>t</sup>*, *α*, *θ*, 1) ≡ *g*<sup>ˆ</sup>(*<sup>t</sup>*, *α*, *<sup>θ</sup>*), for density *g*(*<sup>x</sup>*, *α*, *θ*, 1) ≡ *g*(*<sup>x</sup>*, *α*, *<sup>θ</sup>*), for the distribution function *<sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*, 1) ≡ *<sup>G</sup>*(*<sup>x</sup>*, *α*, *<sup>θ</sup>*), and random variable *<sup>Y</sup>*(*<sup>α</sup>*, *θ*, 1) ≡ *<sup>Y</sup>*(*<sup>α</sup>*, *<sup>θ</sup>*). Everywhere below, for standard strictly stable laws, we will use this notation. To obtain an integral representation, we use the method of inverting the characteristic function of a stable law for the first time proposed by V. Zolotarev in the work [31] and

described in detail in his monograph [32] (see also § 4.4 in [36]). To prove the main theorem, we will need some auxiliary results.

**Property 1** (Property of inversion)**.** *For any admissible set of values of parameters* (*<sup>α</sup>*, *θ*)

$$Y(\mathfrak{a}, -\theta) \stackrel{d}{=} -Y(\mathfrak{a}, \theta). \tag{6}$$

 *θ*)} .

**Proof.** The proof of this property is simple enough. Applying the definition of the characteristic function in Equation (1) and by making the substitution of a variable *t* → <sup>−</sup>*τ*, we obtain

$$\begin{aligned} \mathbb{E}\exp\left\{itY(a,-\theta)\right\} = \exp\{-|t|^a \exp\{i\frac{\pi}{2}a\theta \operatorname{sign} t\}\} &= \\ \exp\{-|t|^a \exp\{-i\frac{\pi}{2}a\theta \operatorname{sign} \tau\}\} &= \mathbb{E}\exp\left\{-i\tau Y(a,\theta)\right\} \end{aligned}$$

It follows directly from here Equation (6).

In terms of the characteristic function *g*<sup>ˆ</sup>(*<sup>t</sup>*, *α*, *<sup>θ</sup>*), probability density functions *g*(*<sup>x</sup>*, *α*, *θ*) and distribution functions *<sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*) of a strictly stable law the property of inversion is written in the form

$$\mathbf{g}(-t, \mathfrak{a}, \boldsymbol{\theta}) = \mathbf{g}(t, \mathfrak{a}, -\boldsymbol{\theta}), \quad \mathbf{g}(-\mathbf{x}, \mathfrak{a}, \boldsymbol{\theta}) = \mathbf{g}(\mathbf{x}, \mathfrak{a}, -\boldsymbol{\theta}), \quad \mathbf{G}(-\mathbf{x}, \mathfrak{a}, \boldsymbol{\theta}) = 1 - \mathbf{G}(\mathbf{x}, \mathfrak{a}, -\boldsymbol{\theta}).\tag{7}$$

The utility of this property consists in the fact that owing to this property it is sufficient to consider the issue of the density representation *g*(*<sup>x</sup>*, *α*, *θ*) or the distribution function *<sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*) only for *x* 0 or for *θ* 0. For negative values of the argumen<sup>t</sup> *x* or the parameter *θ* expressions can be obtained according to the expressions given earlier.

The following property will be useful further.

**Property 2.** *For any two admissible sets of parameters* (*<sup>α</sup>*, *θ*, *λ*) *and* (*<sup>α</sup>*, *θ*, *<sup>λ</sup>*)*, there is such a unambiguously defined real a* > 0*, that <sup>Y</sup>*(*<sup>α</sup>*, *θ*, *λ*) *d*= *aY*(*<sup>α</sup>*, *θ*, *<sup>λ</sup>*)*. For the characteristic function in Equation (1), the value a is connected with parameters in the following way a* = (*λ*/*<sup>λ</sup>*)1/*<sup>α</sup>.*

This property is a full analog of property 2.1 in [32] (see also § 3.7 in [36]) formulated for strictly stable random variables with the characteristic function in Equation (1). This property is proved in the same way to the one which is performed in [32]. In the particular case that is of interest to us *λ* = 1, we obtain

$$Y(\mathfrak{a}, \theta, \lambda) = \lambda^{1/a} Y(\mathfrak{a}, \theta, 1). \tag{8}$$

We now formulate a lemma which makes it possible to perform the inverse Fourier transform of the characteristic function and obtain the density of a strictly stable law.

**Lemma 1.** *The probability density function g*(*<sup>x</sup>*, *α*, *θ*) *for any admissible values of parameters* (*<sup>α</sup>*, *θ*) *and any x can be obtained with the help of the inversion formulas*

$$\log(\mathbf{x}, a, \theta) = \frac{1}{2\pi} \int\_{-\infty}^{\infty} e^{-it\mathbf{x}} \mathcal{G}(t, a, \theta) dt = \begin{cases} \frac{1}{\pi} \mathfrak{R} \int\_{0}^{\infty} e^{it\mathbf{x}} \mathcal{G}(t, a, -\theta) dt, \\\ \frac{1}{\pi} \mathfrak{R} \int\_{0}^{\infty} e^{-it\mathbf{x}} \mathcal{G}(t, a, \theta) dt. \end{cases} \tag{9}$$

**Proof.** Performing the inverse of the characteristic function *g*<sup>ˆ</sup>(*<sup>t</sup>*, *α*, *θ*) we obtain

$$\begin{split} g(\mathbf{x}, \boldsymbol{\mu}, \theta) = \frac{1}{2\pi} \int\_{-\infty}^{\infty} e^{-it\mathbf{x}} \exp\left\{-|t|^{\mu} \exp\{-ia(\boldsymbol{\pi}/2)\theta \operatorname{sign} t\}\right\} dt &= \\ \frac{1}{2\pi} \int\_{-\infty}^{0} \exp\left\{-it\mathbf{x} - |t|^{\mu} \exp\{-ia(\boldsymbol{\pi}/2)\theta \operatorname{sign} t\}\right\} dt + \\ \frac{1}{2\pi} \int\_{0}^{\infty} \exp\left\{-it\mathbf{x} - |t|^{\mu} \exp\{-ia(\boldsymbol{\pi}/2)\theta \operatorname{sign} t\}\right\} dt &= I\_1 + I\_2. \end{split}$$

Let us consider the integral *I*1. By substituting the integration variable in this integral −*t* → *τ*, we obtain

$$\begin{split} l\_1 = \frac{1}{2\pi} \int\_{-\infty}^{0} \exp\left\{-it\mathbf{x} - |t|^a \exp\{-ia(\pi/2)\theta \text{sign}\,t\}\right\} dt \\ &= \frac{1}{2\pi} \int\_{0}^{\infty} \exp\left\{i\tau \mathbf{x} - |\tau|^a \exp\{ia(\pi/2)\theta \text{sign}\,\tau\}\right\} d\tau. \end{split}$$

Now having calculated the sum *I*1 + *I*2, we will obtain

$$\begin{split} l\_{1} + l\_{2} &= \frac{1}{2\pi} \int\_{0}^{\infty} \exp\left\{i\pi x - |\tau|^{\mu} \exp\{in\frac{\pi}{2}\theta \operatorname{sign}\tau\}\right\} d\tau + \\ &\quad \frac{1}{2\pi} \int\_{0}^{\infty} \exp\left\{-itx - |t|^{\mu} \exp\{-ia\frac{\pi}{2}\theta \operatorname{sign} t\}\right\} dt = \\ &\quad \frac{1}{\pi} \int\_{0}^{\infty} \exp\{-t^{\mu} \cos(\frac{\pi}{2}a\theta)\} \cos(tx - t^{\mu} \sin(\frac{\pi}{2}a\theta)) dt = \\ &\quad \frac{1}{\pi} \mathfrak{R} \int\_{0}^{\infty} \exp\{-t^{\mu} \cos(\frac{\pi}{2}a\theta)\} \left(\cos(tx - t^{\mu} \sin(\frac{\pi}{2}a\theta)) + i\sin(tx - t^{\mu} \sin(\frac{\pi}{2}a\theta))\right) dt = \\ &\quad \frac{1}{\pi} \mathfrak{R} \int\_{0}^{\infty} e^{itx} \exp\left\{-t^{\mu} \exp\{i\frac{\pi}{2}a\theta\}\right\} dt = \frac{1}{\pi} \mathfrak{R} \int\_{0}^{\infty} e^{itx} \mathfrak{J}(t, a, -\theta) dt. \end{split} \tag{10}$$

As a result, we obtained the first formula in Equation (9). In order to obtain the second formula in (9) it is necessary to subtract in the penultimate manipulation the imaginary component in Equation (10).

Later, we need analytic continuation of the characteristic function *g*<sup>ˆ</sup>(*<sup>t</sup>*, *α*, *θ*) in the complex plane *z*. We will carry out this analytic continuation in the complex plane *z* with a semiaxis %*z* = *t* > 0 with a cut along the negative part of the real axis arg *z* = <sup>−</sup>*π*. The resulting analytic continuation of the function *g*<sup>ˆ</sup>(*<sup>t</sup>*, *α*, *θ*) with a half-line *t* > 0 will be designated as *<sup>g</sup>*+(*<sup>z</sup>*, *α*, *<sup>θ</sup>*). Using the characteristic function in Equation (1), we obtain

$$\lg^{+}(z,a,\theta) = \exp\left\{-z^{a}\exp\left\{-i\frac{\pi}{2}a\theta\right\}\right\}.\tag{11}$$

The idea of analytic continuation of the integrand in the formula of inversion in Equation (9) in the complex plane *z* and subsequent calculation of the resulting integral Γ exp{*izx*}*g*+(*<sup>z</sup>*, *α*, *θ*)*dz* underlies the method of inverting a characteristic function developed by Zolotarev V.M. in the work [31]. This integral is calculated due to such a change in the integration contour Γ at which its real part does not change (for more details see [32]). To substantiate the change in the integration contour, we need the following lemma.

**Lemma 2.** *For any arbitrarily small ε* > 0 *of any admissible values of parameters α and θ and any x* 0 *the integral is*

$$I(\mathsf{C}\_{\mathsf{R}}) = \int\_{\mathsf{C}\_{\mathsf{R}}} e^{iz\mathsf{x}} g^{+}(z, \mathsf{a}\_{\mathsf{r}} - \theta) dz \to 0,$$

*if:*


*Here, θ*0 = 2*π* arcsin *x and*

$$q\_{\mathcal{O}}(\theta, \mathbf{x}) = \arctan\left(\frac{\cos(\pi \theta / 2)}{\sin(\pi \theta / 2) - \mathbf{x}}\right). \tag{12}$$

**Proof.** Let us consider the integral *I*(*CR*) = *CR eizx <sup>g</sup>*+(*<sup>z</sup>*, *α*, <sup>−</sup>*<sup>θ</sup>*)*dz*. As an integration contour *CR* we will consider contour lines that represent an arc of a circle of radius *R* which has *ϕ*1 arg *z ϕ*2 or

$$\mathcal{C}\_{\mathbb{R}} = \{ z : \ |z| = \mathbb{R}, \ \mathfrak{q}\_1 \lhd\_{\mathbb{R}} \text{arg } z \lhd\_{\mathbb{Q}} \mathfrak{q}\_2 \}. \tag{13}$$

The task is to determine under what conditions imposed on the contour *CR* and the parameters *α* and *θ* limits lim*R*→<sup>0</sup> *I*(*CR*) = 0 and lim*R*→∞ *I*(*CR*) = 0.

For any contour *CR* the inequality

$$|I(\mathbb{C}\_{\mathbb{R}})| \lesssim \int\_{\mathcal{C}\_{\mathbb{R}}} |e^{iz\chi}g^+(z,\alpha,-\theta)| |dz|\_{\prime\prime}$$

is true. Assuming in this expression *z* = *rei<sup>ϕ</sup>* and taking into account that *CR* is an arc of a circle of radius *R* , we obtain

$$\begin{split} |I(\mathbb{C}\_{\mathbb{R}})| &\leqslant \mathcal{R} \int\_{\varrho\_{1}}^{\varrho\_{2}} |e^{ixR\exp\{i\eta\}}g^{+}(Rx^{i\rho},a,-\theta)| d\rho \\ &= \mathcal{R} \int\_{\varrho\_{1}}^{\varrho\_{2}} \left| \exp\left\{ixR e^{i\rho} - R^{a} \exp\left\{ia\left(\varrho + \pi\theta/2\right)\right\} \right\} \right| d\rho = \int\_{\varrho\_{1}}^{\varrho\_{2}} \mathcal{U}(R,\varrho,a,\theta) d\rho, \end{split} \tag{14}$$

where

$$\mathcal{U}(r,\varphi,a,\theta) = \exp\left\{-r\mathbf{x}\sin\varphi - r^\kappa \cos(a(\varphi + \pi\theta/2)) + \ln r\right\}.\tag{15}$$

(1) Let us consider the behavior of this integral at *R* → ∞ and *α* < 1 and *x* > 0. Let us consider the integrand *<sup>U</sup>*(*<sup>R</sup>*, *ϕ*, *α*, *<sup>θ</sup>*). Assuming that *R* → ∞ and taking into account that *α* < 1 we obtain

$$\begin{aligned} \lim\_{R \to \infty} \mathcal{U}(R, \varrho, a, \theta) &= \lim\_{R \to \infty} \exp \{-R(x \sin \varrho + R^{a-1} \cos(a(\varrho + \pi \theta/2))) + \ln R\} \\ &= \lim\_{R \to \infty} \exp \{-Rx \sin \varrho + \ln R\}. \end{aligned}$$

In view of the fact that *x* > 0 and the fact that the linear function grows faster than ln *R*, we obtain lim*R*→∞ exp{−*Rx* sin *ϕ* + ln *R*} = 0, if 0 < *ϕ* < *π* and any *α* < 1 and −1 *θ* 1. If *ϕ* = 0, then

$$\mathcal{U}(\mathbb{R},0,\alpha,\theta) = \exp\{-\mathbb{R}^{\mathbb{A}}\cos(a\pi\theta/2) + \ln \mathbb{R}\}.$$

Taking into account that 0 < *α* < 1 and −1 *θ* 1, we deduce that −*<sup>π</sup>*/2 < *αθπ*/2 < *π*/2 and, therefore, cos(*αθπ*/2) > 0. Assuming that *R* → ∞ in this expression and in view of the fact that *Rα* grows faster than ln *R*, we deduce lim*R*→∞ *<sup>U</sup>*(*<sup>R</sup>*, 0, *α*, *θ*) = 0. Without loss of generality the case *ϕ* = *π* can be excluded from consideration. Hence, *ϕ*1 = 0, and *ϕ*2 we represent in the form *ϕ*2 = *π* − *ε*, where 0 < *ε π*/6 is an arbitrary fixed number. As a result, we deduce

$$\lim\_{R \to \infty} \mathcal{U}(\mathbb{R}, \mathfrak{q}, \mathfrak{a}, \theta) = 0 \quad \text{if} \quad 0 \lessapprox \mathfrak{q} \prec \pi - \mathfrak{e},\tag{16}$$

and integration contour in Equation (13) takes the form

$$\mathcal{C}\_{\mathbb{R}} = \{ z \colon \left| z \right| = \mathcal{R}\_{\prime} \ 0 \leqslant \arg z \leqslant \pi - \varepsilon \}. \tag{17}$$

Now assuming that *R* → ∞ in Equation (14) and using Equation (16) we deduce

$$\lim\_{R \to \infty} |I(\mathbb{C}\_R)| \lesssim \lim\_{R \to \infty} \int\_0^{\pi - \varepsilon} \mathcal{U}(\mathbb{R}, \, \varphi, a, \theta) d\varphi = 0.$$

From here it follows that in the case *α* < 1 and *x* > 0 the integral *I*(*CR*) = 0, at *R* → ∞ where contour integration has the form (17). The first item of the lemma is proved.

(2) Let us consider the case *α* > 1. As it is known that at 1 < *α* 2 the parameter *θ* can vary within the range −(2/*<sup>α</sup>* − 1) *θ* (2/*α* − <sup>1</sup>). We are interested in the conditions under which the integral in Equation (14) will tend to zero at *R* → ∞. Assuming *R* → ∞ in Equation (14), we deduce

$$\lim\_{R \to \infty} |I(\mathbb{C}\_R)| \lesssim \lim\_{R \to \infty} \int\_{\varrho\_1}^{\varrho\_2} \mathbb{U}(R, \varrho, a, \theta) d\varrho. \tag{18}$$

From here it follows that the behavior of this integral at *R* → ∞ is defined by the behavior *<sup>U</sup>*(*<sup>R</sup>*, *ϕ*, *α*, *<sup>θ</sup>*). Applying Equation (15), we deduce

$$\lim\_{R \to \infty} \mathcal{U}(R, \varphi, a, \theta) = \lim\_{R \to \infty} \exp \{-R^{\mu} (R^{1-\mu} \mathbf{x} \sin \varphi + \cos(a(\varphi + \pi \theta/2))) + \ln R\}.\tag{19}$$

Taking into account that 1 < *α* 2, then 1 − *α* < 0. Consequently, at *R* → ∞ we can ignore the summand *R*1−*<sup>α</sup>* sin *ϕ* in comparison with the second summand in these brackets. As a result, we have lim*R*→∞ *<sup>U</sup>*(*<sup>R</sup>*, *ϕ*, *α*, *θ*) = lim*R*→∞ exp{−*R<sup>α</sup>* cos(*α*(*ϕ* + *πθ*/2)) + ln *<sup>R</sup>*}. Taking into account that *Rα* grows faster than ln *R* we deduce that lim*R*→∞ *<sup>U</sup>*(*<sup>R</sup>*, *ϕ*, *α*, *θ*) = 0 if cos(*α*(*ϕ* + *πθ*/2)) > 0, or −*<sup>π</sup>*/2 + 2*kπ* < *α*(*ϕ* + *πθ*/2) < *π*/2 + 2*k<sup>π</sup>*, *k* = 0, ±1, ±2, ... . In this problem we will be interested in the case *k* = 0.

It should be noted that here there is a strict inequality, that is why the cases *α*(*ϕ* + *πθ*/2) = ±*π*/2 will be considered separately. As a result, we have −*<sup>π</sup>*/2 < *α*(*ϕ* + *πθ*/2) < *π*/2. From here it follows that

$$\lim\_{R \to \infty} \mathcal{U}(R, \varphi, a, \theta) = 0, \quad \text{if} \quad -\frac{\pi}{2a} - \frac{\pi \theta}{2} < \varphi < \frac{\pi}{2a} - \frac{\pi \theta}{2}. \tag{20}$$

Now we consider the cases *α*(*ϕ* + *πθ*/2) = ±*π*/2. With this value of the argumen<sup>t</sup> cos(*α*(*ϕ* + *πθ*/2)) = 0 and, therefore, it is impossible to ignore the summand *R*1−*αx* sin *ϕ* in Equation (19) now. We will assume that *x* > 0. In view of the aforesaid Equation (19) takes the form

$$\lim\_{R \to \infty} \mathcal{U}(R, \varphi, \mathfrak{a}, \theta) = \lim\_{R \to \infty} \exp\{-Rx \sin \varphi + \ln R\}, \quad \text{at} \quad \mathfrak{a}(\varphi + \pi \theta/2) = \pm \pi/2. \tag{21}$$

Since *R* grows faster than ln *R*, then we will take interest in the constraints imposed on *ϕ* at which sin *ϕ* > 0. For the case *α*(*ϕ* + *πθ*/2) = *π*/2 we deduce that *ϕ* = *π*2*α* − *πθ*2 . Now we will take into account that the parameter *θ* can take values in a range −(2/*<sup>α</sup>* − 1) *θ* (2/*α* − <sup>1</sup>). Thus, if *θ* = −(2/*<sup>α</sup>* − <sup>1</sup>), then *ϕ* = 3*π*2*α* − *π*2 . If *θ* = 2/*α* − 1, then *ϕ* = − *π*2*α* + *π*2 . As a result, we deduce

$$-\frac{\pi}{2\alpha} + \frac{\pi}{2} \preccurlyeq \varphi \lessapprox \frac{3\pi}{2\alpha} - \frac{\pi}{2} \tag{22}$$

We will take into account now that 1 < *α* 2. Substituting the values *α* = 1 and *α* = 2 in Equation (22), alternately we deduce 0 < *ϕ* < *π*, if *α* → 1 and *π*/4 *ϕ π*/4, if *α* = 2. It should be pointed out that since in the considered case *α* cannot take the value *α* = 1, then in the corresponding inequality there is a strict inequality. As a result we deduce sin *ϕ* > 0, if *α*(*ϕ* + *πθ*/2) = *π*/2 for any 1 < *α* 2 and any −(2/*<sup>α</sup>* − 1) *θ* 2/*α* − 1. Applying this result in Equation (21), we deduce

$$\lim\_{R \to \infty} \mathcal{U}(R, \varphi, a, \theta) = 0, \quad \text{if} \quad a(\varphi + \pi \theta/2) = \pi/2 \quad \text{and} \quad \mathbf{x} > 0. \tag{23}$$

The case *α*(*ϕ* + *πθ*/2) = −*<sup>π</sup>*/2 is considered in the similar way as the previous case. As a result, we obtain −*π* < *ϕ* < 0 if *α* → 1 and −*<sup>π</sup>*/4 *ϕ* −*<sup>π</sup>*/4, if *α* = 2. From here it follows that sin *ϕ* < 0 if *α*(*ϕ* + *πθ*/2) = −*<sup>π</sup>*/2 for any 1 < *α* 2 and any −(2/*<sup>α</sup>* − 1) *θ* 2/*α* − 1. Making use of this result in Equation (21), we obtain

$$\lim\_{R \to \infty} \mathcal{U}(R, \varphi, \mathfrak{a}, \theta) = \infty, \quad \text{if} \quad \mathfrak{a}(\varphi + \pi \theta/2) = -\pi/2. \tag{24}$$

and, consequently, it is necessary to exclude this case from our consideration. Putting together Equations (20) and (23), and taking account of Equation (24) we obtain lim*R*→∞ *<sup>U</sup>*(*<sup>R</sup>*, *ϕ*, *α*, *θ*) = 0 if − *π* 2*α* − *πθ* 2 < *ϕ π*2*α* − *πθ*2 and *x* > 0. In view of the obtained result, the expression in Equation (18) takes the form

$$\lim\_{R \to \infty} |I(\mathbb{C}\_R)| \ll\_{\varepsilon} \lim\_{R \to \infty} \int\_{-\frac{\pi}{2a} - \frac{\pi\theta}{2} + \varepsilon}^{\frac{\pi}{2a} - \frac{\pi\theta}{2}} \mathcal{U}(R, \varphi, \alpha, \theta) d\rho = 0, \alpha$$

where *ε* > 0 is an arbitrary small number. Thus, in the considered case contour integration Equation (13) takes the form *CR* = *z* : |*z*| = *R*, − *π*2*α* − *πθ*2 + *ε ϕ π*2*α* − *πθ*2 , and for this contour lim*R*→∞ *I*(*CR*) = 0. This proves the second item of the lemma.

(3) Now let us consider the case 0 < *α* 2 and *R* → 0. It is known that for a specified range of values of the parameter *α* the parameter *θ* can take the values |*θ*| min(1, 2/*α* − <sup>1</sup>). Assuming that *R* → 0 in Equation (14) we obtain

$$\lim\_{R \to 0} |I(\mathbb{C}\_R)| \lesssim \lim\_{R \to 0} \int\_{\varphi\_1}^{\varphi\_2} \mathcal{U}(\mathbb{R}\_\prime, \varphi, a, \theta) \, d\varphi \,. \tag{25}$$

From here we can see that the behavior of this integral at *R* → 0 will be determined by the behavior *<sup>U</sup>*(*<sup>R</sup>*, *ϕ*, *α*, *θ*) at *R* → 0. Using Equation (15), we obtain

$$\lim\_{R \to 0} \mathcal{U}(R, \varphi, a, \theta) = \lim\_{R \to 0} \exp \left\{-Rx \sin \varphi - R^a \cos(a(\varphi + \pi \theta/2)) + \ln R\right\} = 0 \tag{26}$$

For any values of *ϕ* and any *x*. Choosing *ϕ*1 = −*π* and *ϕ*2 = *π* a contour *CR* will take the form

$$\mathcal{C}\_{\mathbb{R}} = \{ z : \ |z| = \mathbb{R}\_{\prime} - \pi \leqslant \arg z \leqslant \pi \}. \tag{27}$$

Now using Equations (26) and (27) in Equation (25) we will finally obtain lim*R*→<sup>0</sup> |*I*(*CR*)| lim*R*→<sup>0</sup> *π*−*π <sup>U</sup>*(*<sup>R</sup>*, *ϕ*, *α*, *<sup>θ</sup>*)*dϕ* = 0, for any admissible values of parameters *α* and *θ* and for any *x*. It proves the third item of the lemma.

(4) We will consider the case *x* = 0 and *R* → ∞. In this case, the expression in Equation (15) will take the form

$$\mathcal{U}(\mathbb{R}, \boldsymbol{\varrho}, \boldsymbol{u}, \boldsymbol{\theta}) = \exp\left\{-\mathbb{R}^{\mathbb{R}} \cos(\boldsymbol{u}(\boldsymbol{\varrho} + \boldsymbol{\pi}\boldsymbol{\theta}/2)) + \ln \boldsymbol{\mathcal{R}}\right\}.\tag{28}$$

Thus, the integral in Equation (14) will tend to zero at *R* → ∞ if cos(*α*(*ϕ* + *πθ*/2)) > 0. Consequently, *ϕ* must meet the conditions

$$-\frac{\pi}{2\alpha} - \frac{\pi\theta}{2} < \varphi < \frac{\pi}{2\alpha} - \frac{\pi\theta}{2}.\tag{29}$$

Since no additional limitations for parameters *α* and *θ* were introduced here, then this result is true for any admissible values of these parameters. It should be pointed out that there are strict inequalities here. In fact, if *ϕ* = ± *π*2*α* − *πθ*2 , then cos(*α*(*ϕ* + *πθ*/2)) = 0 and in this case *<sup>U</sup>*(*<sup>R</sup>*, *ϕ*, *α*, *θ*) → ∞ at *R* → ∞. Choosing in Equation (14) the values *ϕ*1 = −*<sup>π</sup>*/(<sup>2</sup>*α*) − *πθ*/2 + *ε* and *ϕ*2 = *π*/(<sup>2</sup>*α*) − *πθ*/2 − *ε* as the limits of integration, contour integration *CR* at *x* = 0 takes the form

$$\mathbb{C}\_{R} = \left\{ z \, : \, |z| = R, \, -\frac{\pi}{2a} - \frac{\pi\theta}{2} + \varepsilon \leqslant \arg z \leqslant \frac{\pi}{2a} - \frac{\pi\theta}{2} - \varepsilon \right\},\tag{30}$$

where *ε* is an arbitrary small positive number. Now using Equations (28) and (30) in Equation (14), we obtain

$$\lim\_{R \to \infty} |I(\mathbb{C}\_R)| \lessapprox\_{\mathbb{R}} \lim\_{R \to \infty} \int\_{-\frac{\pi}{2\pi} - \frac{\pi\theta}{2} + \varepsilon}^{\frac{\pi}{2\pi} - \frac{\pi\theta}{2} - \varepsilon} \mathcal{U}(R, \,\!\!\! \rho, \,\!\!\iota \,\!\! \iota \,\!\! \theta = 0)$$

at *x* = 0 and any admissible values of parameters *α* and *θ*. It proves the fourth item of the lemma.

(5) Let us consider the case *α* = 1. In this case the parameter *θ* can vary within the limits −1 *θ* 1. It is necessary to determine under which conditions the integral in Equation (14) tends to zero at *R* → ∞. As in previous cases, assuming that *R* → ∞ in Equation (14), we obtain

$$\lim\_{R \to \infty} |I(\mathbb{C}\_R)| \lesssim \lim\_{R \to \infty} \int\_{\varrho\_1}^{\varrho\_2} \mathcal{U}(R, \,\varrho, 1, \theta) d\rho. \tag{31}$$

From this it follows that these conditions are a consequence of the behavior *<sup>U</sup>*(*<sup>R</sup>*, *ϕ*, 1, *θ*) at *R* → ∞. Applying Equation (15), we obtain

$$\lim\_{R \to \infty} \mathcal{U}(R, \varphi, 1, \theta) = \lim\_{R \to \infty} \exp\left\{-R(\mathbf{x} \sin \varphi + \cos(\varphi + \pi \theta/2)) + \ln R\right\}.\tag{32}$$

In view of the fact that *R* grows faster than ln *R* at *R* → ∞ we obtain that

$$\exp\left\{-\mathsf{R}\left(x\sin\varphi+\cos(\varphi+\pi\theta/2)\right)+\ln\mathsf{R}\right\}\_{\mathbb{R}\to\infty}0$$

if

$$
\propto \sin \varphi + \cos(\varphi + \pi \theta / 2) > 0. \tag{33}
$$

This inequality allows us to determine the conditions imposed on *ϕ*. It should be noted that one should exclude the case *x* = 1, *θ* = 1 from consideration. In fact, substituting these values in Equation (32), we obtain lim*R*→∞ *<sup>U</sup>*(*<sup>R</sup>*, *ϕ*, 1, 1) = ∞ for any *ϕ*. Thus, from Equation (31) we obtain

$$I(\mathbb{C}\_R) \underset{R \to \infty}{\to} \infty, \quad \text{if} \quad \mathfrak{a} = 1, \quad \mathfrak{x} = 1, \quad \theta = 1. \tag{34}$$

Making some transformations the inequality in Equation (33) can be written in the form

$$\tan \varphi > \frac{\cos(\pi \theta / 2)}{\sin(\pi \theta / 2) - x} \tag{35}$$

From this inequality we can see that it is necessary to consider three cases: *x* > 1, 0 *x* < 1, and *x* = 1.

In the case *x* > 1 for any *θ* the difference sin(*πθ*/2) − *x* is negative. Taking into account that the function arctan *x* is a multi-valued function, then choosing the principal branch arctan *x* we deduce *π*/2 *ϕ* > max(−*<sup>π</sup>*/2, *ϕ*0(*<sup>θ</sup>*, *x*)) or

$$
\varphi\_0(\theta, \mathbf{x}) < \varphi \lhd \pi/2, \quad \mathbf{x} > 1, \quad -1 \leqslant \theta \leqslant 1,\tag{36}
$$

where

$$\varphi\_0(\theta, \mathbf{x}) = \arctan\left(\frac{\cos(\pi \theta / 2)}{\sin(\pi \theta / 2) - \mathbf{x}}\right). \tag{37}$$

It should be noted here that this expression is the solution of the equation

$$
\tan \varphi + \cos(\varphi + \pi \theta/2) = 0. \tag{38}
$$

We need to consider the case 0 *x* < 1. Here we need to consider three possible situations: (1) sin(*πθ*/2) − *x* < 0, if *θ* < *θ*0, (2) sin(*πθ*/2) − *x* = 0, if *θ* = *θ*0, (3) sin(*πθ*/2) − *x* > 0, if *θ* > *θ*0. Here *θ*0 = (2/*π*) arcsin *x*. We will introduce the notation as follows *f*(*ϕ*, *θ*, *x*) = *x* sin *ϕ* + cos(*ϕ* + *πθ*/2). In view of this notation, the condition in Equation (33) will take the form *f*(*ϕ*, *θ*, *x*) > 0. In Figure 1 the function graph *f*(*ϕ*, *θ*, *x*) is plotted for 0 *x* < 1 with different values of the parameter *θ* that correspond to three possible situations: *θ* < *θ*0, *θ* = *θ*0, *θ* > *θ*0.

**Figure 1.** The graph of the function *f*(*ϕ*, *θ*, *x*) at *x* < 1 (the curve is plotted at *x* = 0.7) and different values of the parameter *θ*: *θ*0 = (2/*π*) arcsin *x*, *θ*1 < *θ*0, *θ*2 > *θ*0. Heavy dots demonstrate the value of the solution of the Equation *f*(*ϕ*, *θ*, *x*) = 0 for values *θ* = *θ*0, *θ*1, *θ*2.

Let us consider the case *θ* < *θ*0. In Figure 1, a curve *f*(*ϕ*, *θ*1, *x*) corresponds to this case. From this figure one can see that the condition *f*(*ϕ*, *θ*, *x*) > 0 is met for values *ϕ* > *ϕ*0(*<sup>θ</sup>*, *<sup>x</sup>*). We need to remind that *ϕ*0(*<sup>θ</sup>*, *x*) is the solution of the Equation (38). Thus, choosing the principal branch of the function arctan *x* we ge<sup>t</sup> *π*/2 *ϕ* > max(−*<sup>π</sup>*/2, *ϕ*0(*<sup>θ</sup>*, *x*)) or

$$
\varphi\_0(\theta, \mathbf{x}) \preccurlyeq \varphi \lessapprox \pi/2, \quad 0 \preccurlyeq \mathbf{x} \preccurlyeq 1, \quad \theta \preccurlyeq \theta\_0. \tag{39}
$$

This condition is easy to determine from Figure 2. The values of *ϕ* lying above the curve *ϕ*0(*<sup>θ</sup>*, *x*) and not exceeding the value *π*/2 correspond to the condition *f*(*ϕ*, *θ*, *x*) > 0 at *θ* < *θ*0.

**Figure 2.** The graph of the function *ϕ*0(*<sup>θ</sup>*, *x*) at *x* 1. The figure shows the graphs for values *x* = 0, 0.7, 1. Heavy dots designate the values of this function with the value of the parameter *θ* = *θ*0, *θ*1, *θ*2 and *x* = 0.7 corresponding to Figure 1.

In the case of *x* = 0, we ge<sup>t</sup> *θ*0 = 0 and *ϕ*0(*<sup>θ</sup>*, 0) = arctan(cot(*πθ*/2)). Therefore, the condition *θ* < *θ*0 will take the form *θ* < 0. One should pay attention that the argumen<sup>t</sup> *πθ*/2 at −1 *θ* 1 takes the values in the range from −*<sup>π</sup>*/2 to *π*/2. In this range of values cot(*πθ*/2) in the point *θ* = 0 has the point of discontinuity. We will write *ϕ*0(*<sup>θ</sup>*, 0) in the form

$$\varphi\_0(\theta, 0) = \arctan(1/\tan(\pi \theta / 2)). \tag{40}$$

and will make use of the following trigonometric identities

$$\arctan y = \begin{cases} \operatorname{arccot}(1/y) - \pi, & y \lessapprox 0, \\ \operatorname{arccot}(1/y), & y > 0 \end{cases} \tag{41}$$

and

$$
\arctan y + \operatorname{arccot} y = \pi/2.\tag{42}
$$

Now using these two trigonometric identities in Equation (40), we find

$$
\varphi\_0(\theta, \mathbf{x}) = -\frac{\pi}{2} - \frac{\pi \theta}{2}, \quad \theta < 0 \tag{43}
$$

Thus, at *x* = 0 the condition in Equation (39) will take the form −*π*2 − *πθ*2 < *ϕ π*2 , if *α* = 1 and *θ* < 0. The condition obtained is the same as the condition in Equation (29), if to limit the latter above with the value *π*/2 at *θ* < 0.

Now we need to consider the case *θ* = *θ*0. In Figure 2 we can see that when increasing the parameter *θ* the point *ϕ*0(*<sup>θ</sup>*, *x*) will approach the value −*<sup>π</sup>*/2. With such a change of the parameter *θ* the function graph *f*(*ϕ*, *θ*, *x*) in Figure 1 will shift to the left and with the value *θ* = *θ*0 will take the form that is given in Figure 1. As we can see from Figure 2, in this case the function *ϕ*0(*<sup>θ</sup>*, *x*) has a discontinuity in the point *θ* = *θ*0. It is connected with the fact that the principal branch of the function arctan *x* is investigated. In the vicinity of this point we have lim*θ*→*θ*0−<sup>0</sup> *ϕ*0(*<sup>θ</sup>*, *x*) = lim*θ*→*θ*0−<sup>0</sup> arctan cos(*πθ*/2) sin(*πθ*/2)−*<sup>x</sup>* = −*π*2 and lim*θ*→*θ*0+<sup>0</sup> *ϕ*0(*<sup>θ</sup>*, *x*) = lim*θ*→*θ*0+<sup>0</sup> arctan cos(*πθ*/2) sin(*πθ*/2)−*<sup>x</sup>* = *π*2 . From this we can see that the point *θ*0 is the point of discontinuity of the first kind. It is possible to eliminate the discontinuity of the function *ϕ*0(*<sup>θ</sup>*, *x*) in the point *θ*0 by defining the value of this function in the given point. We will select

$$
\rho\_0(\theta\_0, \mathbf{x}) = -\pi/2, \quad 0 \preccurlyeq \mathbf{x} < 1. \tag{44}
$$

Thus, the condition in Equation (33) is met if

$$-\pi/2 < \mathfrak{q} < \pi/2, \quad 0 \preccurlyeq \ge <1, \quad \theta = \theta\_0. \tag{45}$$

It should be noted that there are precisely strict inequalities here. In fact, if to take *ϕ* = ±*π*/2, then in the case considered these values will be the solution of the Equation (38) which will lead to divergence of the integral in Equation (31).

Now we will consider the case *θ* > *θ*0. From Figure 1 one can see that by increasing the parameter *θ* from *θ*0 to 1 the half-period of the function *f*(*ϕ*, *θ*, *x*) satisfying the condition *f*(*ϕ*, *θ*, *x*) > 0 will keep moving to the left. At the same time, the left point which is the solution of an equation *f*(*ϕ*, *θ*, *x*) = 0 will become smaller than −*<sup>π</sup>*/2. Since we take interest in the interval from −*<sup>π</sup>*/2 to *π*/2, then the left bound of the interval will be the value −*<sup>π</sup>*/2. Keeping in mind that the half-period of the function *f*(*ϕ*, *θ*, *x*) is equal to *π*, the right bound of the interval *f*(*ϕ*, *θ*, *x*) > 0 will be the second solution of the equation *f*(*ϕ*, *θ*, *x*) = 0. As a result, the case considered the condition in Equation (33) will be met with values *ϕ* satisfying the inequality

$$-\pi/2 \lessapprox \mathfrak{q} < \mathfrak{q}\_0(\theta, \mathfrak{x}), \quad 0 \lessapprox \mathfrak{x} < 1, \quad \theta > \theta\_0. \tag{46}$$

The graph of the function *ϕ*0(*<sup>θ</sup>*, *x*) is given in Figure 2.

If *x* = 0, we have *θ*0 = 0 and condition *θ* > *θ*0 turns into a condition *θ* > 0. As it was pointed out earlier, in this case it is convenient to represent *ϕ*0(*<sup>θ</sup>*, *x*) in the form of Equation (40). Now using the trigonometric identities in Equations (41) and (42) in Equation (40) we obtain

$$
\varphi\_0(\theta, \mathbf{x}) = \frac{\pi}{2} - \frac{\pi \theta}{2}, \quad \theta > 0. \tag{47}
$$

Thus, at *x* = 0 the condition in Equation (46) takes the form −*π* 2 *ϕ* < *π* 2 − *πθ* 2 , if *x* = 0 and *θ* > 0. The condition obtained is the same as the condition in Equation (29), if to limit the latter below with a value − *π*/2 at *θ* > 0.

If now we put together the expressions in Equations (43), (44), and (47), then *ϕ*0(*<sup>θ</sup>*, 0) takes the form

$$\varphi\_0(\theta,0) = \begin{cases} -\pi/2 - \pi\theta/2, & \theta \ll 0\\ \pi/2 - \pi\theta/2, & \theta > 0. \end{cases}$$

The graph of the function *ϕ*0(*<sup>θ</sup>*, 0) is given in Figure 2. Consider now the case *x* = 1. In this case the inequality in Equation (35) has the form

$$
\tan\phi (1 - \sin(\pi\theta/2)) > -\cos(\pi\theta/2). \tag{48}
$$

From this it follows that 1 − sin(*πθ*/2) > 0, for any *θ* < *θ*0, where *θ*0 = (2/*π*) arcsin 1 = 1. As a result, selecting the principal branch of the function arctan *y*, the inequality in Equation (48) takes the form

$$
\varphi\_0(\theta, \mathbf{x}) < \varphi \lhd \pi/2, \quad \mathbf{x} = 1, \quad -1 \lhd \theta < 1. \tag{49}
$$

This inequality gives a condition under which the inequality will be satisfied in Equation (35) in case *x* = 1. The graph of the function *ϕ*0(*<sup>θ</sup>*, 1) is given in Figure 2.

Thus, the inequalities in Equations (36), (39), (45), (46), and (49) completely define the condition under which the inequality is satisfied in Equation (33). Combining these conditions for Equation (32), we obtain

$$\lim\_{R\to\infty}\mathcal{U}(R,\varphi,1,\theta)=0,\quad\text{if}\quad\begin{cases}\begin{aligned}&\varrho\_{0}(\theta,\mathbf{x})<\mathsf{q}\leqslant\pi/2, &\mathbf{x}>1,&-1\leqslant\theta\leqslant1,\\ &\varrho\_{0}(\theta,\mathbf{x})<\mathsf{q}\leqslant\pi/2, &\mathbf{0}\leqslant\mathbf{x}<1,&\theta<\theta\_{0},\\ &-\pi/2<\mathsf{q}<\pi/2, &\mathbf{0}\leqslant\mathbf{x}<1,&\theta=\theta\_{0},\\ &-\pi/2\leqslant\mathbf{q}<\varphi\_{0}(\theta,\mathbf{x}), & \mathbf{0}\leqslant\mathbf{x}<1,&\theta>\theta\_{0},\\ &\varphi\_{0}(\theta,\mathbf{x})<\mathsf{q}\leqslant\pi/2, &\mathbf{x}=1,&-1\leqslant\theta<1,\end{aligned}\end{cases}\tag{50}$$

where *θ*0 = (2/*π*) arcsin *x*. It should be noted that if *ϕ* = *ϕ*0(*<sup>θ</sup>*, *<sup>x</sup>*), then in this case *x* sin *ϕ* + cos(*ϕ* + *πθ*/2) = 0 and, as a consequence, lim*R*→ ∞ *<sup>U</sup>*(*<sup>R</sup>*, *ϕ*0(*<sup>θ</sup>*, *<sup>x</sup>*), 1, *θ*) = ∞. Therefore, *ϕ* -= *ϕ*0(*<sup>θ</sup>*, *<sup>x</sup>*). In addition, in terms of Equation (34), one should exclude the case *x* = 1, *θ* = 1 from consideration.

Now using Equation (50) in Equation (31), we ultimately obtain

$$I(\mathbb{C}\_{\mathbb{R}}) \underset{R \to \infty}{\to} 0 \text{ if } \begin{cases} \mathbb{C}\_{\mathbb{R}} = \{z : |z| = R, \ \operatorname{\boldsymbol{\rho}}\theta(\theta, \mathbf{x}) + \varepsilon \leqslant \operatorname{\arg} z \leqslant \pi/2\}, & \mathbf{x} > 1, \ -1 \leqslant \theta \leqslant 1, \\\ \mathbb{C}\_{\mathbb{R}} = \{z : |z| = R, \ \operatorname{\boldsymbol{\rho}}\theta(\theta, \mathbf{x}) + \varepsilon \leqslant \operatorname{\arg} z \leqslant \pi/2\}, & \mathbf{0} \leqslant \mathbf{x} < 1, \ -1 \leqslant \theta < \theta\_{0}, \\\ \mathbb{C}\_{\mathbb{R}} = \{z : |z| = R, \ -\pi/2 + \varepsilon \leqslant \arg z \leqslant \pi/2 - \varepsilon\}, & \mathbf{0} \leqslant \mathbf{x} < 1, \ \theta = \theta\_{0}, \\\ \mathbb{C}\_{\mathbb{R}} = \{z : |z| = R, -\pi/2 \leqslant \arg z \leqslant \rho \boldsymbol{\eta}(\theta, \mathbf{x}) - \varepsilon\}, & \mathbf{0} \leqslant \mathbf{x} < 1, \ \theta \leqslant \mathbf{0} \leqslant 1, \\\ \mathbb{C}\_{\mathbb{R}} = \{z : |z| = R, \ \boldsymbol{\rho}\_{0}(\theta, \mathbf{x}) + \varepsilon \leqslant \arg z \leqslant \pi/2\}, & \mathbf{x} = 1, -1 \leqslant \theta < 1. \end{cases}$$

where *ε* is any arbitrary small positive number. This proves items 5, 6, 7, 8, and 9 of the lemma and proves the lemma completely.

**Remark 1.** *The statements of the Lemma 2 are preserved if to substitute contours CR with their parts.*

Using the lemma which was proved one can substantiate the validity of transition from an integral ∞0 exp{*itx*}*g*<sup>ˆ</sup>(*<sup>t</sup>*, *α*, *θ*)*dt* along the real variable *t* to an integral Γ exp{*izx*}*g*+(*<sup>z</sup>*, *α*, *θ*)*dz* in the complex variable *z* along some contour Γ for inversion formulas in Equation (9). We will state this result in the form of a lemma.

**Lemma 3.** *Let us consider the family of contours* {Γ} *in the complex plane z with a cut along a half-line* arg *z* = −*π satisfying the following conditions:*

*1. Every contour starts in the point z* = 0*.*

−

*2. None of the contours* Γ *intersects the lines of the cut.*

*3. Moving from the point z* = 0 *along the contour* Γ *we let it tend to infinity but in such a way that starting from some place all points z* ∈ Γ *have values of arguments within the limits:*

$$0 \leqslant \arg z \leqslant \pi - \varepsilon, \quad \text{if} \quad 0 < \mathfrak{a} < 1, \ -1 \leqslant \theta \leqslant 1, \ \mathfrak{x} > 0,\tag{51}$$

$$-\frac{\pi}{2a} - \frac{\pi\theta}{2} + \varepsilon \leqslant \arg z \leqslant \frac{\pi}{2a} - \frac{\pi\theta}{2}, \quad \text{if} \quad 1 < a \leqslant 2, \ |\theta| \leqslant (2/a - 1), \ \mathbf{x} > 0,\tag{52}$$

$$\text{or} \quad \pi\theta \leqslant \pi\theta, \quad \text{and} \quad \pi \leqslant \pi\theta \qquad \text{and} \quad \pi\dots\pi \leqslant \pi\dots\pi \quad \dots \dots \dots \quad \dots \dots$$

$$\frac{\pi}{2\mathfrak{a}} - \frac{\pi\mathfrak{o}}{2} + \varepsilon \leqslant \arg z \leqslant \frac{\pi}{2\mathfrak{a}} - \frac{\pi\mathfrak{o}}{2} - \varepsilon, \quad \text{if} \quad 0 < \mathfrak{a} \leqslant 2, \ \mathfrak{x} = 0, \ |\theta| \leqslant \min(1, 2/\mathfrak{a} - 1), \tag{53}$$

$$q\_0(\theta, \mathbf{x}) + \varepsilon \leqslant \arg z \leqslant \frac{\varkappa}{2}, \quad \text{if} \quad \mathfrak{a} = 1, \text{ x} > 1, \text{ } -1 \leqslant \theta \leqslant 1,\tag{54}$$

$$-\frac{\cdots}{2} \lessapprox \arg z \lessapprox q\_0(\theta, \mathbf{x}) - \epsilon, \quad \text{if} \quad a = 1, \ 0 \lessapprox \mathbf{x} < 1, \ \theta\_0 < \theta \lessapprox 1,\tag{55}$$

$$-\frac{\pi}{2} + \varepsilon \leqslant \arg z \leqslant \frac{\pi}{2} - \varepsilon, \quad \text{if} \quad \mathfrak{a} = 1, \ 0 \leqslant \mathfrak{x} < 1, \ \theta = \theta\_0. \tag{56}$$

$$\rho\_0(\theta, \mathbf{x}) + \varepsilon \leqslant \arg z \leqslant \frac{\pi}{2}, \quad \text{if} \quad \mathfrak{a} = 1, \ 0 \leqslant \mathbf{x} < 1, \ -1 \leqslant \theta < \theta\_0. \tag{57}$$

$$
\varphi\_0(\theta, \mathbf{x}) + \varepsilon \lessapprox \arg z \lessapprox \frac{\pi}{2}, \quad \text{if} \quad \mathfrak{a} = 1, \; \mathbf{x} = 1, \; -1 \lessgtr \theta < 1,\tag{58}
$$

*where θ*0 = (2/*π*) arcsin *x and ϕ*0(*<sup>θ</sup>*, *x*) *have the form in Equation (12) and ε* > 0 *is any arbitrary small number. Then for any contour of the specified type and any pair of admissible parameters* (*<sup>α</sup>*, *θ*) *and any x* 0 *(with the exception of the point x* = 1, *α* = 1, *θ* = 1*)*

$$\int\_0^\infty e^{it\mathbf{x}} \mathbf{g}(t, \mathfrak{a}\_\prime - \theta) dt = \int\_\Gamma e^{iz\mathbf{x}} \mathbf{g}^+(z, \mathfrak{a}\_\prime - \theta) d\mathfrak{z},\tag{59}$$

*where t is real.*

**Proof.** From constraints imposed on contours Γ one can see that it is possible to divide the whole family of contours {Γ} into two kinds. The contours which start in the point *z* = 0 and tend to infinity without intersecting the positive part of a real semiaxis are referred to the contours of the first kind. The contours intersecting the positive part of a real semiaxis are referred to the contours of the second kind.

We will consider, at first, contours of the first kind. We will introduce the following notation: *zr* is the intersection point of the contour Γ with a circle |*z*| = *r*, *Cr* is an arc of a circle of radius *r* (not crossing the cut) which is formed when moving from the point *z* = *r* to the point *zr* and Γ*<sup>r</sup>*,*<sup>R</sup>* is a part of a contour Γ which is formed when moving from the point *zr* to the point *zR*. We form a closed contour *Gr*,*<sup>R</sup>* = [*r*, *R*] ∪ *CR* ∪ <sup>Γ</sup>¯*r*,*<sup>R</sup>* ∪ *Cr* (see Figure 3). The line means that we go along the contour in the opposite

direction. Since the function *h*(*z*) = *eizx <sup>g</sup>*+(*<sup>z</sup>*, *α*, −*<sup>θ</sup>*) is analytic in the region restricted by the contour *Gr*,*R*, then by using the Cauchy theorem

$$\int\_{G\_{r,\mathbb{R}}} h(z)dz = \int\_{r}^{\mathbb{R}} h(z)dz + \int\_{\mathbb{C}\_{\mathbb{R}}} h(z)dz - \int\_{\Gamma\_{r,\mathbb{R}}} h(z)dz + \int\_{\mathbb{C}\_{r}} h(z)dz = 0.$$

We will assume in this expression that *r* → 0 and *R* → ∞. Using Lemma 2 and Remark 1 we find that *CRh*(*z*)*dz* + *Crh*(*z*)*dz* → 0, at *R* → ∞ and *r* → 0. Therefore, the equality in Equation (59) is true.

**Figure 3.** Auxiliary contour *Gr*,*<sup>R</sup>* (heavy curve).

Now consider the contours of the second kind. These contours are characterized by the feature that they intersect the real axis. Therefore, to prove the lemma, we consider two closed auxiliary contours: the contour *Gr*,*<sup>x</sup>* = [*r*, *x*] ∪ <sup>Γ</sup>¯*r*,*<sup>x</sup>* ∪ *Cr* and contour *Gx*,*<sup>R</sup>* = [*x*, *R*] ∪ *CR* ∪ Γ*<sup>x</sup>*,*<sup>R</sup>* (see Figure 4). Here *x* is the intersection point of a contour Γ with a real axis. Since the function *h*(*z*) is analytic within the regions restricted with the contours Γ*<sup>r</sup>*,*<sup>x</sup>* and Γ*<sup>x</sup>*,*R*, then by using the Cauchy theorem

$$\int\_{\mathcal{G}\_{\mathcal{E}r}} h(z)dz = \int\_{r}^{\mathcal{X}} h(z)dz - \int\_{\Gamma\_{r,\mathcal{X}}} h(z)dz + \int\_{\mathcal{C}\_{r}} h(z)dz = 0,\tag{60}$$

$$\int\_{\mathcal{C}\_{x,\mathbb{R}}} h(z)dz = -\int\_{x}^{\mathbb{R}} h(z)dz + \int\_{\Gamma\_{x,\mathbb{R}}} h(z)dz + \int\_{\mathcal{C}\_{\mathbb{R}}} h(z)dz = 0. \tag{61}$$

Now we will assume in these expressions that *r* → 0 and *R* → ∞. Using Lemma 2 and Remark 1, we find that *Cr h*(*z*)*dz* → 0 at *r* → 0 and *CR h*(*z*)*dz* → 0 at *R* → ∞. Now summing up Equations (60) and (61) in view of this result we obtain

$$\int\_0^\chi h(z)dz + \int\_x^\infty h(z)dz = \int\_{\Gamma\_{0,x}} h(z)dz + \int\_{\Gamma\_{x,\infty}} h(z)dz,$$

hence, the equality in Equation (59) is true.

**Figure 4.** Auxiliary contours *Gr*,*<sup>x</sup>* (heavy curve) and *Gx*,*<sup>R</sup>* (dashed heavy curve).

The following lemma will be useful further

**Lemma 4.** *For any x* 0 *and any* −1 *θ* 1 *the function*

$$\varphi(\theta, \mathbf{x}) = \arctan\left(\frac{\mathbf{x} - \sin\left(\frac{\pi}{2}\theta\right)}{\cos\left(\frac{\pi}{2}\theta\right)}\right) \tag{62}$$

*and function ϕ*0(*<sup>θ</sup>*, *<sup>x</sup>*)*, defined by (12), are connected with relations between each other:*

*1. if x* 1*, then*

$$
\varphi\_0(\theta, \mathbf{x}) = \varphi(\theta, \mathbf{x}) - \pi/2;\tag{63}
$$

*2. if* 0 *x* < 1*, then*

$$\varphi\_0(\theta, \mathbf{x}) = \begin{cases} \begin{array}{ll} \varrho(\theta, \mathbf{x}) - \pi/2, & -1 \leqslant \theta \leqslant \theta\_0, \\ \varrho(\theta, \mathbf{x}) + \pi/2, & \theta\_0 < \theta \leqslant 1, \end{array} \tag{64}$$

*where θ*0 = (2/*π*) arcsin *x.*

**Proof.** Let us consider the case *x* > 1. Using the identity in Equation (42) as well as the identity

$$
\operatorname{arccot}(-y) = \pi - \operatorname{arccot}(y),
\tag{65}
$$

we obtain

$$\rho\_0 q\_0(\theta, \mathbf{x}) = \arctan\left(\frac{\cos(\pi \theta / 2)}{\sin(\pi \theta / 2) - \mathbf{x}}\right) = \frac{\pi}{2} - \operatorname{arccot}\left(\frac{\cos(\pi \theta / 2)}{\sin(\pi \theta / 2) - \mathbf{x}}\right) = -\frac{\pi}{2} + \operatorname{arccot}\left(\frac{\cos(\pi \theta / 2)}{\mathbf{x} - \sin(\pi \theta / 2)}\right).$$

Since *x* > 1, then the argumen<sup>t</sup> cos(*πθ*/2) *<sup>x</sup>*−sin(*πθ*/2) > 0 for any −1 *θ* 1. Now using the identity in Equation (41) for *y* > 0 we ge<sup>t</sup>

$$
\varrho\_0(\theta, \mathbf{x}) = \varrho(\theta, \mathbf{x}) - \pi/2, \quad \mathbf{x} > 1, \quad -1 \leqslant \theta \leqslant 1. \tag{68}
$$

Consider the case *x* = 1. In the same way as in our previous case, using the identities in Equations (42) and (65), we find

$$\varphi\_0(\theta, 1) = \arctan\left(\frac{\cos(\pi \theta / 2)}{\sin(\pi \theta / 2) - 1}\right) = -\frac{\pi}{2} + \operatorname{arccot}\left(\frac{\cos(\pi \theta / 2)}{1 - \sin(\pi \theta / 2)}\right). \tag{67}$$

From this expression one can see that the argumen<sup>t</sup> cos(*πθ*/2) <sup>1</sup>−sin(*πθ*/2) 0 for any −1 *θ* < 1 and it has an indeterminate form 0/0 at *θ* = 1. Evaluating indeterminate forms according to L'Hôpital's rule we ge<sup>t</sup> lim*θ*→1−<sup>0</sup> cos(*πθ*/2) <sup>1</sup>−sin(*πθ*/2) = ∞. Consequently, cos(*πθ*/2)/(<sup>1</sup> − sin(*πθ*/2)) 0, if −1 *θ* 1. Now using in Equation (67) the identity of Equation (41) for *y* > 0, we obtain

$$
\varphi\_0(\theta, 1) = \varphi(\theta, 1) - \pi/2, \quad -1 \lessapprox \theta \lessapprox 1.
$$

Combining now this expression and the expression of Equation (66) we come to Equation (63). Thus, the first item of the lemma is proved.

Now we will consider the case 0 *x* < 1. Using the identities in Equations (42) and (65) we ge<sup>t</sup>

$$q\_0(\theta, \mathbf{x}) = \arctan\left(\frac{\cos(\pi \theta / 2)}{\sin(\pi \theta / 2) - \mathbf{x}}\right) = -\frac{\pi}{2} + \operatorname{arccot}\left(\frac{\cos(\pi \theta / 2)}{\mathbf{x} - \sin(\pi \theta / 2)}\right). \tag{68}$$

Taking into consideration that cos(*πθ*/2) 0 for any *θ* ∈ [−1, 1], we find that the sign of cos(*πθ*/2) *<sup>x</sup>*−sin(*πθ*/2) is defined by the denominator. We have three possible situations: (1) *x* − sin(*πθ*/2) < 0, if *θ* > *θ*0, (2) *x* − sin(*πθ*/2) = 0, if *θ* = *θ*0, (3) *x* − sin(*πθ*/2) > 0, if *θ* < *θ*0. Here *θ*0 = (2/*π*) arcsin *x*. Taking into consideration Equation (44), we obtain cos(*πθ*/2) *<sup>x</sup>*−sin(*πθ*/2) = +<sup>∞</sup>, if *θ* = *θ*0. Thus, cos(*πθ*/2) *<sup>x</sup>*−sin(*πθ*/2) 0, if *θ* > *θ*0, and cos(*πθ*/2) *<sup>x</sup>*−sin(*πθ*/2) 0, if *θ θ*0. Now in Equation (68) applying the identity in Equation (41) we ge<sup>t</sup>

$$\eta\_0(\theta, x) = \begin{cases} -\frac{\pi}{2} + \arctan\left(\frac{x - \sin\left(\pi \theta / 2\right)}{\cos\left(\pi \theta / 2\right)}\right), & \theta \le \theta\_0, \\\ -\frac{\pi}{2} + \pi + \arctan\left(\frac{x - \sin\left(\pi \theta / 2\right)}{\cos\left(\pi \theta / 2\right)}\right), & \theta > \theta\_0. \end{cases}$$

From here it follows Equation (64). The lemma is completely proved.

In Figures 5 and 6, the graphs of the functions *ϕ*0(*<sup>θ</sup>*, *x*) and *ϕ*(*<sup>θ</sup>*, *x*) are given for the cases *x* 1 and 0 *x* < 1 which clearly illustrate the validity of the lemma that has just been proved.

**Figure 5.** The graph of the function *ϕ*(*<sup>θ</sup>*, *x*) and *ϕ*0(*<sup>θ</sup>*, *x*) depending on the parameter *θ* in the case of *x* 1 (The graphs are plotted for *x* = 1.5 and *x* = 1).

**Figure 6.** The graph of the function *ϕ*(*<sup>θ</sup>*, *x*) and *ϕ*0(*<sup>θ</sup>*, *x*) depending on the parameter *θ* in the case *x* < 1 (The graphs are plotted for *x* = 0.7 and *x* = 0).

### **3. Main Results**

Now we can formulate the main theorem which gives an integral representation for the probability density of a stable law *g*(*<sup>x</sup>*, *α*, *θ*) with the characteristic function in Equation (1).

**Theorem 1.** *The distribution density g*(*<sup>x</sup>*, *α*, *θ*) *of a strictly stable law with a characteristic function as in Equation (1) can be represented in the form*

*1. If α* -= 1 *and x* -= 0 *for any values* |*θ*| min(1, 2/*α* − 1)

$$\log(\mathbf{x}, a, \theta) = \frac{\mathfrak{a}}{\pi |\mathfrak{a} - 1|} \int\_{-\pi \theta^\*/2}^{\pi/2} \exp\left\{-|\mathbf{x}|^{a/(a-1)} \mathcal{U}(\mathfrak{q}, a, \theta^\*)\right\} \mathcal{U}(\mathfrak{q}, a, \theta^\*) |\mathbf{x}|^{1/(a-1)} d\mathfrak{q},\tag{69}$$

*where θ*∗ = *θ* sign *x and*

$$\mathcal{U}(\boldsymbol{\varrho},a,\theta) = \left(\frac{\sin\left(a\left(\boldsymbol{\varrho} + \frac{\pi}{2}\theta\right)\right)}{\cos\boldsymbol{\varrho}}\right)^{a/\left(1-a\right)} \frac{\cos\left(\boldsymbol{\varrho}(1-a) - \frac{\pi}{2}a\theta\right)}{\cos\boldsymbol{\varrho}}.\tag{70}$$

*2. If x* = 0*, then for any* 0 < *α* 2 *and* |*θ*| min(1, 2/*α* − 1)

$$\mathcal{S}(0, a, \theta) = \frac{1}{\pi} \cos \left( \pi \theta / 2 \right) \Gamma \left( 1/a + 1 \right) \tag{71}$$

*3. If α* = 1*, then for any* |*θ*| 1 *and any values x*

$$\log(\mathbf{x}, 1, \theta) = \frac{\cos(\pi \theta / 2)}{\pi (\mathbf{x}^2 - 2\mathbf{x}\sin(\pi \theta / 2) + 1)}.\tag{72}$$

**Proof.** To obtain the expression for the probability density we use the inversion formulas in Equation (9). In principle, it makes no difference which formula to use. The result will differ then only by the sign of the parameter *θ*. We will use the first formula in Equation (9), we have

$$\log(\mathbf{x}, \boldsymbol{\alpha}, \theta) = \frac{1}{\pi} \boldsymbol{\Re} \int\_0^\infty e^{it\mathbf{x}} \boldsymbol{\hat{g}}(t, \boldsymbol{\alpha}, -\theta) dt. \tag{73}$$

Without loss of generality, we assume that *x* > 0. The density *g*(*<sup>x</sup>*, *α*, *θ*) for *x* < 0 can be obtained with the inversion property in Equation (7). Next, let us make the substitution of integration variable

*t* → *z*, where *z* is complex-valued. Such a substitution means that we analytically extend the integral to the complex plane. With this analytical continuation is carried out starting with the positive part of a real semiaxis. As a result, we have that *g*<sup>ˆ</sup>(*<sup>t</sup>*, *α*, −*<sup>θ</sup>*) → *<sup>g</sup>*+(*<sup>z</sup>*, *α*, <sup>−</sup>*<sup>θ</sup>*), where *<sup>g</sup>*+(*<sup>z</sup>*, *α*, *θ*) is defined by Equation (11). The improper integral becomes the integral along the contour Γ. We will define the contour Γ in such a way that *eizx <sup>g</sup>*+(*<sup>z</sup>*, *α*, *θ*) = 0, and the contour itself Γ should start in the point *z* = 0 tend to infinity. Since in the inversion formula of Equation (73) the variable is *t* > 0 then from this it follows that arg *z* must lie within the limits − *π*/2 arg *z π*/2. As result, the contour Γ will take the form

$$\Gamma = \{ z : \text{"\'o\'e"} \text{g}^+(z, a, \theta) = 0, \, -\not\!\!f \, \text{\"\'o\'eg} \, z \lessapprox \not\!\!f, \, |z| \gg 0\}, \tag{74}$$

However, the specific type of the contour has to be defined. In view of the foregoing, the expression in Equation (73) will be written in the form

$$\log \left( x, a, \theta \right) = \frac{1}{\pi} \Re \int\_{\Gamma} \exp \{i zx - z^a \exp \{i \frac{\pi}{2} a \theta \} \} dz. \tag{75}$$

As a result, the problem consists in determining the contour form Γ, in proving the validity of transition from Equation (73) to Equation (75) and in calculating this integral.

Let us consider the case *α* -= 1. Representing *z* = *rei<sup>ϕ</sup>* and using this representation in Equation (75) for the intergrand we obtain

$$\exp\{i\mathbf{x}x - z^{a}\exp\{i\frac{\pi}{2}a\theta\}\} = \exp\left\{i\mathbf{x}r^{i\rho} - r^{a}\varepsilon^{ia\rho}\exp\{i\frac{\pi}{2}a\theta\}\right\} = $$

$$\exp\left\{-r\mathbf{x}\sin\varphi - r^{a}\cos\left(a(\varphi + \frac{\pi}{2}\theta)\right)\right\}\exp\left\{i[r\mathbf{x}\cos\varphi - r^{a}\sin\left(a(\varphi + \frac{\pi}{2}\theta)\right)]\right\}. \tag{76}$$

From the condition *eizx <sup>g</sup>*+(*<sup>z</sup>*, *α*, *θ*) = 0 and −*π* 2 arg *z π* 2follows that

$$r x \cos \varphi - r^{\alpha} \sin \left( a \left( \varphi + \frac{\pi}{2} \theta \right) \right) = 0. \tag{77}$$

The solution of this equation is in an explicit form

$$r(\varphi) = \left(\frac{\sin(a(\varphi + \frac{\varphi}{2}\theta))}{x\cos(\varphi)}\right)^{\frac{1}{1-x}}.\tag{78}$$

This expression determines the form of contour integration Γ.

We will determine the admissible region arg *z* ≡ *ϕ*. From the condition |*z*| 0 follows that (sin(*α*(*ϕ* + *π* 2 *<sup>θ</sup>*)))/(*x* cos(*ϕ*)) 0. Taking into account that *x* > 0 we obtain that the condition *r*(*ϕ*) > 0 is met if sin(*α*(*ϕ* + *π* 2 *θ*)) 0 and cos *ϕ* 0. From the condition sin(*α*(*ϕ* + *π* 2 *θ*)) 0 we obtain that −*π* 2 *θ ϕ π α* − *π* 2 *θ*, and from the condition cos *ϕ* 0 we ge<sup>t</sup> that −*π* 2 *ϕ π* 2 . Combining these two inequalities we obtain

$$\max\left(-\frac{\pi}{2}, -\frac{\pi}{2}\theta\right) \leqslant \varphi \leqslant \min\left(\frac{\pi}{2}, \frac{\pi}{\alpha} - \frac{\pi}{2}\theta\right). \tag{79}$$

Taking into consideration that at 0 < *α* < 1 the parameter *θ* takes values from the range −1 *θ* 1 and at 1 < *α* 2 from the range −( 2 *α* − 1) *θ* 2 *α* − 1 we obtain that max - −*π* 2 , −*π* 2 *θ* = −*π* 2 *θ* and min - *π* 2 , *π α* − *π* 2 *θ* = *π* 2 . Using this result in Equation (79) we ge<sup>t</sup> −*πθ*/2 *ϕ π*/2. Taking account of this condition and the expression in Equation (78) we obtain that in the case *α* -= 1 the contour Γ takes the form 

$$\Gamma = \left\{ z \, : \, r(\varphi) = \left( \frac{\sin(a\left(\varphi + \frac{\pi}{2}\theta\right))}{x\cos(\varphi)} \right)^{\frac{1}{1-a}}, -\frac{\pi}{2}\theta \lessdot \varphi \lessdot \frac{\pi}{2}, \, a \neq 1 \right\}. \tag{80}$$

We consider now the contour form Γ with different values of parameters. Replacing the boundary values *ϕ* in the expression for *<sup>r</sup>*(*ϕ*), we obtain

$$\begin{array}{llll}r\left(-\frac{\pi}{2}\theta\right)=0,&r\left(\frac{\pi}{2}\right)=\infty,&\quad\text{at}\quad 0$$

Thus, at 0 < *α* < 1 the contour Γ starts in the point *z* = 0 at *ϕ* = −*π*2 *θ* and tends to infinity at *ϕ* = *π*2 . In the case 1 < *α* 2 the situation is opposite: the contour Γ starts in the point *z* = 0 at *ϕ* = *π*2 and tends to infinity at *ϕ* = −*π*2*θ*.

As we can see the contours Γ described in Equation (80) with different values of the parameter *θ* differ in its type (see Figures 7 and 8). They can be divided into 7 main groups.


**Figure 7.** The type of a contour Γ at *α* ∈ (0, 1) and various values of parameter *θ*. The contours are given for the case *α* = 0.6 and Γ1 - *θ* = −1, Γ2 - *θ* = −0.75, Γ2 - *θ* = −0.5, Γ2 - *θ* = −0.25, Γ 2 - *θ* = 0, Γ3 - *θ* = 0.25, Γ3- *θ* = 0.5, Γ3- *θ* = 0.75, Γ 3- *θ* = 0.98, Γ4- *θ* = 1.

**Figure 8.** The type of a contour Γ at *α* ∈ (1, 2] and with different values of the parameter *θ*. The contours are given for the case *α* = 1.25 and Γ5 - *θ* = −0.6, Γ5 - *θ* = −0.4, Γ5 - *θ* = −0.2, Γ 5 - *θ* = 0, Γ6 - *θ* = 0.2, Γ6 - *θ* = 0.57, Γ7 - *θ* = 0.6.

This indicates that the contours of the first, second and third groups (0 < *α* < 1, −1 *θ* < 1) and contours of the fifth and sixth groups (1 < *α* 2, −(2/*<sup>α</sup>* − 1) *θ* < 2/*α* − 1) satisfy the conditions of the Lemma 3. In fact, the contours of these groups start in the point *z* = 0, do not cross the line of the cut that goes through a half-line arg *z* = −*π* and tend to infinity. The contours of the first, second, third groups tend to infinity at *ϕ* = *π*2 . This means that the contours of these groups satisfy the condition in Equation (51). The contours of the fifth and sixth groups tend to infinity at *ϕ* = −*π*2 *θ*. Thus, these contours satisfy the condition in Equation (52).

Consider the contours of the fourth group Γ4. The contours of this group start in the point *z* = −*<sup>i</sup>*(*α*/*x*)1/(<sup>1</sup>−*<sup>α</sup>*) and tend to infinity at *ϕ* → *π*2 . Consider an auxiliary contour Γ∗ = {*z* : arg *z* = −*<sup>π</sup>*/2, 0 |*z*| (*α*/*x*)1/(<sup>1</sup>−*<sup>α</sup>*)}. With the help of a contour Γ∗ we form a new contour *S* = Γ∗ 5 Γ4. The specific feature of the contour Γ∗ is then that exp{*izx* − *z<sup>α</sup>* exp{*<sup>i</sup> π*2 *αθ*}} = 0 for *z* ∈ Γ<sup>∗</sup>. From this it follows that

$$\frac{1}{2\pi} \Re \int\_{\Gamma^\*} \exp\{izx - z^a \exp\{i\frac{\pi}{2}a\theta\}\} dz = 0.$$

Therefore, for the case 0 < *α* < 1, *θ* = 1 we ge<sup>t</sup>

$$\frac{1}{\pi}\Re\int\_{S}\exp\{izx - z^{\underline{a}}\exp\{i\frac{\pi}{2}a\theta\}\}dz = \frac{1}{\pi}\Re\int\_{\Gamma\_{4}}\exp\{izx - z^{\underline{a}}\exp\{i\frac{\pi}{2}a\theta\}\}dz.\tag{82}$$

However, now the contour *S* = Γ∗ 5 Γ4 completely satisfies the conditions of the Lemma 3: it starts in the point *z* = 0 without crossing the line of the cut and tends to infinity at *ϕ* = *π*/2.

For the contours of the seventh froup Γ7 we do the same. Consider an auxiliary contour Γ∗∗ = {*z* : arg *z* = *π*/2, 0 |*z*| (*x*/*α*)1/(*<sup>α</sup>*−<sup>1</sup>)}. With the help of this contour we form the contour *S*∗ = Γ∗∗ 5 Γ7. Now, the contour *S*∗ completely satisfies the conditions of the Lemma 3. As in the previous case exp{*izx* − *z<sup>α</sup>* exp{*<sup>i</sup> π*2*αθ*}} = 0 for *z* ∈ Γ∗∗. From this it follows

$$\frac{1}{\pi}\mathfrak{R}\int\_{\Gamma^{\*\*}}\exp\{izx - z^{a}\exp\{i\oint a\theta\}\}dz = 0.$$

As a result, for the case 1 < *α* 2, *θ* = 2/*α* − 1 we ge<sup>t</sup>

$$\frac{1}{\pi} \Re \int\_{S^\*} \exp\{izx - z^a \exp\{i\frac{\pi}{2}a\theta\}\} dz = \frac{1}{\pi} \Re \int\_{\Gamma\_\mathcal{I}} \exp\{izx - z^a \exp\{i\frac{\pi}{2}a\theta\}\} dz. \tag{83}$$

Now, applying Lemma 3 and taking account of the equalities in Equations (82) and (83) we find that

$$\frac{1}{2\pi}\mathfrak{R}\int\_{0}^{\infty}\mathfrak{e}^{i\text{tr}}\mathfrak{F}(t,\mathfrak{a},-\theta)dt = \frac{1}{\pi}\int\_{\Gamma}\mathfrak{e}^{i\text{zx}}\mathfrak{g}^{+}(z,\mathfrak{a},-\theta)dz,\tag{84}$$

where the contour Γ is defined by the expression in Equation (80). As a result, an improper integral along the positive part of the real axis in the expression of Equation (73) can be replaced with an integral along the contour Γ. Thus, in the case considered (*α* -= 1) we showed the validity of transition from Equation (73) to Equation (75).

Returning to Equation (75), taking into consideration Equation (76) and representing complex *z* in the form *z* = *reiϕ*, we obtain

$$\begin{split} \log(z,a,\theta) &= \frac{1}{\pi} \Re \int\_{\Gamma} \exp\{i\mathbf{z}\mathbf{x} - z^{a} \exp\{i\frac{\pi}{2}a\theta\}\} dz = \frac{1}{\pi} \int\_{\Gamma} \Re \exp\{i\mathbf{z}\mathbf{x} - z^{a} \exp\{i\frac{\pi}{2}a\theta\}\} d(\theta \mathbf{z}) = \\ &\frac{1}{\pi} \int\_{\Gamma} \exp\{-r\mathbf{x}\sin\varphi - r^{a} \cos(a(\varphi + \frac{\pi}{2}\theta))\} \cos\left(r\mathbf{x}\cos\varphi - r^{a} \sin(a(\varphi + \frac{\pi}{2}\theta))\right) d[r\cos\varphi]. \end{split}$$

We will place here the expression in Equation (78). One should pay attention that in this case the Equation (77) is valid, consequently, we obtain

$$\log(\mathbf{x}, a, \theta) = \frac{1}{\pi} \int\_{\Gamma} \exp\{A(\varphi)\} d[r(\varphi)\cos\varphi],\tag{85}$$

where *<sup>A</sup>*(*ϕ*) = −*<sup>r</sup>*(*ϕ*)*<sup>x</sup>* sin *ϕ* − *r*(*ϕ*)*<sup>α</sup>* cos(*α*(*ϕ* + *π*2 *θ*)) and *r*(*ϕ*) is defined by Equation (78). We transform the function *<sup>A</sup>*(*ϕ*)

$$A(\boldsymbol{\varrho}) = -\mathbf{x} \left( \frac{\sin(a(\boldsymbol{\varrho} + \frac{\pi}{2}\theta))}{\mathbf{x} \cos(\boldsymbol{\varrho})} \right)^{1/(1-a)} \sin \boldsymbol{\varrho} - \left( \frac{\sin(a(\boldsymbol{\varrho} + \frac{\pi}{2}\theta))}{\mathbf{x} \cos(\boldsymbol{\varrho})} \right)^{a/(1-a)} \cos(a(\boldsymbol{\varrho} + \frac{\pi}{2}\theta))$$

$$= -\mathbf{x}^{a/(a-1)} \left( \frac{\sin(a(\boldsymbol{\varrho} + \frac{\pi}{2}\theta))}{\cos(\boldsymbol{\varrho})} \right)^{a/(1-a)} \frac{\cos \left(\boldsymbol{\varrho}(1-a) - \frac{\pi}{2}a\theta\right)}{\cos \boldsymbol{\varrho}} = -\mathbf{x}^{a/(a-1)} \mathcal{U}(\boldsymbol{\varrho}, a, \boldsymbol{\varrho}), \quad \text{(86)}$$

where *<sup>U</sup>*(*ϕ*, *α*, *θ*) has the form of Equation (70).

Now we consider the differential, we have *d*[*r*(*ϕ*) cos *ϕ*] = cos *ϕd*[*r*(*ϕ*)] − *r*(*ϕ*) sin *ϕdϕ*. For the differential *d*[*r*(*ϕ*)] we ge<sup>t</sup>

$$d[r(\boldsymbol{\varrho})] = d\left[ \left( \frac{\sin(a(\boldsymbol{\varrho} + \frac{\pi}{2}\theta))}{\mathbf{x}\cos(\boldsymbol{\varrho})} \right)^{1/(1-a)} \right] = \frac{r(\boldsymbol{\varrho})}{1-a} \left( a \cot\left( a\left(\boldsymbol{\varrho} + \frac{\pi}{2}\theta\right)\right) + \tan\left\{\boldsymbol{\varrho}\right\} d\boldsymbol{\varrho} \cdot \mathbf{x} \right)$$

Using now this result in the expression for *d*[*r*(*ϕ*) cos *ϕ*], we have

$$d[r(\boldsymbol{\varrho})\cos\boldsymbol{\varrho}] = \frac{r(\boldsymbol{\varrho})}{1-a}\cos\boldsymbol{\varrho}\cdot\left(a\cot\left(a\left(\boldsymbol{\varrho}+\frac{\pi}{2}\theta\right)\right)+\tan\,\boldsymbol{\varrho}\right)d\boldsymbol{\varrho}-r(\boldsymbol{\varrho})\sin\boldsymbol{\varrho}d\boldsymbol{\varrho}$$

$$=\frac{ar(\boldsymbol{\varrho})}{1-a}\left(\cos\boldsymbol{\varrho}\cot\left(a\left(\boldsymbol{\varrho}+\frac{\pi}{2}\theta\right)\right)+\sin\,\boldsymbol{\varrho}\right)d\boldsymbol{\varrho}$$

$$=\frac{a}{1-a}\left(\frac{\sin\left(a\left(\boldsymbol{\varrho}+\frac{\pi}{2}\theta\right)\right)}{x\cos\boldsymbol{\varrho}}\right)^{1/\left(1-a\right)}\frac{\cos\left(a\left(\boldsymbol{\varrho}+\frac{\pi}{2}\theta\right)-\boldsymbol{\varrho}\right)}{\sin\left(a\left(\boldsymbol{\varrho}+\frac{\pi}{2}\theta\right)\right)}d\boldsymbol{\varrho}=\frac{a}{1-a}\mathbf{x}^{1/\left(a-1\right)}\mathcal{U}(\boldsymbol{\varrho},a,\theta)d\boldsymbol{\varrho},\quad(87)$$

where *<sup>U</sup>*(*ϕ*, *α*, *θ*) has the form of Equation (70).

Now using Equations (86) and (87) in Equation (85) and also taking into consideration that the motion along the contour Γ, having the form of Equation (80) now described by the parameter change *ϕ*, we obtain

$$\log(\mathbf{x}, a, \theta) = \frac{\mathfrak{a}}{\pi(1-\mathfrak{a})} \int\_{\varphi\_{\min}}^{\varphi\_{\max}} \exp\left\{-\mathbf{x}^{a/(a-1)} \mathcal{U}(\varphi, \mathbf{a}, \theta)\right\} \mathbf{x}^{1/(a-1)} \mathcal{U}(\varphi\_{\mathbf{a}}, \theta) d\mathfrak{a} \tag{88}$$

From the expression in Equation (84) one can see that the integration limits *ϕmin* and *ϕmax* should be selected in such a way that the motion along the contour Γ could correspond to a change from *r* = 0 to *r* = ∞. From the expression in Equation (81) it is clear that in the case 0 < *α* < 1 a change in the angle *ϕ* from −*π*2 *θ* to *π*2 corresponds to the motion along the contour Γ from *r* = 0 to *r* = ∞. Therefore, in this case *ϕmin* = −*π*2 *θ*, *ϕmax* = *π*2 . In the case 1 < *α* 2 a change in the angle *ϕ* from *π*2 to −*π*2 *θ* corresponds to the motion from *r* = 0 to *r* = ∞. Therefore, in this case *ϕmin* = *π*2 , and *ϕmax* = −*π*2 *θ*. Combining these two cases the expression in Equation (88) takes the form

$$\log(\mathbf{x}, a, \theta) = \frac{\mathbf{a}}{\pi |1 - a|} \int\_{-\frac{\pi}{2}\theta}^{\frac{\pi}{2}} \exp\left\{-\mathbf{x}^{a/(a-1)} \mathcal{U}(\varphi, a, \theta)\right\} \mathbf{x}^{1/(a-1)} \mathcal{U}(\varphi, a, \theta) d\varphi, \quad \mathbf{a} \neq 1. \tag{89}$$

It should be pointed out that this formula was obtained on the assumption *x* > 0. The case *x* < 0 is easy to obtain using the property of inversion for the density of probabilities in Equation (7). For this it is enough to replace the parameter *θ* in the formula in Equation (89) with <sup>−</sup>*θ*. It is possible to combine these two cases if in the formula of Equation (89) to substitute the parameter *θ* for the parameter *θ*∗ = *θ* sign *x* and the value *x* is taken in absolute value. As a result, we obtain the expression in Equation (69) valid for any *x* -= 0 and any admissible *α* -= 1.

Now we consider the case 0 < *α* 2 and *x* = 0. In this case the inversion formula of Equation (75) takes the form

$$\log(0, a, \theta) = \frac{1}{\pi} \Re \int\_{\Gamma} \exp\left\{-z^a \exp\left\{\frac{\pi}{2} a\theta\right\}\right\} dz,\tag{90}$$

where the contour Γ is determined by the expression in Equation (74), but it is necessary to determine the specific type of a contour in the case under consideration.

We represent the complex number *z* in the form *z* = *reiϕ*. As a result, the integrand in Equation (90) takes the form

$$\exp\left\{-z^{a}\exp\left\{\frac{\pi}{2}a\theta\right\}\right\}=\exp\left\{-r^{a}\cos\left(a\left(\varphi+\frac{\pi}{2}\theta\right)\right)-ir^{a}\sin\left(a\left(\varphi+\frac{\pi}{2}\theta\right)\right)\right\}.\tag{91}$$

Using the condition *g*+(0, *α*, *θ*) = 0, we obtain

$$r^a \sin\left(a\left(\varphi + \frac{\pi}{2}\theta\right)\right) = 0.\tag{92}$$

This equation has two solutions: *r* = 0, for any *ϕ*, and *ϕ* = −*πθ*/2, *r* 0. It is clear that if *r* = 0, then *g*+(0, *α*, *θ*) = 0 for any value of *ϕ*, for definiteness we will select *ϕ* = −*πθ*/2 if *r* = 0. As a result, the contour of integration of Equation (74) takes the form Γ = {*z* : |*z*| 0, arg *z* = <sup>−</sup>*πθ*/2}, where |*θ*| min(1, 2/*α* − <sup>1</sup>), 0 < *α* 2. From this it is clear that in the case under consideration the family of contours Γ are represented by half-lines starting from the point *z* = 0 at an angle arg *z* = −*πθ*/2. Consequently, the family of contours Γ satisfy items 1 and 2 of Lemma 3 and also the condition in Equation (53) in item 3. Thus, we substantiate the transition from the improper integral of Equation (73) to the contour integral in Equation (90).

Now we put the expression in Equation (91) in Equation (90) and taking into account Equation (92), we obtain

$$\log(0, a, \theta) = \frac{1}{\pi} \int\_{\Gamma} \exp\left\{-r^a \cos\left(a\left(\varphi + \frac{\pi}{2}\theta\right)\right)\right\} \,\mathbb{R}[dz] = \frac{\cos(\pi\theta/2)}{\pi} \int\_0^\infty \exp\left\{-r^a\right\} dr.$$

Here it was taken into account that %[*dz*] = cos(*πθ*/2)*dr* on the contour Γ. The limits of integration were selected in such a way that when moving along the contour *r* could change from 0 to ∞. Since the contour Γ in the case under consideration is a half-line coming out of the point *z* = 0 at an angle arg *z* = −*πθ*/2, then the motion along the contour Γ from 0 to ∞ corresponds to a change *r* from 0 to ∞. Making in this integral a substitution of a variable *r<sup>α</sup>* = *y* and using the definition of the gamma-function <sup>Γ</sup>(*n*) = ∞0 *<sup>t</sup><sup>n</sup>*−1*e*<sup>−</sup>*tdt*, *n* > 0 we obtain the expression Equation (71). Thus, the second item of the theorem is proved.

Now we consider the case *α* = 1. We will make an assumption that *x* 0. In this case, the inversion formula in Equation (75) takes the form

$$\log(\mathbf{x}, 1, \theta) = \frac{1}{\pi} \Re \int\_{\Gamma} \exp\left\{i\mathbf{z}\mathbf{x} - \mathbf{z}\exp\left\{i\frac{\pi}{2}\theta\right\}\right\} d\mathbf{z},\tag{93}$$

where the contour of integration defined by Equation (74) will be written in the form

$$\Gamma = \{ z : \text{"\Sot"} ^{\sharp \text{xy}} \text{g}^+ (z, 1, \theta) = 0, \text{ } ^\pi \text{g} \text{ or } \text{arg} \, z \leqslant \frac{\pi}{2}, \, |z| \geqslant 0 \}, \tag{94}$$

Here it should be noted that analytic continuation of the function *g*<sup>ˆ</sup>(*<sup>t</sup>*, 1, −*<sup>θ</sup>*) from the positive part of the real axis *t* to the complex plane *z* at *α* = 1 is an analytic function and it has the form *<sup>g</sup>*+(*<sup>z</sup>*, 1, −*<sup>θ</sup>*) = exp −*<sup>z</sup>* exp *i π*2 *<sup>θ</sup>*.

As in the previous case, we begin by defining the form of the integration contour Γ. Consider the integrand in Equation (93) and we represent the complex number *z* in the form *z* = *reiϕ*. As a result, we obtain

$$\begin{aligned} \exp\{i\mathbf{z}x - z\exp\{i\overline{\frac{\alpha}{2}}\theta\}\} &= \exp\left\{i\mathbf{x}r^{i\rho} - r\varepsilon^{j\rho}\exp\{i\overline{\frac{\alpha}{2}}\theta\}\right\} = \\ &\exp\{-r\left(\mathbf{x}\sin\varphi + \cos(\left(\varphi + \frac{\pi}{2}\theta\right)\rangle\right)\}\exp\left\{i\mathbf{r}\left(\mathbf{x}\cos\varphi - \sin(\left(\varphi + \frac{\pi}{2}\theta\right)\rangle\right)\right\}. \end{aligned} \tag{95}$$

From this expression we ge<sup>t</sup> that the condition *eizx <sup>g</sup>*+(*<sup>z</sup>*, 1, −*<sup>θ</sup>*) = 0 leads to an equation

$$r\left(\mathbf{x}\cos\varphi-\sin\left(\left(\varphi+\frac{\pi}{2}\theta\right)\right)\right)=0.\tag{96}$$

This equation has two solutions. The first solution is *r* = 0. The second solution we obtain from the equation *x* cos *ϕ* − sin((*ϕ* + *π*2 *θ*)) = 0, *r* 0. Solving it with respect to *ϕ* we ge<sup>t</sup>

$$\varphi(\theta, \mathbf{x}) = \arctan\left(\frac{\mathbf{x} - \sin(\pi \theta / 2)}{\cos(\pi \theta / 2)}\right). \tag{97}$$

As a result, the contour of integration in Equation (94) takes the form

$$\Gamma = \{ z \colon \arg z = \varphi(\theta, \mathbf{x}), \ |z| \ni 0, -1 \leqslant \theta \leqslant 1 \} \tag{98}$$

Thus, the contours Γ with different values *θ* are half-lines coming out of the origin at an angle *ϕ*(*<sup>θ</sup>*, *x*) and tending to infinity.

Next, it is necessary to substantiate the transition from the integral in Equation (73) to the integral along the contour in Equation (93). For this we will use Lemmas 3 and 4. From the definition in Equation (98), it is clear that for all admissible values *θ* and *x* 0 the contour Γ satisfies item 1 and 2 of the Lemma 3. There is only one thing left, to find out if the contour in Equation (98) satisfies item 3 of this lemma.

Consider the case *x* > 1 at first. According to item 3 of the Lemma 3 for the equality in Equation (59) to be valid the contour of Equation (98) must satisfy the condition in Equation (54). Now using the Lemma 4. According to this lemma in the case *x* 1 the functions *ϕ*(*<sup>θ</sup>*, *x*) and *ϕ*0(*<sup>θ</sup>*, *x*) are connected between each other with a ratio (63) from which it directly follows that *ϕ*0(*<sup>θ</sup>*, *x*) < *ϕ*(*<sup>θ</sup>*, *x*) for all *θ* ∈ [−1, 1]. Therefore, at *x* > 1 the condition in Equation (54) is met and we can move from the integral in Equation (73) to the integral in Equation (93). In the case *x* = 1 the contour Γ must meet the condition in Equation (58). Now using the Lemma 4 we obtain that in this case for all *θ* ∈ [−1, 1] the inequality *ϕ*0(*<sup>θ</sup>*, *x*) < *ϕ*(*<sup>θ</sup>*, *x*) is true. Therefore, in this case the condition in Equation (58) is satisfied. Now we consider the case 0 *x* < 1. Similar to previous case, applying the Lemma 4 namely, the formula in Equation (64) we ge<sup>t</sup> that in this

case the contour in Equation (98) satisfies the conditions in Equations (55)–(57) Lemma 3. Thus, the contour in Equation (98) completely satisfies the conditions of the Lemma 3 and therefore, the equality

$$\frac{1}{\pi\pi}\mathfrak{R}\int\_0^\infty e^{it\chi}\mathfrak{J}(t,1,-\theta)dt = \frac{1}{\pi\pi}\int\_\Gamma e^{iz\chi}\mathfrak{J}^+(z,1,-\theta)dz$$

is true. This makes it possible to replace the improper integral in the expression of Equation (73) over the real variable with the integral over the contour Γ. Thus, the possibility for the transition from the formula in Equation (73) to the formula in Equation (75) is substantiated. We need to note that in the case under consideration *α* = 1 the expression in Equation (75) takes the form of Equation (93). It should also be pointed out that that due to Equation (34) the case *α* = 1, *θ* = 1, *x* = 1 is excluded in the Lemma 3. That is why, here, this case should also be excluded from consideration.

Now taking into account Equation (95), the expression of Equation (93) will be written in the following form

$$\log(\mathbf{x},1,\theta) = \frac{1}{\pi} \Re \int\_{\Gamma} \exp\{-r\left(\mathbf{x}\sin\varphi + \cos(\varphi + \frac{\pi}{2}\theta)\right)\} \exp\left\{i\left(\mathbf{x}\cos\varphi - \sin(\varphi + \frac{\pi}{2}\theta)\right)\right\} d\mathbf{z} = \frac{1}{\pi} \Im \int\_{\Gamma} \exp\left\{-r\left(\mathbf{x}\sin\varphi + \cos(\varphi + \frac{\pi}{2}\theta)\right)\right\} \text{Res}\left\{r\left(\mathbf{x}\cos\varphi - \sin(\varphi + \frac{\pi}{2}\theta)\right)\right\} \Re[d\mathbf{z}] = \frac{1}{\pi} \Im \int\_{\Gamma} \exp\left\{-r\left(\mathbf{x}\cos\varphi - \sin(\varphi + \frac{\pi}{2}\theta)\right)\right\} d\mathbf{z}$$

Now using here the definition of Equation (98) we obtain that %[*dz*] = cos(*ϕ*(*<sup>θ</sup>*, *<sup>x</sup>*))*dr*. The motion along the contour Γ should take place in such a way that it would start in the point *z* = 0 in the process of moving it would tend to infinity. Therefore, the motion from *r* = 0 to *r* = ∞ corresponds to such motion. Taking into consideration that on the contour Γ the condition of Equation (96) is met, we obtain

$$\log(\mathbf{x}, 1, \theta) = \frac{1}{\pi} \int\_0^\infty \exp\{-r \left(\mathbf{x} \sin \varphi(\theta, \mathbf{x}) + \cos(\varphi(\theta, \mathbf{x}) + \frac{\pi}{2} \theta) \right) \} \cos(\varphi(\theta, \mathbf{x})) dr$$

As we can see, the integral obtained is easy to calculate. As a result, we obtain

$$g(\mathbf{x}, 1, \theta) = \frac{1}{\pi} \frac{\cos(\varrho(\theta, \mathbf{x}))}{\mathbf{x} \sin \varrho(\theta, \mathbf{x}) + \cos(\varrho(\theta, \mathbf{x}) + \pi \theta/2)}.$$

Using now the definition of the function *ϕ*(*<sup>θ</sup>*, *x*) (97) after simple transformations we ge<sup>t</sup>

$$\log(\mathbf{x}, 1, \theta) = \frac{\cos(\pi \theta / 2)}{\pi \left(\mathbf{x}^2 - 2\mathbf{x}\sin(\pi \theta / 2) + 1\right)}.\tag{99}$$

Recall that consideration was carried out for the case *x* 0. The case *x* < 0 can be obtained using the inversion property for the density of probability of Equation (7). For this it is enough to replace *θ* with <sup>−</sup>*θ*. It is possible to combine these two cases if to introduce a parameter *θ*∗ = *θ* sign *x*. However, we want to note that if to perform this replacement in the expression of Equation (99), then the expression itself will not change *g*(*<sup>x</sup>*, 1, *θ*) → *g*(*<sup>x</sup>*, 1, *θ*∗) ≡ *g*(*<sup>x</sup>*, 1, *<sup>θ</sup>*). Consequently, the formula in Equation (99) is true for any *x*. Thus, the theorem is completely proved.

We will make some remarks on the proved theorem.

**Remark 2.** *The proof of the case α* = 1 *was carried out under the assumption x* 0*. Therefore, the formula in Equation (99) was obtained for the case x* 0*. The generalization of this formula for the case x* < 0 *was carried out using the inversion property g*(−*x*, 1, *θ*) = *g*(*<sup>x</sup>*, 1, <sup>−</sup>*<sup>θ</sup>*)*. It should be noted here that in the process of proving this case the point θ* = 1*, x* = 1 *was excluded from consideration. Therefore, in view of the inversion property,* *the point θ* = <sup>−</sup>1*, x* = −1 *should also be excluded. In these two points the density g*(*<sup>x</sup>*, 1, *θ*) *has a peculiarity. In fact, substituting the value θ* = ±1*, we obtain*

$$g(\mathbf{x}, 1, \pm 1) = \frac{\cos(\pi/2)}{\pi(\mathbf{x} \mp 1)^2} = \begin{cases} 0, & \mathbf{x} \neq \pm 1, \\ 0/0, & \mathbf{x} = \pm 1. \end{cases}$$

*However, an indeterminate form* 0/0 *can be evaluated in the following remark.*

**Remark 3.** *Equations (71) and (72) can be obtained directly from the inversion formulas in Equation (9) without resorting to analytic continuation of the characteristic function g*<sup>ˆ</sup>(*<sup>t</sup>*, *α*, *θ*) *to the complex plane with the subsequent transition from the improper integral over the real variable (73) to the integral along the contour in Equation (75). We first consider the case α* = 1*. Using the inversion formula (the first formula in Equation (9)) we obtain*

$$\log(\mathbf{x}, 1, \theta) = \frac{1}{\pi} \mathfrak{R} \int\_0^\infty e^{it\mathbf{x}} \mathfrak{J}(t, 1, -\theta) dt = \frac{1}{\pi} \mathfrak{R} \int\_0^\infty \exp\left\{ t \left( \mathbf{i}\mathbf{x} - \exp\left\{ i\frac{\pi}{2}\theta \right\} \right) \right\} dt = \frac{1}{\pi} \mathfrak{R} I,\tag{100}$$

*where I* = ∞ 0 exp *t ix* − exp *i π* 2 *θ dt. We determine under which conditions this integral will converge. For the integral I the inequality is valid*

$$|I| \leqslant \int\_0^\infty \left| \exp\left\{ t \left( i\mathbf{x} - \exp\left\{ i\frac{\pi}{2}\theta \right\} \right) \right\} \right| dt \tag{101}$$

*For the integrand we have*

$$\left| \exp \left\{ t \left( \mathbf{ix} - \exp \left\{ i \frac{\pi}{2} \theta \right\} \right) \right\} \right| = \exp \left\{ -t \cos(\pi \theta / 2) \right\} \cos(t(\mathbf{x} - \sin(\pi \theta / 2))). \tag{102}$$

*Since in the case considered (α* = 1*) the parameter θ varies within the limits* −1 *θ* 1*, we obtain* cos(*πθ*/2) > 0, *if* − 1 < *θ* < 1*. From here we get* lim*<sup>t</sup>*→ ∞ exp *t ix* − exp *i π* 2 *θ* = 0*, if* −1 < *θ* < 1*. Thus, the integral in Equation (101) will converge, and, therefore, and the integral I will also converge at* −1 < *θ* < 1*.*

*We consider now the cases θ* = ±1*. If θ* = 1*, then from Equation (102) we obtain*

$$\left| \exp \left\{ t \left( i \mathbf{x} - \exp \{ i \pi/2 \} \right) \right\} \right| = \cos(t(\mathbf{x} - 1)). \tag{103}$$

*Substituting this result in Equation (101) we get*

$$|I| \leqslant \int\_0^\infty \cos(t(\mathbf{x} - 1))dt = \begin{cases} 0, & \mathbf{x} \neq 1 \\ \rightsquigarrow, & \mathbf{x} = 1. \end{cases} \tag{104}$$

*Similarly, if θ* = <sup>−</sup>1*, we have*

$$|I| \leqslant \int\_0^\infty \cos(t(\mathbf{x} + 1))dt = \begin{cases} 0, & \mathbf{x} \neq -1 \\ \infty, & \mathbf{x} = -1. \end{cases} \tag{105}$$

*Thus, at θ* = ±1 *the integral I* = 0 *for all x* -= ±1*, and integral I will diverge in the points x* = 1, *θ* = 1 *and x* = −1, *θ* = <sup>−</sup>1*.*

> *Returning to Equation (100) and by calculating the integral directly, we get*

$$\begin{split} \log(\mathbf{x}, 1, \theta) &= \frac{1}{\pi} \mathfrak{R} \int\_0^\infty \exp\left\{ t \left( \mathbf{ix} - \exp\left\{ i \frac{\pi}{2} \theta \right\} \right) \right\} dt = -\frac{1}{\pi} \mathfrak{R} \frac{1}{i\mathbf{x} - \exp\{ i \pi \theta / 2 \}} \\ &= -\frac{1}{\pi} \mathfrak{R} \frac{-\mathbf{ix} - \exp\{ -i \pi \theta / 2 \}}{(i\mathbf{x} - \exp\{ i \pi \theta / 2 \})(-i\mathbf{x} - \exp\{ -i \pi \theta / 2 \})} = \frac{1}{\pi} \frac{\cos(\pi \theta / 2)}{\mathbf{x}^2 - 2 \mathbf{x} \sin(\pi \theta / 2) + 1}. \end{split} \tag{106}$$

*From the formula obtained it is clear that if θ* = ±1

$$g(\mathbf{x}, 1, \pm 1) = \frac{\cos(\pi/2)}{\pi(\mathbf{x} \mp 1)^2} = \begin{cases} 0, & \mathbf{x} \neq \pm 1 \\ \quad \infty, & \mathbf{x} = \pm 1. \end{cases} \tag{107}$$

*Thus the behavior of the formula in Equation (106) coincides with the behavior of the integral I in Equations (104) and (105) in the cases θ* = ±1*. Therefore, the formula in Equation (106) is true for any* −1 *θ* 1*. The expression in Equation (107) means that the density g*(*<sup>x</sup>*, 1, ±1) *is a degenerate distribution in the point x* = ±1*. In other words,*

$$g(\mathbf{x}, \mathbf{1}, \pm \mathbf{1}) = \delta(\mathbf{x} \mp \mathbf{1}).\tag{108}$$

Thus, the obtained expression in Equation (106) completely coincides with the one previously obtained in the Theorem 1 the density in Equation (72), and the conclusion presented in this remark is an alternative way of deducing this density.

**Remark 4.** *By a similar method, one can obtain the density value at x* = 0*. Using the first formula in Equation (9) and making a substitution of the integration variable tα* = *τ we get*

$$\begin{split} g(0, a, \theta) = \frac{1}{\pi} \mathfrak{R} \int\_0^\infty \hat{g}(t, a, -\theta) dt &= \frac{1}{\pi} \mathfrak{R} \int\_0^\infty \exp\left\{-t^a \exp\left\{i\frac{\pi}{2} a\theta\right\}\right\} dt \\ &= \frac{1}{a\pi} \mathfrak{R} \int\_0^\infty \exp\left\{-\tau \cos\left(\frac{\pi}{2} a\theta\right) - i\tau \sin\left(\frac{\pi}{2} a\theta\right)\right\} \tau^{1/a - 1} d\tau \\ &= \frac{1}{a\pi} \int\_0^\infty \exp\left\{-\tau \cos\left(\frac{\pi}{2} a\theta\right)\right\} \cos\left(y \sin\left(\frac{\pi}{2} a\theta\right)\right) \tau^{1/a - 1} d\tau. \end{split} \tag{109}$$

*To calculate the integral obtained it is necessary to use the formula (see [45], Section 1.5. the Equation (35))*

$$\int\_0^\infty t^{a-1} e^{-\varepsilon t \cos \beta} \cos(\varepsilon t \sin \beta) dt = \Gamma(a) \varepsilon^{-a} \cos(a\beta), \quad \varepsilon > 0, \ \Re a > 0, \ -\pi/2 < \beta < \pi/2.$$

*It is clear that the integral in Equation (109) completely satisfies the conditions of this integral. Consequently, using it in Equation (109) we get g*(0, *α*, *θ*) = 1*π* <sup>Γ</sup>(1/*α* + 1) cos(*πθ*/2)*. The formula obtained coincides completely with the formula in Equation (71).*

We will make another useful remark.

**Remark 5.** *The Theorem 1 formulates an integral representation for the density of a standard strictly stable law. However, it is useful to have a formula that allows one to convert the density of a standard strictly stable law to the density of a strictly stable law with arbitrary λ. The Property 2 and, in particular, the formula in Equation (8) allows one to obtain such a formula.*

> *In fact, in terms of characteristic functions, the formula in Equation (8) will be written as*

$$\mathcal{G}(t,\mathfrak{a},\theta,\lambda) = \mathcal{G}\left(\lambda^{1/\mathfrak{a}}t,\mathfrak{a},\theta\right). \tag{110}$$

*The density is obtained using the inverse Fourier transform of the characteristic function g*(*<sup>x</sup>*, *α*, *θ*, *λ*) = 1 2*π* ∞−∞ *eitx g*<sup>ˆ</sup>(*<sup>t</sup>*, *α*, *θ*, *λ*)*dt. Using now the relation in Equation (110) and changing the integration variable λ*1/*<sup>α</sup>t* = *τ we obtain the relation for densities.*

$$\mathbb{S}(\mathbf{x}, \mathfrak{a}, \theta, \lambda) = \lambda^{-1/a} \mathbb{S}(\mathbf{x}\lambda^{-1/a}, \mathfrak{a}, \theta). \tag{111}$$

Thus, the Theorem 1 defines an integral representation for the probability density of a standard strictly stable law with the characteristic function in Equation (1). Using this representation, we can obtain an integral representation for the distribution function of the standard strictly stable law. We formulate this result as a corollary to the Theorem 1.

**Corollary 1.** *The distribution function of a stable law <sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*) *with the characteristic function in Equation (1) can be represented as*

*1. If α* -= 1*, then for any* |*θ*| min(1, 2/*α* − 1) *and x* -= 0

$$G(\mathbf{x}, \mathbf{a}, \theta) = \frac{1}{2} (1 - \text{sign}\,(\mathbf{x})) + \text{sign}\,(\mathbf{x}) G^{(+)} (|\mathbf{x}|, \mathbf{a}, \theta^\*), \tag{112}$$

*where θ*∗ = *θ* sign (*x*)*,*

$$\begin{split} G^{(+)}\left(\mathbf{x},a,\theta\right) &= 1 - \frac{(1+\theta)}{4}(1+\operatorname{sign}\left(1-a\right)) \\ &\quad + \frac{\operatorname{sign}\left(1-a\right)}{\pi} \int\_{-\pi\theta/2}^{\pi/2} \exp\left\{-\mathbf{x}^{a/(a-1)}\mathcal{U}(\boldsymbol{\varphi},a,\theta)\right\} d\boldsymbol{\varphi}, \quad \mathbf{x}>0, \end{split} \tag{113}$$

*and <sup>U</sup>*(*ϕ*, *α*, *θ*) *is defined by Equation (70).*

*2. If α* = 1*, then for any* −1 *θ* 1 *and any x*

$$G(\mathbf{x}, 1, \theta) = \frac{1}{2} + \frac{1}{\pi} \arctan\left(\frac{\mathbf{x} - \sin(\pi \theta / 2)}{\cos(\pi \theta / 2)}\right). \tag{114}$$

*3. If x* = 0*, then for any admissible α and θ*

$$G(0, a, \theta) = (1 - \theta) / 2. \tag{115}$$

**Proof.** We consider the case *α* -= 1. It is necessary to obtain the distribution function of a stable law with the characteristic function in Equation (1). For the density of distribution a stable law, one should choose the expression in Equation (89) which defines the density *g*(*<sup>x</sup>*, *α*, *θ*) for *x* > 0. In view of this, we write the distribution function in the form

$$G^{(+)}({\bf x}, {\bf a}, \theta) = 1 - \int\_{{\bf x}}^{\infty} \g(\xi', {\bf a}, \theta) d\xi', \quad {\bf x} > 0.$$

> Now substituting here Equation (89) we find

$$\begin{split} G^{(+)}(x,a,\theta) &= 1 - \frac{\pi}{\pi|1-a|} \int\_{-\pi\theta/2}^{\pi/2} \mathcal{U}(\rho,a,\theta) d\rho \int\_{x}^{\infty} \xi^{1/(a-1)} \exp\left\{-\xi^{\mu/(a-1)}\mathcal{U}(\rho,a,\theta)\right\} d\xi \\ &= \left\{1 + \frac{\text{sign}\,(a-1)}{\pi} \int\_{-\pi\theta/2}^{\pi/2} \left(1 - \exp\left\{-\mathbf{x}^{a/(a-1)}\mathcal{U}(\rho,a,\theta)\right\}\right) d\rho, & a < 1, \, x > 0, \\ 1 - \frac{\text{sign}\,(a-1)}{\pi} \int\_{-\pi\theta/2}^{\pi/2} \exp\left\{-\mathbf{x}^{a/(a-1)}\mathcal{U}(\rho,a,\theta)\right\} d\rho, & a > 1, \, x > 0. \end{split}$$
 
$$= \left\{1 - \frac{(1+\theta)}{2} + \frac{1}{\pi} \int\_{-\pi\theta/2}^{\pi/2} \exp\left\{-\mathbf{x}^{a/(a-1)}\mathcal{U}(\rho,a,\theta)\right\} d\rho, & a < 1, \, x > 0, \\ 1 - \frac{1}{\pi} \int\_{-\pi\theta/2}^{\pi/2} \exp\left\{-\mathbf{x}^{a/(a-1)}\mathcal{U}(\rho,a,\theta)\right\} d\rho, & a > 1, \, x > 0. \right\}$$

If we combine the cases *α* < 1 and *α* > 1 we obtain

$$\mathcal{G}^{(+)}({\bf x},a,\theta) = 1 - \frac{(1+\theta)}{4}(1 + \text{sign}\,(1-a)) + \frac{\text{sign}\,(1-a)}{\pi} \int\_{-\pi\theta/2}^{\pi/2} \exp\left\{-\mathbf{x}^{a/(a-1)}\mathcal{U}(\boldsymbol{\varphi},a,\theta)\right\}d\boldsymbol{\rho}\dots$$

This formula defines the distribution function of stable law for the case *x* > 0 and *α* -= 1. The case *x* < 0 is reduced to the case *x* > 0 using the property of an inversion, namely, the formula in Equation (7) for *<sup>G</sup>*(*<sup>x</sup>*, *α*, *<sup>θ</sup>*). As a result, we ge<sup>t</sup> *<sup>G</sup>*(−)(−*x*, *α*, *θ*) = 1 − *<sup>G</sup>*(+)(*<sup>x</sup>*, *α*, <sup>−</sup>*<sup>θ</sup>*), *x* > 0. This formula gives the distribution function for negative *x*. Combining the formulas for *<sup>G</sup>*(+)(*<sup>x</sup>*, *α*, *θ*) and *<sup>G</sup>*(−)(−*x*, *α*, *<sup>θ</sup>*), we obtain the formula in Equation (112) which is true for any *x* -= 0 and *α* -= 1.

We now consider the case *α* = 1. According to the definition *<sup>G</sup>*(*<sup>x</sup>*, 1, *θ*) = *x*−∞ *g*(*ξ*, 1, *<sup>θ</sup>*)*dξ*. Substituting the density (72) here and replacing the variable *ξ* − sin(*πθ*/2) = *τ* we obtain

$$\begin{split} G(\mathbf{x},1,\boldsymbol{\theta}) &= \frac{\cos(\pi\boldsymbol{\theta}/2)}{\pi} \int\_{-\infty}^{\infty} \frac{d\boldsymbol{\xi}}{\overline{\boldsymbol{\xi}}^{2} - 2\underline{\boldsymbol{x}}\sin(\pi\boldsymbol{\theta}/2) + 1} \\ &= \frac{\cos(\pi\boldsymbol{\theta}/2)}{\pi} \int\_{-\infty}^{\infty - \sin(\pi\boldsymbol{\theta}/2)} \frac{d\boldsymbol{\tau}}{\overline{\boldsymbol{\tau}}^{2} + \cos^{2}(\pi\boldsymbol{\theta}/2)} = \frac{1}{2} + \frac{1}{\pi} \arctan\left(\frac{\mathbf{x} - \sin(\pi\boldsymbol{\theta}/2)}{\cos(\pi\boldsymbol{\theta}/2)}\right). \end{split}$$

Thus, the second item of the corollary has been proved.

To calculate *<sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*) at *x* = 0 we use the formula in Equation (113). Performing the passage to the limit *x* → 0 in this expression we ge<sup>t</sup> *G*(0, *α*, *θ*)=(1 − *θ*)/2. Since the formula in Equation (113) is valid for *α* -= 1, then the result obtained is valid only for *α* -= 1. To calculate *G*(0, 1, *θ*) it is necessary to use (114). Substituting the value *x* = 0 in (114) we obtain *G*(0, 1, *θ*)=(1 − *θ*)/2. Thus, the expression (115) is true for any admissible *α* and *θ*.

We make some remarks on the proved corollary.

**Remark 6.** *In the Remark 3 it is emphasized that in the case α* = 1, *θ* = ±1 *the density g*(*<sup>x</sup>*, 1, ±1) *is a degenerate distribution in the point x* = ±1 *and has the form of Equation (108). Consequently, for the indicated parameter values, the distribution function <sup>G</sup>*(*<sup>x</sup>*, 1, ±1) *will have the form of the Heaviside function <sup>G</sup>*(*<sup>x</sup>*, 1, ±1) = *<sup>H</sup>*(*x* ∓ 1)*. This is directly seen from the form of the distribution function at α* = 1*. Indeed, substituting the values θ* = ±1 *in Equation (114), we obtain*

$$G(\mathbf{x}, 1, \pm 1) = \frac{1}{2} + \frac{1}{\pi} \arctan\left(\frac{\mathbf{x} \mp \sin(\pi/2)}{\cos(\pi/2)}\right) = \begin{cases} 1, & \mathbf{x} > \pm 1, \\ 0, & \mathbf{x} < \pm 1 \end{cases} = H(\mathbf{x} \mp 1).$$

**Remark 7.** *The proved corollary gives an integral representation for the distribution function of a standard strictly stable law. In order to get the distribution function of a strictly stable law for an arbitrary λ it is necessary to use the Remark 5. By definition <sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*, *λ*) = *x* −∞ *g*(*y*, *α*, *θ*, *λ*)*dy. Using now the relation in Equation (111) and changing the variable xλ*−1/*<sup>α</sup>* = *τ we arrive at the relation*

$$G(\mathfrak{x}, \mathfrak{a}, \theta, \lambda) = G\left(\mathfrak{x}\lambda^{-1/\alpha}, \mathfrak{a}, \theta\right).$$

It should be noted the density in Equation (72) and distribution function in Equation (114) is not new. These formulas were deduced by V.M. Zolotarev (see § 2.3 in [32]). As we can see, this distribution is the stable law for *α* = 1 and any −1 *θ* 1 and it is expressed in terms of elementary functions.

#### **4. The Calculation of the Density and Distribution Function of a Stable Law**

Integral representations for the probability density and distribution function of a stable law with the characteristic function in Equation (1) were obtained in the Theorem 1 and Corollary 1. These integral representations in Equations (69) and (112) express the probability density and the distribution function in terms of a definite integral. That is why, using the methods of numerical integration, it is possible to calculate the values of these integrals without much difficulty.

In this paper, to calculate definite integrals in Equations (69) and (112) we used the adaptive Gaussian–Kronrod numerical integration algorithm for 31 points. To implement the program for calculating the functions *g*(*<sup>x</sup>*, *α*, *θ*) and *<sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*) we used the implementation of this algorithm in the gsl library (GNU Scientific Library) of version 1.8 [46]. The calculation results for the functions *g*(*<sup>x</sup>*, *α*, *θ*) and *<sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*) by the Equations (69) and (112) are given in Figures 9–14.

The figures show the probability density and distribution function for the values of the characteristic exponent *α* = 0.3, 0.6, 0.9, 1.2, 1.5, 1.8 and specified values of the parameter *θ* and *λ* = 1. It should be noted that the admissible region of the parameter *θ* is determined by the inequality |*θ*| min(1, 2/*α* − <sup>1</sup>). Thus, −1 *θ* 1, if 0 < *α* 1, and −(2/*<sup>α</sup>* − 1) *θ* 2/*α* − 1, if 1 < *α* 2. It is clear that if *α* > 1 then the admissible region of the parameter *θ* narrows and at *α* = 2 the parameter *θ* may take a single value *θ* = 0.

Let us analyze the results presented in more detail. We first consider the case *α* < 1. The results related to this case are given in Figures 9–11. From these figures it can be seen that when *θ* = 1 the probability density is concentrated on the positive semiaxis. Thus, *g*(*<sup>x</sup>*, *α*, 1) = 0, *<sup>G</sup>*(*<sup>x</sup>*, *α*, 1) = 0, if *x* < 0, *α* < 1. Similarly, for the case *θ* = −1 and *α* < 1, we obtain that the negative semiaxis is the area of concentration of the probability density. Consequently, *g*(*<sup>x</sup>*, *α*, −<sup>1</sup>) = 0, *<sup>G</sup>*(*<sup>x</sup>*, *α*, −<sup>1</sup>) = 1, if *x* > 0, *α* < 1. This result is in complete agreemen<sup>t</sup> with remarks 3 and 4 on page 79 of theorem 2.2.3 from the book by [32].

In the introduction, it was noted that in the work [31] (see also [32]) integral representations were obtained for the probability density and distribution function of a stable law with the characteristic function in Equation (5). In order to avoid any confusion the parameters of a stable law with the characteristic function in Equation (5) will be designated as (*α* , *β*, *λ* , *<sup>γ</sup>*). For the parameters of a strictly stable law with the characteristic function in Equation (1) we keep the notation (*<sup>α</sup>*, *θ*, *<sup>λ</sup>*). The parameters (*α* , *β*, *λ* , *γ*) are related to the parameters (*<sup>α</sup>*, *θ*, *λ*) by the relations (see [32,36]): *α* = *α*,

$$\begin{aligned} \theta &= \beta K(\alpha')/\alpha', & \lambda &= \lambda', & \text{if} \quad \alpha \neq 1, \\ \theta &= (2/\pi) \arctan(2\gamma/\pi), & \lambda &= \lambda' \left(\pi^2/4 + \gamma^2\right)^{1/2}, & \text{if} \quad \alpha = 1. \end{aligned}$$

In the case of *α* -= 1 from the relation for the parameters *θ* and *β* we ge<sup>t</sup> *θ* = *β*, if *α* < 1, and *θ* = *β*(1 − 2/*α*), if *α* > 1.

*Mathematics* **2020**, *8*, 775

**Figure 9.** The probability density function *g*(*<sup>x</sup>*, *α*, *θ*) (**on the left**) and cumulative distribution function *<sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*) (**on the right**) of a standard strictly stable law with *α* = 0.3 and specified values *θ*.

**Figure 10.** The probability density function *g*(*<sup>x</sup>*, *α*, *θ*) (**on the left**) and cumulative distribution function *<sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*) (**on the right**) of a standard strictly stable law with *α* = 0.6 and specified values *θ*.

**Figure 11.** The probability density function *g*(*<sup>x</sup>*, *α*, *θ*) (**on the left**) and cumulative distribution function *<sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*) (**on the right**) of a standard strictly stable law with *α* = 0.9 and specified values *θ*.

*Mathematics* **2020**, *8*, 775

**Figure 12.** The probability density function *g*(*<sup>x</sup>*, *α*, *θ*) (**on the left**) ) and cumulative distribution function *<sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*) (**on the right**) of a standard strictly stable law with *α* = 1.2 and specified values *θ*.

**Figure 13.** The probability density function *g*(*<sup>x</sup>*, *α*, *θ*) (**on the left**) ) and cumulative distribution function *<sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*) (**on the right**) of a standard strictly stable law with *α* = 1.5 and specified values *θ*.

**Figure 14.** The probability density function *g*(*<sup>x</sup>*, *α*, *θ*) (**on the left**) and cumulative distribution function *<sup>G</sup>*(*<sup>x</sup>*, *α*, *θ*) (**on the right**) of a standard strictly stable law with *α* = 1.8 and specified values *θ*.

Thus, at *α* < 1 the parameters *θ* and *β* coincide. Therefore, the corresponding properties for the probability density and the distribution function of stable laws with characteristic functions in Equations (1) and (5) coincide. The situation slightly changes if *α* > 1. It can be seen from the above relation that, firstly, the admitted region of the parameter *θ* narrows in comparison with the admitted region of the parameter *β*, secondly, the parameter *θ* changes its sign to the opposite with respect to the parameter *β*. Since the parameter *θ* has the meaning of an asymmetry parameter, then a change of the sign of this parameter when passing from *α* < 1 to *α* > 1 will affect the form of probability density. This is clearly seen in Figures 11 and 12 for densities with extreme values of the parameter *θ*. Comparing densities for the values *α* = 0.9, *θ* = 1 given in Figure 11 and density for the values *α* = 1.2, *θ* = 0.66 given in Figure 12 one can see that these densities are turned into different directions. This fact is a consequence of the fact that the parameter *θ* changed the sign compared to the sign of the parameter *β*. Thus, for the same sign of the parameter *θ* of density *g*(*<sup>x</sup>*, *α*, *θ*) at *α* < 1 and *α* > 1 will be turned into different directions. The reason for this behavior is related to the selected parameter system (*<sup>α</sup>*, *θ*, *λ*) of the characteristic function in Equation (1). As it was mentioned in the book by Zolotarev V.M. (see page 19 in [32]), distributions from the class of strictly stable laws are continuous in the totality of their parameters in the entire range of their admissible values it is with this choice of parameters.

It should be pointed out that the Theorem 1 and Corollary 1 formulate expressions for probability density and distribution functions for strict stable laws with a scale parameter *λ* = 1. To obtain the density and distribution function with an arbitrary value of the scale parameter *λ* it is necessary to use Remarks 5 and 7.
