Direct Methods For Singular Integral Equations and Non-Homogeneous Parabolic PDEs

A. Aghili


In this article, the author presented some applications of the Laplace, 2, and Post-Widder transforms for solving fractional singular integral equation, impulsive differential equation and systems of differential equations. Finally, analytic solution for a non-homogeneous partial differential equation with non-constant coefficients is given. The obtained results reveal that the integral transform method is an effective and convenient tool.

MSC. 44A10; 44A20; 45E10; 35C15; 35C05.

Keywords. Laplace transform; 2-transform; modified Bessel function; Post-Widder transform.

Author address: Department of Applied Mathematics, University of Guilan, Rasht, Iran, P.O. Box 1841, e-mail:,

1 Introduction

The use of integral transforms in applications is quite extensive. In applied probability, as inventory and risk theory, queueing theory, the Laplace transform is a valuable tool for finding underlying probability density functions. As the Laplace transforms, the 2-transform is used in a variety of applications, the most common usage of the 2-transform is in the solution of the singular integral equations and initial value problems. Our interest in this transform stems from the potential applications to boundary value problems. To the best of our knowledge, the properties of the 2-transform have not been studied in any detail.

For solving partial differential equations, two methods, have been more extensively used the Laplace type integral transformations on the one hand and separation of variables on the other hand. New methods have also been proposed, the first integral method, the (GG)-expansion method and many more. The main purpose of this work has been to employ the integral transform method for studying certain mathematical models. It is worth mentioning that the integral transform methods are mostly suitable for the solution of linear differential and partial differential equations. Finally, this article presents the exact solution of non-homogeneous partial differential equation with non-constant coefficients, which is solved by direct application of the 2- transform. The primary advantage of this approach is that it solves boundary value problems characterized by partial differential equation without classical differential equation theory. We confine ourselves here to a few non-trivial examples which illustrate the method and lead to some interesting new results needed in the paper.

Definition 1.1 The Laplace transform of function f(t) is given by [123]

{f(t);s} =0estf(t)dt = F(s). (1)

If {f(t);s} = F(s), then 1{F(s)} is as follows

f(t) = 1 2πici∞c+i∞estF(s)ds, (2)

where F(s) is analytic in the region Re (s) > c.

The above complex integral is known as Bromwich integral [4]. The existence of the Laplace transform will depend on the function f(t) and the parameter s.

Lemma 1.2 The following identities hold true.

  1. 1 (s+λ) = [ 1 πt λeλ2t Erfc (λt) ],
  2. 1 s(s+λ) = [eλ2t Erfc (λt) ],

Proof. See [13].  _

Singular integral equations arise frequently in the mathematical modeling of continuum phenomena, and many a time cannot be treated by known analytical techniques. Though certain problems had received the attention of aerodynamicists long ago, by contrast, extensive development of theory and methods for the approximate numerical solution is of recent vintage. Many physical problems dealing with radiative transfer, neutron transport, dispersal of aerosol like particles, fluid flow and waveguides can be reduced to singular integral equations. Comprehensive accounts of techniques for numerical solution of integral equations can be found in the monograph by Prossdorf and Silberman [5].

Theorem 1.3 Let us consider fractional singular integro-differential equation

D0,tc,αϕ(t) = f(t) + λt+ϕ(ξ),0 < t < 1

ϕ(0) = u0,0+ϕ(ξ) = k,0 < α < 1. Then the above fractional singular integro-differential equation has the following solution. ϕ(t) = u0 n=0n=+(λ)nt(α+1)n Γ(1+(α+1)n) + n=0n=+(λ)n0tf(t η) η(α+1)n Γ(1+(α+1)n) λk n=0n=+(λ)nt(α+1)(1+n)1 Γ(1+(α+1)n) .

Note 1 This kind of singular integral equation is not considered in the literature.

Solution 1 Taking the Laplace transform of the above fractional singular integral equation term wise, leads to

sαΦ(s) sα1u 0 = F(s) + λΦ(s)Φ(0) s = F(s) + λΦ(s)k s .

After solving the transformed equation, we obtain

Φ(s) = sF(s) λ+sα+1 + u0sαλk λ+sα+1 ,


Φ(s) = n=0+(λ)n [ F(s) sn(α+1)+α + u0 s(α+1)n+1 λk s(α+1)(n+1) ] .

Taking the inverse Laplace transform term-wise, we get

ϕ(t) = n=0+(λ)n [0tf(t ξ) ξn(α+1)+α1 Γ(n(α+1)+α) + u0t(α+1)n Γ(n(α+1)+1) λkt(α+1)(n+1)1 Γ((α+1)(n+1)) ],

0 < t < 1. It is easy to verify that ϕ(0) = u0.

In applied mathematics, engineering and mathematical physics, Bessel functions are associated most commonly with the partial differential equations of the wave or diffusion equations in cylindrical or spherical coordinates. No other special functions have received such a detailed treatment as have the Bessel functions [6].

Lemma 1.4 By using an appropriate integral representation for the modified Bessel functions of the second kind of order ν, Kν(s), we have the following

1 {Kν(s) sν } = π Γ (ν+1 2 )2ν(t2 1)ν1 2 . (3)

Proof. In view of the Definition 1.1 , taking the inverse Laplace transform of the given Kν(s) sν , we obtain

h(t) = 1 2πiai∞a+i∞est (Kν(s) sν )ds, (4)

by using the following integral representation for Kν(s)

Kν(s) sν = π Γ(ν+1 2)2ν0es cosh t sinh 2νtdt. (5)

By inserting relation (5 ) in (4 ), we get

h(t) = 1 2πiai∞a+i∞est ( π Γ(ν+1 2)2ν0es cosh r sinh 2νrdr )ds, (6)

and changing the order of integration in relation (6 ) leads to

h(t) = π Γ(ν+1 2)2ν0sinh 2νr ( 1 2πiai∞a+i∞es(tcosh r)ds )dr. (7)

The inner integral is δ(t cosh r), therefore

h(t) = π Γ(ν+1 2)2ν0δ(t cosh r)sinh 2νrdr. (8)

Making the change of variable t cosh r = u, and considerable algebra and elimination process, we obtain

h(t) = π Γ(ν+1 2)2νt1δ(u)((tu)21)ν (tu)2 1 du = π Γ(ν+1 2)2ν(t2 1)ν1 2 . (9)

For the special case ν = 0, we get the following relation

1{K 0(s)} = (t2 1)1 2 . (10)


Example 1.5 Consider the following generalized Abel singular integral equation of the second kind.

ϕ(t) = eλ2 t tν+βk + 1 Γ(α)0t ϕ(ξ) (tξ)1α,ϕ(0) = 0,0 < α,ν < 1. (11)

This type of integral equation arises in the theory of wave propagation over a flat surface. Such integral equations occur rather frequently in mathematical physics and possess very interesting properties.

Note 2 The above mentioned singular integral equation can be written in terms of the Riemann-Liouville fractional integral as below

ϕ(t) = eλ2 t tν+βk + I0,tRL,αϕ(t),0 < α,ν < 1. (12)

Solution 2 By taking the Laplace transform of the given integral equation, after simplifying we arrive at

Φ(s) = ( eλ2 t tν+βk ) + Φ(s) sα . (13)

Solving transformed equation leads to

Φ(s) = ( eλ2 t tν+βk )(1 1 sα)1 = ( eλ2 t tν+βk ) + ( eλ2 t tν+βk )( 1 sα1). (14)

Upon taking the inverse Laplace transform, we get

ϕ(t) = eλ2 t tν+βk +0t e λ2 tη (tη)ν+βk1( 1 sα1). (15)

In the special case α = 0.5, we have the following relation

ϕ(t) = eλ2 t tν+βk +0t e λ2 tη (tη)ν+βk ( 1 πη + eη Erfc (η) ). (16)

Note 3 We may check that ϕ(0) = 0.

In recent years, fractional calculus appeared as an important tool to deal with anomalous diffusion processes. A more physical approach of anomalous diffusion processes has several applications in many fields such as diffusion in porous media or long range correlation of DNA sequence [7]. The closed form solution of the time fractional impulsive heat equation has been presented. At this stage we use the joint transform method to obtain a solution of a time fractional impulsive heat equation. The joint transform method provides an effective procedure for exact solution of a wide class of systems representing real physical problems.

Problem 1.6 Let us consider the following non-homogeneous time fractional impulsive heat equation,

D0,tc,0.5u = λu xx + ηδ(t a)δ(x b)t > 0, < x < +, (17)

with initial and boundary conditions

u(x,0) = 2πδ(x),lim |x|+|u(x,t)| < M0,

where a,b,η are constants and λ,M0 are positive constants.

Note 4 Fractional derivative is in the Caputo sense.

Solution 3 The joint Laplace-Fourier transform of function u(x,t) is defined as follows

[ [u(x,t);x w ];t s ] =0+est [ 1 2π+eiwxu(x,t)dx ]dt = U(w,s).

Let us take the joint Laplace-Fourier transform of the above equation (17 ) term wise and using boundary conditions, we get

sU(w,s) 1 s = λw2U(w,s) + ηebseiwa.

Solving the transformed equation leads to

U(w,s) = 1 s(s+λw2) + ηeiaw[ ebs s+λw2].

At this point taking the inverse joint Laplace-Fourier transforms and by virtue of the Lemma 1.2 we obtain

u(x,t) = 1 2π [+eixw[etλ2w4 Erfc (λtw2)]dw +
+η+ei(xa)w [ 1 π(tb) λw2eλ2(tb)w4 Erfc (λt bw2) ]dw ].

2 The 2-transform

In the literature, we have significant generalizations of the integral transforms and new uses of the transformation method in engineering, applied mathematics and physics applications. The 2-transformation was first introduced by Yurekli [5] and denoted as follows

2{f(t);s} =0tes2t2 f(t)dt. (18)

In the absence of methods for the inversion of the 2-transform, recently the authors [8], established a simple formula to invert the 2-transform of a desired function. We present certain new inversion techniques for the 2-transform and an application of generalized product theorem for solving singular integral equations and boundary value problems are given.

Lemma 2.1 The following identity holds true

2{tηδ(atm λ);s} = 1 am(λ a)ηm+2 m es2(λ a) 2 m ,λ,η > 0,k,m > 1.

Solution 4 By definition of the 2-transform, we have

2{tηδ(atm λ);s} =0+tη+1es2t2 δ(atm λ)dt.

Making a change of variable atm λ = ξ, then we have

2 {tηδ(atm λ);s } =λ+(ξ+λ a )ηm+2 m es2(ξ+λ a ) 2 m δ(ξ) am,

using elementary property of Dirac-delta function, yields

2 {tηδ(atm λ);s } = 1 am(λ a)ηm+2 m es2(λ a) 2 m .

Consider the special case η = 0,m = 1,a = 1, we get

2 [δ(t λ);s ] = λeλ2s2 .

3 Elementary Properties of the 2-Transform

Here, we will derive a relation between the 2-transform of the derivative of the function and the 2-transform of the function itself. We recall a useful lemma about the 2 -transform of the δ-derivatives.

Lemma 3.1 If f,f,,f(n1) are all continuous and of exponential order exp (c2t2) as t for some real constant c and piecewise continuous derivative f(n) on the interval t 0

  1. For n = 1,2, then

    2 {δtnf(t);s } = 2ns2n 2{f(t);s} 2n1s2(n1)f(0+) (19) 2n2s2(n2)(δ tf)(0+) (δ tn1f)(0+).
  2. For n = 1,2,

    2 {t2nf(t);s } = (1)n 2n δsn 2 {f(t);s }, (20)

where the differential operators δt,δt2, are defined as below

δt = 1 t d dt,δt2 = δ tδt = 1 t2 d2 dt2 1 t3 d dt.

Proof. See [9].  _

4 Inversion Formula for the 2-Transform and Efros’ Theorem

Lemma 4.1 Let us assume that F(s) is analytic function (s = 0 is not a branch point) except at finite number of poles each of which lies to the left hand side of the vertical line Re (s) = c and if F(s) 0 as s through the left plane Re (s) c, and

2{f(t);s} =0texp (s2t2)f(t)dt = F(s), (21)


21{F(s)} = 1 2πici∞c+i∞2F(s)est2 ds = k=1m [Res{2F(s)est2 },s = sk ]. (22) 

Proof. See [8].  _

Example 4.2 Solving the following Cauchy’s problem attached to a second order impulsive differential equation with non-constant coefficients

1 t2y′′ t 1 t3y + 4λ2y(t) = tβδ(t ξ),y(0+) = (δ ty)(0+) = 0. (23)

Solution 5 By taking the 2-transform of the above equation term wise, we get

2 (δt2y(t) ) + 4λ2 2(y(t)) = 2(tβδ(t ξ)). (24)

Let us assume that 2(y(t)) = Y (s), then after evaluation of the 2-transform each term, we arrive at

4s4Y (s) + 4λ2Y (s) = ξβ+1eξ2s2 , (25)


Y (s) = ξβ+1eξ2s2 4s4+4λ2 , (26)

using inversion formula for the 2-transform, we have

y(t) = 1 2πici∞c+i∞ (ξβ+1eξ2s 2(s2+λ2) )est2 ds, (27)

finally, using the second part of the Lemma 3.1 leads to the following solution:

y(t) = ξβ+1 1 2λ sin (λ(t2 ξ2)). (28)

Lemma 4.3 (Efros’ Theorem for 2-Transforms)

Let 2(f(t)) = F(s) and assuming Φ(s), q(s) are analytic and such that, 2(Φ(t,τ)) = Φ(s)τeτ2q2(s) , then the following relation holds true,

2 {0f(τ)Φ(t,τ)dτ } = F(q(s))Φ(s). (29)

Proof. By definition of the 2-transform

2 {0f(τ)ϕ(t,τ) } =0tes2t2 (0f(τ)ϕ(t,τ) )dt, (30)

and changing the order of integration we arrive at

0f(τ) (0tes2t2 ϕ(t,τ)dt ) = Φ(s)0f(τ)τeτ2q2(s) = Φ(s)F(q(s)). (31)


In the sequel we will show that the 2-transform is suitable for solving singular integral equation with trigonometric kernel.

Example 4.4 By means of the above Lemma 4.3 , we may solve the singular integral equation with trigonometric kernel

0f(ξ)sin () = H(t λ), (32)


f(t) = 2 π cos (λt) t . (33)

Solution 6 Applying the 2-transform followed by the generalized product theorem and using the fact that

2{sin (τt)} = π 4s3τeτ2 4s2 , (34)


2 [H(t λ,t s ] = 1 2s2eλ2s2 , (35)


F( 1 2s)π 4s3 = 1 2s2eλ2s2 , (36)


F(s) = 1 πseλ2 4s2 , (37)

using inversion formula for the 2-transform leads to

f(t) = 21 [ 1 πseλ2 4s2 ;s t ] = 2 π cos (λt) t . (38)

5 Solution to the System of Second Order Differential Equations With Non Constant Coefficients Via The 2-Transform

Theorem 5.1 We may consider the system of non-homogeneous second order differential equations in general form as follows

1 t2X (t) 1 t3X (t) = AX (t) + B (39)

where A and B are coefficient and constants matrices of type (n × n) and (n × 1) respectively and X (t), X (t) X (t) are column vectors. Thus, the above system of equations has the following solution.

X (t) = 21 {(4s4I A)1(2s2X (0) + δ tX (0) + 1 2s2B)}.

Proof. For the solution of the above system, first, we take 2-transform of the above system, we get,

2 [ 1 t2X 1 t3X (t) ] = A 2 [X (t) ] + B2[1]


4s4 2 [X (t) ] 2s2X (0) δ tX (0) = A2 [X (t) ] + 1 2s2B

after simplifying the above relation, we obtain

(4s4I A) 2 [X (t) ] = 2s2X (0) + δ tX (0) + 1 2s2B


2 [X (t) ] = (4s4I A)1(2s2X (0) + δ tX (0) + 1 2s2B)


X (t) = 21 {(4s4I A)1(2s2X (0) + δ tX (0) + 1 2s2B)}.


In the sequel, we give certain illustrative examples and lemmas related to the 2, Post-Widder transforms, and inversion formula for the Post-Widder transform.

6 Illustrative Lemmas and Examples

Lemma 6.1 By using an integral representation for the modified Bessel functions of the second kind of order ν, Kν(s), the following identity holds true

21 {Kν(s2) s2ν } = π Γ(ν+1 2)2ν(t4 1)ν1 2 . (40)

Solution 7 By applying Lemma 4.1 and taking the inverse 2-transform of the given Kν(s2) s2ν , we arrive at

g(t) = 1 2πiai∞a+i∞est2 (Kν(s) sν )ds, (41)

at this point, we use an integral representation for Kν(s)

Kν(s) sν = π Γ(ν+1 2)2ν0es cosh t sinh 2νtdt. (42)

By inserting relation (3 ) in (2 ), we get

g(t) = 1 2πiai∞a+i∞est2 ( π Γ(ν+1 2)2ν0es cosh r sinh 2νrdr )ds, (43)

in relation (4 ), we may change the order of integration to obtain

g(t) = π Γ(ν+1 2)2ν0sinh 2νr ( 1 2πiai∞a+i∞es(t2cosh r)ds )dr, (44)

the inner integral is δ(t2 cosh r), therefore

g(t) = π Γ(ν+1 2)2ν0δ(t2 cosh r)sinh 2νrdr, (45)

let us introduce a change of variable t2 cosh r = u, and considerable algebra and elimination process, we obtain

g(t) = π Γ(ν+1 2)2νt21δ(u)((t2u)21)ν (t2 u)2 1 du = π Γ(ν+1 2)2ν(t4 1)ν1 2 . (46)

In view of the definition of the 2-transform we have

π Γ(ν+1 2)2ν0+es2t2 t(t4 1)ν1 2 dt = Kν(s2) s2ν .

Let us consider the special case ν = 0, we get the following relation

21{K 0(s2)} = (t4 1)1 2 . (47)

Lemma 6.2 The following singular integral equation of Post-Widder type

0+ (u) u2+s2du = c s2a2k, (48)

has solution as below

ϕ(u) = 4c sin (π k ) πu2+a2k. (49)

Proof. Using the inverse Post-Widder transform (second iteration of the 2-transform) we have

ϕ(u) = 1 2πici∞c+i∞esu2 ( 1 2πici∞c+i∞ eps2 pa2kdp)ssds, (50)

introducing the new variable w = p a2 leads to

ϕ(u) = 1 2πici∞c+i∞2esu2 ( 1 2πiδi∞δ+i∞2ews2+a2s2 wk dw)ssds = (51)

= 1 2πici∞c+i∞4ces(a2+u2) (s1 k1 Γ(1 k) ) ds,


ϕ(u) = 4c Γ(1 k) 1 2πici∞c+i∞e(a2+u2)s 1 s11k ds, (52)

therefore, the final solution is as below

ϕ(u) = 4c Γ(1 k) 1 Γ(11 k). 1 u2+a2k = 4c sin (π k ) πu2+a2k. (53)


Lemma 6.3 Let us show the following Post-Widder type singular integral equation

0+ (u) u2+s2du = sn s2kλ2, (54)

has a solution as below

ϕ(u) = 2 π [u 2 n sin (πn) λu4k ]. (55)

Proof. Let us consider the following inversion formula for the Post-Widder transform [1]

𝒫1{F(s)} = 1 πi {F(u2e) F(u2e) }. (56)

In view of the above inversion formula for the Post-Widder transform, we have that

𝒫1 { 2sn s2kλ } = 1 πi { u2en (u2e)2kλ u2en (u2e)2kλ } = 2 π [u 2 n sin (πn) λu4k ].

Since the obtained solution satisfies integral equation, we get the following interesting integral identity,

4 sin (π n) π 0+ u1+ 2n (u2+s2)(λu4k)du = sn s2kλ2. (57)

Let us take s = 1, then after simplifying we get the following integral

0+ u1+ 2n (u2+1)(λu4k)du = π 4 sin (π n). (58)


Example 6.4 Let us consider the following singular integral equation with trigonometric kernel

0+(x)cos ξxdx = eξ Erfc (ξ). (59)

Solution 8 Taking the Laplace transform of both sides of the integral equation with respect to variable ξ, we arrive at

{0+(x)cos ξxdx } = 1 (s1)s, (60)

or, equivalently

0+ x x2+s2 ϕ(x)dx = 1 s(s1)s, (61)

the left hand side of the above relation can be written as Widder potential transform of ϕ(x). We have

𝒫{ϕ(x);s} = 1 s(s1)s, (62)


ϕ(x) = 1 πi ( 1 x2e(x2e1)x2 e 1 x2e(x2e1)x2 e ), (63)

after simplifying, we get

ϕ(x) = 2 πx3(x2+1), (64)

we deduce that

0+ 2 cos ξx πx2(x2+1)dx = eξ Erfc (ξ). (65)

Finally, from the above integral and inversion formula for the Fourier-cosine transforms we arrive at

0+eξ cos Erfc (ξ) = 1 x2(x2+1). (66)

Through an application of Leibnitz’s rule, by differentiating the above integral with respect to x under the integral sign and after simplifying we arrive at

0+ξeξ sin Erfc (ξ) = 4x2+2 x4(x2+1)2. (67)

At this point upon using inverse Fourier-sine transform we have

2 π0+sin (ξx) 4x2+2 x4(x2+1)2dx = ξeξ Erfc (ξ) (68)

In special case if we choose x = 1 in the above integral (68 ), we have the following

0+ξeξ sin ξErfc (ξ) = 3 2. (69)

Note 5 In the next section we will briefly illustrate the method of residues as it applies to both the 2-transform and its inversion formula.

7 Main Results. The 2-Transform For Non-Homogeneous Parabolic PDEs.

The second order PDEs, with non-constant coefficients have a number of applications in electrical and mechanical engineering, medical sciences and economics. The heat equation plays an important role in a number of fields of science. In this section, we will study the application of such PDEs. Note that separation of variables will only work if both the partial differential equation and the boundary conditions are linear and homogeneous.

Problem 7.1 Let us consider the following non-homogeneous parabolic PDE

1 2tut = urr + 1 rur λ2u k0 < r < 1,t > 0, (70)

with initial and boundary conditions

u(1,t) = u(r,0) = 0,|u(r,t)| < M0.

where k > 0,M0 > 0.

Solution 9 By taking the 2-transform of the above equation, we have

(s2 + λ2)U(r,s) 0.5u(r,0) U rr(r,s) 1 rUr(r,s) = k 2s2, (71)


Urr + 1 rUr (s2 + λ2)U = k 2s2U(1,s) = 0,|U(r,s)| < M. (72)

The general solution of the transformed equation is given in terms of the modified Bessel functions of order zero as follows

U(r,s) = c1I0(rs2 + λ2) + c2K0(rs2 + λ2) k 2s2(s2+λ2), (73)

since K0(sr) is unbounded as r 0, we have to choose c2 = 0, thus

U(r,s) = c1I0(rs2 + λ2) k 2s2(s2+λ2), (74)

from U(1,s) = 0, we find c1 = k 2s2(s2+λ2)I0(s2 +λ2), therefore

U(r,s) = k 2s2(s2+λ2) + kI0(rs2 +λ2) 2s2(s2+λ2)I0(s2 +λ2). (75)

Using complex inversion formula for the 2-transform, we have

u(r,t) = k λ2(1 eλ2t2 ) + k 2πici∞c+ est2I0(s+λ2r) s(s+λ2)I0(s+λ2)ds, (76)

the integrand in the above integral has simple poles at at s + λ2 = 0,s = 0 and s + λ2 = αn2,n = 1,2,3,. hence, the residue of integrand at s = λ2 is

lim sλ2(s + λ2) est2I0(s+λ2r) s(s+λ2)I0(s+λ2) = eλ2t2 λ2 , (77)

and the residue of integrand at s = 0 is

lim s0(s) I0(s+λ2r) s(s+λ2)I0(s+λ2) = I0(λr) λ2I0(λ), (78)

the residue of integrand at s = λ2 αn2 is

lim sλ2αn2(s + αn2 + λ2) est2I0(s+λ2r) s(s+λ2)I0(s+λ2), (79)

the integrand in the above integral has simple poles at s + λ2 = αn2,n = 1,2,3, and also at s + λ2 = 0 where ηn are simple zeros of the modified Bessel function of order zero I0 , as s + λ2 = η1,η2,,ηn,, hence, we have

lim sαn2λ2(s + λ2 + α n2) est2I0(s+λ2r) (s+λ2)I0(s+λ2) = = (lim sαn2λ2 s+λ2+αn2 I0(s+λ2) ) (lim sαn2λ2est2 I0(s+λ2r) s+λ2 ) = (lim sαn2λ2 1 I0(s+λ2) 1 2s+λ2 ) (eαn2t2 I0(iαnr) iαn2 ) = 2eαn2t2 I0(iαnr) iαnI1(iαn)


u(r,t) = 2k λ2 + 2kI0(λr) λ2I0(λ) 2 n=1n=eαn2t2I0(iαnr) iαnI1(iαn) . (80)

In view of the properties of the Bessel function, we have

I0(t) = J0(it) = J0(it),I1(t) = iJ1(it),J1(t) = J1(t),

therefore, we obtain the exact solution as below

u(r,t) = 2kI0(λr) λ2I0(λ) 2k λ2 + 2 n=1+eαn2t2J0(αnr) αnJ1(αn) . (81)

Note 6 We can check that u(1,t) = 0 and |u(r,t)| < M0.

8 Conclusion

The main goal of the present paper is to extend the application of the 2-transform to derive an analytic solution of boundary value problems. We have presented certain methods of solution for singular integral equations and boundary value problem using the Laplace and 2-transform. Certain non-trivial examples are also provided. The formulations presented in this article are simple and can be extended to other problems in the field of integral transforms.

Acknowledgements 1 The author would like to express his sincere thanks to the referee for helpful comments and suggestions that lead to a significant improvement of the paper.


[1]    A. Aghili, Solution to time fractional non-homogeneous first order PDE with non-constant coefficients, Tbilisi Math. J., 12 (2019) no. 4, pp. 149–155.

[2]    A. Aghili, Complete solution for the time fractional diffusion problem with mixed boundary conditions by operational method, Appl. Math. Nonlin. Sci., 6 (2020) no. 1, pp. 1–12.

[3]    A. Apelblat, Laplace transforms and their applications, Nova Science Publishers, Inc., New York, 2012.

[4]    H.J. Glaeske, A.P. Prudnikov, K.A. Skornik, Operational calculus and related topics, Chapman and Hall/CRC, 2006.

[5]    S. Prossdorf, B. Silberman, Numerical analysis for integral and related operator equations, Academie Verlag, Berlin, 1991.

[6]    A. Aghili, Special functions, integral transforms with applications, Tbilisi Math. J., 12 (2019) no. 1, pp. 33–44.

[7]    I. Podlubny, Fractional differential equations, Academic Press, New York, 1999.

[8]    A. Aghili, A. Ansari, Solving partial fractional differential equations using the A-transform, Asian-Eur. J. Math., 3 (2010) no. 2, pp. 209–220.

[9]    O. Yurekli, I. Sadek, A Parseval-Goldstein type theorem on the Widder potential transform and its applications, Int. J. Math. Math. Sci., 14 (1991), pp. 517–524.