arxiv: v1 [math.ca] 16 Aug 2017

Transcription

1 Some remarks on the theorems of Wright and Braaksma on the Wright function p Ψ q (z) arxiv: v [math.ca] 6 Aug 207 R. B. Paris University of Abertay Dundee, Dundee DD HG, UK Abstract We carry out a numerical investigation of the asymptotic expansion of the so-called Wright function pψ q(z) (a generalised hypergeometric function) in the case when exponentially small terms are present. This situation is covered by two theorems of Wright and Braaksma. We demonstrate that a more precise understanding of the behaviour of pψ q(z) is obtained by taking into account the Stokes phenomenon.. Introduction We consider the Wright function (a generalised hypergeometric function) defined by ( ) (α,a pψ q (z) p Ψ ),...,(α p,a p ) q (β,b ),...,(β q,b q ) ;z = g(n) = n=0 g(n) zn n!, (.) p Γ(α rn+a r ) q Γ(β rn+b r ), (.2) where p and q are nonnegative integers, the parameters α r and β r are real and positive and a r and b r are arbitrary complex numbers. We also assume that the α r and a r are subject to the restriction α r n+a r 0,, 2,... (n = 0,,2,... ; r p) (.3) so that no gamma function in the numerator in (.) is singular. In the special case α r = β r =, the function p Ψ q (z) reduces to a multiple of the ordinary hypergeometric function pψ q (z) = p Γ(a ( ) r) a,...,a q Γ(b r) p F p q ;z ; b,...,b q see, for example, [3, p. 40]. This paper is partly based on the internal report [7].

2 2 R. B. Paris We introduce the parameters associated with g(n) given by κ = + ϑ = q β r p a r p α r, h = p α αr r q β βr r, q b r + 2 (q p), ϑ = ϑ. (.4) If it is supposed that α r and β r are such that κ > 0 then p Ψ q (z) is uniformly and absolutely convergent for all finite z. If κ = 0, the sum in (.) has a finite radius of convergence equal to h, whereas for κ < 0 the sum is divergent for all nonzero values of z. The parameter κ will be found to play a critical role in the asymptotic theory of p Ψ q (z) by determining the sectors in the z-plane in which its behaviour is either exponentially large, algebraic or exponentially small in character as z. The determination of the asymptotic expansion of p Ψ q (z) for z and finite values of the parameters has a long history; for details, see [0, 2.3]. Detailed investigations were carried out by Wright [6, 7] and by Braaksma [2] for a more general class of integral functions than (.). We present a summary of their results related to the asymptotic expansion of p Ψ q (z) for large z in Section 2. Our purpose here is to consider two of the expansion theorems involving the presence of exponentially small expansions valid in certain sectors of the z-plane. We demonstrate by numerical computation that a more precise understanding of the asymptotic structure of p Ψ q (z) can be achieved by taking into account the Stokes phenomenon. 2. Standard asymptotic theory for z We first state the standard asymptotic expansion of the integral function p Ψ q (z) as z for κ > 0 and finite values of the parameters given in [7] and [2]; see also [, 2.3]. To present this expansion we introduce the exponential expansion E p,q (z) and the algebraic expansion H p,q (z) associated with p Ψ q (z). The exponential expansion E p,q (z) can be obtained from the Ford-Newsom theorem [3, 4]. A simpler derivation of this result in the case p Ψ q (z) based on the Abel-Plana form of the well-known Euler-Maclaurin summation formula is given in [0, pp ]. We have the formal asymptotic sum E p,q (z) := Z ϑ e Z A j Z j, Z = κ(hz) /κ, (2.) j=0 where the coefficients A j are those appearing in the inverse factorial expansion of g(s)/s! given by { M g(s) } Γ(+s) = κ(hκκ ) s A j Γ(κs+ϑ +j) + ρ M (s) Γ(κs+ϑ. (2.2) +M) j=0 Empty sums and products are to be interpreted as zero and unity, respectively.

3 Remarks on the theorems of Wright and Braaksma 3 Here g(s) is defined in (.2) with n replaced by s, M is a positive integer and ρ M (s) = O() for s in arg s < π. The leading coefficient A 0 is specified by p q A 0 = (2π) 2 (p q) κ 2 ϑ β 2 br r. (2.3) α ar 2 r The coefficients A j are independent of s and depend only on the parameters p, q, α r, β r, a r and b r. An algorithm for their evaluation is described in the appendix. The algebraic expansion H p,q (z) follows from the Mellin-Barnes integral representation [, 2.4] pψ q (z) = i Γ(s)g( s)(ze πi ) s ds, arg( z) < π( 2πi 2κ), (2.4) i where the path of integration is indented near s = 0 to separate 2 the poles of Γ(s) from those of g( s) situated at s = (a r +k)/α r, k = 0,,2,... ( r p). (2.5) In general there will be p such sequences of simple poles though, depending on the values of α r and a r, some of these poles could be multiple poles or even ordinary points if any of the Γ(β r s+b r ) are singular there. Displacement of the contour to the right over the poles of g( s) then yields the algebraic expansion of p Ψ q (z) valid in the sector in (2.4). If it is assumed that the parameters are such that the poles in (2.5) are all simple we obtain the algebraic expansion given by H p,q (z), where k=0 H p,q (z) := p m= α m z am/αm S p,q (z;m) (2.6) and S p,q (z;m) denotes the formal asymptotic sum ( ) k ( ) k p +am S p,q (z;m) := Γ Γ(a r α r (k +a m )/α m ) k! q Γ(b, (2.7) r β r (k +a m )/α m ) z k/αm α m with the prime indicating the omission of the term corresponding to r = m in the product. This expression in (2.6) consists of (at most) p expansions each with the leading behaviour z am/αm ( m p). When the parameters α r and a r are such that some of the poles are of higher order, the expansion (2.7) is invalid and the residues must then be evaluated according to the multiplicity of the poles concerned; this will lead to terms involving log z in the algebraic expansion. The three main expansion theorems are as follows. Throughout we let ǫ denote an arbitrarily small positive quantity. Theorem. If 0 < κ < 2, then E p,q (z)+h p,q (ze πi ) in arg z 2 πκ pψ q (z) (2.8) H p,q (ze πi ) in 2 πκ+ǫ arg z π as z. The upper or lower sign in H p,q (ze πi ) is chosen according as arg z > 0 or arg z < 0, respectively. 2 This is always possible when the condition (.3) is satisfied.

4 4 R. B. Paris It is seen that the z-plane is divided into two sectors, with a common vertex at z = 0, by the rays arg z = ± 2 πκ. In the sector arg z < 2 πκ, the asymptotic character of pψ q (z) is exponentially large whereas in the complementary sector 2πκ < arg z π, the dominant expansion of p Ψ q (z) is algebraic in character. On the rays arg z = ± 2πκ the exponential expansion is oscillatory and is of a comparable magnitude to H p,q (ze πi ). Theorem 2. If κ = 2 then pψ q (z) E p,q (z)+e p,q (ze 2πi )+H p,q (ze πi ) (2.9) as z in the sector arg z π. The upper or lower signs are chosen according as arg z > 0 or arg z < 0, respectively. The rays arg z = ± 2 πκ now coincide with the negative real axis. It follows that pψ q (z) is exponentially large in character as z except in the neighbourhood of the negative real axis, where the algebraic expansion becomes asymptotically significant. Theorem 3. When κ > 2 we have 3 pψ q (z) N n= N E p,q (ze 2πin )+H p,q (ze πiz ) (2.0) as z in the sector arg z π. The integer N is chosen such that it is the smallest integer satisfying 2N + > 2κ and the upper or lower is chosen according as arg z > 0 or arg z < 0, respectively. In this case the asymptotic behaviour of p Ψ q (z) is exponentially large for all values of arg z and, consequently, the algebraic expansion may be neglected. The sums E p,q (ze 2πin ) are exponentially large (or oscillatory) as z for values of arg z satisfying arg z+2πn 2 πκṫhe division of the z-plane into regions where pψ q (z) possesses exponentially large or algebraic behaviour for large z is illustrated in Fig.. When 0 < κ < 2, the exponential expansion E p,q (z) is still present in the sectors 2πκ < arg z < min{π,πκ}, where it is subdominant. The rays arg z = ±πκ (0 < κ < ), where E p,q (z) is maximally subdominant with respect to H p,q (ze πi ), are called Stokes lines. 4 As these rays are crossed (in the sense of increasing arg z ) the exponential expansion switches off according to Berry s now familiar error-function smoothing law []; see [8] for details. The rays arg z = ± 2πκ, where E p,q (z) is oscillatory and comparable to H p,q (ze πi ), are called anti-stokes lines. In view of the above interpretation of the Stokes phenomenon a more precise version of Theorem is as follows: Theorem 4. When 0 < κ 2, then E p,q (z)+h p,q (ze πi ) in arg z min{π ǫ,πκ ǫ} pψ q (z) H p,q (ze πi ) in πκ+ǫ arg z π (0 < κ < ) E p,q (z)+e p,q (ze 2πi )+H p,q (ze πi ) in arg z π ( < κ 2) (2.) 3 In [6], the expansion was given in terms of the two dominant expansions only, viz. E p,q(z) and E p,q(ze 2πi ), corresponding to n = 0 and n = ± in (2.0). 4 The positive real axis arg z = 0 is also a Stokes line where the algebraic expansion is maximally subdominant.

5 Remarks on the theorems of Wright and Braaksma 5 Stokes line Θ ΠΚ EXPONENTIALLY SMALL + ALGEBRAIC Θ ΠΚ 2 anti - Stokes line ALGEBRAIC EXPONENTIALLY LARGE + ALGEBRAIC 0 Κ Stokes line Θ ΠΚ EXPONENTIALLY SMALL + ALGEBRAIC anti Stokes line Θ ΠΚ 2 Figure : The exponentially large and algebraic sectors associated with pψ q(z) in the complex z-plane with θ = arg z when 0 < κ <. The Stokes and anti-stokes lines are indicated. as z. The upper or lower signs are chosen according as arg z > 0 or arg z < 0, respectively. We omit the expansion on the Stokes lines arg z = ±πκ; the details in the case p =, q 0 are discussed in [9]. The expansions in (2.a) and (2.8a) were given by Wright [6, 7] in the sector arg z min{π, 3 2πκ ǫ} as he did not take into account the Stokes phenomenon. Since E p,q (z) is exponentially small in 2πκ < arg z π, then in the sense of Poincaré, the expansion E p,q (z) can be neglected and there is no inconsistency between Theorems and 4. Similarly, E p,q (ze 2πi ) is exponentially small compared to E p,q (z) in 0 arg z < π and there is no inconsistency between the expansions in (2.8a) and (2.c) when < κ < 2. However, in the vicinity of arg z = π, these last two expansions are of comparable magnitude and, for real parameters, they combine to generate a real result on this ray. A similar remark applies to E p,q (ze 2πi ) in π < arg z 0. The following theorem was given by Braaksma [2, p. 33]. Theorem 5. If p = 0, so that g(s) has no poles and κ >, then H 0,q (z) 0. When < κ < 2, we have the expansion 0Ψ q (z) E 0,q (z)+e 0,q (ze 2πi ) (2.2) as z in the sector arg z π The upper or lower sign is chosen according as arg z > 0 or arg z < 0, respectively. The dominant expansion 0 Ψ q (z) E p,q (z) holds in the reduced sector arg z π ǫ. It can be seen that (2.2) agrees with (2.c) when H p,q (z) 0. Braaksma gave the result (2.2) valid in a sector straddling the negative real axis given by π δ arg z π + δ, where 0 < δ < 2 π( 2 κ). It is our purpose here to examine Theorems 4 and 5 in more detail by means of a series of examples. We carry out a numerical investigation to show that (2.c) is valid when < κ < 2 and, when 0 < κ <, that the exponential expansion E p,q (z) in Theorem 4 switches off (as arg z increases) across the Stokes lines arg z = ±πκ, where E p,q (z) is

6 6 R. B. Paris maximally subdoiminant with respect to H p,q (ze πi ). Similarly in Theorem 5, we show that when < κ < 2theexpansionsE p,q (ze 2πi )switchoffacrossthestokeslinesarg z = ±π( 2 κ), where they are maximally subdominant with respect to E p,q(z). Thus, although the expansions in (2.a) and (2.2) are valid asymptotic descriptions, more accurate evaluation will result from taking into account the Stokes phenomenon as the above-mentioned rays are crossed. 3. Numerical examples Example 3. Our first example is the Mittag-Leffler function E a,b (z) defined by E a,b (z) = n=0 z n Γ(an+b), where we consider a > 0. This corresponds to a case of Ψ (z) with the parameters κ = a, h = a a, ϑ = b and g(s) = Γ(+s)/Γ(as+b). Then from (2.) (2.3), we have Z = z /a, A 0 = /a with A j = 0 for j. The exponential and algebraic expansions are from (2.), (2.6) and (2.7) given by E, (z) = a z( b)/a exp[z /a ], H, (ze πi ) = k= z k Γ(b ak). Then, from Theorems2, 3 and 4 weobtain the followingasymptotic expansions 5 as z. (i) When 0 < a < a z( b)/a exp[z /a z k ] ( arg z πa ǫ) Γ(b ak) E a,b (z) k= z k (3.) (πa+ǫ arg z π); Γ(b ak) k= (ii) when < a < 2 a z( b)/a exp[z /a z k ] Γ(b ak) k= E a,b (z) a z( b)/a exp[z /a ]+ a (ze 2πi ) ( b)/a exp[(ze 2πi ) /a ] z k Γ(b ak) k= (iii) when a = 2 ( arg z π ǫ) ( arg z π); (3.2) E a,b (z) a z( b)/a exp[z /a ]+ a (ze 2πi ) ( b)/a exp[(ze 2πi ) /a ] 5 When a = we have E,b (z) = z b e z P(b,z), where P(α,z) = γ(α,z)/γ(α) is the normalised incomplete gamma function. It then follows from [5, (8.2.5), (8..2)] that the expansion of E,b (z) is given by (3.a) as z in arg z π.

7 Remarks on the theorems of Wright and Braaksma 7 (iv) when a > 2 z k Γ(b ak) k= ( arg z π); (3.3) E a,b (z) a N n= N (ze 2πin ) ( b)/a exp[z /a e 2πin/a ] k= z k Γ(b ak) ( arg z π), (3.4) where N is the smallest integer 6 satisfying 2N+ > 2a. The upper or lower signs are taken according as arg z > 0 or arg z < 0, respectively. When 0 < a <, it is established in [6] (see also [5]) that the exponential term a exp[z /a ] in (3.a)is multiplied by the approximate factor involving the error function 2 + [ (πa θ) 2 erf a ] z 2 as z in the neighbourhood of the Stokes lines θ = arg z = ±πa, respectively, where it is maximally subdominant. This shows that the above exponential term indeed switches off in the familiar manner [] as one crosses the Stokes lines in the sense of increasing θ and that consequently the expansion in (3.a) is valid in arg z πa ǫ. On the negative real axis we put z = x, with x > 0. From (3.2), we have when < a < 2 E a,b ( x) a (xeπi ) ( b)/a exp[x /a e πi/a ]+ a (xe πi ) ( b)/a exp[(x /a e πi/a ] = F a,b (x) k= as x +, where F a,b (x) = 2 [ a x( b)/a exp x /a cos π ] a k= ( x) k Γ(b ak) ( x) k Γ(b ak), (3.5) [ cos x /a sin π ] a +π( b). (3.6) a The presence of the additional exponential expansion E, (ze 2πi ) in (3.2) is seen to be essential in order to obtain a real result 7 (when b is real) on the negative z-axis. Example 3.2 Our second example is the function F (z) = n=0 Γ( 2 n+a) Γ(n+b) z n n! (κ = 3 2 ), (3.7) 6 The more refined treatment of E a, (z) discussed in [, Section 5..4] has the integer N satisfying N < a < N +. The additional exponential expansions present in (3.4) with this choice of N are, however, 2 exponentially small for arg z π. 7 We remarkthat the result(3.5) can also be deduced by use of the identity E a,b ( x) = 2E 2a,b (x 2 ) E a,b (x) combined with the expansions of E a,b (z) for z +.

8 8 R. B. Paris where a and b are finite parameters, which corresponds to a case of Ψ (z). The exponential expansion is E, (z) = Z ϑ e Z A j Z j, Z = 3 2 (hz)2/3, where, from (2.3), j=0 A 0 = ( 2 3 )ϑ+/2 ( 2 )a /2 and ϑ = a b, h = 2 /2. An algorithm for the computation of the normalised coefficients c j = A j /A 0 is described in the appendix. In our computationswehaveemployed0 j 40; the first ten coefficients c j for F (z) are listed in Table for the particular case a = 4 and b = 3 4. From (2.6), the algebraic expansion is H, (ze πi ) = 2 k=0 ( ) k Γ(2k +2a) k!γ(b 2a 2k) (ze πi ) 2k 2a. j c j j c j Table : The normalised coefficients c j for j 0 (with c 0 = ) for the sum (3.7) when a = 4 and b = 3 4. It is clearly sufficient for real parameters to consider values of z satisfying 0 arg z π and this we do throughout this section. From Theorem 4, we obtain F (z) = E, (z)+e, (ze 2πi )+H, (ze πi ) as z in 0 arg z π, from which we see that F (z) is exponentially large in the sector arg z < 3π/4. We have computed F (z) for a value of z and varying θ = arg z in the range 0.7π θ π. In Table 2 we show the absolute values of R (z) F (z) E opt, (z) Hopt, (ze πi ) compared with E, (ze 2πi ) (which was computed for 0 j 5), where the superscript opt denotes that both the asymptotic sums E, (z) and H, (ze πi ) are truncated at their respective optimal truncation points. The results clearly confirm that (i) the exponential expansion E, (z) is present in the algebraicsector 3 4π < arg z π and (ii) the subdominant expansion E, (ze 2πi ) is present in (at least) the sector 0.7π θ π. It was not possible to penetrate very far into the exponentially large sector arg z < 3 4π, since the error in

9 Remarks on the theorems of Wright and Braaksma 9 θ/π R (z) E, (ze 2πi ) Table 2: Values of the absolute error R (z) in the computation of F (z) using an optimal truncation of both E, (z) and H, (ze πi ) compared with E, (ze 2πi ) as a function of θ for z = 00e iθ, a = 4 and b = 3 4. the computation of E, (z) even at optimal truncation swamps the algebraic and subdominant exponential expansions. Such a computation would require a hyperasymptotic evaluation of the dominant expansion on the lines of that described for the generalised Bessel function 0 Ψ (z) in Wong and Zhao [4]. Example 3.3 Consider the function F 2 (z) = n=0 Γ( 2 3 n+a)zn Γ( 3 n+b)n! (κ = 2 3 ). According to Theorem 4, the expansion of F 2 (z) for large z is F 2 (z) E, (z)+h, (ze πi ) (0 arg z 2 3 π ǫ). The algebraic expansion is, from (2.6), given by H, (ze πi ) = 3 ( ) k Γ( 3 2 k a) 2 k!γ(b 2 a ) 3(k+a)/2 2 k)(ze πi k=0 and the exponential expansion E, (z) is obtained from (2.) with the parameters ϑ = a b, h = ( 2 3 3( )2 3 ) 3 and A 0 = κ 2 ϑ ( 2 3 2( )a 3 ) 2 b. The coefficients A j are obtained as indicated in Example 3.2. The function F 2 (z) is exponentially large in the sector arg z < 3π, whereas in the sector 3 π < arg z π the algebraic expansion H,(ze πi ) is dominant. The expansion E, (z) is maximally subdominant with respect to H, (ze πi ) on the ray arg z = πκ = 2 3 π. Consequently, as arg z increases, the exponential expansion E, (z) should switch off across the Stokes line arg z = 2 3 π, to leave the algebraic expansion H,(ze πi ) in the sector 2 3 π < arg z π. To demonstrate this, we define the Stokes multiplier S(θ) by F 2 (z) = H opt (ze πi )+A 0 Z ϑ e Z S(θ). In Table 3 we show the absolute values of R 2 (z) := F 2 (z) H opt, (ze πi ) and of the leading term of E, (z) as a function of θ = arg z. We also show the values 8 of Re(S(θ)) in the 8 The Stokes multiplier S(θ) has a small imaginary part that we do not show.

10 0 R. B. Paris neighbourhood of the Stokes line arg z = 2 3 π for the case z = 0eiθ and a = 3, b = 4. It is seen that the Stokes multiplier has the value when θ = 2π (before the transition commences) and 0 when θ = 3 4π (after the transition is almost completed). θ/π R 2 (z) A 0 Z ϑ e Z Re(S(θ)) Table 3: Values of the absolute error in R 2 (z) F 2 (z) H opt, (ze πi ) in the computation of F 2 (z) using an optimal truncation of the algebraic expansion compared with the leading term of E, (z) as a function of θ for z = 0e iθ, a = 3 and b =. The final column shows the real part of the computed Stokes multiplier 4 S(θ) for transition across the ray arg z = 2 3 π. Example 3.4 Our final example is the function of the type 0 Ψ 2 (z) given by F 3 (z) = n=0 z n n!γ(cn+a)γ(cn+b) (κ = +2c), (3.8) where 0 < c 2. Since p = 0, the algebraic expansion H 0,2(z) 0. From Theorem 5 we obtain the asymptotic expansion F 3 (z) E 0,2 (z)+e 0,2 (ze 2πi ) ( arg z π), where the associated parameters are ϑ = a b, h = c 2c and A 0 = cϑ κ ϑ /2 2π. The function F 3 (z) is exponentially large in the sector arg z < 2π( + 2c). The other expansion E 0,2 (ze 2πi ) is subdominant in the upper half-plane but combines with E 0,2 (z) on the negative real axis to produce (for real a and b) a real expansion. Since the exponential factors associated with E 0,2 (z) and E 0,2 (ze 2πi ) are exp[ Z e iθ/κ ] and exp[ Z e i(θ 2π)/κ ], where θ = arg z and we recall that Z is defined in (2.), the greatest difference between these factors occurs when ( ) ( ) θ θ 2π sin = sin ; κ κ

11 Remarks on the theorems of Wright and Braaksma that is, when θ = 2π(2 κ). Consequently, as arg z increases in the upper half-plane, we expect that the expansion E 0,2 (ze 2πi ) should switch on across the Stokes line arg z = 2 π(2 κ); similarconsiderationsapply to E 0,2(ze 2πi ) and the Stokes line arg z = 2 π(2 κ) in the lower half-plane. To demonstrate the correctness of this claim, we choose c = 0 (so that κ = 6 5 ) and a = 4, b = 3 4. The function F 3(z) is therefore exponentially large in the sector arg z < 3 5 π and the Stokes line in the upper half-plane is arg z = 2 5π. We have chosen a b to have a half-integer value for a very specific reason. The more detailed treatment in [8] shows that there is a third (subdominant) exponential series present in the expansion of F 3 (z) given by 2cosπ(a b)x ϑ e X A j ( X) j, X = κ(hze πi ) /κ. j=0 Our present choice of a and b therefore eliminates this third expansion and enables us to deal with a case comprising only two exponential expansions. In Table 4, we show for z = 20 and varying θ = arg z the values of F 3 (z) E opt 0,2 (z) and E 0,2 (ze 2πi ) together with the real part of the Stokes multiplier S(θ) defined by F 3 (z) = E opt 0,2 (z)+a 0(Ze 2πi/κ ) ϑ exp[ze 2πi/κ ]S(θ). The results clearlydemonstratethe switching-onofthe subdominant expansione 0,2 (ze 2πi ) across the Stokes line arg z = 2 5π as arg z increases in the upper half-plane. θ/π R 3 (z) E 0,2 (ze 2πi ) Re(S(θ)) Table 4: Values of the absolute error in R 3 (z) F 3 (z) E opt 0,2 (z) in the computation of F 3(z) using an optimal truncation of E 0,2 (z) compared with E 0,2 (ze 2πi ) as a function of θ for z = 20e iθ, a = 4 and b = 3. The final column shows the real part of the computed Stokes multiplier S(θ) for transition across 4 the ray arg z = 2 5 π. Appendix: An algorithm for the computation of the coefficients c j = A j /A 0 We describe an algorithm for the computation of the normalised coefficients c j = A j /A 0 appearing in the exponential expansion E p,q (z) in (2.). Methods of computing these coefficients by recursion in the case α r = β r = have been given by Riney [2] and Wright [8];

12 2 R. B. Paris see [, Section 2.2.2] for details. Here we describe an algebraic method for arbitrary α r > 0 and β r > 0. The inverse factorial expansion (2.2) can be re-written as g(s)γ(κs+ϑ ) Γ(+s) { M = κa 0 (hκ κ ) s j=0 c j (κs+ϑ ) j + O() (κs+ϑ ) M } (A.) for s uniformly in arg s π ǫ, where g(s) is defined in (.2) with n replaced by s. Introduction of the scaled gamma function Γ (z) = Γ(z)(2π) 2e z z 2 z leads to the representation Γ(αs+a) = (2π) 2 e αs (αs) αs+a 2 e(αs;a)γ (αs+a), where ( e(αs;a) := e a + a ) αs+a 2 ( = exp [(αs+a αs 2 )log + a ) ] a. αs as where Then, after some routine algebra we find that the left-hand side of (A.) can be written Υ p,q (s) := g(s)γ(κs+ϑ ) Γ(+s) p Γ (α r s+a r ) q Γ (β r s+b r ) = κa 0 (hκ κ ) s R p,q (s)υ p,q (s), (A.2) Γ (κs+ϑ p ) Γ (+s), R p,q(s) := e(α rs;a r ) e(κs;ϑ ) q e(β rs;b r ) e(s;). Substitution of (A.2) in (A.) then yields the inverse factorial expansion in the form R p,q (s)υ p,q (s) = M j=0 c j (κs+ϑ ) j + O() (κs+ϑ ) M (A.3) as s in arg s π ǫ. We nowexpand R p,q (s) andυ p,q (s) fors + makinguse ofthe well-knownexpansion (see, for example, [, p. 7]) Γ (z) ( ) k γ k z k k=0 where γ k are the Stirling coefficients, with ( z ; arg z π ǫ), Then we find whence γ 0 =, γ = 2, γ 2 = 288, γ 3 = , γ 4 = ,.... Γ (αs+a) = γ αs +O(s 2 ), R p,q (s) = + A 2s +O(s 2 ), e(αs;a) = + a(a ) 2αs +O(s 2 ), Υ p,q (s) = + B 2s +O(s 2 ),

13 Remarks on the theorems of Wright and Braaksma 3 where we have defined the quantities A and B by A = p a r (a r ) α r q b r (b r ) ϑ p β r κ ( ϑ), B = α r Upon equating coefficients of s in (A.3) we then obtain q β r + κ. c = 2 κ(a+ 6 B). (A.4) The higher coefficients are obtained by continuation of this expansion process in inverse powers of s. We write the product on the left-hand side of (A.3) as an expansion in inverse powers of κs in the form M C j R p,q (s)υ p,q (s) = + (κs) j +O(s M ) j= (A.5) as s +, where the coefficients C j are determined with the aid of Mathematica. From the expansion of the ratio of two gamma functions in [5, (5..3)] we obtain (κs+ϑ = { M ( ) k } (j) k ) j (κs) j (κs) k B ( j) k (ϑ )+O(s M ), k! j=0 where B (σ) k (x) are the generalised Bernoulli polynomials defined by ( ) σ t e t e xt = k=0 B (σ) k (x) t k k! ( t < 2π). Here we have σ = j 0 and B (σ) 0 (x) =. Then the right-hand side of (A.3) as s + becomes with M + j= M c j (κs+ϑ +O(s M ) = + ) j j= M c j ( ) k (j) k (κs) j (κs) k B ( j) k (ϑ )+O(s M ) k! k=0 M D j = + (κs) j +O(s M ) j= j ( ) j D j = ( ) k c k j k B (k j+) k (ϑ ), k=0 (A.6) where we have made the change in index j + k j and used triangular summation (see [3, p. 58]). Substituting (A.5) and (A.6) into (A.3) and equating the coefficients of like powers of κs, we then find C j = D j for j M, whence j ( c j = C j ( ) k j k k= ) c j k B (k j+) k (ϑ ).

14 4 R. B. Paris Thus we find c = C, c 2 = C 2 c B (0) (ϑ ), c 3 = C 3 2c 2 B ( ) (ϑ )+c B (0) 2 (ϑ ),... and so on, from which the coefficients c j can be obtained recursively. With the aid of Mathematica this procedure is found to work well in specific cases when the various parameters have numerical values, where up to a maximum of 00 coefficients have been so calculated. References [] M.V. Berry, Uniform asymptotic smoothing of Stokes s discontinuities, Proc. Roy. Soc. London A422 (989) 7 2. [2] B.L.J. Braaksma, Asymptotic expansions and analytic continuations for a class of Barnes-integrals, Compos. Math. 5 (963) [3] W.B. Ford, The Asymptotic Developments of Functions Defined by Maclaurin Series, University of Michigan, Science Series, Vol. II, 936. [4] C.V. Newsom, On the character of certain entire functions in distant portions of the plane, Amer. J. Math. 60 (938) [5] F.W.J. Olver, D.W. Lozier, R.F. Boisvert and C.W. Clark (eds.), NIST Handbook of Mathematical Functions, Cambridge University Press, Cambridge, 200. [6] R.B. Paris, Exponential asymptotics of the Mittag-Leffler function, Proc. Roy. Soc. London 458A (2002) [7] R.B. Paris, Some remarks on the theorems of Wright and Braaksma concerning the asymptotics of the generalised hypergeometric functions, Technical Report MS 09:0, Abertay University, [8] R.B. Paris, Exponentially small expansions in the asymptotics of the Wright function, J. Comput. Appl. Math. 234 (200) [9] R.B. Paris, Exponentially small expansions of the Wright function on the Stokes lines, Lithuanian Math. J. 54 (204) [0] R.B. Paris and A.D. Wood, Asymptotics of High Order Differential Equations, Pitman Research Notes in Mathematics, 29, Longman Scientific and Technical, Harlow, 986. [] R.B. Paris and D. Kaminski, Asymptotics and Mellin-Barnes Integrals, Cambridge University Press, Cambridge, 200. [2] T.D. Riney, On the coefficients in asymptotic factorial expansions, Proc. Amer. Math. Soc. 7 (956) [3] L.J. Slater, Generalized Hypergeometric Functions, Cambridge University Press, Cambridge, 966. [4] R. Wong and Y.-Q. Zhao, Smoothing of Stokes s discontinuity for the generalized Bessel Proc. Roy. Soc. London A455 (999) [5] R. Wong and Y.-Q. Zhao, Exponential asymptotics of the Mittag-Leffler function, Constr. Approx. 8 (2002)

15 Remarks on the theorems of Wright and Braaksma 5 [6] E.M. Wright, The asymptotic expansion of the generalized hypergeometric function, J. Lond. Math. Soc. (Ser. 2) 0 (935) [7] E.M. Wright, The asymptotic expansion of the generalized hypergeometric function, Proc. Lond. Math. Soc. (Ser. 2) 46 (940) [8] E.M. Wright, A recursion formula for the coefficients in an asymptotic expansion, Proc. Glasgow Math. Assoc. 4 (958) 38 4.