Arxiv:1904.11370V1 [Math.GM] 19 Apr 2019

NEW INTEGRAL TRANSFORM: SHEHU TRANSFORM A GENERALIZATION OF SUMUDU AND LAPLACE TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS

SHEHU MAITAMA∗, WEIDONG ZHAO

Abstract. In this paper, we introduce a Laplace-type integral transform called the Shehu transform which is a generalization of the Laplace and the Sumudu integral transforms for solving differential equations in the time domain. The proposed integral transform is successfully derived from the clas- sical Fourier integral transform and is applied to both ordinary and partial differential equations to show its simplicity, efficiency, and the high accuracy.

1. Introduction Historically, the origin of the integral transforms can be traced back to the work of P. S. Laplace in 1780s and Joseph Fourier in 1822. In recent years, differential and integral equations have been solved using many integral transforms ([1]-[11]). The Laplace transform, and Fourier integral transforms are the most commonly used in the literature. The Fourier integral transform [12] was named after the French mathematician Joseph Fourier. Mathematically, Fourier integral transform is defined as: 1 Z ∞ z[f(t)] = f(ω) = √ exp (−iωt) f(t)dt. (1.1) 2π −∞ The Fourier transform have many applications in physics and engineering processes [13]. The Laplace integral transform is similar with the Fourier transform and is defined as: Z ∞ £[f(t)] = F (s) = exp (−st) f(t)dt. (1.2) −∞ The Laplace transform is highly efficient for solving some class of ordinary and partial differential equations [14]. By replacing the variable iω with the variable s in Equ.(1.1), the well-known Fourier transform will become a Laplace transform and the vice-versa. The only difference between the Laplace transform, and the arXiv:1904.11370v1 [math.GM] 19 Apr 2019 Fourier transform is that the Laplace transform can be defined for both stable and unstable system while the Fourier transform can only be defined on a stable system. In mathematical literature, the discrete-time equivalent of the Laplace transform called z-transform [15] converts a discrete-time signal into a complex frequency-domain representation. The basic idea of the z-transform was known to Laplace and later it was re-introduced by the Jewish-Polish mathematician Witold Hurewicz to treat a sampled-data control systems used with radar in 1947 ([16]- [17]). In mathematics and signal processing, the bilateral or two-sided z-transform

2010 Mathematics Subject Classification. 44A10, 44A15, 44A20, 44A30, 44A35. Key words and phrases. Shehu transform; Fourier integral transform; Laplace transform; natural transform; Sumudu transform; ordinary and partial differential equations. 1 NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 2 of a discrete-time signal x[n] is the normal power series X(z) which is defined as: ∞ X X(z) = Z{x[n]} = x[n]z−n, (1.3) n=−∞ where n is an integer and z is in general a complex number [18]. The multiplicative version of the two-sided Laplace transform called the Mellin integral transform is defined as [19]: Z ∞ M[f(s); s] = f ∗(s) = xs−1f(x)dx. (1.4) 0 The Mellin integral transform is similar with the Laplace transform and Fourier transform and is widely applied in computer science and number theory due to its invariant property [20, 21]. In railway engineering, the Laplace-Carson transform [22] which is a Laplace-type integral transform named after Pierre Simon Laplace and John Renshaw Carson is defined as:

Z ∞ ˆ fC (p) = p exp (−pt) f(t)dt, t ≥ 0. (1.5) 0 The Laplace-Carson integral transform have many applications in physics and engineering and can easily be converted into a Mellin deconvolution problem, see ([23],[24]). In mathematics, the Hankel’s integral transform [25] which is similar to the Fourier transform was first introduced by the German mathematician Hermann Hankel and was widely used in physical science and engineering [26]. The Hankel’s transform is defined as: Z ∞ Fv(s) = Hv[f(r)] = rf(r)Jv(sr)dr, r ≥ 0, (1.6) 0 1 where Jv is the Bessel function of the first kind of order v with v ≥ − 2 . In 1993, Watugala introduced a Laplace-like integral transform called the Sumudu integral transform [27]. In recent years, Sumudu transform has been applied to many real-life problems because of its scale and unit preserving properties ([28]- [31]). The mathematical definition of the Sumudu transform is given by: 1 Z ∞ −t S[f(t)](u) = G(u) = exp f(t)dt, (1.7) u 0 u provided the integral exists for some u. Based on the basic idea of the Laplace and the Sumudu integral transform, the Elzaki transform was proposed in 2011. The Elzaki transform is closely related with the Laplace transform, Sumudu transform, and the natural transform. Elzaki transform is defined as [32]: Z ∞ −t E[f(t)] = T (u) = u exp f(t)dt, (1.8) 0 u provided the integral exists for some u. The natural transform [33] which is similar to Laplace and Sumudu integral transform was introduced in 2008. In recent years, natural transform was successfully applied to many applications (see [34, 35]). The natural transform is defined by the following integral: Z ∞ + 1 −st N [f(t)](s, u) = R(s, u) = exp f(t)dt, s > 0, u > 0, (1.9) u 0 u NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 3 provided the integral exists for some variables u and s. Recently, a new integral transform called the M-transform which is also similar to natural transform is introduced by Srivastava et al. in 2015. Mathematically speaking, M-transform is closely connected with the well-known Laplace transform and the Sumudu integral transform. M-transform was successfully applied to first order initial-boundary value problem (see Srivastava et al. [36]). The M-transform is defined as: Z ∞ exp (−ut) f(vt) Mρ,m[f(t)](u, v) = m m ρ dt, (1.10) 0 (t + v )

(ρ ∈ C; <(ρ) ≥ 0, m ∈ Z+ = 1, 2, 3, ···) , where both u ∈ C and v ∈ R+ are the M-transform variables. In 2013, Atangana and Kilicman introduced a novel integral transform called the Abdon-Kilicman integral transform [37] for solving some differential equations with some kind of singularities. The novel integral transform is defined as: Z ∞ n Mn(s) = Mn[f(x)](s) = x exp (−xs) f(x)dx. (1.11) 0 The Atangana-Kilicman integral becomes Laplace transform when n = 0. Recently, a Laplace-type integral transform called the Yang transform ([38]-[40]) for solving steady heat transfer problems was introduced in 2016. The integral transform is defined as: Z ∞ −τ Y [φ(τ)] = φ(ω) = exp φ(τ)dτ, (1.12) 0 ω provided the integral exists for some ω. Due to the rapid development in the physical science and engineering models, there are many other integral transforms in the literature. However, most of the existing integral transforms have some limitations and cannot be used directly to solved nonlinear problems or many complex mathematical models. As a result, many authors became highly interested to come up with the alternative approach for solving many real-life problems. In 2016, Atangana and Alkaltani introduced a new double integral equation and their properties based on the Laplace transform and decomposition method. The double integral transform was successfully applied to second order partial differential equation with singularity called the two- dimensional Mboctara equation [41]. Recently, Eltayeb applied double Laplace decomposition method to nonlinear partial differential equations [42]. In 2017, Bel- gacem el at. extended the applications of the natural and the Sumudu transforms to fractional diffusion and Stokes fluid flow realms [43]. Motivated by the above-mentioned researches, in this paper we proposed a Laplace-type integral transform called Shehu transform for solving both ordinary and partial differential equations. The Laplace-type integral transform converges to Laplace transform when u = 1, and to Yang integral transform when s = 1. The proposed integral transform is successfully applied to both ordinary and partial differential equations. All the results obtained in the applications section can easily be verified using the Laplace or Fourier integral transforms. Throughout this paper, the Shehu transform is denoted by an operator S[.]. NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 4

2. Main result Deﬁnition 1. The Shehu transform of the function v(t) of exponential order is deﬁned over the set of functions, n |t| i o A = v(t): ∃ N, η1, η2 > 0, |v(t)| < N exp , if t ∈ (−1) × [0, ∞) , ηi by the following integral Z ∞ −st S [v(t)] = V (s, u) = exp v(t)dt 0 u Z α −st = lim exp v(t)dt; s > 0, u > 0. (2.1) α→∞ 0 u It converges if the limit of the integral exists, and diverges if not. The inverse Shehu transform is given by −1 S [V (s, u)] = v(t), for t ≥ 0. (2.2) Equivalently

Z α+i∞ −1 1 1 st v(t) = S [V (s, u)] = exp V (s, u) ds, (2.3) 2πi α−i∞ u u where s and u are the Shehu transform variables, and α is a real constant and the integral in Equ.(2.3) is taken along s = α in the complex plane s = x + iy. Theorem 1. The sufficient condition for the existence of Shehu transform. If the function v(t) is piecewise continues in every finite interval 0 ≤ t ≤ β and of exponential order α for t > β. Then its Shehu transform V (s, u) exists. Proof. For any positive number β, we algebraically deduce Z ∞ st Z β st Z ∞ st exp − v(t)dt = exp − v(t)dt + exp − v(t)dt. (2.4) 0 u 0 u β u Since the function v(t) is piecewise continues in every finite interval 0 ≤ t ≤ β, then the first integral on the right hand side exists. Besides, the second integral on the right hand side exists, since the function v(t) is of exponential order α for t > β. To verify this claim, we consider the following case Z ∞ Z ∞ st st exp − v(t)dt ≤ exp − v(t) dt β u β u Z ∞ st ≤ exp − |v(t)| dt 0 u Z ∞ st ≤ exp − N exp (αt) dt β u Z ∞ (s − αu)t = N exp − dt β u uN (s − αu)tγ = − lim exp − γ→∞ (s − αu) u 0 uN = . s − αu The proof is complete. NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 5

Property 1. Linearity property of Shehu transform. Let the functions αv(t) and βw(t) be in set A, then (αv(t) + βw(t)) ∈A, where α and β are nonzero arbitrary constants, and

S [αv(t) + βw(t)] = αS [v(t)] + βS [w(t)] . (2.5) Proof. Using the Deﬁnition 1 of Shehu transform, we get Z ∞ −st S [αv(t) + βw(t)] = exp (αv(t) + βw(t))dt 0 u Z ∞ −st Z ∞ −st = exp (αv(t)) dt + exp (βw(t))dt 0 u 0 u Z ∞ −st Z ∞ −st = α exp v(t)dt + β exp w(t)dt 0 u 0 u Z ∞ Z ∞ = αu exp (−st) v(ut)dt + βu exp (−st) w(ut)dt 0 0 = αS [v(t)] + βS [w(t)] . The proof is complete. In particular, using the Deﬁnition 1 and Property 1, we obtain

S [3 cos(t) + 5 sin(2t)] = 3S [cos(t)] + 5S [sin(2t)]

3us 5u2 = + , s2 + u2 s2 + (2u)2 see entries of table 1. Property 2. Change of scale property of Shehu transform. Let the function v(βt) be in set A, where β is an arbitrary constant. Then u s [v(βt)] = V , u . (2.6) S β β Proof. Using the Deﬁnition 1 of Shehu transform, we deduce Z ∞ −st S [v(βt)] = exp v(βt)dt (2.7) 0 u η dη Substituting η = βt which implies t = β and dt = β in Equ.(2.7) yields 1 Z ∞ −sη S [v(βt)] = exp v(η)dη β 0 uβ 1 Z ∞ −st = exp v (t) dt β 0 uβ u Z ∞ −st = exp v(ut)dt β 0 β u s = V , u . β β

The proof is complete. NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 6

Theorem 2. Derivative of Shehu transform. If the function v(n)(t) is the nth derivative of the function v(t) ∈ A with respect to 0t0, then its Shehu transform is deﬁned by n−1 h i sn X s n−(k+1) v(n)(t) = V (s, u) − v(k)(0). (2.8) S un u k=0 When n=1, 2, and 3 in Equ. (2.8) above, we obtain the following derivatives with respect to t. s [v0(t)] = V (s, u) − v(0). (2.9) S u

s2 s [v00(t)] = V (s, u) − v(0) − v0(0). (2.10) S u2 u

s3 s2 s [v000(t)] = V (s, u) − v(0) − v0(0) − v00(0). (2.11) S u3 u2 u Proof. Now suppose Equ. (2.8) is true for n = k. Then using Equ. (2.9) and the induction hypothesis, we deduce h i s h i (v(k)(t))0 = v(k)(t) − v(k)(0) S uS " k−1 # s sk X s k−(i+1) = [v(t)] − v(i)(0) − v(k)(0) u uk S u i=0 k s k+1 X s k−i = [v(t)] − v(i)(0), u S u i=0 which implies that Equ. (2.8) holds for n = k + 1. By induction hypothesis the proof is complete. The following important properties are obtain using the Leibniz’s rule ∂v(x, t) Z ∞ −st ∂v(x, t) ∂ Z ∞ −st S = exp dt = exp v(x, t) dt ∂x 0 u ∂x ∂x 0 u ∂ ∂v(x, t) d = [V (x, s, u)] ⇒ = [V (x, s, u)] , ∂x S ∂x dx

∂2v(x, t) Z ∞ −st ∂2v(x, t) ∂2 Z ∞ −st S 2 = exp 2 dt = 2 exp v(x, t) dt ∂x 0 u ∂x ∂x 0 u ∂2 ∂2v(x, t) d2 = [V (x, s, u)] ⇒ = [V (x, s, u)] , ∂x2 S ∂x2 dx2 and ∂nv(x, t) Z ∞ −st ∂nv(x, t) ∂n Z ∞ −st S n = exp n dt = n exp v(x, t)dt ∂x 0 u ∂x ∂x 0 u ∂n ∂nv(x, t) dn = [V (x, s, u)] ⇒ = [V (x, s, u)] . ∂xn S ∂xn dxn NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 7

3. Some useful results of Shehu transform

Property 3. Let the function v(t) = 1 be in set A. Then its Shehu transform is given by u [1] = . (3.1) S s Proof. Using Equ.(2.1), we deduce Z ∞ −st u −stγ u [1] = exp dt = − lim exp = . S γ→∞ 0 u s u 0 s This ends the proof. Property 4. Let the function v(t) = t be in set A. Then its Shehu transform is given by u2 [t] = . (3.2) S s2 Proof. Using the Deﬁnition 1 of the Shehu transform and integration by parts, we get Z ∞ −st u −stγ u Z ∞ −st [t] = t exp dt = lim t exp + exp dt S γ→∞ 0 u s u 0 s 0 u u2 −stγ u2 = − lim exp = . 2 γ→∞ 2 s u 0 s Thus the proof ends. tn Property 5. Let the function v(t) = n! n = 0, 1, 2.. be in set A. Then its Shehu transform is given by tn un+1 = . (3.3) S n! s Proof. From the Deﬁnition 1 of the Shehu transform and integration by parts, we deduce Z ∞ Z ∞ n n −st u n−1 −st S [t ] = t exp dt = n t exp dt 0 u s 0 u 2 Z ∞ u n−2 −st = 2 n(n − 1) t exp dt s 0 u 3 Z ∞ u n−3 −st = 3 n(n − 1)(n − 2) t exp dt s 0 u 4 Z ∞ u n−4 −st = 4 n(n − 1)(n − 2)(n − 3) t exp dt s 0 u 5 Z ∞ n+1 u n−5 −st u = 5 n(n − 1)(n − 2)(n − 3)(n − 4) t exp dt = ··· = n! . s 0 u s The proof is completed. tn Property 6. Let the function v(t) = Γ(n+1) n = 0, 1, 2, ··· be in set A. Then its Shehu transform is given by tn un+1 = . (3.4) S Γ(n + 1) s NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 8

The proof of property 6 follows immediately from the previous property 5. Property 7. Let the function v(t) = exp(αt) be in A. Then its Shehu transform is given by u [exp(αt)] = . (3.5) S s − αu Proof. Using Equ.(2.1), we get Z ∞ (s − αu)t S [exp(αt)] = exp − dt 0 u u (s − αu)tγ u = − lim exp − = . γ→∞ s − αu u 0 s − αu This ends the proof. Property 8. Let the function v(t) = t exp(αt) be in set A. Then its Shehu transform is given by u2 [t exp(αt)] = . (3.6) S (s − αu)2 Proof. Using the Deﬁnition 1 of the Shehu transform and integration by parts, we get Z ∞ (s − αu)t S [t exp(αt)] = t exp − dt 0 u u (s − αu)tγ = − lim t exp − γ→∞ s − αu u 0 u Z ∞ (s − αu)t + exp − dt s − αu 0 u u2 (s − αu)tγ u2 = − lim exp − = . 2 γ→∞ 2 (s − αu) u 0 (s − αu) The proof is complete. tn exp(αt) Property 9. Let the function v(t) = n! n = 0, 1, 2, ... be in set A. Then its Shehu transform is given by tn exp(αt) un+1 = . (3.7) S n! (s − αu)n+1 Proof. Using the Deﬁnition 1 of the Shehu transform and integration by parts, we deduce Z ∞ n n (s − αu)t S [t exp(αt)] = t exp − dt 0 u un Z ∞ (s − αu)t = tn−1 exp − dt (s − αu) 0 u 2 Z ∞ u n(n − 1) n−2 (s − αu)t n! = 2 t exp − dt = ··· = n+1 . (s − αu) 0 u (s − αu) Thus the proof is complete. NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 9

tn Property 10. Let the function v(t) = Γ(n+1) exp(αt) n = 0, 1, 2, ... be in set A. Then its Shehu transform is given by tn exp(αt) un+1 = . (3.8) S Γ(n + 1) (s − αu)n+1

The proof of Property 10 follows as a direct consequence of Property 9. Property 11. Let the function v(t) = sin(αt) be in set A. Then its Shehu transform is given by αu2 [sin(αt)] = . (3.9) S s2 + α2u2 Proof. Using the Deﬁnition 1 of the Shehu transform and integration by parts, we get Z ∞ st S [sin(αt)] = exp − sin(αt)dt 0 u u st γ uα Z ∞ st = − lim exp − sin(αt) + exp − cos(αt)dt γ→∞ s u 0 s 0 u αu2 st γ α2u2 Z ∞ st = − lim exp − cos(αt) − exp − sin(αt)dt 2 γ→∞ 2 s u 0 s 0 u αu2 α2u2 Z ∞ st = 2 − 2 exp − sin(αt)dt. s s 0 u

Simplifying the required integrals complete the proof of Property 11. Property 12. Let the function v(t) = cos(αt) be in set A. Then its Shehu transform is given by us [cos(αt)] = . (3.10) S s2 + α2u2 Proof. Using the Deﬁnition 1 of the Shehu transform and integration by parts, we deduce Z ∞ st S [cos(αt)] = exp − cos(αt)dt 0 u u st γ αu Z ∞ st = − lim exp − cos(αt) − exp − sin(αt)dt γ→∞ s u 0 s 0 u u αu2 st γ α2u2 Z ∞ st = − lim exp − sin(αt) − exp − cos(αt)dt 2 γ→∞ 2 s s u 0 s 0 u u α2u2 Z ∞ st = − 2 exp − cos(αt)dt. s s 0 u

Simplifying the required integrals complete the proof of Property 12. sinh(αt) Property 13. Let the function v(t) = α be in set A. Then its Shehu transform is given by sinh(αt) u2 = . (3.11) S α s2 − α2u2 NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 10

Proof. From the Deﬁnition 1 of the Shehu transform and integration by parts, we get

Z ∞ st S [sinh(αt)] = exp − sinh(αt)dt 0 u u st γ uα Z ∞ st = − lim exp − sinh(αt) + exp − cosh(αt)dt γ→∞ s u 0 s 0 u αu2 st γ α2u2 Z ∞ st = − lim exp − cos(αt) + exp − sinh(αt)dt 2 γ→∞ 2 s u 0 s 0 u αu2 α2u2 Z ∞ st = 2 + 2 exp − sinh(αt)dt. s s 0 u

Simplifying the required integrals complete the proof of Property 13.

Property 14. Let the function v(t) = cosh(αt) be in set A. Then its Shehu transform is given by

us [cosh(αt)] = . (3.12) S s2 − α2u2

Proof. Applying the Deﬁnition 1 of the Shehu transform and integration by parts, we get

Z ∞ st S [cosh(αt)] = exp − cosh(αt)dt 0 u u st γ αu Z ∞ st = − lim exp − cos(αt) + exp − sinh(αt)dt γ→∞ s u 0 s 0 u u αu2 st γ α2u2 Z ∞ st = − lim exp − sinh(αt) + exp − cos(αt)dt γ→∞ 2 s s u 0 s 0 u u α2u2 Z ∞ st = + 2 exp − cosh(αt)dt. s s 0 u

Collecting the required integrals complete the proof of Property 14.

exp(βt) sin(αt) Property 15. Let the function α be in set A. Then its Shehu transform is given by

exp (βt) sin(αt) u2 = . (3.13) S α (s − βu)2 + α2u2 NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 11

Proof. Using the Deﬁnition 1 of the Shehu transform and integration by parts, we deduce Z ∞ (s − βu) S [exp (βt) sin(αt)] = exp − t sin(αt)dt 0 u −u (s − βu) γ = lim exp − t sin(αt)dt γ→∞ (s − βu) u 0 uα Z ∞ (s − βu) + exp − t cos(αt)dt s − βu 0 u u2α (s − βu) γ = − lim exp − t cos(αt)dt 2 γ→∞ (s − βu) u 0 αu2α2 Z ∞ (s − βu) − 2 exp − t sin(αt)dt (s − βu) 0 u u2α u2α2 Z ∞ (s − βu) = 2 − 2 exp − t sin(αt)dt. (s − βu) (s − βu) 0 u Simplifying the required integrals complete the proof of property 15. This ends the proof. Property 16. Let the function exp (βt) cos(αt) be in set A. Then its Shehu transform is given by u(s − αu) [exp (βt) cos(αt)] = . (3.14) S (s − βu)2 + α2u2 Proof. Applying the Deﬁnition 1 of the Shehu transform and integration by parts, we get Z ∞ (s − βu) S [exp (βt) cos(αt)] = exp − t cos(αt)dt 0 u u (s − βu) γ = − lim exp − t cos(αt) γ→∞ s − βu u 0 αu Z ∞ (s − βu) + exp − t sin(αt)dt s − βu 0 u u αu Z ∞ (s − βu) = + exp − t sin(αt)dt s − βu s − βu 0 u u αu2 (s − βu) γ = + lim exp − t sin(αt) 2 γ→∞ s − βu (s − αu) u 0 α2u2 Z ∞ (s − βu) − 2 exp − t cos(αt)dt (s − βu) 0 u u α2u2 Z ∞ (s − βu) = − 2 exp − t cos(αt)dt. s − βu (s − βu) 0 u Simplifying the required integrals complete the proof of property 16. exp(βt)−exp(αt) Property 17. Let the function β−α be in set A. Then its Shehu transform is given by exp (αt) u2 = . (3.15) S β − α (s − αu)(s − βu) NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 12

Proof. Using the deﬁnition of Shehu transform, we get exp (αt) u Z ∞ −st S = exp (exp(βt) − exp (αt)) dt β − α β − α 0 u Z ∞ Z ∞ 1 (βu−s) t 1 (αu − s)t = e u dt − exp dt β − α 0 β − α 0 u u (s − βu)tγ = lim exp − γ→∞ (β − α)(βu − s) u 0 u (s − βu)tγ − lim exp − γ→∞ (β − α)(αu − s) u 0 u u = − + (β − α)(βu − s) (β − α)(αu − s) −u(αu − s) + u(βu − s) u2 = = . (β − α)(αu − s)(βu − s) (s − αu)(s − βu) The proof is complete. β exp(βt)−α exp(αt) Property 18. Let the function β−α be in set A. Then its Shehu transform is given by: β exp (βt) − α exp (αt) us = . (3.16) S β − α (s − αu)(s − βu) Proof: Using the deﬁnition of Shehu transform, we get β exp (βt) − α exp (αt) 1 Z ∞ −st S = exp (β exp (βt) − α exp (αt)) dt β − α β − α 0 u β Z ∞ (βu − s) α Z ∞ (αu − s)t = exp t dt − exp dt β − α 0 u β − α 0 u uβ (s − βu)tγ = lim exp − γ→∞ (β − α)(βu − s) u 0 uα (s − αu)tγ − lim exp − γ→∞ (β − α)(αu − s) u 0 uβ uα = − + (β − α)(βu − s) (β − α)(αu − s) −uβ(αu − s) + uα(βu − s) us = = . (β − α)(αu − s)(βu − s) (s − αu)(s − βu) This ends the proof. More properties of the Shehu transform and their converges to the natural transform, the Sumudu transform, and the Laplace transform are presented in table 1. The comprehensive summary of Shehu transform properties are presented in table 2.

4. Applications In this section, the applications of the proposed transform are presented. The simplicity, eﬃciency and high accuracy of the Shehu transform are clearly illus- trated. NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 13

Example 1. Consider the following first order ordinary differential equation dv(t) + v(t) = 0, (4.1) dt subject to the initial condition v(0) = 1. (4.2) Applying the Shehu transform on both sides of Equ. (4.1), we get s V (s, u) − v(0) + V (s, u) = 0. (4.3) u Substituting the given initial condition and simplifying, we deduce u V (s, u) = (4.4) s + u Taking the inverse Shehu transform of Equ. (4.4), yields v(t) = exp(−t). (4.5) Example 2. Consider the following second order ordinary differential equation d2v(t) dv(t) + = 1 (4.6) dt2 dt subject to the initial conditions dv(0) v(0) = 0, = 0. (4.7) dt Applying the Shehu transform on both sides of Equ. (4.6), we obtain s2 s s u V (s, u) − v(0) − v0(0) + V (s, u) − v(0) = . (4.8) u2 u u s Substituting the given initial conditions and simplifying, we deduce u u2 u V (s, u) = − + + . (4.9) s s2 s + u Taking the inverse Shehu transform of Equ. (4.9), we get v(t) = −1 + t + exp(−t). (4.10) Example 3. Consider the following second nonhomogeneous order ordinary differential equation d2v(t) dv(t) − 3 + 2v(t) = exp(3t) (4.11) dt2 dt subject to the initial conditions dv(0) v(0) = 1, = 0. (4.12) dt Applying the Shehu transform on both sides of Equ. (4.11), yields s2 s s u V (s, u) − v(0) − v0(0) − 3 V (s, u) − v(0) + 2V (s, u) = . (4.13) u2 u u s − 3u Substituting the given initial conditions and simplifying, we obtain 5 u u 1 u V (s, u) = − 2 + . (4.14) 2 (s − u) s − 2u 2 (s − 3u) Taking the inverse Shehu transform of Equ. (4.14), we get 5 1 v(t) = exp(t) − 2 exp(2t) + exp(3t). (4.15) 2 2 NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 14

Example 4. Consider the following ordinary diﬀerential equation d2v(t) dv(t) + 2 + 5v(t) = exp(−t) sin(t) (4.16) dt2 dt subject to the initial conditions dv(0) v(0) = 0, = 1. (4.17) dt Applying the Shehu transform on both sides of Equ. (4.16), we get s2 s s s u2 V (s, u) − v(0) − v0(0) + 2 V (s, u) − v(0) + 5V (s, u) = u2 u u u (s + u)2 + u2 (4.18) Substituting the given initial conditions and simplifying, we get 1 u2 2 u2 V (s, u) = + (4.19) 3 ((s + u)2 + u2) 3 ((s + u)2 + (2u)2) Taking the inverse Shehu transform of Equ. (4.19), we get 1 2 v(t) = exp(−t) sin(t) + exp(−t) sin(2t). (4.20) 3 3 Example 5. Consider the following homogeneous partial diﬀerential equation ∂v(x, t) ∂2v(x, t) = (4.21) ∂t ∂x2 subject to the boundary and initial conditions v(0, t) = 0, v(1, t) = 0, v(x, 0) = 3sin(2πx). (4.22) Applying the Shehu transform on both sides of Equ. (4.21), we get s d2V (x, s, u) V (x, s, u) − v(x, 0) = . (4.23) u dx2 Substituting the given initial condition and simplifying, we get d2V (x, s, u) s − V (x, s, u) = −3sin(2πx). (4.24) dx2 u The general solution of Equ. (4.24) can be written as

V (x, s, u) = Vh(x, s, u) + Vp(x, s, u), (4.25) where Vh(x, s, u) is the solution of the homogeneous part which is given by r s r s V (x, s, u) = α exp x + α exp − x , (4.26) h 1 u 2 u and Vp(x, s, u) is the solution of the nonhomogeneous part which is given by

Vp(x, s, u) = β1sin(2πx) + β2cos(2πx). (4.27) Applying the boundary conditions on Equ. (4.26), we get r s r s α + α = 0 and α exp + α exp − = 0 ⇒ V (x, s, u) = 0, 1 2 1 u 2 u h NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 15 since α1 = α2 = 0. Using the method of undetermined coeﬃcients on the nonhomogeneous part, we get 3u V (x, s, u) = sin(2πx), (4.28) p s + 4π2u 3u Since, β1 = s+4π2u , and β2 = 0. Then Equ. (4.25) will become 3u V (x, s, u) = sin(2πx), (4.29) s + 4π2u Taking the inverse Shehu transform of Equ. (4.29), we get v(x, t) = 3 exp(−4π2t) sin(2πx). (4.30)

Figure 1. 3D and 2D surfaces of the analytical solution of Equ. (4.21) in the ranges −1 < x < 1, and −1 < t < 1, when t = 1. Example 6. Consider the following nonhomogeneous partial diﬀerential equation ∂2v(x, t) ∂2v(x, t) = β2 + sin(πx) (4.31) ∂t2 ∂x2 subject to the boundary and initial conditions ∂v(x, 0) v(0, t) = 0, v(1, t) = 0, v(x, 0) = 0, = 0, β2 = 1. (4.32) ∂t Applying the Shehu transform on both sides of Equ. (4.31), we get s2 s d2V (x, s, u) u V (x, s, u) − v(x, 0) − v0(x, 0) = + sin(πx). (4.33) u2 u dx2 s Substituting the given initial condition and simplifying, we get d2V (x, s, u) s2 u − V (x, s, u) = − sin(πx). (4.34) dx2 u2 s The general solution of Equ. (4.34) can be written as

V (x, s, u) = Vh(x, s, u) + Vp(x, s, u), (4.35) where Vh(x, s, u) is the solution of the homogeneous part which is given by s s V (x, s, u) = λ exp x + λ exp − x , (4.36) h 1 u 2 u and Vp(x, s, u) is the solution of the nonhomogeneous part which is given by

Vp(x, s, u) = η1 sin(πx) + η2 cos(πx). (4.37) NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 16

Applying the boundary conditions on Equ. (4.36), we deduce s s λ + λ = 0 and λ exp + λ exp − = 0 ⇒ V (x, s, u) = 0, 1 2 1 u 2 u h since λ1 = λ2 = 0. Using the method of undetermined coeﬃcients on the nonhomogeneous part, we get 1 u us V (x, s, u) = − sin(πx), (4.38) p π2 s s2 + u2π2 since, u3 1 u us η = = − , and η = 0. 1 s(s2 + u2π2) π2 s s2 + u2π2 2 Then Equ. (4.35) will becomes 1 u us V (x, s, u) = − sin(πx). (4.39) π2 s s2 + u2π2 Taking the inverse Shehu transform of Equ. (4.39), we get 1 v(x, t) = (1 − cos(πt)) sin(πx). (4.40) π2

Figure 2. 3D and 2D surfaces of the analytical solution of Equ. (4.31) in the ranges −1 < x < 1, and −1 < t < 1.

5. Conclusion We introduced an efficient Laplace-type integral transform called the Shehu transform for solving both ordinary and partial differential equations. We presented its existence and inverse transform. We presented some useful properties and their applications for solving ordinary and partial differential equations. We provide a comprehensive list of the Laplace transform, Sumudu transform, and the natural transform in table 1 to show their mutual relationship with the Shehu transform. Finally, based on the mathematical formulations, simplicity and the findings of the proposed integral transform, we conclude that it is highly efficient because of the following advantages: • It is a generalization of the Laplace and the Sumudu integral transforms. • Its visualization is easier than the Sumudu transform, the natural transform, and the Elzaki transform. NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 17

• The Laplace-type integral transform become Laplace transform when the variable u = 1 and the Yang integral transform when the variable s = 1. • It can easily be applied directly to some class of ordinary and the partial diﬀerential equations as demonstrated in the application section. • For advanced research in physical science and engineering, the proposed integral transform can be considered a stepping-stone to the Sumudu transform, the natural transform, the Elzaki transform, and the Laplace transform.

6. Acknowledgements The authors are highly grateful to the editor’s and the anonymous referees’ for their useful comments and suggestions in this paper. This research is partially supported by the National Natural Science Foundations of China under Grants No. (11571206). The ﬁrst author also acknowledges the ﬁnancial support of China Scholarship Council (CSC) in Shandong University with grand (CSC No: 2017GXZ025381).

References

[1] H.A. Agwa, F.M. Ali, A. Kilicman, A new integral transform on time scales and its applications, Advances in Difference Equations, 60(2012), 1–14. [2] C. Ahrendt, The Laplace transform on time scales, Pan. Am. Math. J., 19(2009), 1–36. [3] A. Atangana, A note on the triple Laplace transform and its applications to some kind of third-order differential equation, Abstract and Applied Analysis, 2013(2013), Article ID 769102, 1–10. [4] H.M. Srivastava, A.K. Golmankhaneh, D. Baleanu, X.Y. Yang, Local fractional Sumudu transform with applications to IVPs on Cantor sets, Abstr Appl Anal, 2014(2014), Article ID 620529, 1–7. [5] G. Dattoli, M. R. Martinelli, P. E. Ricci, On new families of integral transforms for the solution of partial differential equations, Integral Transforms and Special Functions, 8(2005), 661–667. [6] H. Bulut, H.M. Baskonus, and F.B.M. Belgacem, The analytical solution of some fractional ordinary differential equations by the Sumudu transform method, Abstract and Applied Anal- ysis, 2013(2013), Article ID 203875, 1–6. [7] S.Weerakoon, The Sumudu transform and the Laplace transform: reply, International Jour- nal of Mathematical Education in Science and Technology. 28(1997), 159–160. [8] D. Albayrak, S.D. Purohit, and F. U¸car,Certain inversion and representation formulas for q-sumudu transforms, Hacettepe Journal of Mathematics and Statistics, 43(2014), 699–713. [9] S. Weerakoon, Application of Sumudu transform to partial differential equations, Interna- tional Journal of Mathematical Education in Science and Technology, 25(1994), 277–283. [10] X.Y. Yang, Y. yang, C. Cattani, and M. Zhu, A new technique for solving 1-D Burgers equation, Thermal Science, 21(2017), S129–S136. [11] A. Kilicman, H. Eltayeb, On new integral transform and differential equations, J. Math. Probl. Eng, 2010(2010), Article ID: 463579, 1–13. NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 18

[12] S. Bochner, K. Chandrasekharan, Fourier transforms, Princeton University Press, Princeton, NJ, USA, (1949). [13] R. N. Bracewell, The Fourier transform and its applications, McGraw-Hill, Boston, Mass, USA, 3rd edition, (2000). [14] R. Murray, Spiegel. Theory and problems of Laplace transform. New York, USA: Schaum’s Outline Series, McGraw–Hill, (1965). [15] L. Debnath, D. Bhatta.. Integral transform and their applications. CRC Press, New York, NY, USA (2010). [16] B. Davies, Integral transforms and their applications, Texts in Applied Mathematics, Springer, New York, NY, USA, (2002). [17] E.I. Jury, Theory and applications of the z-transform Method, John Wiley and Sons, New York, NY,USA, (1964). [18] K. Liu, R.J. Hu, C. Cattani, G.N. Xie, X.J. Yang, and Y. Zhao, Local fractional z-transforms with applications to signals on Cantor sets, Abstract and Applied Analysis, 2014(2014),Ar- ticle ID: 638648, 1–6. [19] P.M. Morse, H. Feshbach, Methods of theoretical physics, McGraw-Hill, New York, (1953), 484–485. [20] P. Flajolet, X. Gourdon, and P. Dumas, Mellin transforms and asymptotics: harmonic sums, Theoretical Computer Science, 144(1995), 3–58. [21] C. Donolato, Analytical and numerical inversion of the Laplace-Carson transform by a differential method, Computer Physics Communications, 145(2002), 298–309. [22] A.M. Makarov, Application of the Laplace-Carson method of integral transformation to the theory of unsteady visco-plastic flows, J. Engrg. Phys. Thermophys 19(1970), 94–99. [23] E. Sjntoft, A straightforward deconvolution method for use in small computers, Nucl. In- strum. Methods, 163(1979), 519–522. [24] A.S. Vasudeva Murthy, A note on the differential inversion method of Hohlfield et al., SIAM J. Appl. Math., 55(1995), 712–722. [25] I.N. Sneddon, The Use of integral transform, McGraw-Hill, New York, (1972). [26] K. Xie, Y. Wanga, K. Wang, and X. Cai, Application of Hankel transforms to boundary value problems of water flow due to a circular source, Applied Mathematics and Computation, 216(2010), 1469–1477. [27] Watugala GK. Sumudu transform–a new integral transform to solve differential equations and control engineering problems. Math. Engg. in Indust., 6(1998), 319–329. [28] M.A. Asiru, Sumudu transform and solution of integral equations of convolution type, Int. J. of Math. Edu. Sci. and Tech., 33(2002), 944–949. [29] F.B.M. Belgacem, S.L. Kalla, A.A. Karaballi, Analytical investigations of the Sumudu transform and applications to integral production equations, Math. Probl. in Engg., 3(2003), 103-118. [30] F.B.M. Belgacem, A.A. Karaballi, Sumudu transform fundamental properties, investigations and applications. J. of Appl. Math. and Stoch. Anal., 2006(2006), Article ID 91083, 1–23. [31] H. ELtayeh and A. kilicman, On Some applications of a new integral transform, Int. Journal of Math. Analysis, 4(2010), 123–132. NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 19

[32] T.M. Elzaki. The new integral transform ”Elzaki transform”. Glob. J. of Pur. and Appl. Math., 7(2011), 57–64. [33] Z.H. Khan, W.A. Khan, N-transform-properties and applications. NUST J. of Engg. Sci., 1(2008), 127–133. [34] F.B.M. Belgacem, R. Silambarasan, Theory of natural transform. Math. in Engg. Sci., and Aeros., 3(2012), 99–124. [35] F.B.M. Belgacem, R. Silambarasan, Advances in the natural transform. AIP Conference Proceedings; 1493 January 2012; USA: American Institute of Physics. (2012), 106–110. [36] H.M. Srivastava, Minjie Luo, R.K.Raina, A new integral transform and its applications, Acta Mathematica Scientia, 35(2015), 1386–1400. [37] A. Atangana, A. Kilicman, A novel integral operator transform and its application to some FODE and FPDE with some kind of singularities, Mathematical Problems in Engineering, 2013(2013), Article ID: 531984, 1–7. [38] X.J. Yang, A new integral transform method for solving steady heat-transfer problem, Ther- mal Science, 20(2016), S639–S642. [39] X.J. Yang, A new integral transform operator for solving the heat-diffusion problem, Applied Mathematics Letters, 64(2017), 193–197. [40] Y.X. Jun, F. Gao , A new technology for solving diffusion and heat equations, Thermal Science, 21(2017), 133–140. [41] A. Atangana, B.S.T. Alkaltani, A novel double integral transform and its applications, Jour- nal of Nonlinear Science and Applications, 9(2016), 424–434. [42] H. Eltayeb, A note on double Laplace decomposition method and nonlinear partial differential equations, New Trends in Mathematical Sciences, 5(2017), 156–164. [43] F.B.M. Belgacem, R. Silambarasan, H. Zakia, T. Mekkaoui, New and extended applications of the natural and Sumudu transforms: Fractional diffusion and Stokes fluid flow realms. Advances Real and Complex Analysis with Applications, Publisher: Birkhuser, Singapore, (2017). 107–120.

Appendix Table 1: Here we present a comprehensive list of the Shehu transform of some special functions and their relationship with the natural transform N [v(t)], the Sumudu transform S [v(t)], and the Laplace transform.

S.NO. v(t) S [v(t)] N [v(t)] S [v(t)] £ [v(t)]

u 1 1 1 1 s s 1 s

u2 u 1 2 t s2 s2 u s2

u 1 1 1 3 exp(α(t)) s−αu s−αu 1−αu s−α

sin(αt) u2 u u 1 4 α s2+α2u2 s2+α2u2 1+α2u2 s2+α2

us s 1 s 5 cos(αt) s2+α2u2 s2+α2u2 1+α2u2 s2+α2 NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 20

us s 1 s 6 cosh αt s2−u2 s2−u2 1−u2 s2−1

tn u n+1 un n 1 7 n! n = 0, 1, 2... s sn+1 u sn+1

tn u n+1 un n 1 8 Γ(n+1) n = 0, 1, 2, ... s sn+1 u sn+1

us s 1 1 9 cos(t) s2+u2 s2+u2 1+u2 s2+1

u2 u u 1 10 sin(t) s2+u2 s2+u2 1+u2 s2+1

sinh(αt) u2 αu αu2 α 11 α s2−α2u2 s2−α2u2 1−α2u2 s2−α2

us s 1 s 12 cosh(αt) s2−α2u2 s2−α2u2 1−α2u2 s2−α2

u(s−βu) s−βu 1−βu s−β 13 exp(βt)cosh(αt) (s−βu)2−α2u2 (s−βu)2−α2u2 (s−βu)2−α2u2 (s−βu)2−α2

exp(βt) sinh(αt) u2 u u 1 14 α (s−βu)2−α2u2 (s−βu)2−α2u2 (1−βu)2−α2u2 (s−β)2−α2

t sin(αt) u3s u2s u3 s 15 2α (s2+α2u2)2 (s2+α2u2)2 (1+α2u2)2 (s2+α2)2

u2(s2−α2u2)2 u(s2−α2u2)2 u(1−α2u2)2 (s2−α2)2 16 t cos(αt) (s2+α2u2)2 (s2+α2u2)2 (1+α2u2)2 (s2+α2)2

sin(αt)+αt cos(αt) u2s2 us2 u s2 17 2α (s2+α2u2)2 (s2+α2u2)2 (1+α2u2)2 (s2+α2)2

αt sin(αt) us3 s3 1 s3 18 cos αt − 2 (s2+α2u2)2 (s2+α2u2)2 (1+α2u2)2 (s2+α2)2

sin(αt)−αt cos(αt) u4 u3 u3 1 19 2α3 (s2+α2u2)2 (s2+α2u2)2 (1+α2u2)2 (s2+α2)2

u2 u u2 1 20 t sinh(αt) + t cosh(αt) (s−αu)2 (s−αu)2 (1−αu)2 (s−α)2

t sinh(αt) u3s u2s u2 s 21 2α (s2−α2u2)2 (s2−α2u2)2 (1−α2u2)2 (s2−α2)2 √ u −1 αu 1 −1 αu −1 2 1 −1 α 22 Si(αt) (Sine integral) s tan s s tan s tan u α s tan s

u s2+α2 1 s2+α2u2 1 α2u2+1 1 s2+α2 23 Ci(αt) (Cosine integral) − 2s log α2 − 2s log α2u2 − 2 log α2u2 − 2s log α2

u αu−s 1 αu−s αu−1 1 α−s 24 Ei(αt) (Exp. integral) − s log αu − s log αu log αu − s log α

(3−α2t2) sin(αt)−3αt cos(αt) u6 u5 u5 1 25 8α5 (s2+α2u2)3 (s2+α2u2)3 (1+α2u2)3 (s2+α2)3

(3−α2t2) sin(αt)+5αt cos(αt) u2s4 us4 u s4 26 8α (s2+α2u2)3 (s2+α2u2)3 (1+α2u2)3 (s2+α2)3

(8−α2t2) cos(αt)−7αt sin(αt) us5 s5 1 s5 27 8 (s2+α2u2)3 (s2+α2u2)3 (1+α2u2)3 (s2+α2)3

t2 sin(αt) u4(3s2−α2u2) u3(3s2−α2u3) u3(−3+α2u2) (3s2−α2) 28 2α (s2+α2u2)3 (s2+α2u2)3 (1+α2u2)3 (s2+α2)3

t2 cos(αt) u3(s3−3α2u2s) u2(s3−3α2u2s) u2(1−3α2u2) (s3−3α2s) 29 2 (s2+α2u2)3 (s2+α2u2)3 (1+α2u2)3 (s2+α2)3

t3 sin(αt) su5(s−αu)2 su4(s−αu)2 u4(1−αu)2 s(s−α)2 30 24α (s2+α2u2)4 (s2+α2u2)4 (1+α2u2)4 (s2+α2)4

exp(αt)−exp(βt) u2 u u 1 31 α−β α 6= β (s−αu)(s−βu) (s−αu)(s−βu) (1−βu)(1−αu) (s−β)(s−α) NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 21

α exp(αt)−β exp(βt) us s 1 s 32 α−β α 6= β (s−βu)(s−αu) (s−βu)(s−αu) (1−βu)(1−αu) (s−β)(s−α)

u 1 1 1 33 I0(αt) √ √ √ √ s2−α2u2 s2−α2u2 1−α2u2 s2−α2

−αs 1 −αs 1 −α 34 δ(t − α) u exp u u exp u u exp u exp(−αs)

u 1 1 1 35 J0(αt) √ √ √ √ s2+α2u2 s2+α2u2 1+α2u2 s2+α2

Table 2: General properties of Shehu transform

S.NO. Deﬁnition/Property S [v(t)] = V (s, u) Transforms

R ∞ −st 1 Deﬁnition S [v(t)] 0 exp u v(t) dt; s > 0, u > 0

−1 1 R α+i∞ 1 st 2 Inverse v(t) = S [V (s, u)] 2πi α−i∞ u exp u V (s, u) ds

3 Linearity S [αv(t) + βw(t)] αS [v(t)] + βS [w(t)]

u s 4 Change of scale S [v(αt)] α V α , u

0 s 5 Derivatives S [v (t)] u V (s, u) − v(0)

00 s2 s 0 S [v (t)] u2 V (s, u) − u v(0) − v (0)

000 s3 s2 s 0 00 S [v (t)] u3 V (s, u) − u2 v(0) − u v (0) − v (0) . .

h (n) i sn Pn−1 s n−(k+1) S v (t) un V (s, u) − k=0 u NEW INTEGRAL TRANSFORM FOR SOLVING DIFFERENTIAL EQUATIONS 22

Table 3: Summary of some integral transform and their deﬁnitions

S.No. Integral transform Deﬁnition

R ∞ 1 Laplace transform £[f(t)] = F (s) = 0 exp (−st) f(t)dt

∞ 2 Fourier transform [f(t)] = f(ω) = √1 R exp (−iωt) f(t)dt z 2π −∞

∗ R ∞ s−1 3 Melling transform M[f(s); s] = f (s) = 0 x f(x)dx R ∞ 4 Hankel’s transform Fv(s) = Hv[f(r)] = 0 rf(r)Jv(sr)dr, r ≥ 0

1 R ∞ −t 5 Sumudu transform S[f(t)] = G(u) = u 0 exp u f(t)dt ˆ R ∞ 6 Laplace-Carson transform fC (p) = p 0 exp (−pt) f(t)dt, t ≥ 0

R ∞ n 7 Atangana-Kilicman transform Mn(s) = Mn[f(x)](s) = 0 x exp (−xs) f(x)dx

R ∞ −t 8 El-zaki transform E[f(t)] = T (u) = u 0 exp u f(t)dt

R ∞ −τ 9 Yang transform Y [φ(τ)] = φ(ω) = 0 exp ω φ(τ)dτ

+ 1 R ∞ −st 10 natural transform N [f(t)] = R(s, u) = u 0 exp u f(t)dt, s > 0, u > 0

P∞ −n 11 z-transform X(z) = Z{x[n]} = n=−∞ x[n]z , n ∈ Z, z ∈ C

R ∞ exp(−ut)f(vt) 12 M-transform Mρ,m[f(t)](u, v) = 0 (tm+vm)ρ dt, ρ ∈ C, <(ρ) ≥ 0, m ∈ Z+

School of Mathematics, Shandong University, Jinan, Shandong 250100, P.R. China

∗Corresponding author: [email protected]