CC BY
253Probl. Anal. Issues Anal. Vol. 7(25), No. 1, 2018, pp. 41-58
DOI: 10.15393/j3.art.2018.4350
41
UDC 519.71, 517.977, 517.44
A. A. Hamüüd, K. P. Ghadle
THE APPROXIMATE SOLUTIONS OF FRACTIONAL VOLTERRA - FREDHOLM INTEGRO-DIFFERENTIAL EQUATIONS BY USING ANALYTICAL TECHNIQUES
Abstract. This paper demonstrates a study on some significant latest innovations in the approximation techniques to find the approximate solutions of Caputo fractional Volterra - Fredholm integro-differential equations. To apply this, the study uses Ado-mian decomposition method and modified Laplace Adomian decomposition method. A wider applicability of these techniques is based on their reliability and reduction in the size of the computational work. This study provides analytical approximate to determine the behavior of the solution. It proves the existence and uniqueness results and convergence of the solution. In addition, it brings an example to examine the validity and applicability of the proposed techniques.
Key words: Adomian decomposition method, Laplace transform, Volterra - Fredholm integro-differential equation, Caputo fractional derivative
2010 Mathematical Subject Classification: 26A33, 49M27
1. Introduction. In this paper, we consider the Caputo fractional Volterra - Fredholm integro-differential equations of the form:
X 1
c Day(x) = g(x)+ a(x)y(x) + J Ki (x,t)Fi (y(t))dt + j K2(x,t)F2(y(t))dt,
0 0
(1)
with the initial conditions
y
« (0) = St, i = 0,1, ••• ,n - 1, (2)
©Petrozavodsk State University, 2018
where n — 1 < a < n, n E N, y : [0,1] —> R is the continuous function to be determined, g, a : [0,1] —> R and K : [0,1] x [0,1] —> R,i = 1, 2, are continuous functions. Fi : R —> R,i = 1, 2, are Lipschitz continuous functions. The fractional derivative cDa is understood here in the Caputo sense.
The fractional integro-differential equations have attracted much more interest of mathematicians and physicists because they provide an efficiency for the description of many practical dynamical problems arising in engineering and scientific disciplines such as, physics, biology, electrochemistry, chemistry, economy, electromagnetic, control theory and viscoelasticity [3, 19, 20]. In recent years, numerous papers have been concentrating on the development of numerical and analytical techniques for fractional integro-differential equations. There are many approaches for seeking exact solutions of linear and nonlinear equations, for example, Al-Samadi and Gumah [3] applied the homotopy analysis method for fractional SEIR epidemic model, Jafarian et al. [14] applied successfully artificial neural networks approach to solve fractional order Volterra integro-differential equations. Also, some comparisons had been made between this iterative approach and another traditional technique. The obtained results reveal that this method is very effective, Momani [19] and Qaralleh [20] applied Adomian polynomials to solve fractional integro-differential equations and systems of fractional integro-differential equations, Ka-dem and Kilicman [15] utilized the HPM and VIM methods for integro-differential equation of fractional order with initial-boundary conditions, Yang [21] used the hybrid of block pulse function and Chebyshev polynomials to solve nonlinear Fredholm fractional integro-differential equations, Yang and Hou [22] applied the Laplace decomposition method to solve the fractional integro-differential equations, Mittal and Nigam [18] utilized the Adomian decomposition method to approximate solutions of fractional integro-differential equations, and Ma and Huang [17] applied hybrid collocation method to study integro-differential equations of fractional order. Moreover, properties of the fractional integro-differential equations have been studied by several authors [4, 12, 16].
The main objective of the present paper is to study the behavior of the solution that can be formally determined by such analytical approximation methods as the Adomian decomposition method and the modified Laplace Adomian decomposition method. Moreover, we prove the convergence and uniqueness solution of the fractional Volterra - Fredholm integro-differential equation.
The rest of the paper is organized as follows. In section 2 some preliminaries and basic definitions related to fractional calculus and Laplace transform are recalled. In section 3 Adomian decomposition method is constructed for solving fractional Volterra-Fredholm integro-differential equations. In section 4 modified Laplace Adomian decomposition method is constructed for solving Volterra-Fredholm integro-differential equations of fractional order. In section 5 the convergence and uniqueness of the solution are proved. In section 6 the analytical example is presented to illustrate the accuracy of these methods. Finally, a brief conclusion is given in section 7.
2. Preliminaries. The mathematical definitions of fractional derivative and fractional integration are the subject of several different approaches. The most frequently used definitions of the fractional calculus involve the Riemann - Liouville fractional derivative, Caputo derivative, Riesz derivative and Granewald-Letnikov fractional derivative [6, 8, 16, 18, 22]. In our work the Caputo's definition of fractional derivative has been used.
Definition 1. The Riemann - Liouville fractional integral of order a > 0 of a function f is defined as
x
Jaf (x) = - t)a-1 f (t)dt, x > 0, a G R+,
r(a) J
0
J°f (x) = f (x), where R+ is the set of positive real numbers.
Definition 2. The fractional derivative of f (x) in the Caputo sense is defined by
cDa f (x) = Jn-a Dnf (x) =
x
1 C dn f (t)
1 1 (x - t)n-a-1 f^ dt, n - 1 <a<n,
^ ' J-Ln. ' '
=
r(n — a) J dt
dnf (x)
dxn
a=n,
where the parameter a is the order of the derivative and is allowed to be real or even complex.
In this paper only real and positive a will be considered. Hence, we have the following properties:
• JaJv f = Ja+v f, a,v> 0.
• Jaxp = F(( + ^ x xa+p, a > 0, ( > -1, x > 0.
r(a + ( + 1)
n— 1 k
• JaDa f (x) = f (x) - £ f (k)(0+) —, x > 0, n - 1 < a < n.
k=0 '
Definition 3. The Riemann - Liouville fractional derivative of order a> 0 is normally defined as
Daf (x) = DnJn—af (x), n - 1 < a < n. (3)
Definition 4. The Laplace transform of a function f (x), is defined as
L[f (x)] = F(s) = J f (x)e—sxdx, 0
where x > 0, s can be either real or complex.
Definition 5. Given two functions f and g, we define, for any x > 0,
x
(f * g)(x) = y f (t)g(x - t)dt.
0
The function f * g is called the convolution of f and g.
Theorem 1. The Laplace transform of the Caputo derivative is defined by
n— 1
L[cDaf (x)] = saF(s) sa—k—1 f(k) (0), n - 1 < a < n. (4)
k=0
3. Adomian Decomposition Method. Consider the equation (1) where cDa is the operator defined as (3). Operating with Ja on both sides of the equation (1), we get [1, 2, 5, 7, 18]
n 1 k k
y(x) = £ y(k) (0+)-TT + Ja(g(x))+
+ Ja (a(x)y(x) + J Ki(x,t)Fi(y(t))dt + j K2(x,t)F2(y(t))d?j.
0 0
Adomian decomposition method defines the solution y(x) by the series
y — Y1 yn
n=0
and the nonlinear functions F1 and F2 are decomposed as
Fi = ^ An, F2 = ^ Bn,
n=0
2 — Bn-,
n=0
where An, Bn are the Adomian polynomials given by
An — ,
n!
Bn — —. n!
dn dфn
' dn
dфп
n
Fi(J2 Фг У*
i=0 n
Фг y
i=0
ф=0
ф=0
(5)
(6)
(7)
(8)
The Adomian polynomials were introduced in [1, 2, 10, 11, 13] as:
Ao = Fi (yo), Ai = yi Fi (yo),
A2 — У2 F' (y0 ) + 2 У2 F'' (У0 ),
Лз —
1
УзFi (y0) + yiy2Fi (y0) + 3у 1 Fi (У0)-
and
B0 B
B2
F2 (У0),
У iF2 (У0)-
1
y2F2 Ы + 2У2F2 (У0)-
B3 = y3F2 (yo) + yiV2F2 (y0) + 1 y3F2'' M,
The components y0,y1,y2, ■ ■ ■ are determined recursively by
n — 1 k
yo = £ y(k) (0+) J" g(x), k=0 '
x 1
yk+1 = J" (a(x)y(x) + J Ki (x,t)Ak dt + J K2(x,t)Bk d?j.
0 0
Having defined the components y0,y1,y2, ■ ■ ■, the solution y in a series form defined by (5) follows immediately. It is important to note that the decomposition method suggests that the 0th component y0 is to be defined by the initial conditions and the function g(x) as described above. The other components, namely y1 ,y2, ■ ■ ■, are derived recurrently.
4. Modified Laplace Adomian Decomposition Method. We
consider the fractional Volterra - Fredholm integro-differential equation (1). We apply the Laplace transform to both sides [8, 9, 22]:
L [cD"y(x)] = L[g(x)] + L[a(x)y(x)] +
+ L
K1(x,s)F1 (y(s))ds + K2(x,s)F2 (y(s))ds
Using the differentiation property of Laplace transform (4), we get saL[y(x)] - yo = L[g(x)] + L[a(x)y(x)] +
x 1
+ L J Ki(x,s)Fi(y(s))ds + J K2(x,s)F2(y(s))ds
oo
m-1
where y0 = xa-k-iy(k (0). Thus, the given equation is equivalent to
k=0
L[y(x)] = yo + saL[g(x)] + L [a(x)y(x)] +
1
x
+ — L
sa
K1(x,s)F1 (y(s))ds + K2(x,s)F2(y(s))ds
(9)
Substituting (5) and (6) into (9), we have yo , 1
LIZ yn (x) = sa + ¡aL [g(x)] + sa L[a(x)y(x)] + sa
n=0
x
x 1
/00 » oo
Ki(x, s) ^ Ands + / K2 (x, s) ^ Bnds
n n=0 n n=0
(10)
Matching both sides of (10) yields the following iterative algorithm:
£[yo(x)]
yo +1 L[»(x)]-
L[yi(x)] =
L[y2(x)] =
J-L
sa
a(x)y(x)+ K1(x,s)A0 ds + K2 (x,s)B0 ds
0
0
a(x)y(x) + J K1(x,s)A1 ds + J K2 (x,s)B1 ds 00
L[yn+1(x)] =
a(x)y(x)+ K1(x,s)An ds + K2(x,s)Bnds
The solution y(x) is defined by the series
00
y(x) = ^2 yn(x).
n=o
(11)
5. Uniqueness and Convergence. In this section we present an existence and uniqueness results of equation (1) with the initial condition (2) and prove it. Before starting and proving the main results, we introduce the following hypotheses:
1
x
1
x
1
x
s
1
x
a
s
(A1) There exist two constants LFl ,LF2 > 0 such that, for any y1 and y2 E C(J, R)
|Fi(yi(x)) - Fi(y2(x))|< Lfi |yi - y2|
and
|F2(yi(x)) - F2 (y2(x))| < Lf2 |yi - y2 | .
(A2) There exist two functions K\,K% E C(D, R+), where C(D, R+) is the set of all positive functions, continuous on D = {(x,t) E R x R : 0 < t < x < 1}, such that
x i
= sup / |Ki(x, t)| dt < K*2 = sup / |K2(x, t)| dt < xe[o,i]0 xe[o,i]0
(A3) The two functions a, g : J ^ R are continuous.
Lemma 1. If y0 E C (J, R), then y(x) E C (J, R+) is a solution of the problem (1) - (2) iff y satisfies
x x
y(x) = -0)/(x - s)a-ia(s)y(s)ds + —0) J(x - s)a-ig(s)ds +
oo
x _ s
1 ' ' 1
+ T(a)J(x - s)
00
1
+ j K2(s,T)F2(v(t))dT
0
K1(s,T )F1(y(T ))dT +
ds + y0
for x E ■]■
Our first result is based on the Banach contraction principle. Theorem 2. Assume that (A1), (A2) and (A3) hold. If
( ||a|L + KtLFl + K|Lf2 ) ^ . (12)
V r(a +1) )< 1 (12)
then there exists a unique solution y(x) E C(J) of (1) - (2).
Proof. By Lemma 1 we know that y is a solution of (1) - (2) iff it satisfies
x x
y(x) = -— (x - s)a-ia(s)y(s)ds + -—■ (x - s)a-ig(s)ds + r(a) J r(a) J
o o
x
+ PR / (x -
0 0
1
+ [ K2(s,r)F2(y(r))dT
0
Ki(s,T )Fi(y(T ))dT +
ds + yo.
Let the operator T : C(J, R) ^ C(J, R) be defined by
x x
(Ty)(x) = r(ay/(x - s)a-1a(s)y(s)ds + J(x - s)a-1 g(s)ds + 00
x
+ W)S(x - s)"1
00
1
+ y K2 (s,T )F2(y(T ))dT 0
Ki (s,T )Fi(y(T ))dT +
ds + yo.
We can see that, if y G C(J, R) is a fixed point of T, then y is a solution
of (1)-(2).
Now we prove that T has a fixed point y in C (J, R). For that, let y1,y2 G C (J, R) be such that for any x G [0,1]
x x
y1 (x) = -— (x - s)a-1a(s)y1(s)ds + -—■ (x - s)a-1g(s)ds + r(a) J r(a) J
00
x
+ r(a) /(X - s)"-1
00
1
+ J K2(s,T )F2(yi (t ))dT 0
Ki(s,T )Fi(yi (t ))dT +
ds + yo,
x x
V2(x) = -— (x - s)"-1a(s)V2(s)ds + -—■ (x — s)a-1g(s)ds +
r(a) J r(a) J
00
x
+ pm /(x - s)"-1
00
k1(s,t )f1(y2 (t ))dT +
+ J K2(s,T )F2(V2 (T ))dT
0
Consequently, we get
|(Ty1 )(x) - (TV2)(x)|<
ds + V0 ■
x
1 f 1
< (x — s)"-1 |a(s)||V1(s) - V2(s)| ds + x
r(a) J r(a)
0
( f K^T^^y (t)) - f1(y2(t))| dT + ^
x J (x — s)
0
a— 1
0
1
ds
+ / |K2 (s,T)||F2 (V1(T)) — F2 (V2(t))| dT ,
0
a Kt L
< |V1 (x) — V2w + rO+b |V1(x) — V2(x)| +
+r(aaLiFi)|y1(x)—y2(x)|=
("""" +Ka+1+ K|Lp2) |V1 (x) — V2(x)|-
+ Kt Lfi + K|Lf2 r(a + 1)
From the inequality (12) we have
||TV1 — TV2^ < Hv1— V2^ ■
So, T is a contraction map. By the Banach contraction principle, we conclude that T has an unique fixed point y in C (J, R) □
Theorem 3. Suppose that (A1) - (A3) and (12) hold. If the series solution y(x) = Vi(x) and ||y:i.||œ < ^, obtained by the m-order deforma-
i=0
tion, is convergent, then it converges to the exact solution of the fractional Volterra -Fredholm integro-differential equation (1) - (2).
1
x
Proof. Denote as (C[0,1], || • ||^) the Banach space of all continuous functions on J, with |yi(x)| < to for all x in J.
Frist we define the sequence of partial sums sn. Let sn and sm be
n
arbitrary partial sums with n > m. We are to prove that sn = ^^yi(x)
i=o
is a Cauchy sequence in this Banach space:
||sn - sm= max \sn - sm\ = max
VxeJ VxeJ
n m
J2yi(x) -Y1 yi(x) i=0 i=0
= max
WxeJ
i=m+1
a(t)yi (t)+
t 1 + J K1 (t, s)Ai (s)ds + J K2 (t,s)Bi (s)ds 00
dt
max
WxeJ
1
r(a)
(x - t)
a1
n1
a(t) V yi(t)+
0
i=m
„ n—1 1 n—1
+ K1 (t,s)J2 Ai(s)ds + K2(t,s)J2 Bi(s)d,
i=m
n— 1
2(t,s) V Bi(s)ds
dt
From (5) and (6), we have
n1
"Y^Ai = F1 (sn—1) - F1 (sm—1),
n—1
^ Bi = F2 (sn —1) - F2 (sm —1),
i=m n1
^2yi = y(sn— 1) - y(sm — 1)-
x
So,
||sn sm||TO — max
VxEJ
r(a)
(x - t)
a—i
a(t)(y(sn-i) - y(sm-i)) +
+ J Ki(t, s)(Fi(sn-i) - Fi(sm-i))ds + o
+ / K2(t,s)(F2(sn-i) - F2(sm-i))ds
dt
<
max
WxeJ
( 1 x
r(a)
lx - tl
a— i
|a(t)||y(sn-i) - y(sm — i)| +
+ J Ki(t,s)||(Fi(sn—i) - Fi(sm—i))|ds +
o
+ / |K2 (t,s)||(F2(sn—i) - F2(sm —i))|ds
dt I <
<
1
r(a + 1)
||a(t)|| oo||sn — i sm — i H oo + k{ lfl x
X ||sn —i - sm — Joo + K2LF2 y sn— i sm — i Hoo
(
||a||oo + KILfi + KILF2 r(a + 1)
)
||sn— i sm — i y oo
5ysn — i sm — i Hoo)
where
5 =
(
||a||oo + Kl Lfi + KZLFV r(a + 1)
-)
Let n = m + 1, then
||sn - sm ||oo < 5|sm - sm — i||oo < . . . < 5m|si - so|c
x
1
t
i
t
i
^n - sm \\« < £ \\sk - sk-i \\« <
k=m+i
< (5m + Sm+i + ... + 5n-i)l\si - so\\« <
< 5m(1 + 5 + 52 + • • • + 5n-m-i)|si - so\U <
< S"' {) ii* u«.
Since 0 < 5 < 1, we have (1 - 5n-m) < 1, then
5m
\sn - sm\« < --Ayi
1 - 0
As |yi(x)| < to, then \\sn - sm\\« —> 0 when m —> to.
We conclude that sn is a Cauchy sequence in C[0,1], so y = lim yn. Then, the series converges and the proof is completed. □
_ n ■
6. Illustrative Example. In this section, we present the analytical techniques based on the Adomian decomposition method and the modified Laplace Adomian decomposition method to solve Caputo fractional Volterra - Fredholm integro-differential equation.
Example. Consider the following Caputo fractional Volterra - Fredholm integro-differential equation.
t i 6t2.25 t2et f f
cD075[y(t)] = - -5-y(t) + J etsy(s)ds + J (4 - s-3)y(s)ds, (13)
oo
with the initial condition y(0) =0 and the exact solution y(t) = t3. Firstly, we apply the Adomian decomposition method.
Applying the operator J075 to both sides of (13),
y(t) = E ^ F + J0 75 [t225] - 5J0 75 tety(t)] +
dtk k! r(3.25)
+ J
0.75
etsy(s)ds + I (4 - s 3)y(s)ds
t
i
Then, yo(t) =
TO— 1
£ dy(°) tk + 6 j0.7512.25 =
k=0
dtk k! r(3.25)
6
0 + ___r(9/4+1) t(9/4+3/4) = t3
r(3.25)r(9/4 +3/4 + 1) ' ''
y1(t) = J
0.75
-t2 etyo(t) 5
= - 5 J075 \t2etyo (t)] + J0 = -5 J0.75 [t2etyo(t)] + J0
+ J e sy0(s)ds + J (4 - s à)y0(s)ds
00 t 1
+ /(4 - s"3 )s3ds
00
1 et t5 + 0 5
-5 J0.75 [t2ety0(t)] + 1 J0.75 [ett2y0(t)] = 0,
yn(t) = 0.
Therefore, the obtained solution is
y(t) = t3.
Secondly, the modified Laplace Adomian decomposition method.
We apply the Laplace transform to both sides of (13).
L [cD°.75y(t)] = L
y(t)
+L
6t
2.25
r(3.25)
+
+
e sy(s)ds + / (4 - s 3)y(s)ds
Using the property of Laplace transform and the initial condition, we get
s4 L [y(t)] = L
y(t)
+
6t
2.25
r(3.25)
+
1
t
1
t
+ C
t 1
Je'+ /(4 - s"3
0 0
C
y(t)
= — L
S 4
y(t)
+ 4t L
S 4
6t
2.25
_r(3.25)_
+ L
S 4
t 1
i / (A „-3n
J e sy(s)ds + J (4 - s 3)y(s)ds 00
Substituting (5) and (6) into the above equation, we obtain - - - - 6t2.25
C
J^yn(t)
Ln=0
1
s 4
= — L
r(3.25)
+ 4 L
s 4
(-£
v / n
yn(t)
+
+ L
s 4
/(X> „ —
e's £ Ands + / (4 - s"3) £ Bnds
n n=0 n n=0
By matching both sides of above equation, we have the following relation:
L [y0 (t)] L [yi (t)]
4- L s 4 r 6t2.25 1
r(3.25)
1L s 4 r-t2 e'y0 (t)
5
t i + J e'sA0ds + J(4 - s"3)B0ds
00
L [yn+i(t)] =
1
_3
s 4
—L
-t2e yn(t) 5
+ / e'sAnds + / (4 - s 3)Bnds
Applying inverse Laplace transform to above equations, we get yo (t) = t3,
yi (t) = L
-i
4 L
s 4
-t e5y0(t) + / e's4ds + /(4 - s"3)s3ds
= 0
i
t
t
i
yn (t) = 0.
Therefore, the obtained solution is y(t) = t3.
7. Conclusion. In this paper, Adomian decomposition method and modified Laplace Adomian decomposition method have been successfully applied to find the approximate solution of the fractional Volterra-Fredholm integro-differential equation. The reliability of the methods and reduction in the size of the computational work give these methods a wider applicability. The methods are very powerful and efficient in finding analytical as well as numerical solutions for wide classes of linear and nonlinear fractional Volterra - Fredholm integro-differential equations. They provide more realistic series solutions that converge very rapidly in real physical problems. Finally, the behavior of the solution can be formally determined by analytical approximate. The proposed methods can be applied to other nonlinear fractional differential equations, systems of differential and integral equations.
References
[1] Abbaoui K., Cherruault Y. Convergence of Adomian's method applied to nonlinear equations. Math. Comput. Model., 1994, vol. 20, no. 9, pp. 6973. DOI: 10.1016/0895-7177(94)00163-4.
[2] Adomian G. A review of the decomposition method in applied mathematics. J. Math. Anal. Appl., 1988, vol. 135, no. 2, pp. 501-544. DOI: 10.1016/0022-247X(88)90170-9.
[3] AL-Smadi M., Gumah G. On the homotopy analysis method for fractional SEIR epidemic model. Res. J. Appl. Sci. Eng. Technol., 2014, vol. 18, no. 7, pp. 3809-3820.
[4] Alhendi F., Shammakh W., Al-Badrani H. Numerical .solutions for quadratic integro-differential equations of fractional orders. Open J. Appl. Sci., 2017, vol. 7, pp. 157-170. DOI: 10.4236/ojapps.2017.74014.
[5] Araghi M., Behzadi S. Solving nonlinear Volterra - Fredholm integro-differential equations using the modified, Adomian decomposition method. Comput. Methods Appl. Math., 2009, vol. 9, no. 4, pp. 321-331. DOI: 10.2478/cmam-2009-0020.
[6] Daftardar-Gejji V., Jafari H. Solving a multi-order fractional differential equation using Adomian decomposition. Appl. Math. Comput., 2007, vol. 189, no. 1, pp. 541-548. DOI: 10.1016/j.amc.2006.11.129.
[7] Ghadle K., Hamoud A. Study of the approximate solution of fuzzy Volterra - Fredholm integral equations by using (ADM), Elixir Appl. Math., 2016, vol. 98, pp. 42567-42573.
[8] Hamoud A., Ghadle K. The reliable modified of Laplace Adomian decomposition method to solve nonlinear interval Volterra - Fredholm integral equations. Korean J. Math., 2017, vol. 25, no. 3, pp. 323-334. DOI: 10.11568/kjm.2017.25.3.323.
[9] Hamoud A., Ghadle K. The combined modified Laplace with Adomian decomposition method for solving the nonlinear Volterra - Fredholm integro-differential equations. J. Korean Soc. Ind. Appl. Math., 2017, vol. 21, no. 1, pp. 17-28. DOI: 10.12941/jksiam.2017.21.017.
[10] Hamoud A., Ghadle K. Modified Adomian decomposition method for -solving fuzzy Volterra - Fredholm integral equations. J. Indian Math. Soc., 2018, vol. 85, no. (1-2), pp. 01-17. DOI: 10.18311/jims/2018/16260.
[11] Hamoud A., Ghadle K. A study of some reliable methods for solving fuzzy Volterra - Fredholm integral equations. Acta Universitatis Apulensis, 2018, vol. 53, pp. 65-92. DOI: 10.17114/j.aua.2018.53.06.
[12] Hamoud A., Ghadle K. Recent advances on reliable methods for solving Volterra - Fredholm integral and integro-differential equations. Asian Journal of Mathematics and Computer Research, 2018, vol. 24, pp. 128-157.
[13] Hamoud A., Ghadle K. Existence and uniqueness results for fractional Volterra - Fredholm integro-differential equations. Int. J. Open Problems Compt. Math., 2018, vol. 11, no. 3, pp. 16-30. To Appear.
[14] Jafarian A., Rostami F., Golmankhaneh A., Baleanu D. Using ANNs approach for solving fractional order Volterra integro-differential equations. Int. J. Comput. Intell. Syst., 2017, vol. 10, pp. 470-480. DOI: 10.2991/ij-cis.2017.10.1.32.
[15] Kadem A., Kilicman A. The approximate solution of fractional Fredholm integro-differential equations by variational iteration a,nd, homotopy perturbation methods. Abstr. Appl. Anal., 2012, vol. 2012, pp. 1-10. DOI: 10.1155/2012/486193.
[16] Kilbas A., Srivastava H., Trujillo J. Theory and applications of fractional differential equations. North-Holland Math. Stud., 2006.
[17] Ma X., Huang C. Numerical solution of fractional integro-differential equations by a hybrid collocation method. Appl. Math. Comput., 2013, vol. 219, no. 12, pp. 6750-6760. DOI: 10.1016/j.amc.2012.12.072.
[18] Mittal R., Nigam R. Solution of fractional integro-differential equations by Adomian decomposition method. Int. J. Appl. Math. Mech., 2008, vol. 4, no. 2, pp. 87-94.
[19] Momani S., Noor M. Numerical methods for four-order fractional integro-differential equations. Appl. Math. Comput., 2006, vol. 182, pp. 754-760. DOI: 10.1016/j.amc.2006.04.041.
[20] Momani S., Qaralleh R. An efficient method for solving systems of fractional integro-differential equations. Comput. Math. Appl., 2006, vol. 52, pp. 459-470. DOI: 10.1016/j.camwa.2006.02.011.
[21] Yang C. Numerical solution of nonlinear Fredholm integro-differential equations of fractional order by using hybrid of block-pulse functions and Chebyshev polynomials. Math. Probl. Eng., 2011, vol. 2011, pp. 1-11. DOI: 10.1155/2011/341989.
[22] Yang C., Hou J. Numerical solution of integro-differential equations of fractional order by Laplace decomposition method. Wseas Trans. Math., 2013, vol. 12, no. 12, pp. 1173-1183.
Received January 19, 2018.
In revised form, May 12, 2018.
Accepted May 14, 2018.
Published online May 25, 2018.
Ahmed A. Hamoud
Department of Mathematics,
Dr. Babasaheb Ambedkar Marathwada University,
Aurangabad, 431004, India;
Department of Mathematics,
Taiz University, Taiz, Yemen.
E-mail: drahmedselwi985@gmail.com
Kirtiwant P. Ghadle Department of Mathematics,
Dr. Babasaheb Ambedkar Marathwada University, Aurangabad, 431004, India. E-mail: drkp.ghadle@gmail.com
