![Arxiv:2012.07068V3 [Math.CA] 22 Feb 2021 Lxne Apelblat Alexander 40 Ersea 40,Ire.Eal Apelblat@Bgu.Ac.I Email: Israel](http://data.docslib.org/img/4c4e1ed234235722b1502996726bdb2e-1.webp)
arXiv:2012.07068v3 [math.CA] 22 Feb 2021 lxne Apelblat Alexander 40 erSea 40,Ire.Eal [email protected] Email: Israel. b 84105, Sheva, Beer 84105 a Le riso h nes alc rnfr of transform Laplace inverse the of erties Mittag-Le the fun functions, hyperbolic and fun trigonometric elementary logarithmic, include ponential, integrals determined Mittag-Le of the a integrands identities with integral integrals The convolution the of transform. Laplace for inverse der identities the integral to by and possible integrals was finite it integrals, calculus, operational and Wlodarski, Abstract: M ujc lsicto : Classification Subject AMS 10.3390/sym13020354 DOI: pages. 15 354, art Paolo (2021), Editor 13 Vol Symmetry guest of with Issue Polynomials”, Special and the Functions ”Special in published been has paper This integrals. convolution and r presented. are Keywords: ∗ 0000-0003-4858-7309 ORCID: [email protected] Email: Italy. Bologna, 40126 eateto hmclEgneig e uinUniversity Gurion Ben Engineering, Chemical of Department eateto hsc n srnm,Uiest fBologna of University Astronomy, and Physics of Department ucinrpeetdb h nes Laplace inverse the by represented function orsodn author. Corresponding ffl rfntos otrafntos oie eslfunction Bessel modified functions; Volterra functions; er plcto fteErstermt the to theorem Efros the of Application fo hoe;ivreLpaetasom;Wih functio Wright transforms; Laplace inverse theorem; Efros sn pca aeo h fo hoe hc a eie by derived was which theorem Efros the of case special a Using rnfr of transform ffl rfntosadteVler ucin.Sm prop- Some functions. Volterra the and functions er 2Fbur 2021 February 22 a 63,3C0 31,3A5 44A20 34A25, 33E12, 33C10, 26A33, n rnec Mainardi Francesco and 1 s s − − µ µ ffl exp( exp( radVler ucin.The functions. Volterra and er − s ν with ) − s ucinrepresented function .it ν tos n h error the and ctions) l ) µ emil nterms in mainly re tos(oe,ex- (power, ctions v ayinfinite many ive ≥ ;fiie infinite finite, s; i rei 46, Irnerio Via , n 0 and 0 mloRicci, Emilio fteNegev, the of (MDPI) b ∗ s Mittag- ns; < ν < 1 1 Introduction
Inversions of the Laplace transforms of exponential functions
c+i 1 ∞ f (t) = 1 F(s) = est F(s)ds ; c> 0, ν,µ L− { } 2πi ν Zc i (1) e s − ∞ F(s) = − ; 0 <ν< 1 ; µ 0, sµ ≥ were during the 1945 - 1970 period in a focus of attention of a number of well- known mathematicians like Humbert [1], Pollard [2], Wlodarski [3], Mikusin- ski [4-7], Wintner [8], Ragab [9] and Stankovicˇ [10]. In the case µ = 0, Mikusin- ski was able to obtain the inverse Laplace transform in terms of integral repre- sentations [7]
∞ st sν fν,0(t) = e− fν,0(t)dt = e− ; 0 <ν< 1 ; t> 0 L Z0 1 ∞ ut uν cos(πν) ν fν,0(t) = e− e− sin[u sin(πν)]du (2) π Z0 2 ∞ ν πν u cos( 2 ) ν πν fν,0(t) = e− cos[u cos( 2 )] cos(ut)du π Z0 and also as the finite trigonometric integral
π ν ξ fν,0(t) = ξe− du ; 0 <ν< 1 ; t> 0 π(1 ν)t Z0 − ν/(1 ν) (3) 1 sin(νu) − sin[(1 ν)u] ξ = − ν/(1 ν) t − sinu ! sinu
It was established that the functions fν,0(t) can be expressed in terms of expo- nential and parabolic cylinder functions when ν = 1/2 [9,11] and by help of the Airy functions and their first derivatives for ν = 1/3 and ν = 2/3 [9,12]. For ν = 1/4, the solution was deduced by Barkai [13] in 2001 using Mathematica as a sum of three generalized hypergeometric functions, but the numerical result was uncertain, presumably for a bag in the computing program . In 2010-2012 Gorska and Penson [14]-[15] were able to represent fν,0(t) in terms of Mejer G functions. Earlier, in 1958 Ragab [9] expressed fν,0(t) in terms of MacRobert E functions for ν = 1/n. In 1952 Wlodarski [3] showed that when the generalized product Efros the- orem [16] is applied to the Laplace transform of fν,µ(t) which is given in (1), then it is possible to derive the following formula
g(t) = G(s) ; 0 <ν< 1 ; µ 0 L{ } ≥ 1 G(sν ) ∞ 1 µ = g(u) F(s) du = L− s L− { } (4) n o Z0 usν ∞ 1 e− ∞ g(u) − sµ du = g(u)fν,µ(u)du Z0 L Z0
2 It was established in our recent paper [17] that for µ =0and µ = 1 ν this func- − tional expression can be written in terms of specific Wright functions Wν,µ( t), sometimes referred to as the Mainardi functions −
Fν(t) = W ν,0( t) ; 0 <ν< 1 − − (5) Mν(t) = W ν,1 ν( t) − − − in the following way
1 G(sν ) ∞ u du ν − 1 ν = g(u)Fν tν u ; 0 <ν< 1 L s − Z0 (6) n νo ν 1 G(s ) ∞ u t − 1 ν = g(u)Mν tν du L s − n o Z0
This follows from the fact that the functions Fν(t) and Mν(t) satisfy
ν 1 λ 1 λ e λs = = − L Fν ν L ν Mν ν 1 ν λν t t t s − (7) 0 <ν< 1 ; λ> 0
It was also illustrated by us in [17], that in many cases, by using standard tables of the Laplace transforms [18-21], the left hand side Laplace transforms in (6) can be inverted and numerous infinite integrals, finite integrals and integral identities for the functions Fν(t) and Mν(t) can be derived. Let us recall that the Wright functions [22-23], considered initially as a some kind generalization of the Bessel functions, are defined as an entire functions of the argument z C and parameters λ> 1 and µ C by ∈ − ∈ ∞ zk W (z) = ; λ> 1, µ C. (8) λ,µ k! Γ (λk + µ) − ∈ Xk=0 It is usual to distinguish them in two kinds, the first kind with λ 0 and the second kind with λ = ν and ν (0,1), see e.g [11]. ≥ Restricting our attention− to positive∈ argument t> 0, the Mainardi functions turn out to be particular Wright functions of the second kind expressed by the following series
∞ ( t)k 1 ∞ ( 1)k +1tk F (t) = − = − Γ (νk + 1) sin(πνk), ν k!Γ ( νk) π k! kX=1 − kX=1 k k 1 ∞ ( t) 1 ∞ ( t) − (9) M (t) = − = − Γ (νk) sin(πνk), ν k!Γ ( ν(k + 1) + 1) π (k 1)! kX=0 − kX=1 −
Fν(t) = νtMν(t).
3 The interest in the Wright functions of the second kind comes from the fact that they play an important role in solution of the linear partial differential equations of fractional order which describe a wide spectrum of phenomena including probability distributions, anomalous diffusion and diffusive waves [24-33]. In the section 2 of this paper, by using the complex inversion formula in (1) and the Bromwich contour, the integral representation of fν,µ(t) with 0 <ν< 1 and 0 µ< 1 is derived and some basic properties of this inverse transform are established.≤ The cases µ 1 can be dealt as well, see Appendix for details. The next two sections≥ 3 and 4 are devoted to evaluation of infinite, fi- nite and convolution integrals by using the Efros theorem in the Wlodarski form. The integrands of these integrals or integral identities include the ele- mentary functions (power, exponential, logarithmic, trigonometric and hyper- bolic functions) and the special functions (the error function, Mittag-Leffler functions and the Volterra functions). The last section provides concluding remarks. In derivations of these integrals direct and inverse Laplace transforms which are taken from tables of transforms [18-21] are always presented in mathe- matical expressions. All mathematical operations and manipulations with el- ementary and special functions, integrals and transforms are formal and their validity is assured by considering the restrictions usually imposed in the oper- ational calculus.
2 Integral representations of the inverse Laplace µ ν transform of s− exp( s ) − Before infinite integrals in (6) will be evaluated, it is of interest to derive the function fν,µ(t) by performing the complex integration from (1). In investi- gated case, the branch point of the integrand exists and is located at the origin s = 0 and therefore the equivalent Bromwich contour is plotted in Figure 1. The closed contour of integration ABCDEFA includes the line AB, the arcs BC and FA of a circle of radius R with center at origin, the arc DE of a circle of radius r 0 with center at→ origin, ∞ and two parallel lines CD and EF. The the cut along→ the negative axis ensures that F(s) is a single-valued function. However, according to the Cauchy lemma the integrals along the arcs BC and FA vanish as R . As a consequence→ ∞ there are only three contributions coming from integrals on the CD and EF lines and from the small circle round origin,D E, so we have
sν e− F(s) = µ ; 0 <ν< 1 ; µ 0, s ≥ π 1 ∞ ut πi πi 1 treiθ iθ iθ fν,µ(t) = e− F(ue− ) F(ue ) du + lim e F(re )re dr. 2πi Z0 − 2πi r 0Z π h i → − (10)
4 Figure 1: The equivalent Bromwich contour.
However, the last trigonometric integral vanishes for µ< 1 and therefore the final result of the complex integration from (10) is
ut uν cos(πν) 1 ∞ e− − ν fν,µ(t) = µ sin[u sin(πν)+ πµ]du π Z0 u (11) t> 0 ; 0 <ν< 1 ; µ< 1
The same integral representation has been derived in 1970 by Stankovic˙ [10] for the Wright function, but on the negative values of the argument t and which is presented here in our notation
ut uν cos(πν) µ 1 1 1 ∞ e− − ν t − Wν,µ ν = µ sin[u sin(πν)+ πµ]du −t π Z0 u (12) t> 0 ; 0 <ν< 1 ; µ< 1 and comparing (11) with (12) we have that the inverse exponential functions can be expressed in terms of the Wright functions, in agreement with the sur- vey analysis by Mainardi and Consiglio [33] where also plots are presented,
µ 1 1 fν,µ(t) = t − Wν,µ ν ; 0 <ν< 1 ; µ 0,. (13) −t ≥ In particular, for µ = 1 ν the Wright functions are reduced to the Mainardi − functions Fν(t) and Mν(t) [11] 1 1 1 1 fν,1 ν(t) = Fν ν = ν Mν ν ; 0 <ν< 1. (14) − ν t t t
5 For ν = 1/2, µ =0 and ν = µ = 1/2, the integrals in (11) become the Laplace transforms of trigonometric functions [18-21]
1/4t 1 ∞ e f (t) = e ut sin(u1/2)du = − 1/2,0 − 3/2 π Z0 2 √πt ut ut 1/4t 1 ∞ e π 1 ∞ e e f (t) = − sin(u1/2 + )du = − cos(u1/2)du = − 1/2,1/2 1/2 1/2 π Z0 u 2 π Z0 u √π (15) Since [11,12] 1 1 2 F = K (16) 1/3 1/3 1/2 1/3 t 3πt √27t ! it follows from (11) and (14) that
ut u1/3/2 1 ∞ e √3 √3 f (t) = − − √3cos u1/3 sin u1/3 du = 1/3,2/3 2π 2/3 2 − 2 Z0 u " ! !# (17) 1 λ 2 K1/3 π r t √27t ! Differentiation of the integral (11) with respect to the argument t gives
dnf (t) n ut uν cos(πν) ν,µ ( 1) ∞ e− − ν n = fν,µ n(t) = − µ n sin[u sin(πν)+ πµ]du dt − π Z0 u − (18) t> 0 ; 0 <ν< 1 ; µ< 1 ; n = 0,1,2,3,... and using rules of the operational calculus we have the initial and final values of the function from [33]
1 µ sν fν,µ(t +0) = lim [sF(s)] = s − e− =0 ; 0 <ν< 1 ; µ< 1 → s →∞ (19) 1 µ sν fν,µ(t ) = lim[sF(s)] = s − e− = 0 → ∞ s 0 → This is in an agreement with the finding of Pollard [2] that for µ = 0, the inverse Laplace transform fν,0 is positive almost everywhere, but for µ< 0 Stankovic˙ [10] postulated that in a some interval this function is negative and has at least one zero. Expanding the exponent in F(s) into series it is possible to obtain the behaviour of the function for large values of the argument t
ν e s 1 sν s2ν 1 F(s) =Ł 1 − = 1 1 + ... L− − sµ L− sµ − 1! 2! − { } ( ) ( " #) (20) 1 1 1 1 fν,µ(t ) + ... → ∞ ∼ t1 µ (Γ (µ) − Γ (µ ν)tν 2Γ (µ 2ν)t2ν − ) − − − The integral of the fν,µ can be derived from directly from
t sν 1 e− fν,µ(u)du = − 1+µ = fν,µ+1(t) (21) Z0 L (s )
6 In order to obtain the recurrence relations the following operational rule can be applied
sν sν sν d d e− e− e− L tfν,µ 1(t) = [F(s)] = = (µ 1) µ + ν µ ν (22) − −ds −ds (sµ 1 ) − s s n o − − and the inverse transforms in (22) are
tfν,µ 1(t) = (µ 1)fν,µ(t)+ νfν,µ ν(t). (23) − − − However, from (18) with n = 1, we have
dfν,µ(t) = fν,µ 1(t) (24) dt − and therefore dfν,µ(t) t = (µ 1)fν,µ(t)+ νfν,µ ν(t). (25) dt − − 3 Integrals of the inverse Laplace transform of µ ν s− exp(s ) with elementary functions
In the first example the Wlodarski integral formula (4) is applied to the power function g(t) = tλ,
g(t) = tλ ; λ> 0 Γ (λ + 1) Γ (λ + 1) G(s) = ; G(sν) = sλ+1 s(λ+1)ν (26) Γ Γ Γ (λ+1)ν +µ 1 1 1 (λ + 1) 1 (λ + 1) (λ + 1)t − − = − = L (sµ • s(λ+1)ν ) L (s(λ+1)ν +µ ) Γ [(λ + 1)ν + µ)] and therefore we have Γ (λ+1)ν +µ 1 ∞ λ (λ + 1)t − u fν,µ(u)du = Γ ; 0 <ν< 1 ; µ< 1 ; λ> 0. (27) Z0 [(λ + 1)ν + µ)] The above results can be extended to functions g(t) which are defined at finite intervals and to step functions jumping at integral values of variable t. Let start with 1 ; 0
7 which leads to the finite integral
λ tν +µ 1 f (u)du = − λ[(µ 1)/ν +1]f ( t ) ν,µ Γ − ν,ν +µ λ1/ν Z0 (ν + µ) − (29) 0 <ν< 1 ; µ< 1 ; λ> 0
In the next example, from
1 t ; 0
1 2ν +µ 1 ν +µ 1 t − t − (1 u)fν,µ(u)du = Γ Γ + fν,2ν +µ(t) (31) Z0 − (2ν + µ) − (ν + µ) Using 0 ; 0
t ∞( )µ ( ) = [(1/ν +1)µ 1/ν +1] (33) t λ fν,µ u du λ − fν,ν +νµ+µ 1/ν . Zλ − λ In the next two examples of the case of exponential functions are considered
αt g(t) = e− ; α> 0 1 1 G(s) = ; G(sν) = s + α sν + α (34) ν 1 ν +µ 2 1 1 1 1 1 s − t − ν − = − = ⋆ Eν ( αt ) L sµ • sν + α L (sν +µ 1 • sν + α ) Γ (ν + µ 1) − − −
8 From (4) and (34) the convolution integral with the Mittag-Leffler function
∞ tk E (t) = , ν Γ (kν + 1) kX=0 (35) sν 1 E ( αtν) = − , L{ ν − } sν + α is derived ν +µ 2 ∞ αu t − ν e− fν,µ(u)du = Γ ⋆ Eν( αt ) Z0 (ν + µ 1) − (36) − 0 <ν< 1 ; µ< 1 ; α> 0. By changing variable of integration x = t (cosθ)2, all convolution integrals can be expressed as in terms of finite trigonometric integrals. The shifted increas- ing and decreasing exponential functions are considered in the following two examples. From
0 ; 0
∞ (u λ) (ν +µ 1)/ν ν t 1 e− − fν,µ(u)du = λ − Eν( t ) ⋆fν,2ν +µ 1( ). (38) − − − λ1/ν Zλ h i Similarly from
0 ; 0
∞ t e (u λ) f (u)du = λ(ν +µ 2)/ν E ( tν) ⋆f ( ). (40) − − ν,µ − ν ν,ν +µ 1 1/ν Zλ − − λ
9 The logarithmic functions are the next group of elementary functions to be considered. In the simplest case from
g(t)=lnt ; C = eγ ln(Cs) ln(Csν) G(s) = ; G(sν) = − s − sν 1 ln(Csν) 1 ν ln(Cs) + (ν 1)γ 1 = 1 = (41) − µ − ν − ν +µ 1 − − L (s • s ) L (s − • s ) (ν 1)γtν +µ 1 νtν +µ 2 − − + − lnt Γ (ν + µ) Γ (ν + µ 1) ∗ − where γ is the Euler constant, it follows that
ν +µ 1 ν +µ 2 ∞ (ν 1)γt − νt − lnufν,µ(u)du = − + ⋆ lnt . (42) Z Γ (ν + µ) Γ (ν + µ 1) 0 − In the more general case
λ 1 g(t) = t − lnt ; λ> 0 ; 0 <ν< 1 Γ (λ) Γ (λ)[ψ(λ) ν lns] G(s) = [ψ(λ) lns ]; G(sν) = − sλ − sλν 1 [ψ(λ) ν lns] ψ(λ)+ γν ν ln(Cs) Γ (λ) 1 − = Γ (λ) 1 = − µ λν − (λν +µ) (λν +µ 1) L (s • s ) L ( s − s − • s ) Γ (λ)[ψ(λ)+ γν]t(λν +µ 1) Γ (λ)νt(λν +µ 2) − + − ⋆ lntν Γ (λν + µ) Γ (λν + µ 1) − (43) we have Γ (λν +µ 1) Γ (λν +µ 2) ∞ λ 1 (λ)[ψ(λ)+ γν] t − (λ)νt − ν u − lnufν,µ(u)du = Γ + Γ ⋆ lnt Z0 (λν + µ) (λν + µ 1) − (44) Trigonometric and hyperbolic functions is the last groups of elementary functions to be considered. From g(t) = sin(λt) ; 0 <ν< 1 λ λ G(s) = ; G(sν) = s2 + λ2 s2ν + λ2 1 1 1 s2ν 1 λt2ν +µ 2 λ 1 = λ 1 − = − ⋆ E ( λ2 t2ν), − µ 2ν 2 − 2ν +µ 1 2ν 2 Γ 2ν L s • s + λ L (s − • s + λ ) (2ν + µ 1) − − (45) the following integral identity is derived
2ν +µ 2 ∞ λt − 2 2ν sin(λu)fν,µ(u)du = ⋆ E2ν ( λ t ). (46) Z Γ (2ν + µ 1) − 0 −
10 Similarly as in (45), for the hyperbolic sine function, the change is only in the sign
g(t) = sinh(λt) ; 0 <ν< 1 λ λ G(s) = ; G(sν) = s2 λ2 s2ν λ2 1 − 1 − 1 s2ν 1 λt2ν +µ 2 λ 1 = λ 1 − = − ⋆ E (λ2t2ν ) − µ 2ν 2 − 2ν +µ 1 2ν 2 Γ 2ν L s • s λ L (s − • s λ ) (2ν + µ 1) − − − (47) and therefore we have 2ν +µ 2 ∞ λt − 2 2ν sinh(λu)fν,µ(u)du = E2ν (λ t ). (48) Z Γ (2ν + µ 1) ∗ 0 − In the case of cosine function, from
g(t) = cos(λt) ; 0 <ν< 1 s sν G(s) = ; G(sν) = s2 + λ2 s2ν + λ2 1 sν 1 s2ν 1 tν +µ 2 1 = 1 − = − ⋆ E ( λ2t2ν ), − µ 2ν 2 − ν +µ 1 2ν 2 Γ 2ν L (s • s + λ ) L (s − • s + λ ) (ν + µ 1) − − (49) we have ν +µ 2 ∞ t − 2 2ν cos(λu)fν,µ(u)du = ⋆ E2ν( λ t ). (50) Z Γ (ν + µ 1) − 0 − Similarly as in (47) and (48)
g(t) = cosh(λt) ; 0 <ν< 1 s sν G(s) = ; G(sν) = s2 λ2 s2ν λ2 1 − sν 1− s2ν 1 tν +µ 2 1 = 1 − = − ( 2 2ν) − µ 2ν 2 − ν +µ 1 2ν 2 Γ ⋆ E2ν λ t , L (s • s λ ) L (s − • s λ ) (ν + µ 1) − − − (51) for the hyperbolic cosine function we have
ν +µ 2 ∞ t − 2 2ν cosh(λu)fν,µ(u)du = ⋆ E2ν (λ t ). (52) Z Γ (ν + µ 1) 0 − In the case of the product of trigonometric and hyperbolic sine functions the direct and inverse Laplace transforms are
g(t) = sin(λt)sinh(λt) ; 0 <ν< 1 2λ2s 2λ2sν G(s) = ; G(sν) = s4 + 4λ4 s4ν + 4λ4 1 sν 1 s4ν 1 (53) 2 2 1 = 2 2 1 − = λ − µ 4ν 4 λ − 3ν +µ 1 4ν 4 L (s • s + 4λ ) L (s − • (s + 4λ )) 2λ2 (3ν +µ 2) 4 4ν Γ (3ν +µ 1) [t − ⋆ E4ν( 4λ t )] − −
11 and therefore 2 ∞ 2λ (3ν +µ 2) 4 4ν sin(λu)sinh(λu)fν,µ(u)du = Γ [t − ⋆ E4ν( 4λ t )]. Z0 (3ν + µ 1) − − (54) For the product of trigonometric and hyperbolic sine functions we have
g(t) = cos(λt)cosh(λt) ; 0 <ν< 1 s3 s3ν G(s) = ; G(sν) = s4 + 4λ4 s4ν + 4λ4 1 s3ν 1 s4ν 1 tν +µ 2 1 = 1 − = − ⋆ E ( 4λ4t4ν), − µ 4 4 − ν +µ 1 4ν 4 Γ 4ν L (s • s + 4λ ) L (s − • s + 4λ ) (ν + µ 1) − − (55) and therefore ν +µ 2 ∞ t − 4 4ν cos(λu)(cosh(λu)fν,µ(u)du = ⋆ E4ν ( 4λ t ). (56) Z Γ (ν + µ 1) − 0 −
4 Integrals of the inverse Laplace transform of µ ν s− exp( s ) with the Mitag-Leffler, Error and − Volterra functions.
The Laplace transform of the two parameter Mittag-Leffler function is [11,34,35]
sα β tβ 1 E ( λtα) = − , L − α,β ± sα λ n o ∞ zk ∓ (57) E (z) = , α,β Γ (kα + β) kX=0 and this permits to obtain
g(t) = tβ 1E ( λtα) ; 0 <ν< 1 − α,β ± sα β s(α β)ν G(s) = − ; G(sν) = − sα λ sαν λ (58) (α ∓β)ν αν (βν∓+µ) 1 s − 1 s − (βν +µ 1) α = − = t − E + ( λt ), (sµ • sαν λ) L ( sαν λ ) αν,βν µ ± ∓ ∓ and ∞ β 1 α (βν +µ 1) α u − Eα,β( λu )fν,µ(u)du = t − Eαν,βν +µ( λt ). (59) Z0 ± ± Thus, in both sides of expressions (59) appear the Mittag-Leffler functions. Evidently, for β = 1 they are reduced to the classical Mittag-Leffler functions.
∞ α (ν +µ 1) α Eα,1( λu )fν,µ(u)du = t − Eαν,ν +µ( λt ), (60) Z0 ± ±
12 and for β = α and β = α + 1 we have
∞ α 1 α (αν +µ 1) α u − Eα,α( λu )fν,µ(u)du = t − Eαν,αν +µ( λt ), Z0 ± ± (61) ∞ α α [(α +1)ν +µ 1] α u Eα,α +1( λu )fν,µ(u)du = t − Eαν,(α +1)ν +µ( λt ). Z0 ± ± If β = 1/2, in the integrands of (61) the Mittag-Leffler functions are expressed then by the error functions [35]
2 E ( z) = ez [1 erf(z)] ; z = λuα, 1/2 ± ± 1 z2 E1/2,1/2( z) = ze [1 erf(z)] , ± ( √z ± ± ) (62) ez E1/2,3/2(z) = erf(z). √z
If β is positive integer, the Mittag-Leffler functions are expressed by elemen- tary functions [36]. In particular cases, with ν = 1/2,µ = 0. ν = µ = 1/2 and ν = 1/3,µ = 2/3 the explicit form of the inverse transforms of exponential func- tions is known (see (15) and (17)). The Laplace transform of the error function is
λ g(t)=erf ; 0 <ν< 1 2t1/2 1/2 ν/2 1 e λs 1 e λs G(s) = − − ; G(sν) = − − s sν ν/2 2/ν ν/2 1 1 e λs 1 λ2µ/ν e (λ s) (63) 1 − − = 1 − = − µ ν − µ 2/ν µ L s • s L s − (λ s) µ 1 ν/2 t − 2(µ 1)/ν t λ − fν/2,µ , Γ (µ) − λ2/ν ! which yields
µ 1 ν/2 ∞ λ t t erf f (u)du = − λ2(µ 1)/νf . (64) 1/2 ν,µ Γ − ν/2,µ 2/ν Z0 2u (µ) − λ ! The Volterra functions are defined by the following integrals [36]
u ∞ t ν(t) = Γ du, Z0 (u + 1) u +α ∞ t ν(t,α) = Γ du, (65) Z0 (u + α + 1) β u +α ∞ u t µ((t,β,α) = Γ Γ du, Z0 (β + 1) (u + α + 1)
13 and their Laplace transforms are
1 L ν(λt) = ; λ> 0 s { } s ln λ λα ν(λt,α) = s (66) L{ } sα +1 ln λ λα µ((λt,β,α) = . L s β +1 sα +1 ln λ The logarithmic functions in the Laplace transforms permit to express the in- tegrals of the Volterra functions with inverse Laplace transform of exponential function in terms of convolution integrals. From
g(t) = ν(λt); λ> 0 ; 0 <ν< 1 1 1 G(s) = ; G(sν) = s sν s ln sν ln λ λ ! (67) 1 1 1 1 1 = 1 = − µ ν 1/ν ν − ν +µ 1 1/ν L (s • s ln[(s/λ ) ]) L (s − • νs ln(s/λ ) ) tν +µ 2 − ν(λ1/ν t) νΓ (ν + µ 1) ∗ − it follows that
ν +µ 2 ∞ t − 1/ν ν(λu)fν,µ(u)du = ⋆ ν(λ t). (68) Z νΓ (ν + µ 1) 0 − Similarly from
g(t) = ν(λt,ρ); λ,ρ > 0 ρ ρ λ ν λ G(s) = ; G(s ) = ν ; 0 <ν< 1 ρ+1 s s s ln s(ρ+1)ν ln λ λ ! 1 1 1 1 λρ 1 = λρ 1 = − µ (ρ+1)ν 1/ν ν − µ 1 (ρ+1)ν +1 1/ν L (s • s ln[(s/λ ) ]) L (s − • νs ln(s/λ )) tµ 2 − ⋆,ν[λ1/νt,(ρ + 1)ν] Γ (µ 1)νλ − (69) we have µ 2 ∞ t − 1/ν ν(λu,ρ)fν,µ du = ⋆ ν[λ t,(ρ + 1)ν]. (70) Z Γ (µ 1)νλ 0 −
14 Finally, the Laplace transform of the generalized Volterra function is
g(t) = µ(λt,ξ,ρ); λ,ρ,ξ > 0; ρ ρ λ ν λ G(s) = ; G(s ) = ν ; 0 <ν< 1, ρ+1 s +1 s s [ln ]ξ s(ρ+1)ν[ln ]ξ +1 λ λ ! ρ 1 1 1 λ − = (71) L (sµ • s(ρ+1)ν ln[(s/λ1/ν )ν] ξ +1 ) λρ 1 { 1 } 1 = ξ +1 − µ 1 (ρ+1)ν +1 1/ν ξ +1 ν L (s − • s [ln(s/λ )] ) tµ 2 − ⋆ µ[λ1/ν t,ξ,(ρ + 1)ν] Γ (µ 1)λνξ +1 − and (65) yields
µ 2 ∞ t − 1/ν µ(λu,ξ,ρ)fν,µ du = ⋆ µ[λ t,ξ,(ρ + 1)ν]. (72) Z Γ (µ 1)λνξ +1 0 − The number of similar convolution integrals can be significantly enlarged if the Volterra functions are multiplied by tn with n = 1,2,3,... , then their Laplace transforms should be differentiated n times. The result of such dif- ferentiations are linear combinations of these functions [33,36]. For example with n =1 and λ = 1 we have d 1 ρ + 1 1 tν(t,ρ) = 1 = 1 + = − ρ+1 − ρ+2 ρ+2 2 L (−ds s lns ) L (s lns s (lns) ) (73) (ρ + 1)ν[t,(ρ + 1)ν]+ µ[t,1,(ρ + 1)ν].
Conclusions
By applying the Efros theorem in the form established by Wlodarski it was possible to derive a number of infinite integrals, finite integrals and integral identities with the function which represent the Laplace inverse transform of µ ν s− exp( s ) with 0 <ν< 1 and 0 µ< 1. The extension to the cases µ 1 is dealt in− Appendix ≤ ≥ Derived by us integrals include in integrands elementary functions (power, exponential, logarithmic, trigonometric and hyperbolic functions) and the er- ror functions, the Mittag-Leffler functions and the Volterra functions. Many results appear in form of the convolution integrals. Performing the inversion by the complex integration, it was possible to show that the inverse Laplace inverse transform, which means the original function, can be also expressed in terms of the Wright functions and for par- ticular values of parameters by the Mainardi functions. Using rules of operational calculus some properties of the inverse Laplace transform were derived.
15 Acknowledgments
The research work of F.M. has been carried out in the framework of the activi- ties of the National Group of Mathematical Physics (GNFM, INdAM). Both the authors would like to acknowledge the unknown reviewers for their constric- tive comments.
Appendix: Additional properties of the inverse of µ ν the Laplace transform of s− exp( s ) − The restriction posed on the second parameter of the function fν,µ(t), i.e. µ< 1 µ ν (in the inverse of the Laplace transform s− exp( s )) can be removed by using rules of operational calculus. In the general case,− for a non-negative param- eter ρ 0, this inverse transform can always be expressed as the convolution integral≥ which includes the corresponding power function:
sν sν λ 1 1 e− 1 1 e− t − f (t) = − = − = ⋆f (t) ν,ρ L ( sρ ) L (sλ • sµ ) Γ (λ) ν,µ (A 1) − ρ = λ + µ In the particular case µ = 1, it reduces to the simple integral, because 1/s ex- presses the integration operation: t 1 1 sν fν,1(t) = − e− = fν,0(u)du. (A 2) L s Z0 − The product of two identical or different Laplace inverse transforms can be derived in the form of convolution integrals from
sν 2sν (21/ν s)ν e sν e e 1 − e− 1 − 2µ/ν 1 − − µ = − = 2 − L ( sµ • s ) L ( s2µ ) L (21/ν s)2µ (A 3) sν 2sν (21/ν s)ν e sν e e − 1 − e− 1 − (µ +ρ)/ν 1 − ρ = = 2 − µ s − µ+ρ − 1/ν (µ+ρ) L ( s • ) L ( s ) L (2 s) which gives t f (t) ⋆ (t) = f (t ξ)f (ξ)dξ = 2(2µ 1)/νf t ν,µ ν,µ ν,µ − ν,µ − ν,2µ 21/ν Z0 (A 4) t − f (t) f (t) = f (t ξ)f (ξ)dξ = 2(µ+ρ 1)/νf t ν,µ ∗ ν,ρ ν,µ − ν,ρ − ν,µ+ρ 21/ν Z0 These results can be generalized to n-fold integrals when exists the factor 1/sn, with n = 1,2,3,.... In the case of the product of identical Laplace inverse trans- forms it is reduced to the following convolution integral
t n 1 1 1 (t u) − − F(s) = − f (u)du. (A 5) L sn Z (n 1)! − 0 −
16 In similar way, if the inverse transform is in a more general form
αsν βsν e− e− fν,µ(t);α = ; fν,ρ(t);β = ; α,β > 0 L L( sµ ) L L( sµ ) n αsνo βsν (αn +β)sν o 1 e− e− 1 e− − µ ρ = − µ+ρ = (A 6) L ( s • s ) L ( s ) − 1/ν ν e [(α +β) s] (α + β)(µ+ρ)/ν 1 − − 1/ν µ+ρ L [(α + β) s] we have t fν,µ(t;α) ⋆fν,ρ(t;β) = fν,µ(t ξ;α)fν,ρ(ξ;β)dξ = Z0 − (A 7) − (α + β)(µ+ρ 1)/νf t ;(α + β) − ν,µ+ρ ( + )1/ν α β The differentiation of Laplace inverse transforms with respect to the parame- ters ν and µ gives
ν sν sν ∂fν,µ(t) s e lns e ln(sC) sν − − lnCe− γ = = + µ ν ; C = e µ µ ν 1 s − L( ∂ν ) − s −s − − • s sν (A 8) ∂fν,µ(t) µe − = − L( ∂µ ) − sµ+1 which yields the following inverses of (A-8)
∂fν,µ(t) = lnt⋆fν,µ ν 1(t)+ γfν,µ ν(t) ∂ν − − − t (A 9) ∂fν,µ(t) − ∂µ = µfν,µ+1(t) = µ fν,µ(ξ)dξ, L − − Z0 where γ denotes the Euler constant.
References
[1] Humbert, P. Nouvelle correspondances symboliques. Bull. Soc. Math., France, 69 (1945) 121–129. [2] Pollard, H. The representation of exp( xλ) as a Laplace integral. Bull. Amer. Math. Soc., 52 (1946) 908–910. − [3] Wlodarski, L. Sur une formule de Eftros. Studia Math., 13 (1952) 183–187. [4] Mikusinski, J. Sur les fonctions exponentielles du calcul operatoire.´ Studia Math., 12 (1951) 208–224. [5] Mikusinski, J. Sur la croissance de la function operationelle´ exp( sαλ). Bull. Acad. Polon. 4 (1953) 423–425. − [6] Mikusinski, J.. Sur la function dont la transformee´ de Laplace est exp( sαλ). Bull. Acad. Polon. 6 (1958) 691–693.. −
17 [7] Mikusinski, J. On the function whose Laplace transform is exp( sλ). Studia Math., 18 (1959) 195–198. − [8] Wintner, A. Cauchy’s stable distributions and an ”explicit formula” of Mellin. Amer. J. Math., 78 (1956) 819–861. [9] Ragab, F.M. The inverse Laplace transform of a exponential function. Comm. Pure Appl. Math., 11 (1958) 115–127. [10] Stankovic,ˇ B. On the function of E.M. Wright. Publications de L’Institute Mathematique.´ Nouvelle serie,´ 10 (1970) 113–124. [11] Mainardi, F.. Fractional Calculus and Waves in Linear Viscoelasticity, Im- perial College Press, London, 2010. [12] Hanyga, A. Multidimensional solutions of time-fractional diffusion-wave equations. Proc. Royal Soc. A., London, 458 (2002) 933–957. [13] Barkai, E. Fractional Fokker-Planck equation, solution, and application. Phys. Rev. E. 63 (2001) 046118/1–17. [14] Penson, K.A., Gorska,´ K. Exact and explicit probability densities for one- sided Levy´ stable distribution. Phys. Rev. Letters, 105 (2010) 210604/1–9. [15] Gorska,´ K., Penson, K.A. Levy´ stable distributions via associated integral transform. J. Math. Phys., 53 (2012), 053302/1–10. [16] Efros, A.M. Some applications of operational calculus in analysis. Mat. Sbornik, 42 (1935) 699–705. [17] Apelblat, A., Mainardi, F. Application of the Efros theorem to the Wright functions of the second kind and other results. Lecture Notes of TICMI (Tbilisi Int. Center of Mathematics and Informatics) , 21 (2020), 9–28. Special issue edited by P.E. Ricci, dedicated to Prof. George Jaiani on the occasion of his 75th birthday, Tbilisi University Press, ISSN 1512-0511. [18] Erdelyi,´ A., Magnus, W., Oberhettinger, F., Tricomi, F.G. Tables of Integral Transforms. McGraw-Hill, New York, 1954 [19] Roberts, G.E., Kaufman, H. Tables of Laplace Transforms. W.B. Saunders Co., Philadelphia, 1966. [20] Hladik, J. La Transformation de Laplace a Plusieurs Variables. Masson et Cie Editeurs.´ Paris, 1969. [21] Oberhettinger, F., Badii, L. Tables of Laplace Transforms. Springer-Verlag, Berlin, 1973. [22] Wright, E.M. On the coefficients of power series having exponential sin- gularities. J. London Math. Soc., 8 (1933) 71–79. [23] Wright, E.M. The generalized Bessel function of order greater than one. Quart. J. Math. Oxford, 11 (1940) 36–48. [24] Montroll, E.W., Bendler, J.T. On Levy´ (or stable) distributions and the Williams-Watts model of dielectric relaxation. J. Stat. Phys., 34 (1984) 129– 192. [25] Zolotarev, V.M. One-dimensional Stable Distributions. Translated from Russian by H.H. McFaden. Amer. Math. Soc., Providence, Rhode Island, 1986. [26] Mainardi, F. On the initial value problem for the fractional diffusion-wave equation, in S. Rionero and T. Ruggeri (Editors), 7th Conference on Waves and Stability in Continuous Media (WASCOM 1993), World Scientific, Singapore, 1994, pp. 246–251.
18 [27] Saichev, A.I., Zaslavsky, G.M. Fractional kinetic equations :solutions and applications. Chaos, 7 (1997) 753–764. [28] Mainardi, F., Tomirotti, M.. Seismitic pulse propagation with constant Q and stable probability distributions. Ann. Geofisica, 40 (1997) 1311–1328. [29] Gorenflo, R., Luchko, Yu., Mainardi, F.. Analytical properties and applica- tions of the Wright functions. Fract. Calc. Appl. Anal., 2 (1999) 383–414. [30] F. Mainardi, Yu. Luchko, G. Pagnini. The fundamental solution of the space-time fractional diffusion equation. Fract. Calc. Appl. Anal., 4 (2001) 153–192. [31] Mainardi, F., Mura, A., Pagnini, G., Gorenflo, R. Fractional relaxation and time- fractional diffusion of distributed order. Proc. of the 2nd IFAC Workshop on Fractional Differentiation and its Applications. Porto Portugal, 2006. [33] Garg,, M., Rao, A. Fractional extensions of some boundary value problems in oil strata. Proc. India Acad. Sci. (Math. Sci.) 117 (2007) 267–281. [33] Mainardi, F., Consiglio, A. The Wright functions of the second kind in Mathematical Physics. Mathematics (MDPI), 8 No 6 (2020), 884/1–26. DOI: 10.3390/MATH8060884; E-print arXiv:2007.02098 [math.GM] [34] Apelblat, A. Laplace Transforms and Their Applications. Nova Science Publishers, Inc., New York, 2012. [35] Apelblat, A., Differentiation of the Mittag-Leffler functions with respect to parameters in the Laplace transform approach. Mathematics (MDPI), 8 (2020), 657/1–29. DOI:10.3390/math8050657 [36] Gorenflo, R., Kilbas, A.A., Mainardi, F., Rogosin, S.V. Mittag-Leffler Func- tions, Related Topics and Applications. 2nd Edition, Springer, Heidelberg, 2020. [37] Apelblat, A. Volterra Functions. Nova Science Publishers, Inc., New York, 2008.
19 This figure "Bromwich2-contour.jpg" is available in "jpg" format from:
http://arxiv.org/ps/2012.07068v3