Mathematical Theory and Modeling
ISSN 2224-5804 (Paper)
ISSN 2225-0522 (Online)
Vol.2, No.4, 2012
www.iiste.org
Solution of Linear and Nonlinear Partial Differential
Equations Using Mixture of Elzaki Transform and the
Projected Differential Transform Method
Tarig. M. Elzaki1* & Eman M. A. Hilal2
1. Mathematics Department, Faculty of Sciences and Arts-Alkamil, King Abdulaziz University,
Jeddah-Saudi Arabia.
Mathematics Department, Faculty of Sciences, Sudan University of Sciences and Technology-Sudan.
2. Mathematics Department, Faculty of Sciences for Girles King Abdulaziz University
Jeddah-Saudi Arabia
* E-mail of the corresponding author: Tarig.alzaki@gmail.com
Abstract
The aim of this study is to solve some linear and nonlinear partial differential equations using the new
integral transform "Elzaki transform" and projected differential transform method. The nonlinear terms can
be handled by using of projected differential transform method; this method is more efficient and easy to
handle such partial differential equations in comparison to other methods. The results show the efficiency
and validation of this method.
Keywords: Elzaki transform, projected differential transform method, nonlinear partial differential
equations.
1. Introduction
Many problems of physical interest are described by linear and nonlinear partial differential equations with
initial or boundary conditions, these problems are fundamental importance in science and technology
especially in engineering. Some valuable contributions have already been made to showing differential
equations such as Laplace transform method [lslam, Yasir Khan, Naeem Faraz and Francis Austin (2010),
Lokenath Debnath and Bhatta (2006)], Sumudu transform [A.Kilicman and H.E.Gadain (2009), (2010)],
differential transform method etc. Elzaki transform method is very effective tool for solve linear partial
differential equations [Tarig Elzaki & Salih Elzaki (2011)]. In this study we use the projected differential
transform method [Nuran Guzel and Muhammet Nurulay (2008), Shin- Hsiang Chang , Ling Chang (2008)
] to decompose the nonlinear terms, this means that we can use both Elzaki transform and projected
differential transform methods to solve many nonlinear partial differential equations.
1.1. Elzaki Transform:
The basic definitions of modified of Sumudu transform or Elzaki transform is defined as follows, Elzaki
transform of the function
f (t ) is
50
Mathematical Theory and Modeling
ISSN 2224-5804 (Paper)
ISSN 2225-0522 (Online)
Vol.2, No.4, 2012
∞
−
t
E [ f (t ) ] = v ∫ f (t ) e v dt ,
www.iiste.org
t >0
0
Tarig M. Elzaki and Sailh M. Elzaki in [2011], showed the modified of Sumudu transform [Kilicman &
ELtayeb (2010)] or Elzaki transform was applied to partial differential equations, ordinary differential
equations, system of ordinary and partial differential equations and integral equations.
In this paper, we combined Elzaki transform and projected differential transform methods to solve linear
and nonlinear partial differential equations.
To obtain Elzaki transform of partial derivative we use integration by parts, and then we have:
∂f (x , t ) 1
= T (x ,v ) −vf (x , 0)
E
∂t v
∂ 2 f (x , t ) 1
∂f (x , 0)
= 2 T (x ,v ) − f (x , 0) −v
E
2
∂t
∂t
v
∂f (x , t ) d
=
E
[T (x ,v )]
∂x dx
∂ 2f (x , t ) d 2
E
= 2 [T (x ,v )]
2
∂x
dx
Proof:
To obtain ELzaki transform of partial derivatives we use integration by parts as follows:
p
p
p −t
−t
−t
∞
∂f −vt
∂f
∂f
v
v
Ε ( x ,t ) = ∫ v
e dt = lim ∫ve
dt = lim ve f (x , t ) − ∫ e v f (x , t )dt
p →∞
p →∞
∂t
∂t
0 ∂t
0 0
0
=
T ( x ,v )
v
−vf ( x , 0 )
We assume that f is piecewise continuous and it is of exponential order.
Now
∂f
Ε
∂x
∞
∞
−t
−t
∂f ( x , t )
∂
v
v
=
=
ve
dt
ve
f ( x , t ) dt , using the Leibnitz rule to find:
∫
∫
∂
∂
x
x
0
0
∂f d
T ( x ,v )
Ε =
∂x dx
By the same method we find:
To find:
Let
∂ 2f d 2
Ε 2 = 2 T ( x ,v )
∂x dx
∂2 f
Ε 2 ( x, t )
∂t
∂f
= g , then we have:
∂t
∂g ( x , t )
∂ 2f
g ( x , t ) −vg x , 0
Ε 2 ( x , t ) = Ε
( )
=Ε
v
∂t
∂t
51
Mathematical Theory and Modeling
ISSN 2224-5804 (Paper)
ISSN 2225-0522 (Online)
Vol.2, No.4, 2012
www.iiste.org
∂ 2f
1
∂f
Ε 2 ( x , t ) = 2 T ( x , v ) − f ( x , 0 ) −v
( x , 0)
∂t
∂t
v
We can easily extend this result to the nth partial derivative by using mathematical induction.
1.2. Projected Differential Transform Methods:
In this section we introduce the projected differential transform method which is modified of the
differential transform method.
Definition:
The basic definition of projected differential transform method of function
f (x 1 , x 2 ,LL , x n ) is
defined as
1 ∂ k f (x 1 , x 2 ,LL , x n )
f (x 1 , x 2 ,LL , x n −1 , k ) =
∂x nk
k !
x n =0
Such that
(1)
f (x 1 , x 2 ,LL , x n ) is the original function and f (x 1 , x 2 ,LL , x n −1 , k ) is projected
transform function.
And the differential inverse transform of f
(x 1 , x 2 ,LL , x n −1 , k ) is defined as
∞
f (x 1 , x 2 ,LL , x n ) = ∑ f (x 1 , x 2 ,LL , x n )(x − x 0 ) k
k =o
The fundamental theorems of the projected differential transform are
52
(2)
Mathematical Theory and Modeling
ISSN 2224-5804 (Paper)
ISSN 2225-0522 (Online)
Vol.2, No.4, 2012
www.iiste.org
Theorems:
(1)
If z ( x1 , x2 ,......, xn ) = u ( x1 , x2 ,......, xn ) ± v ( x1 , x2 ,......, xn )
Then z ( x1 , x2 ,......, xn −1 , k ) = u ( x1 , x2 ,......, xn −1 , k ) ± v ( x1 , x2 ,......, xn −1 , k )
( 2)
Ifz ( x1 , x2 ,......, xn ) = c u ( x1 , x2 ,......, xn )
( 3)
Ifz ( x1 , x2 ,......, xn ) =
Then z ( x1 , x2 ,......, xn −1 , k ) = cu ( x1 , x2 ,......, xn −1 , k )
du ( x1 , x2 ,......, xn )
dxn
Then z ( x1 , x2 ,......, xn −1 , k ) = ( k + 1) u ( x1 , x2 ,....., xn −1 , k + 1)
( 4 ) Ifz ( x1 , x2 ,......, xn ) =
d n u ( x1 , x2 ,......, xn )
dxnn
Then z ( x1 , x2 ,......, xn−1 , k ) =
( 5) If
( k + n )! u
k!
( x1 , x2 ,....., xn−1 , k + n )
z ( x1 , x2 ,......, xn ) = u ( x1 , x2 ,......, xn ) v ( x1 , x2 ,......, xn )
k
Then z ( x1 , x2 ,......, xn −1 , k ) = ∑ u ( x1 , x2 ,......, xn −1 , m ) v ( x1 , x2 ,......, xn −1 , k − m )
m =0
( 6 ) If z ( x1 , x2 ,......, xn ) = u1 ( x1 , x2 ,......, xn ) u2 ( x1 , x2 ,......, xn ) ..... un ( x1 , x2 ,......, xn ) Then
z ( x1 , x2 ,......, xn −1 , k ) =
k
k n−1
∑ ∑
k3
......
kn −1 = 0 k n− 2 = 0
k2
∑ ∑ u ( x , x ,......, x
k 2 = 0 k1 = 0
1
1
2
n −1
, k1 ) u2 ( x1 , x2 ,......, xn −1 , k2 − k1 )
×.....un −1 ( x1 , x2 ,......, xn −1 , kn −1 − kn −2 ) un ( x1 , x2 ,......, xn −1 , k − kn −1 )
(7)
If z ( x1 , x2 ,......, xn ) = x1q1 x2q2 ........xnqn
1 k = qn
z ( x1 , x2 ,......, xn −1 , k ) = δ ( x1 , x2 ,......, xn−1 , qn − k ) =
0 k ≠ qn
Note that c is a constant and n is a nonnegative integer.
Then
2. Applications:
Consider a general nonlinear non-homogenous partial differential equation with initial conditions of the
form:
Du (x , t ) + Ru (x , t ) + Nu (x , t ) = g (x , t )
u (x , 0) = h (x ) ,
u t (x , 0) = f (x )
53
(3)
Mathematical Theory and Modeling
ISSN 2224-5804 (Paper)
ISSN 2225-0522 (Online)
Vol.2, No.4, 2012
www.iiste.org
Where D is linear differential operator of order two, R is linear differential operator of less
order than D , N is the general nonlinear differential operator and g ( x , t ) is the source term.
Taking Elzaki transform on both sides of equation (3), to get:
E [Du (x , t )] + E [Ru (x , t )] + E [Nu (x , t )] = E [ g (x , t )]
(4)
Using the differentiation property of Elzaki transforms and above initial conditions, we have:
E [u (x , t )] = v 2 E [ g (x , t )] + v 2 h (x ) + v 3f (x ) −v 2 E [Ru (x , t ) + Nu (x , t )]
(5)
Appling the inverse Elzaki transform on both sides of equation (5), to find:
u (x , t ) = G (x , t ) − E −1 {v 2 E [ Ru (x , t ) + Nu (x , t )]}
(6)
Where G ( x , t ) represents the term arising from the source term and the prescribed initial conditions.
Now, we apply the projected differential transform method.
u (x , m + 1) = −E −1 {v 2 E [ A m + B m ]} , u (x , 0) = G (x , t )
Where
(7)
A m , B m are the projected differential transform of Ru (x , t ) , Nu (x , t ) .
From equation (7), we have:
u (x , 0) = G (x , t )
,
u (x ,1) = − E −1 {v 2 E [ A 0 + B 0 ]}
u (x , 2) = −E −1 {v 2 E [ A1 + B 1 ]}
,
u (x ,3) = − E −1 {v 2 E [ A 2 + B 2 ]}
.....
Then the solution of equation (3) is
u (x , t ) = u (x , 0) + u (x ,1) + u (x , 2) + ........
To illustrate the capability and simplicity of the method, some examples for different linear and nonlinear
partial differential equations will be discussed.
Example 2.1:
Consider the simple first order partial differential equation
∂y
∂y
=2
+y
∂x
∂t
, y (x , 0) = 6e −3x
(8)
Taking Elzaki transform of (8), leads to
E [ y (x , t )] = 6v 2e −3 x +
v
E [Am − B m ]
2
Take the inverse Elzaki transform to find,
y (x , m + 1) = E −1 {
v
E [ Am − B m ]
2
}
54
,
y (x , 0) = 6e −3x
(9)
Mathematical Theory and Modeling
ISSN 2224-5804 (Paper)
ISSN 2225-0522 (Online)
Vol.2, No.4, 2012
Where
Am =
www.iiste.org
∂y (x , m )
, B m = y (x , m ) are the projected differential transform of
∂x
∂y (x , t )
, y (x , t ) .
∂x
The standard Elzaki transform defines the solution
y (x , t ) by the series
∞
y (x , t ) = ∑ y (x , m )
m =0
From equation (9) we find that:
y (x , 0) = 6e −3x
A 0 = −18e −3x , B 0 = 6e −3x , y (x ,1) = E =1 −12v 3e −3x = −12te −3x
A1 = 36te −3x , B 1 = −12e −3x , y (x , 2) = E =1 24v 4e −3x = 12t 2e −3x
..................................................................
y (x ,3) = −8t 3e −3x
The solution in a series form is given by
y (x , t ) = 6e −3x − 12te −3x + 12t 2e −3x + ........ = 6e −3x e −2t = 6e − (3x + 2t )
Example 2.2:
Consider the following linear second order partial differential equation
u xx + u tt = 0 , u (x , 0) = 0, u t (x , 0) = cos x
(10)
To find the solution we take Elzaki transform of equation (10) and making use of the conditions to find,
E [u (x , t ) ] = v 3 cos x − v 2 E [ A m ] ,
Am =
∂ 2u (x , m )
∂x 2
Take the inverse Elzaki transform we get:
u (x , m + 1) = −E −1 {v 2 E [ A m ]} , u (x , 0) = 0 , u (x ,1) = t cos x
By using equation (11), we find that:
t3
cos x
3!
t 3
t 5
t3
A 2 = − cos x ⇒ u (x ,3) = E −1 v 2 E cos x = cos x
3!
3!
5!
.
.
.
A1 = −t cos x ⇒ u (x , 2) = E −1 {v 2 E [t cos x ]} =
55
(11)
Mathematical Theory and Modeling
ISSN 2224-5804 (Paper)
ISSN 2225-0522 (Online)
Vol.2, No.4, 2012
Then the solution is
www.iiste.org
t3 t5
u (x , t ) = cos x t + + + ..... = cos x sinh t
3! 5!
Example 2.3:
Consider the following second order nonlinear partial differential equation
∂u ∂u
∂ 2u
=
+
u
∂t ∂x
∂x 2
2
,
u (x .0) = x 2
(12)
To find the solution take Elzaki transform of (12) and using the condition we get:
E [u (x , t ) ] = v 2 x 2 + vE [ A m + B m ]
h
∂u (x , m ) ∂u (x , h − m )
∂ 2u (h − m )
,
, B m = ∑ u (x , m )
∂x
∂x
∂x 2
m =0
m =0
h
Where A m
=∑
are
projected
∂ 2u
∂u
differential transform of
, u
∂x 2
∂x
2
Take the inverse Elzaki transform to get:
u (x , m + 1) = E −1 {vE [ A m + B m ]} , u (x , 0) = x 2
(13)
From equation (13) we find that:
A 0 = 4x 2 , B 0 = 2x 2 ⇒ u (x ,1) = E −1 6x 2v 3 = 6x 2t
A1 = 48x 2t , B 1 = 24x 2t ⇒ u (x , 2) = E −1 72x 2v 4 = 36x 2t 2
.
.
.
Then the solution of equation (12) is
u (x , t ) = x 2 + 6x 2t + 36x 2t 2 + ..... = x 2 (1 − 6t ) −1 =
x2
1 − 6t
Which is an exactly the same solution obtained by the Adomian decomposition method.
Example 2.4:
Let us consider the nonlinear partial differential equation
2
∂u
∂u
2 ∂ u
= 2u
+
u
∂t
∂x 2
∂x
2
u (x , 0) =
,
x +1
2
Taking Elzaki transform of equation (14) and using the condition leads to:
56
(14)
Mathematical Theory and Modeling
ISSN 2224-5804 (Paper)
ISSN 2225-0522 (Online)
Vol.2, No.4, 2012
www.iiste.org
x +1
E [u (x , t ) ] = v 2
+ vE [ 2A m + B m ]
2
Taking the inverse Elzaki transform to find:
u (x , m + 1) = E −1 {vE [ 2A m + B m ]} , u (x , 0) =
h
k
A m = ∑ ∑ u (x , m )
k =0 m = 0
Where
h
x +1
2
(15)
∂u (x , h − m ) ∂u (x , h − k )
∂x
∂x
k
B m = ∑ ∑ u (x , m ) u (x , h − m )
k =0 m =0
∂ 2u (x , h − k )
∂x 2
From equation (15) we have:
u (x , 0) =
x +1
x +1
, u (x ,1) =
t ,
2
2
u (x , 2) =
3(x + 1) 2
t ,...........
8
Then the solution of equation (14) is
u (x , t ) =
1
x +1
x +1
−
(1 − t ) 2 =
2
2 1−t
3. Conclusion
The solution of linear and nonlinear partial differential equations can be obtained using Elzaki transform
and projected differential transform method without any discretization of the variables. The results for all
examples can be obtained in series form, and all calculations in the method are very easy. This technique is
useful to solve linear and nonlinear partial differential equations.
References
[1] S. lslam, Yasir Khan, Naeem Faraz and Francis Austin (2010), Numerical Solution of Logistic
Differential Equations by using the Laplace Decomposition Method, World Applied Sciences Journal 8
(9):1100-1105.
[2] Nuran
Guzel and Muhammet Nurulay (2008), Solution of Shiff Systems By using Differential
Transform Method, Dunlupinar universities Fen Bilimleri Enstitusu Dergisi, ISSN 1302-3055, PP. 49-59.
[3] Shin- Hsiang Chang , I. Ling Chang (2008), A new algorithm for calculating one – dimensional
differential transform of nonlinear functions, Applied Mathematics and Computation 195, 799-808.
[4] Tarig M. Elzaki (2011), The New Integral Transform “Elzaki Transform” Global Journal of Pure and
Applied Mathematics, ISSN 0973-1768, Number 1, pp. 57-64.
[5] Tarig M. Elzaki & Salih M. Elzaki (2011), Application of New Transform “Elzaki Transform” to Partial
Differential Equations, Global Journal of Pure and Applied Mathematics, ISSN 0973-1768,Number 1, pp.
65-70.
[6] Tarig M. Elzaki & Salih M. Elzaki (2011), On the Connections Between Laplace and Elzaki transforms,
Advances in Theoretical and Applied Mathematics, ISSN 0973-4554 Volume 6, Number 1, pp. 1-11.
[7] Tarig M. Elzaki & Salih M. Elzaki (2011), On the Elzaki Transform and Ordinary Differential Equation
With Variable Coefficients, Advances in Theoretical and Applied Mathematics. ISSN 0973-4554 Volume 6,
57
Mathematical Theory and Modeling
ISSN 2224-5804 (Paper)
ISSN 2225-0522 (Online)
Vol.2, No.4, 2012
www.iiste.org
Number 1, pp. 13-18.
[8] Lokenath Debnath and D. Bhatta (2006). Integral transform and their Application second Edition,
Chapman & Hall /CRC.
[9] A.Kilicman and H.E.Gadain (2009). An application of double Laplace transform and Sumudu transform,
Lobachevskii J. Math.30 (3), pp.214-223.
[10] J. Zhang (2007), A Sumudu based algorithm m for solving differential equations, Comp.
Sci. J.
Moldova 15(3), pp – 303-313.
[11] Hassan Eltayeb and Adem
differential
Equations, Applied
kilicman (2010), A Note on
Mathematical
the
Sumudu
Transforms
and
Sciences, VOL, 4, no.22,1089-1098
[12] Kilicman A. & H. ELtayeb (2010). A note on Integral transform and Partial Differential Equation,
Applied Mathematical Sciences, 4(3), PP.109-118.
[13] Hassan ELtayeh and Adem kilicman (2010), on Some Applications of a new Integral Transform, Int.
Journal of Math. Analysis, Vol, 4, no.3, 123-132.
[14] N.H. Sweilam, M.M. Khader (2009). Exact Solutions of some capled nonlinear partial differential
equations using the homotopy perturbation method. Computers and Mathematics with Applications 58
2134-2141.
[15] P.R. Sharma and Giriraj Methi (2011). Applications of Homotopy Perturbation method to Partial
differential equations. Asian Journal of Mathematics and Statistics 4 (3): 140-150.
[16] M.A. Jafari, A. Aminataei (2010). Improved Homotopy Perturbation Method. International
Mathematical Forum, 5, no, 32, 1567-1579.
[17] Jagdev Singh, Devendra, Sushila (2011). Homotopy Perturbation Sumudu Transform Method for
Nonlinear Equations. Adv. Theor. Appl. Mech., Vol. 4, no. 4, 165-175.
58
Mathematical Theory and Modeling
ISSN 2224-5804 (Paper)
ISSN 2225-0522 (Online)
Vol.2, No.4, 2012
www.iiste.org
Table 1. Elzaki transform of some Functions
f (t )
E [ f (t )] = T (u )
1
v2
t
v3
tn
n ! v n +2
e at
v2
1 − av
sinat
av 3
1 + a 2v 2
cosat
v2
1 + a 2v 2
59
International Journals Call for Paper
The IISTE, a U.S. publisher, is currently hosting the academic journals listed below. The peer review process of the following journals
usually takes LESS THAN 14 business days and IISTE usually publishes a qualified article within 30 days. Authors should
send their full paper to the following email address. More information can be found in the IISTE website : www.iiste.org
Business, Economics, Finance and Management
PAPER SUBMISSION EMAIL
European Journal of Business and Management
EJBM@iiste.org
Research Journal of Finance and Accounting
RJFA@iiste.org
Journal of Economics and Sustainable Development
JESD@iiste.org
Information and Knowledge Management
IKM@iiste.org
Developing Country Studies
DCS@iiste.org
Industrial Engineering Letters
IEL@iiste.org
Physical Sciences, Mathematics and Chemistry
PAPER SUBMISSION EMAIL
Journal of Natural Sciences Research
JNSR@iiste.org
Chemistry and Materials Research
CMR@iiste.org
Mathematical Theory and Modeling
MTM@iiste.org
Advances in Physics Theories and Applications
APTA@iiste.org
Chemical and Process Engineering Research
CPER@iiste.org
Engineering, Technology and Systems
PAPER SUBMISSION EMAIL
Computer Engineering and Intelligent Systems
CEIS@iiste.org
Innovative Systems Design and Engineering
ISDE@iiste.org
Journal of Energy Technologies and Policy
JETP@iiste.org
Information and Knowledge Management
IKM@iiste.org
Control Theory and Informatics
CTI@iiste.org
Journal of Information Engineering and Applications
JIEA@iiste.org
Industrial Engineering Letters
IEL@iiste.org
Network and Complex Systems
NCS@iiste.org
Environment, Civil, Materials Sciences
PAPER SUBMISSION EMAIL
Journal of Environment and Earth Science
JEES@iiste.org
Civil and Environmental Research
CER@iiste.org
Journal of Natural Sciences Research
JNSR@iiste.org
Civil and Environmental Research
CER@iiste.org
Life Science, Food and Medical Sciences
PAPER SUBMISSION EMAIL
Journal of Natural Sciences Research
JNSR@iiste.org
Journal of Biology, Agriculture and Healthcare
JBAH@iiste.org
Food Science and Quality Management
FSQM@iiste.org
Chemistry and Materials Research
CMR@iiste.org
Education, and other Social Sciences
PAPER SUBMISSION EMAIL
Journal of Education and Practice
JEP@iiste.org
Journal of Law, Policy and Globalization
JLPG@iiste.org
New Media and Mass Communication
NMMC@iiste.org
Journal of Energy Technologies and Policy
JETP@iiste.org
Historical Research Letter
HRL@iiste.org
Public Policy and Administration Research
PPAR@iiste.org
International Affairs and Global Strategy
IAGS@iiste.org
Research on Humanities and Social Sciences
RHSS@iiste.org
Developing Country Studies
DCS@iiste.org
Arts and Design Studies
ADS@iiste.org
Global knowledge sharing:
EBSCO, Index Copernicus, Ulrich's
Periodicals Directory, JournalTOCS, PKP
Open Archives Harvester, Bielefeld
Academic Search Engine, Elektronische
Zeitschriftenbibliothek EZB, Open J-Gate,
OCLC WorldCat, Universe Digtial Library ,
NewJour, Google Scholar.
IISTE is member of CrossRef. All journals
have high IC Impact Factor Values (ICV).