Next Article in Journal
One Note for Fractionation and Increase for Mixed-Level Designs When the Levels Are Not Multiple
Next Article in Special Issue
Fractional Calculus in Russia at the End of XIX Century
Previous Article in Journal
Simultaneous Contour Method of a Trigonometric Integral by Prudnikov et al.
Previous Article in Special Issue
Series in Le Roy Type Functions: A Set of Results in the Complex Plane—A Survey
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Asymptotic Expansion of a Function Introduced by L.L. Karasheva

Division of Computing and Mathematics, Abertay University, Dundee DD1 1HG, UK
Mathematics 2021, 9(12), 1454; https://doi.org/10.3390/math9121454
Submission received: 30 April 2021 / Revised: 15 June 2021 / Accepted: 17 June 2021 / Published: 21 June 2021
(This article belongs to the Special Issue Special Functions with Applications to Mathematical Physics)

Abstract

:
The asymptotic expansion for x ± of the entire function F n , σ ( x ; μ ) = k = 0 sin ( n γ k ) sin γ k x k k ! Γ ( μ σ k ) , γ k = ( k + 1 ) π 2 n for μ > 0 , 0 < σ < 1 and n = 1 , 2 , is considered. In the special case σ = α / ( 2 n ) , with 0 < α < 1 , this function was recently introduced by L.L. Karasheva (J. Math. Sciences, 250 (2020) 753–759) as a solution of a fractional-order partial differential equation. By expressing F n , σ ( x ; μ ) as a finite sum of Wright functions, we employ the standard asymptotics of integral functions of hypergeometric type to determine its asymptotic expansion. This was found to depend critically on the parameter σ (and to a lesser extent on the integer n). Numerical results are presented to illustrate the accuracy of the different expansions obtained.
MSC:
33C70; 34E05; 41A30; 41A60

1. Introduction

In a recent paper, L.L. Karasheva [1] introduced the entire function
Θ n , α ( x ; μ ) : = k = 0 sin ( n γ k ) sin γ k x k k ! Γ ( μ α k 2 n ) , γ k : = ( k + 1 ) π 2 n ,
where μ > 0 , 0 < α < 1 and n = 1 , 2 , and, throughout, x is a real variable. This function is of interest as it is involved in the fundamental solution of the differential equation
α u t α + ( 1 ) n 2 n u x 2 n = f ( x , t )
for positive integer n, where the derivative with respect to t is the fractional derivative of the order α . In the simplest case n = 1 , we have Θ 1 , α ( x ; μ ) = ϕ ( σ , μ ; x ) , σ : = α / ( 2 n ) , where ϕ ( σ , μ ; x ) is the Wright function
ϕ ( σ , μ ; x ) : = k = 0 x k k ! Γ ( μ σ k ) ( σ < 1 ) ,
which finds application as a fundamental solution of the diffusion-wave equation [2]. Under the above assumptions on n and α it follows that the parameter σ associated with (1) satisfies 0 < σ < 1 2 .
In this study, however, we shall allow the parameter σ to satisfy 0 < σ < 1 and consider the function
F n , σ ( x ; μ ) : = k = 0 sin ( n γ k ) sin γ k x k k ! Γ ( μ σ k ) ( 0 < σ < 1 ) ,
which coincides with Θ n , α ( x ; μ ) when σ = α / ( 2 n ) . From the well-known expansion
sin ( n γ k ) sin γ k = r = 0 n 1 e i γ k ( 2 r n + 1 ) = r = 0 n 1 e i ( k + 1 ) ω r ,
where
ω r : = ( n 2 r 1 ) π 2 n ( 0 r n 1 ) ,
it follows that (3) can be expressed as a finite sum of Wright functions defined in (2) with rotated arguments (compare [1], Equation (4))
F n , σ ( x ; μ ) = r = 0 n 1 e i ω r ϕ ( σ , μ ; x e i ω r ) .
We note that the extreme values of ω r satisfy ω 0 = ω n 1 = ( n 1 ) π / ( 2 n ) , whence | ω r | < 1 2 π for 0 r n 1 .
We use the representation in (5), with the values of ω r in (4), to determine the asymptotic expansion of F n , σ ( x ; μ ) for x ± by application of the asymptotic theory of the Wright function. A summary of the expansion of ϕ ( σ , μ ; z ) for large | z | is given in Section 3. The expansions of F n , σ ( x ; μ ) for x ± are given in Section 4 and Section 5, where they are shown to depend critically on the parameter σ (and to a lesser extent on the integer n). A concluding section presents our numerical results confirming the accuracy of the different expansions obtained.

2. An Alternative Representation of F n , σ ( x ; μ )

The Wright function appearing in (2) can be written alternatively as
ϕ ( σ , μ ; x ) = 1 π k = 0 x k k ! Γ ( 1 μ + σ k ) sin π ( μ σ k )
= 1 2 π e π i ϑ Ψ ( x e π i σ ) + e π i ϑ Ψ ( x e π i σ )
upon use of the reflection formula for the gamma function, where ϑ : = 1 2 μ . The associated Wright function Ψ ( z ) is defined by
Ψ ( z ) : = k = 0 z k k ! Γ ( σ k + δ ) ( 0 < σ < 1 , δ = 1 μ ) ,
which is valid for | z | < . Hence, we obtain the representation
F n , σ ( x ; μ ) = 1 2 π r = 0 n 1 e i ω r Υ r ( σ ; x ) ,
where
Υ r ( σ ; x ) : = e π i ϑ Ψ ( x e π i σ i ω r ) + e π i ϑ Ψ ( x e π i σ i ω r ) .
If we now exploit the symmetry of the ω r in (4) (and the fact that x is a real variable), we observe that the values of ω r for 0 r N 1 , where N = n / 2 , satisfy
{ ω 0 , ω 1 , , ω N 1 } = ( n 1 ) π 2 n , ( n 3 ) π 2 n , , π 2 n ϵ n , ϵ n = 1 ( n even ) 2 ( n odd ) .
Then, we can write
F n , σ ( x ; μ ) = 1 π r = 0 N 1 e i ω r Υ r ( σ ; x ) + Δ n e π i ϑ Ψ ( x e π i σ ) ,
where
Δ n = 0 ( n even ) 1 ( n odd ) .
The form (8) involves half the number of Wright functions Ψ ( z ) and will be used to determine the asymptotic expansion of F n , σ ( x ; μ ) as x ± in Section 4 and Section 5.

3. The Asymptotic Expansion of Ψ ( z ) for | z |

We first present the large- | z | asymptotics of the function Ψ ( z ) in (6) based on the presentation described in ([3], Section 4); see also ([4], Section 4.2), ([5], §2.3). We introduce the following parameters:
κ = 1 σ , h = σ σ , ϑ = δ 1 2 , δ = 1 μ ,
together with the associated (formal) exponential and algebraic expansions
E ( z ) : = Z ϑ e Z j = 0 A j ( σ ) Z j , H ( z ) : = 1 σ k = 0 ( 1 ) k k ! Γ k + δ σ z ( k + δ ) / σ ,
where (The dependence of the coefficients A j ( σ ) on the parameter δ is not indicated.)
Z : = κ ( h z ) 1 / κ , A 0 ( σ ) = ( 2 π / κ ) 1 / 2 ( σ / κ ) ϑ .
Then, since 0 < κ < 1 , we obtain from ([5], p. 57) the large-z expansion
Ψ ( z ) E ( z ) + H ( z e π i ) ( | arg z | 1 2 π κ ) H ( z e π i ) ( 1 2 π κ < | arg z | π ) ,
where the upper or lower signs are chosen according as arg z > 0 or arg z < 0 , respectively.
The expansion E ( z ) is exponentially large as | z | in the sector | arg z | < 1 2 π κ , and oscillatory (multiplied by the algebraic factor z ϑ / κ ) on the anti-Stokes lines arg z = ± 1 2 π κ . In the adjacent sectors 1 2 π κ < | arg z | < π κ , the expansion E ( z ) continues to be present, but is exponentially small reaching maximal subdominance relative to the algebraic expansion on the Stokes lines (On these rays, E ( z ) undergoes a Stokes phenomenon where it switches off in a smooth manner (see [6], p. 67).) arg z = ± π κ . In our treatment of F n , σ ( x ; μ ) , we will not be concerned with exponentially small contributions, except in one special case when x where the expansion of F n , σ ( x ; μ ) is exponentially small.
The first few normalised coefficients c j = A j ( σ ) / A 0 ( σ ) are [3,4]:
c 0 = 1 , c 1 = 1 24 σ { 2 + 7 σ + 2 σ 2 12 δ ( 1 + σ ) + 12 δ 2 } , c 2 = 1 1152 σ 2 { 4 + 172 σ + 417 σ 2 + 172 σ 3 + 4 σ 4 24 δ ( 6 + 41 σ + 41 σ 2 + 6 σ 3 ) + 120 δ 2 ( 4 + 11 σ + 4 σ 2 ) 480 δ 3 ( 1 + σ ) + 144 δ 4 } , c 3 = 1 414 , 720 σ 3 { ( 1112 + 9636 σ + 163 , 734 σ 2 + 336 , 347 σ 3 + 163 , 734 σ 4 + 9636 σ 5 1112 σ 6 ) δ ( 3600 + 220 , 320 σ + 929 , 700 σ 2 + 929 , 700 σ 3 + 220 , 320 σ 4 + 3600 σ 5 ) + δ 2 ( 65 , 520 + 715 , 680 σ + 1 , 440 , 180 σ 2 + 715 , 680 σ 3 + 65 , 520 σ 4 ) δ 3 ( 161 , 280 + 816 , 480 σ + 816 , 480 σ 2 + 161 , 280 σ 3 ) + δ 4 ( 151 , 200 + 378 , 000 σ + 151 , 200 σ 2 ) 60 , 480 δ 5 ( 1 + σ ) + 8640 δ 6 } .
In addition to the Stokes lines arg z = ± π κ , where E ( z ) is maximally subdominant relative to the algebraic expansion, the positive real axis is also a Stokes line. Here, the algebraic expansion is maximally subdominant relative to E ( z ) . As the positive real axis is crossed from the upper to the lower half plane the factor e π i appearing in H ( z e π i ) changes to e π i , and vice versa. The details of this transition will not be considered here; see ([5], p. 248) for the case of the confluent hypergeometric function 1 F 1 ( a ; b ; z ) .

4. The Asymptotic Expansion of F n , σ ( x ; μ ) for x +

4.1. Asymptotic Character as a Function of σ

Let us denote the arguments of the Ψ functions appearing in (8) by
z r ± = x exp [ i ϕ r ± ] , ϕ r ± = ± π σ ω r .
The representation of the asymptotic structure of the functions Ψ ( z r ± ) is illustrated in Figure 1 for different values of σ . The figures show the rays arg z = ± π σ and the anti-Stokes lines (dashed lines) arg z = ± 1 2 π κ . In the case σ = 2 3 , the exponentially large sector is | arg z | < 1 6 π , and it is seen from Figure a that the arguments z r ± for 0 r N 1 and x e ± π i σ all lie in the domain where Ψ ( z ) has an algebraic expansion; this conclusion applies a fortiori when 2 3 < σ < 1 . When σ = 1 2 , the exponentially large sector is | arg z | < 1 4 π ; when n = 2 , we have ω 0 = 1 4 π so that z 0 + is situated on the boundary of the exponentially large sector.
Other values of n 3 will have some z r + inside this sector, whereas the z r are in the algebraic sector for n 2 . Similarly, the case σ = 1 3 , where the rays arg z = ± π σ and arg z = ± 1 2 π κ coincide, has all the z r + situated in the exponentially large sector, with the z r situated in the algebraic domain. Finally, when σ = 1 6 , the exponentially large sector | arg z | < 5 12 π encloses the rays arg z = ± π σ with the result that all the z r + lie in the exponentially large sector, whereas the z r lie in the algebraic domain (except when n = 2 when z 0 lies on the lower boundary of the exponentially large sector).
To summarise, we have the following asymptotic character of F n , σ ( x ; μ ) when x + as a function of the parameter σ :
0 < σ < 1 2 Exp . large + Algebraic ( for n 2 ) 1 2 σ < 2 3 Exp , large ( dependent on n ) + Algebraic 2 3 σ < 1 Algrbraic ( for n 2 ) .

4.2. Asymptotic Expansion

From (8) and (10), we have the algebraic expansion associated with F n , σ ( x ; μ ) given by
H ( x ) = 1 σ k = 0 x K k ! Γ ( 1 K ) θ n , k , K : = k + δ σ ,
where, with appropriate choices of the factors e ± π i in H ( z ) ,
θ n , k = ( 1 ) k sin π K { r = 0 N 1 e π i ϑ i ω r ( e π i σ i ω r · e π i ) K + e π i ϑ i ω r ( e π i σ i ω r · e π i ) K + Δ n e π i ϑ ( e π i σ · e π i ) K } = ( 1 ) k sin π K r = 0 N 1 e ( K 1 ) i ω r ( e π i ( ϑ + κ K ) + e π i ( ϑ + κ K ) ) + Δ n e π i ( ϑ + κ K ) = 2 r = 0 N 1 e ( K 1 ) i ω r + Δ n ,
as cos π ( ϑ + κ K ) = cos π ( K k 1 2 ) = ( 1 ) k sin π K .
For the exponential component, we introduce the quantities
X = κ ( h x ) 1 / κ , Φ r ± = ± π ϑ κ ω r 1 + ϑ κ
and the formal asymptotic sum
S ( X e i Ω ) : = j = 0 A j ( σ ) ( X e i Ω / κ ) j .
Then, from (8) and (10), we have the exponential expansion in the form
E ( x ) = X ϑ π { r = 0 N 1 exp [ X e i ϕ r + / κ + i Φ r + ] S ( X e i ϕ r + ) + exp [ X e i ϕ r / κ + i Φ r ] S ( X e i ϕ r )
+ Δ n exp [ X e π i σ / κ + π i ϑ / κ ] S ( X e π i σ ) } .
It is important to stress that only the exponential terms with | ϕ r ± | 1 2 π κ , that is those with
| ± π σ ω r | 1 2 π κ ,
are to be retained in E ( x ) in (19). In addition, it is seen by inspection of Figure 1 that the second term involving S ( X e i ϕ r ) does not contribute to E ( x ) when 1 3 σ < 1 , since, for this range of σ , the ray arg z = π σ lies outside (or, when σ = 1 3 , on the lower boundary of) the exponentially large sector | arg z | < 1 2 π κ . Thus, when 1 2 σ < 2 3 , the exponential expansion is significant if π σ ω 0 1 2 π κ ; that is, if n n 0 = 1 / ( 2 3 σ ) .
In summary, we have the following theorem.
Theorem 1.
The following expansion holds for x + :
F n , σ ( x ; μ ) E ( x ) + H ( x ) ( 0 < σ < 1 2 ; n 2 ) E ( x ) + H ( x ) ( 1 2 σ < 2 3 ; n n 0 ) H ( x ) ( 1 2 σ < 2 3 ; n < n 0 ) H ( x ) ( 2 3 σ < 1 ; n 2 ) ,
where n 0 = 1 / ( 2 3 σ ) and the exponential and algebraic expansions E ( x ) and H ( x ) are defined in (15) and (19).

4.3. Karasheva’s Estimate for | Θ n , α ( x ; μ ) |

When σ = α / ( 2 n ) < 1 2 , we see from Theorem 1 that the dominant exponential expansion as x + corresponds to r = 0 , yielding
Θ n , α ( x ; μ ) A 0 ( σ ) X ϑ π exp [ X e i ( π σ ω 0 ) / κ + i Φ 0 + ]
= A 0 ( σ ) X ϑ π exp [ X cos ( π σ ω 0 ) / κ ) ] cos [ X sin ( π σ ω 0 ) / κ ) + Φ 0 + ] ,
where
π σ ω 0 κ = 2 n π σ ( n 1 ) π 2 n α = ( α + 1 n ) π 2 n α .
Thus, we have the leading order estimate
Θ n , α ( x ; μ ) A 0 ( σ ) X ϑ π exp X cos ( n 1 α ) π 2 n α cos X sin ( n 1 α ) π 2 n α Φ 0 +
as x + . When expressed in our notation, Karasheva’s estimate for | Θ n , α ( x ; μ ) | in ([1], §8) agrees with (20) (when the second cosine term is replaced by 1), except that she did not give the value of the multiplicative constant A 0 ( σ ) / π given in (11). However, the presentation of her result as an upper bound is not evident due to the presence of possibly less dominant exponential expansions and also the subdominant algebraic expansion.

5. The Expansion of F n , σ ( x ; μ ) for x

To examine the case of negative x, we replace x by e π i x , with x > 0 , and use the fact that Ψ ( z e 2 π i ) = Ψ ( z ) to find, from (8), that
F n , σ ( x ; μ ) = 1 π r = 0 N 1 e i ω r Υ r ( κ ; x ) + Δ n e π i ϑ Ψ ( x e π i κ ) .
The rays arg z = ± π σ in Figure 1 are now replaced by the Stokes lines arg z = ± π κ . The Stokes and anti-Stokes lines arg z = ± 1 2 π κ are illustrated in Figure 2 when 0 < σ < 1 2 and 1 2 < σ < 1 . In the sectors 1 2 π κ < | arg z | < π κ , we recall that the exponential expansion E ( z ) is still present but is exponentially small as | z | .
For the algebraic component of the expansion two cases arise when the argument π κ ω r of the second Ψ function in Υ r ( κ ; x ) is either (i) positive or (ii) negative. In case (i), the algebraic expansion H ( z ) does not encounter a Stokes phenomenon as its argument does not cross arg z = 0 , whereas in case (ii), a Stokes phenomenon arises for those values of r that make π κ ω r < 0 . In case (i), the algebraic component contains the factor inside the sum over r in (21)
e π i ϑ ( e π i κ i ω r · e π i ) K + e π i ϑ ( e π i κ i ω r · e π i ) K
= e i ω r K ( e π i ( ϑ σ K ) + e π i ( ϑ σ K ) ) = 2 e i ω r K cos π ( k + 1 2 ) 0
upon recalling the definition of K in (15) and noting that δ ϑ = 1 2 . Similarly, the final term involves the factor e π i ϑ ( e π i κ · e π i ) K = cos π ( ϑ σ K ) = 0 . Thus, the algebraic contribution to F n , σ ( x ; μ ) vanishes in case (i).
For case (ii) to apply, we require that π κ ω 0 < 0 ; that is, n > n * = 1 / ( 2 σ 1 ) . Suppose that π κ ω r < 0 for 0 r r 0 . Then, the algebraic component resulting from the terms with r r 0 becomes
1 π σ k = 0 ( 1 ) k Γ ( K ) k ! x K r = 0 r 0 e ( K 1 ) i ω r e π i ϑ ( e π i κ · e π i ) K + e π i ϑ ( e π i κ · e π i ) K
= 2 π σ k = 0 ( 1 ) k Γ ( K ) k ! x K r = 0 r 0 e ( K 1 ) i ω r π i K cos π ( ϑ σ K + π K ) ,
where, in the second term in round braces, we have taken account of the Stokes phenomenon (the first term and that multiplied by Δ n are unaffected). Some routine algebra then produces the algebraic contribution
H ^ ( x ) : = 2 σ k = 0 x K k ! Γ ( 1 K ) θ ^ n , k , θ ^ n , k : = r = 0 r 0 cos π K ( K 1 ) ω r
when n > n * and H ^ ( x ) 0 when n < n * . (We avoid here consideration of the algebraic contribution when π κ ω r = 0 , that is, on the Stokes line arg z = 0 .)
Reference to Figure 2 shows that there is no exponential contribution to F n , σ ( x ; μ ) from the terms Ψ ( x e π i κ ) and Ψ ( x e π i κ i ω r ) . From (10) and (21), we find the exponential expansion results from the terms Ψ ( x e π i κ i ω r ) , which is given by
E ^ ( x ) : = X ϑ π r = 0 N 1 exp [ X e i ω r / κ i Φ ] S ( X e i ω r / κ ) ,
where X and the asymptotic sum S are defined in (17) and (18) with Φ : = ω r ( 1 + ϑ / κ ) . For σ < 1 2 (when the algebraic expansion vanishes), the expansion of F n , σ ( x ; μ ) will be exponentially small provided π κ ω 0 > 1 2 π κ ; that is, when n < 1 / σ . If n = 1 / σ , there is an exponentially oscillatory contribution, and when n > 1 / σ , the expansion is exponentially large.
To summarise, we have the theorem:
Theorem 2.
The following expansion holds for x + :
F n , σ ( x ; μ ) E ^ ( x ) ( 0 < σ 1 2 ) E ^ ( x ) + H ^ ( x ) ( 1 2 < σ < 1 ) ,
where the exponential expansion E ^ ( x ) is defined in (23). This last expansion is exponentially small as x when 0 < σ < 1 2 and n < 1 / σ . The algebraic expansion H ^ ( x ) is given by
H ^ ( x ) : = 2 σ k = 0 x K k ! Γ ( 1 K ) θ ^ n , k ( n > n * ) , 0 ( n < n * ) ,
where n * = 1 / ( 2 σ 1 ) and K, θ ^ n , k are specified in (15) and (22).

6. Numerical Results

In this section, we describe numerical calculations that support the expansions given in Theorems 1 and 2. The function F n , σ ( x ; μ ) was evaluated using the expression in terms of Wright functions (valid for real x)
F n , σ ( x ; μ ) = 2 r = 0 N 1 e i ω r ϕ ( σ , μ ; x e i ω r ) + Δ n ϕ ( σ , μ ; x ) , N = n / 2 ,
which follows from (5) and the symmetry of ω r .
In Table 1, we present the results of numerical calculations for x + compared with the expansions given in Theorem 1. We choose four representative values of σ that focus on the different cases of Theorem 1 and n = 2 , 3 and 4. The numerical value of F n , σ ( x ; μ ) was obtained by high-precision evaluation of (25). The exponential expansion E ( x ) was computed with the truncation index j = 3 and the algebraic expansion H ( x ) was optimally truncated (that is, at or near its smallest term).
The first case σ = 1 3 has an exponentially large expansion with a subdominant algebraic contribution for all three values of n. The second case σ = 1 2 corresponds to n 0 = 2 ; when n = 2 , E ( x ) is oscillatory and makes a similar contribution as H ( x ) , whereas when n = 3 and 4, E ( x ) is exponentially large. The third case σ = 5 9 corresponds to n 0 = 3 ; when n = 2 , there is no exponential contribution, whereas when n = 3 , E ( x ) is oscillatory and thus makes a similar contribution as H ( x ) ; when n = 4 , E ( x ) is exponentially large. Finally, when σ = 2 3 , the expansion of F n , σ ( x ; μ ) is purely algebraic in character.
In Table 2, we present illustrative examples of Theorem 2 when x . The first case, σ = 1 4 ( κ = 3 4 ), has an expansion that is exponential in character; for n < 1 / σ = 4 , E ^ ( x ) is exponentially small, whereas for n = 4 , the argument π κ ω 0 = 3 8 π lies on the upper boundary of the exponentially large sector | arg z | < 3 8 π , and thus E ^ ( x ) is oscillatory. For n 5 , E ^ ( x ) becomes exponentially large as x . In the second case, σ = 2 5 ( κ = 3 5 ), E ^ ( x ) is exponentially small for n = 2 and exponentially large for n 3 .
In the third case, σ = κ = 1 2 , E ^ ( x ) is oscillatory for n = 2 and exponentially large for n 3 . Finally, when σ = 3 4 ( κ = 1 4 ), the function F n , σ ( x ; μ ) is exponentially large for n = 2 , 3 and n 5 . However, for n = 4 , the two values ω 0 = 3 8 π and ω 1 = 1 8 π yield arguments π κ ω r ( r = 0 , 1 ) situated on both boundaries of the exponentially large sector | arg z | < 1 8 π . In this case E ^ ( x ) is oscillatory and, since n * = 2 , there is, in addition, an algebraic contribution H ^ ( x ) .

7. Concluding Remarks

We employed the standard asymptotics of the Wright function Ψ ( z ) defined in (6) to determine the asymptotic expansion of F n , σ ( x ; μ ) for x ± . We found that this behaviour depended critically on the parameter σ . The numerical results presented in Table 1 and Table 2 demonstrate that the asymptotic forms of F n , σ ( x ; μ ) stated in Theorems 1 and 2 agreed well with the numerically computed values of F n , σ ( ± x ; μ ) . In particular, we showed that, when σ < 1 2 , the expansion of F n , σ ( x ; μ ) exponentially decays as x .

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declares no conflict of interest.

References

  1. Karasheva, L.L. On properties of an entire function that is a generalization of the Wright function. J. Math. Sci. 2020, 250, 753–759. [Google Scholar] [CrossRef]
  2. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity; Imperial College Press: London, UK, 2010. [Google Scholar]
  3. Paris, R.B. The asymptotics of the generalised Bessel function. Math. Aeterna 2017, 7, 381–406. [Google Scholar]
  4. Paris, R.B. Asymptotics of the special functions of fractional calculus. In Handbook of Fractional Calculus with Applications; Kochubei, A., Luchko, Y., Eds.; De Gruyter: Berlin, Germany, 2019; Volume 1, pp. 297–325. [Google Scholar]
  5. Paris, R.B.; Kaminski, D. Asymptotics and Mellin-Barnes Integrals; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  6. Olver, F.W.J.; Lozier, D.W.; Boisvert, R.F.; Clark, C.W. (Eds.) NIST Handbook of Mathematical Functions; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
Figure 1. Diagrams representing the rays arg z = ± π σ and the boundaries of the exponentially large sector (shown by dashed rays) | arg z | < 1 2 π κ , κ = 1 σ for (a) σ = 2 / 3 , (b) σ = 1 / 2 , (c) σ = 1 / 3 , and (d) σ = 1 / 6 . Outside the exponentially large sector, the expansion of Ψ ( z ) is algebraic in character. The circular quadrants represent the range of the arguments arg z = ± π σ ω r for 0 r n / 2 1 , with n 2 and the arrow-head corresponds to n = . When σ = 1 / 3 , the rays arg z = ± π σ and arg z = ± 1 2 π κ coincide.
Figure 1. Diagrams representing the rays arg z = ± π σ and the boundaries of the exponentially large sector (shown by dashed rays) | arg z | < 1 2 π κ , κ = 1 σ for (a) σ = 2 / 3 , (b) σ = 1 / 2 , (c) σ = 1 / 3 , and (d) σ = 1 / 6 . Outside the exponentially large sector, the expansion of Ψ ( z ) is algebraic in character. The circular quadrants represent the range of the arguments arg z = ± π σ ω r for 0 r n / 2 1 , with n 2 and the arrow-head corresponds to n = . When σ = 1 / 3 , the rays arg z = ± π σ and arg z = ± 1 2 π κ coincide.
Mathematics 09 01454 g001aMathematics 09 01454 g001b
Figure 2. Diagrams representing the rays arg z = ± π κ and the boundaries of the exponentially large sector (shown by dashed rays) | arg z | < 1 2 π κ , κ = 1 σ for (a) 0 < σ < 1 2 and (b) 1 2 < σ < 1 . The circular quadrants represent the range of the arguments arg z = ± π κ ω r for 0 r N 1 with the arrow-head corresponding to n = . The ± signs in (b) denote the signs to be chosen in H ( z ) on either side of the Stokes line arg z = 0 .
Figure 2. Diagrams representing the rays arg z = ± π κ and the boundaries of the exponentially large sector (shown by dashed rays) | arg z | < 1 2 π κ , κ = 1 σ for (a) 0 < σ < 1 2 and (b) 1 2 < σ < 1 . The circular quadrants represent the range of the arguments arg z = ± π κ ω r for 0 r N 1 with the arrow-head corresponding to n = . The ± signs in (b) denote the signs to be chosen in H ( z ) on either side of the Stokes line arg z = 0 .
Mathematics 09 01454 g002
Table 1. The values of the exponential and algebraic expansions compared with F n , σ ( x ; μ ) for large x > 0 for different values of σ and n when μ = 3 / 4 and x = 8 .
Table 1. The values of the exponential and algebraic expansions compared with F n , σ ( x ; μ ) for large x > 0 for different values of σ and n when μ = 3 / 4 and x = 8 .
σ n = 2 n = 3 n = 4
1/3 E ( x ) 1.81418881× 10 2 1.08294258× 10 3 3.08231679× 10 3
H ( x ) + 0.34241316 + 0.17280892 + 0.34497729
E ( x ) + H ( x ) 1.81076468× 10 2 1.08276977× 10 3 3.08197181× 10 3
F n , σ ( x ; μ ) 1.80709370× 10 2 1.08284759× 10 3 3.08254767× 10 3
1/2 E ( x ) + 0.06317153 + 1.15957937× 10 3 4.47945373× 10 4
H ( x ) + 0.74012019 + 1.09449277 + 1.45169481
E ( x ) + H ( x ) + 0.80329172 + 1.16067387× 10 3 4.47930856× 10 4
F n , σ ( x ; μ ) + 0.80329527 + 1.16069221× 10 3 4.47921506× 10 4
5/9 E ( x ) 0.14805870 + 2.77243091× 10 2
H ( x ) + 0.79825166 + 1.17615555 + 1.55857242
E ( x ) + H ( x ) + 0.79825166 + 1.02809685 + 2.78801663× 10 2
F n , σ ( x ; μ ) + 0.79825119 + 1.02809649 + 2.78801134× 10 2
2/3 H ( x ) + 0.84046066 + 1.23266920 + 1.63072031
F n , σ ( x ; μ ) + 0.84046066 + 1.23266920 + 1.63072031
Table 2. The values of the exponential and algebraic expansions compared with F n , σ ( x ; μ ) for large x < 0 for different values of σ and n when μ = 3 / 4 and | x | = 8 (for σ = 1 / 4 , 1 / 2 , 2 / 5 ), | x | = 5 (for σ = 3 / 4 ).
Table 2. The values of the exponential and algebraic expansions compared with F n , σ ( x ; μ ) for large x < 0 for different values of σ and n when μ = 3 / 4 and | x | = 8 (for σ = 1 / 4 , 1 / 2 , 2 / 5 ), | x | = 5 (for σ = 3 / 4 ).
σ n = 2 n = 3 n = 4
1/4 E ^ ( x ) + 1.59003829× 10 2 + 1.77442984× 10 1 + 6.49578248× 10 1
F n , σ ( x ; μ ) + 1.59003416× 10 2 + 1.77011100× 10 1 + 6.49580223× 10 1
2/5 E ^ ( x ) 4.18901636× 10 2 3.79446870× 10 0 3.02428770× 10 1
F n , σ ( x ; μ ) 4.18889220× 10 2 3.79475882× 10 0 3.02402120× 10 1
1/2 E ^ ( x ) 0.56022532 + 1.23070020× 10 3 1.28808653× 10 4
F n , σ ( x ; μ ) 0.56023534 + 1.23066913× 10 3 1.28803505× 10 4
3/4 E ^ ( x ) + 1.81213632× 10 28 + 7.55354383× 10 13 0.84956415
H ^ ( x ) 1.93112636× 10 1 0.28756658
E ^ ( x ) + H ^ ( x ) + 1.81213632× 10 28 + 7.55354383× 10 13 1.13713072
F n , σ ( x ; μ ) + 1.81213650× 10 28 + 7.55354314× 10 13 1.13713081
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Paris, R. The Asymptotic Expansion of a Function Introduced by L.L. Karasheva. Mathematics 2021, 9, 1454. https://doi.org/10.3390/math9121454

AMA Style

Paris R. The Asymptotic Expansion of a Function Introduced by L.L. Karasheva. Mathematics. 2021; 9(12):1454. https://doi.org/10.3390/math9121454

Chicago/Turabian Style

Paris, Richard. 2021. "The Asymptotic Expansion of a Function Introduced by L.L. Karasheva" Mathematics 9, no. 12: 1454. https://doi.org/10.3390/math9121454

APA Style

Paris, R. (2021). The Asymptotic Expansion of a Function Introduced by L.L. Karasheva. Mathematics, 9(12), 1454. https://doi.org/10.3390/math9121454

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop