# Sumudu Transform Method for Analytical Solutions of Fractional Type Ordinary Differential Equations

Authors: Seyma Tuluce Demiray; Hasan Bulut; Fethi Bin Muhammad Belgacem
Journal: Mathematical Problems in Engineering
 (2015)

---

## Abstract

We make use of the so-called Sumudu transform method (STM), a type of ordinary differential equations with both integer and noninteger order derivative. Firstly, we give the properties of STM, and then we directly apply it to fractional type ordinary differential equations, both homogeneous and inhomogeneous ones. We obtain exact solutions of fractional type ordinary differential equations, both homogeneous and inhomogeneous, by using STM. We present some numerical simulations of the obtained solutions and exhibit two-dimensional graphics by means of Mathematica tools. The method used here is highly efficient, powerful, and confidential tool in terms of finding exact solutions.


---

## Body

## 1. Introduction

Computation and analysis of solutions for nonlinear fractional differential equations now span a half-century or more and play a crucial role in several theoretical and applied sciences such as, but certainly not limited to, theoretical biology and ecology, solid state physics, viscoelasticity, fiber optics, signal processing and electric control theory, stochastic based finance, and hemo-, hydro-, and thermodynamics [1–7]. In particular, modified Kudryashov method [8], generalized Kudryashov method (GKM) [9], adomian decomposition method (ADM) [10], and homotopy decomposition method (HDM) [11] have all been used to seek solutions to fractional differential equations.The Sumudu transform method (STM) was initiated in 1993, by Watugala who used it to solve engineering control problems [12, 13]. Subsequently, he expanded the Sumudu transform to two variables [14]. The first application of the STM to partial differential equations and the establishment of the inverse formula were done by Weerakoon [15, 16]. Sumudu transform based solutions to convolution type integral equations and discrete dynamic systems were later obtained by Asiru [17–19]. Following in this trend, Belgacem et al. established applications of the STM [20, 21]. Most of the STM applications in the first decade of the 21st century concentrated on integer type differential equations, and although Belgacem’s practical interests for applying to fractional differential equations started as early as 2008 in various local presentations and conference talks, concrete results in this direction started appearing only in the second decade, with various teams of research [22], including the Bulut et al. team [23]. In fact, the Bulut et al. teams of research have been most active in the integer type differential equations area using the STM [24, 25]. Also, Atangana and Kilicman have tackled Sumudu transform method for nonlinear fractional heat-like equations [26].In this work, our aim is to exhibit exact solutions of some homogeneous and nonhomogeneous fractional ordinary differential equations by using the STM. For paper layout, we first recall basic features of fractional calculus, in Section2. In Section 3, we remind the reader of the properties of the STM. In Section 4, we apply the STM to obtain new exact solutions of both homogeneous and inhomogeneous fractional ordinary differential equations.


## 2. Basic Features of Fractional Calculus

In this section, we primarily introduce main features of fractional calculus following notations [1–7], making consensus among them. Fractional derivatives and integrals are most commonly introduced and used in the sense of Abel-Riemann (A-R), described by (1)
D
α
f
t
=
1
Γ
m
-
α
d
d
t
m
∫
0
t
f
τ
t
-
τ
α
-
m
+
1
d
τ
,
k
k
k
k
k
k
k
k
k
k
k
k
k
m
-
1
<
α
<
m
,
d
m
d
t
m
f
t
,
k
k
k
k
k
α
=
m
,where m
∈
Z
+
,
α
∈
R
+. D
α is derivative operator and (2)
D
-
α
f
t
=
1
Γ
α
∫
0
t
t
-
τ
α
-
1
f
τ
d
τ
,
0
<
α
≤
1
.Also, according to A-R, an integral of fractional order is described by performing the integration operator J
α in the following format: (3)
J
α
f
t
=
1
Γ
α
∫
0
t
t
-
τ
α
-
1
f
τ
d
τ
,
t
>
0
,
α
>
0
.Following Podlubny [3] set-up, we have (4)
J
α
t
n
=
Γ
1
+
n
Γ
1
+
n
+
α
t
n
+
α
,
D
α
t
n
=
Γ
1
+
n
Γ
1
+
n
-
α
t
n
-
α
.Another important and basic definition was first established by Caputo [1], given by (5)
C
D
α
f
t
=
1
Γ
m
-
α
∫
0
t
f
m
τ
t
-
τ
α
-
m
+
1
d
τ
,
k
k
k
k
k
k
k
k
k
k
m
-
1
<
α
<
m
,
d
m
d
t
m
f
t
,
α
=
m
.A basic property of the Caputo fractional derivative is that (see for instance [22]) (6)
J
α
D
C
α
f
t
=
f
t
-
∑
k
=
0
∞
f
k
0
+
t
k
k
!
.


## 3. Basic Features of the Sumudu Transform Method (STM)

The Sumudu transform is obtained over the set of functions [20] as (7)
A
=
f
t
∣
∃
M
,
τ
1
,
τ
2
>
0
,
f
t
<
M
e
t
/
τ
i
,
k
k
k
k
k
k
k
i
f
t
∈
-
1
j
×
0
,
∞by (8)
F
u
=
S
f
t
=
∫
0
∞
f
u
t
e
-
t
d
t
,
u
∈
-
τ
1
,
τ
2
.Theorem 1.
IfF
u is Sumudu transform of f
t, then the Sumudu transform of the derivatives with integer order is as follows [20–23]:(9)
S
d
n
f
t
d
t
n
=
u
-
n
F
u
-
∑
k
=
0
n
-
1
u
k
d
n
f
t
d
t
n
t
=
0
.Proof.
The Sumudu transform of the first derivative off
t, f
′
t
=
d
f
t
/
d
t, is given by (10)
S
d
f
t
d
t
=
∫
0
∞
d
f
u
t
d
t
e
-
t
d
t
=
l
i
m
p
→
∞
∫
0
p
d
f
u
t
d
t
e
-
t
d
t
=
l
i
m
p
→
∞
1
u
e
-
t
/
u
f
t
0
p
+
1
u
2
∫
0
p
e
-
t
/
u
f
t
d
t
=
l
i
m
p
→
∞
1
u
e
-
t
/
u
f
t
0
p
k
k
k
k
k
k
k
k
+
1
u
1
u
∫
0
p
e
-
t
/
u
f
t
d
t
=
l
i
m
p
→
∞
-
1
u
f
0
+
1
u
1
u
∫
0
p
e
-
t
/
u
f
t
d
t
=
-
1
u
f
0
+
1
u
F
u
.Proceeding in the same manner, we get the Sumudu transform of the second order derivative as (11)
S
d
2
f
t
d
t
2
=
1
u
2
F
u
-
f
0
-
u
d
f
t
d
t
t
=
0
.To finish the proof of Theorem 1, we proceed by induction in the same way only to reach the general formula for the Sumudu transform of any integer n
-order derivative [23]; namely, (12)
S
d
n
f
t
d
t
n
=
u
-
n
F
u
-
∑
k
=
0
n
-
1
u
k
d
n
f
t
d
t
n
t
=
0
.Theorem 2.
IfF
u is the Sumudu transform of f
t, then the Sumudu transform of the Riemann-Liouville fractional derivative is given by (see [22])(13)
S
D
α
f
t
=
u
-
α
F
u
-
∑
k
=
1
n
u
α
-
k
D
α
-
k
f
t
t
=
0
,
K
K
K
K
K
K
K
K
K
K
K
K
K
-
1
<
n
-
1
≤
α
<
n
.At the moment, we harness the properties developments of the STM and utilize them for finding exact solutions of fractional ordinary differential equations.We consider the general linear fractional ordinary differential equation (FODE) as follows:(14)
∂
α
y
t
∂
t
α
=
∂
2
y
t
∂
t
2
+
∂
y
t
∂
t
+
y
t
+
c
,subject to the initial condition(15)
y
0
=
f
0
.When we get Sumudu transform of (14) taking into consideration (12) and (13), we obtain Sumudu transform of (14) as follows:(16)
S
∂
α
y
t
∂
t
α
=
S
∂
2
y
t
∂
t
2
+
S
∂
y
t
∂
t
+
S
y
t
+
S
c
,
u
-
α
F
u
-
∑
k
=
1
n
u
α
-
k
D
α
-
k
f
t
t
=
0
=
1
u
2
F
u
-
f
0
-
u
∂
f
t
∂
t
t
=
0
+
1
u
F
u
-
f
0
+
F
u
+
c
,
F
u
-
∑
k
=
1
n
u
α
-
k
D
α
-
k
f
t
t
=
0
=
u
α
-
2
F
u
-
f
0
-
u
∂
f
t
∂
t
t
=
0
+
u
α
-
1
F
u
-
f
0
+
u
α
F
u
+
c
u
α
,
F
u
=
u
α
-
2
F
u
-
u
α
-
2
f
0
+
∑
k
=
1
n
u
α
-
k
D
α
-
k
f
t
t
=
0
-
u
α
-
1
∂
f
t
∂
t
t
=
0
+
u
α
-
1
F
u
-
u
α
-
1
f
0
+
u
α
F
u
+
c
u
α
,
F
u
-
u
α
-
2
F
u
-
u
α
-
1
F
u
-
u
α
F
u
=
-
u
α
-
2
f
0
+
∑
k
=
1
n
u
α
-
k
D
α
-
k
f
t
t
=
0
-
u
α
-
1
∂
f
t
∂
t
t
=
0
-
u
α
-
1
f
0
+
c
u
α
,
F
u
=
1
1
-
u
α
-
2
-
u
α
-
1
-
u
α
×
g
u
-
u
α
-
1
f
0
-
u
α
-
2
f
0
+
c
u
α
,where g
u is described by ∑
k
=
1
n
u
α
-
k
D
α
-
k
f
t
t
=
0
-
u
α
-
1
∂
f
t
/
∂
t
t
=
0.When we get inverse Sumudu transform of (16) by using inverse transform table in [21, 22], we obtain solution of (14) by using STM in the following manner: (17)
y
t
=
S
-
1
1
1
-
u
α
-
2
-
u
α
-
1
-
u
α
k
k
k
k
k
k
×
g
u
-
u
α
-
1
f
0
-
u
α
-
2
f
0
+
c
u
α
1
1
-
u
α
-
2
-
u
α
-
1
-
u
α
.


## 4. STM Implementations to Homogeneous and Inhomogeneous FODEs

In this Section, we implement the STM to homogeneous and inhomogeneous fractional ordinary differential equations (HFODEs and IHFODEs) in the following three examples.Example 3.
Initially, we consider the inhomogeneous fractional ordinary differential equation [4] as(18)
D
2
y
t
+
D
3
/
2
y
t
+
y
t
=
t
+
1
,
D
0
t
-
1
/
2
y
t
t
=
0
=
0
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
i
k
k
k
t
>
0
,subject to the initial conditions(19)
y
0
=
y
′
0
=
1
.In order to find exact solution of (18), we take the Sumudu transform of both sides of (18) as follows: (20)
S
D
2
y
t
+
S
D
3
/
2
y
t
+
S
y
t
=
S
t
+
S
1
,
1
u
2
F
u
-
f
0
-
u
d
f
t
d
t
t
=
0
+
F
u
u
3
/
2
-
D
1
/
2
f
t
u
t
=
0
+
F
u
=
u
+
1
,
1
u
2
F
u
-
1
-
u
+
F
u
u
3
/
2
-
1
u
F
u
u
1
/
2
-
D
-
1
/
2
f
t
u
t
=
0
+
F
u
=
u
+
1
,
1
u
2
F
u
-
1
-
u
+
D
-
1
/
2
f
t
u
t
=
0
+
F
u
=
u
+
1
,
1
+
1
u
2
F
u
=
1
+
u
+
1
u
+
1
u
2
⟹
1
+
u
2
F
u
=
1
+
u
+
u
2
+
u
3
,
F
u
=
1
+
u
+
u
2
+
u
3
1
+
u
2
=
1
+
u
.Upon Sumudu inverting (20) which is linear (also see inverse transform table in [21]), we obtain exact solution as (21)
y
t
=
1
+
t
.Remark 4.
The exact solution of (18) found here by using Sumudu transform method agrees with the solution in [4] for corresponding values of tuned parameters.Example 5.
Secondly, we consider homogeneous fractional ordinary differential equation as follows [6]:(22)
d
y
t
d
t
+
d
1
/
2
y
t
d
t
1
/
2
-
2
y
t
=
0
,
d
1
/
2
y
0
+
d
t
1
/
2
=
c
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
i
k
k
k
t
>
0
,with initial condition (23)
y
0
+
=
0
.We first apply the Sumudu transform of both sides of (22) to get(24)
S
d
y
t
d
t
+
S
d
1
/
2
y
t
d
t
1
/
2
-
2
S
y
t
=
0
,which yields(25)
1
u
F
u
-
f
0
+
F
u
u
1
/
2
-
D
-
1
/
2
f
t
u
t
=
0
-
2
F
u
=
0
,
F
u
u
+
F
u
u
1
/
2
-
c
u
-
2
F
u
=
0
,
1
+
u
1
/
2
-
2
u
F
u
=
c
⟹
F
u
=
c
1
+
u
1
/
2
-
2
u
k
k
k
k
k
k
k
k
k
k
k
=
2
c
3
1
+
2
u
+
c
3
1
-
u
,
F
u
=
2
c
3
1
-
2
u
1
-
4
u
+
c
3
1
+
u
1
-
u
,
F
u
=
2
c
3
1
1
-
4
u
-
2
u
1
-
4
u
+
c
3
1
1
-
u
-
u
1
-
u
.
When we take inverse Sumudu transform of (25) by using inverse transform table in [21], we obtain exact solution of (22) by using STM as follows: (26)
y
t
=
2
c
3
e
4
t
-
e
4
t
e
r
f
2
t
+
c
3
e
t
+
e
t
e
r
f
t
.Remark 6.
To our knowledge, the solution of (22), obtained by using the STM, is new and does not figure in the published literature.Example 7.
Thirdly, we investigate the inhomogeneous fractional ordinary differential equation [6] as(27)
D
0
t
1
/
2
y
t
-
y
t
=
-
1
,
D
0
t
-
1
/
2
y
t
t
=
0
=
0
,
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
k
i
k
k
k
t
>
0
.In order to obtain the exact solution of (27) by using the STM, we apply the Sumudu transform to both sides of (27) as follows: (28)
S
D
1
/
2
y
t
-
S
y
t
+
S
1
=
0
,
S
D
1
/
2
y
t
-
F
u
=
-
S
1
,
F
u
u
1
/
2
-
D
-
1
/
2
f
t
u
t
=
0
-
F
u
=
-
S
1
,
F
u
u
1
/
2
-
F
u
=
-
1
⟹
1
u
1
/
2
-
1
F
u
=
-
1
⟹
1
-
u
1
/
2
F
u
=
-
u
1
/
2
,
F
u
=
u
u
-
1
=
u
u
+
1
u
-
1
=
u
+
u
u
-
1
=
u
-
1
+
1
+
u
u
-
1
=
1
-
1
1
-
u
-
u
1
-
u
.When we get the inverse Sumudu transform of (28) by using the table in [21], we find exact solution of (27) by the STM as follows:(29)
y
t
=
1
-
e
t
-
e
t
e
r
f
t
.Remark 8.
Again, the exact solution of (27) exhibited here, found through implementing the Sumudu transform method, is to our knowledge new and not trackable in previous literature.
We plot solution (21) of (18) in Figure 1, solution (26) of (22) in Figure 2, solution (29) of (27) in Figure 3, which shows the dynamics of solutions with suitable parametric choices.Figure 1
Two-dimensional graphic for exact solution of (18), when 0
<
t
<
5, via the STM.Figure 2
Two-dimensional graphic for exact solution of (22), for 0
<
t
<
5 and c
=
3, via the STM.Figure 3
Two-dimensional graphic for exact solution of (27), for 0
<
t
<
5, via the STM.


## 5. Conclusions

In this paper, we considered homogeneous and inhomogeneous fractional ordinary differential equations (HFODEs and IHFODEs) and treated them by applying the Sumudu transform method and obtained their exact solutions. To date and to our knowledge, the exact solutions obtained in the last two examples are new. We also exhibited two-dimensional graphics for the obtained solutions by means of programming language Mathematica. Accordingly, we can conclude that not only does the Sumudu transform method play an important role in treating homogeneous and inhomogeneous fractional ordinary differential equations, but also it is highly effective with regard to yielding exact solutions. We expect the STM to be equally successful when treating even more complex applications than presented involving FODEs.



---

*Source: 131690-2015-06-29.xml*