Next Article in Journal
A New Approach to the Solution of Non-Linear Integral Equations via Various FBe-Contractions
Previous Article in Journal
A New Flexible Sigmoidal Growth Model
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Analytical Numerical Method for Solving Fuzzy Fractional Volterra Integro-Differential Equations

by
Mohammad Alaroud
1,
Mohammed Al-Smadi
2,
Rokiah Rozita Ahmad
1,* and
Ummul Khair Salma Din
1
1
Center for Modelling and Data Science, Faculty of Science and Technology, Universiti Kebangsaan Malaysia, 43600 Bangi, Selangor DE, Malaysia
2
Department of Applied Science, Ajloun College, Al-Balqa Applied University, Ajloun 26816, Jordan
*
Author to whom correspondence should be addressed.
Submission received: 9 January 2019 / Revised: 26 January 2019 / Accepted: 30 January 2019 / Published: 12 February 2019

Abstract

:
The modeling of fuzzy fractional integro-differential equations is a very significant matter in engineering and applied sciences. This paper presents a novel treatment algorithm based on utilizing the fractional residual power series (FRPS) method to study and interpret the approximated solutions for a class of fuzzy fractional Volterra integro-differential equations of order 0 < β 1 which are subject to appropriate symmetric triangular fuzzy conditions under strongly generalized differentiability. The proposed algorithm relies upon the residual error concept and on the formula of generalized Taylor. The FRPS algorithm provides approximated solutions in parametric form with rapidly convergent fractional power series without linearization, limitation on the problem’s nature, and sort of classification or perturbation. The fuzzy fractional derivatives are described via the Caputo fuzzy H -differentiable. The ability, effectiveness, and simplicity of the proposed technique are demonstrated by testing two applications. Graphical and numerical results reveal the symmetry between the lower and upper r -cut representations of the fuzzy solution and satisfy the convex symmetric triangular fuzzy number. Notably, the symmetric fuzzy solutions on a focus of their core and support refer to a sense of proportion, harmony, and balance. The obtained results reveal that the FRPS scheme is simple, straightforward, accurate and convenient to solve different forms of fuzzy fractional differential equations.

1. Introduction

Fuzzy theory of fractional differential and integro-differential equations is a new and important branch of fuzzy mathematics. It has ample applications due to the fact that many practical problems in industrial engineering, computer science, physics, artificial intelligence, and operations research may be converted to uncertain process problems of fractional order. The topic of fuzzy fractional integro-differential equations (FFIDEs) has gained the attention of researchers in recent times because it is considered a powerful tool by which to present vague parameters and to handle with their dynamical systems in natural fuzzy environments. Indeed, it has a great significance in the fuzzy analysis theory and its applications in fuzzy control models, artificial intelligence, quantum optics, measure theory, and atmosphere, etc. [1,2,3,4]. In several cases, information about the real-life problems included is always pervaded with uncertainty. This uncertainty results from several factors, such as measurement errors, deficient data or if the constraint conditions are determined. So, it is necessary to have some mathematical tools to understand this uncertainty. Consequently, to form a convenient and applicable algorithm is important to achieve a mathematical structure that would suitably process FFIDEs and solve them.
Most FFIDE problems cannot be solved analytically, and hence, finding good approximate solutions using numerical methods will be very valuable. Recently, numerous scholars have devoted their interest to studying and investigating solutions to FFIDEs utilizing different numerical and analytical techniques; these solutions include the Fuzzy Laplace transforms technique [5], two-dimensional Legendre wavelet technique [6], Adomian decomposition technique [7], variational iteration technique [8], and the fitted reproducing kernel Hilbert space technique [9]. Besides that, other scholars have shown an interest in the existence and uniqueness of the solutions for FDEs. (For more details, see [10,11,12,13,14].) Applications of the RPS scheme are extensive and numerous, especially for the simulation of models of integer order. Such approaches have recently been developed as a powerful numerical tool for treating different problems of arbitrary order. From a methodological point of view, the Volterra integro-differential equations considered in this study under uncertainty concerned with fractional derivatives are simply generalizations of classical forms, provided that β = 1 , which have various uses in physics and engineering, including in dynamical systems, nonlinear propagation of a traveling wave, damped nonlinear string, electronics, and telecommunications, etc. For more details about the integro-differential equations, refer to [15,16,17,18].
The present work aims to expand the applications of the fractional residual power series (FRPS) technique to determine approximated solutions for a class of fuzzy fractional Volterra integro-differential equations (FFVIDEs) subject to certain fuzzy initial conditions in the following form:
D a + β y ( x ) + a x φ ( x , t ) y ( t ) d t = g ( x ) ,   x a ,   0 < β 1 ,
with the fuzzy initial conditions
y ( a ) = y 0
where λ is a parameter, g : [ a , b ] F is a continuous fuzzy-valued function,   φ ( x , t ) is a continuous crisp arbitrary kernel function, D a + β is the fuzzy fractional derivative in Caputo sense, y 0 F and y ( x ) is unknown analytical fuzzy function to be determined. Throughout this paper, F stands for the set of all fuzzy numbers on the real line .
The FRPS technique developed in [19] is considered an easy and applicable tool to create a power series solution for strongly linear and nonlinear equations without being linearized, discretized, or exposed to perturbation [20,21,22,23,24,25,26,27]. This technique is featured by the following characteristics: firstly, the method provides the solutions in Taylor expansions; therefore, the exact solutions will be available when the solutions are polynomials. Secondly, the solutions along with their derivatives can be applied for each arbitrary point in the given interval. Thirdly, the method does not require modifications while converting from the lower to the higher order. Consequently, this method can be easily and directly applied to the given system by selecting an appropriate value for the initial guesses/approximations. Fourthly, the FRPS technique needs minor computational requirements with less time and more accuracy.
The remainder of the current paper is structured as follows: In the next section, some essential notations, definitions, and results relating to fuzzy calculus and fuzzy fractional calculus are shown. The process to solve fuzzy FVIDEs is discussed in Section 3. In Section 4, the construction of solutions using the FRPS algorithm is presented. Numerical examples will be performed to show the capability, potentiality and simplicity of the method. The last section of the current paper is a conclusion.

2. Overview of Fuzzy Calculus and Fuzzy Fractional Calculus

In the current section, the essential notations, definitions and the basic results relating to fuzzy calculus and fuzzy fractional calculus are presented. In general, a fuzzy number ρ is a fuzzy subset of with normal, convex, and upper semi-continuous membership function of bounded support. For more details, refer to [13,28,29,30,31,32,33,34].
Definition 1.
[28]A fuzzy number ρ is a mapping such that ρ : [ 0 , 1 ] with the following properties:
  • ρ is fuzzy convex, that is, ρ ( λ s + ( 1 λ ) t ) m i n { ρ ( s ) , ρ ( t ) } for all s , t , λ [ 0 , 1 ] .
  • ρ is normal, that is, s * for which ρ ( s * ) = 1 .
  • ρ is upper-semi continuous, that is, ρ ( s * ) l i m s s * + ρ ( s ) for any s * .
  • s u p p ( ρ ) = { s : ρ ( s ) > 0 } is the support of ρ , and { s : ρ ( s ) > 0 } ¯ is compact, where { } ¯ denotes the closure of a subset.
The parametric form or the r -cut representation of a fuzzy number ρ F , is defined as: [ ρ ] r = { s : ρ ( s ) r } , if r ( 0 , 1 ] , and [ ρ ] r = s u p p ( ρ ) ¯ , if r = 0 . Clearly, the parametric form of ρ is closed and bounded interval [ ρ 1 r , ρ 2 r ] in which ρ 1 r is the lower r -cut representation of ρ , and ρ 2 r is the upper r -cut representation of ρ . An equivalent parametric definition is also given in [13] as follows:
Definition 2.
[28]A fuzzy number ρ in parametric form is a pair ( ρ 1 ,   ρ 2 ) of functions ρ 1 ( r ) ,   ρ 2 ( r ) , for r [ 0 , 1 ] , which satisfy the following requirements:
  • ρ 1 ( r ) is a bounded non-decreasing left continuous for each r ( 0 , 1 ] , and right continuous at r = 0 .
  • ρ 2 ( r ) is a bounded non-increasing left continuous for each r ( 0 , 1 ] , and right continuous at r = 0 .
  • ρ 1 ( r ) ρ 2 ( r ) , for each r [ 0 , 1 ] .
The metric structure of F is defined by the Hausdorff distance mapping D H : F × F + { 0 } such that D H ( ρ , ω ) = sup 0   r 1 max { | ρ 1 r ω 1 r | , | ρ 2 r ω 2 r | } , for arbitrary fuzzy numbers ρ and ω . The metric ( F ,   D H ) has been proved in [30] as complete metric space.
The next definition explains the concept of strongly generalized differentiability, which was introduced by Bede and Gal in [31], to obtain solutions which have a decreasing length of their support, and hence, a decreasing level of uncertainty.
Definition 3.
Let y : ( a , b ) F and for fixed x 0 ( a , b ) . y is called a strongly generalized differentiable at x 0 , if there is an element y ( x 0 ) F such that either:
  • The H-differences y ( x 0 + ξ ) y ( x 0 ) ,   y ( x 0 ) y ( x 0 ξ ) exist, for each ξ > 0 sufficiently tends to 0 and lim ξ 0 + y ( x 0 + ξ   ) y ( x 0 ) ξ = y ( x 0 ) = lim ξ 0 + y ( x 0 ) y ( x 0 ξ ) ξ ,
  • The H-differences y ( x 0 ) y ( x 0 + ξ ) , y ( x 0 ξ ) y ( x 0 ) exist, for each ξ > 0 sufficiently tends to 0 and   lim ξ 0 + y ( x 0 ) y ( x 0 + ξ ) ξ = y ( x 0 ) = lim ξ 0 + y ( x 0 ξ ) y ( x 0 ) ξ .
where the limit here is taken in the complete metric space ( F ,   D H ) .
Definition 4.
Let y : [ a , b ] F .Then the function y is continuous at x 0 [ a , b ] if for every ϵ > 0 ,   δ = δ ( x 0 , ϵ ) > 0 such that   D H ( y ( x ) , y ( x 0 ) ) < ϵ , for each x [ a , b ] , whenever   | x x 0 | < δ .
The concept of fuzzy integral functions was first presented in 1982 by Dubois and Prade. Since then, many approaches have been proposed, such as Riemann integral approach by Goetschel and Voxman [29] and Lebesgue-type approach by Kaleva [28]. In this article, we will use the Riemann integral approach, which is characterized by dealing with fuzzy-valued function integration.
Definition 5.
[29] Let y : F be a continuous fuzzy-valued function. For each partition = { x 0 * , x 1 * , , x n * } of and for arbitrary point η i [ x i 1 * , x i * ] , 1 i n . Let Ψ = i = 1 n y i ( η i ) ( x i * x i 1 * ) and Δ = max 1 i n | x i * x i 1 * | . Then, the definite integral of y ( x ) over is defined by y ( x ) d x = lim Δ 0 Ψ provided that the limit exists in the metric space ( F , D H ) .
The theorems below were proposed and studied by Chalco-Cano and Roman-Flores in [32] to assist us to convert the fuzzy fractional differential equations (FFDEs) into a system of ordinary fractional differential equations (OFDEs), ignoring the fuzzy setting approach.
Theorem 1.
Suppose that y : [ a , b ] F , where [ y ( x ) ] r = [ y 1 r ( x ) ,   y 2 r ( x ) ] , r [ 0 , 1 ] ,
  • If y is (1)-differentiable on [ a , b ] , then y 1 r and y 2 r are two differentiable functions on [ a , b ] , and [ D 1 1 y ( x ) ] r = [ y 1 r ( x ) ,   y 2 r ( x ) ] ,
  • If y is (2)-differentiable on [ a , b ] , then y 1 r and y 2 r are two differentiable functions on [ a , b ] , and [ D 2 1 y ( x ) ] r = [ y 2 r ( x ) ,   y 1 r ( x ) ] .
Theorem 2.
Let y : [ a , b ] F be continuous fuzzy–valued function and put [ y ( x ) ] r = [ y 1 r ( x ) ,   y 2 r ( x ) ] for each r [ 0 , 1 ] , then a b y ( x ) d x exist, belong to F , y 1 r ( x )   and y 2 r ( x ) are integrable functions on [ a , b ] and [ a b y ( x ) d x ] r = [ a b y 1 r ( x ) d x ,   a b y 2 r ( x ) d x ] , r [ 0 , 1 ] .
Definition 6.
Let y : [ a , b ]   F and y C F [ a , b ] L F [ a , b ] . One can say y is Caputo fuzzy H -differentiable at x when D a + β y ( x ) = 1 Γ ( 1 β ) a x y ( ζ ) ( x ζ ) β   d ζ   exists, where 0 < β 1 . As well, we say that y is Caputo [ ( 1 ) β ] differentiable when y is (1)-differentiable and y is Caputo [ ( 2 ) β ] differentiable when y is (2)-differentiable, where C F [ a , b ] and L F [ a , b ] stand for the space of all continuous and Lebesque integrable fuzzy-valued functions on [ a , b ] , respectively.
Theorem 3.
Let 0 < β 1 and y C F [ a , b ] . Then, for each r [ 0 , 1 ] , the Caputo fuzzy fractional derivative exists on ( a , b ) such that
[ D a + β y ( x ) ] r = [ 1 Γ ( 1 β ) a x y 1 r ( ζ ) ( x ζ ) β   d ζ   ,   1 Γ ( 1 β ) a x y 2 r ( ζ ) ( x ζ ) β   d ζ ]   for   ( 1 ) - differentiable , [ D a + β y ( x ) ] r = [ 1 Γ ( 1 β ) a x y 2 r ( ζ ) ( x ζ ) β   d ζ   ,   1 Γ ( 1 β ) a x y 1 r ( ζ ) ( x ζ ) β   d ζ ]   for   ( 2 ) - differentiable .

3. Formulation of Fuzzy Fractional Volterra IDEs

In the current section, the fuzzy fractional Volterra integro-differential equations (FFVIDEs) will be thoroughly studied using three main concepts, namely: strongly generalized differentiability, Caputo’s H-differentiability and Riemann integrability, where the FFVIDEs are converted into corresponding system of crisp ordinary fractional Volterra integro-differential equations (OFVIDEs) for every differentiability type. The (FFVIDEs) will be appropriately solved if the initial value is a fuzzy number and the solution is a fuzzy function, and then the integral, as well as the derivative, should be individually considered as a fuzzy integral and fuzzy derivative. The formulation of the problem is considered the most significant part of this process. It is the determination process of the parametric form of g , the integrability type’s selection, the differentiability type’s selection, as well as the separation kernel of ϕ . Then, FFVIDEs (1) and (2) are framed as a set of crisp OFVIDEs; as a result, a new form of discretization will be used. Besides that, we assume that without loss of generality   ϕ ( x , t , y ( t ) ) = φ ( x , t ) y ( t ) , so we can apply the FRPS technique to solve FFVIDEs (1) and (2).
Now, for all r ( 0 , 1 ] and a x b , let the parametric form of the fuzzy functions g ( x ) and y ( x ) , as [ g ( x ) ] r = [ g 1 r ( x ) , g 2 r ( x ) ] and [ y ( x ) ] r = [ y 1 r ( x ) , y 2 r ( x ) ] , respectively, as well as [ y ( a ) ] r = [ y 1 r ( a ) , y 2 r ( a ) ] = [ y 0 , 1 r , y 0 , 2 r ] . Hence, one can write the FFVIDEs (1) and (2) in the r -cut representation as follows:
[ D a + β y ( x ) ] r + λ a x [ ϕ ( x , t , y ( t ) ) ] r d t = [ g ( x ) ] r ,
where the r -cut representation of the function, [ ϕ ( x , t , y ( t ) ) ] r is given by:   [ ϕ ( x , t , y ( t ) ) ] r = [ ϕ 1 r ( x , t , y 1 r ( t ) , y 2 r ( t ) ) , ϕ 2 r ( x , t , y 1 r ( t ) , y 2 r ( t ) ) ] , in which
ϕ 1 r ( x , t , y 1 r ( t ) , y 2 r ( t ) ) = { φ ( x , t )   y 1 r ( t ) ,   φ ( x , t ) 0 φ ( x , t )   y 2 r ( t ) ,   φ ( x , t ) < 0 ϕ 2 r ( x , t , y 1 r ( t ) , y 2 r ( t ) ) = { φ ( x , t )   y 2 r ( t ) ,   φ ( x , t ) 0 φ ( x , t )   y 1 r ( t ) ,   φ ( x , t ) < 0
The ( n ) -solution of fuzzy FVIDEs (1) and (2) is a function y : [ a , b ] F that has Caputo [ ( n ) β ]-differentiable and satisfies the FFVIDEs (1) and (2). To compute it, we perform the next algorithm.
Algorithm 1: To obtain the ( n ) -solution of the FFVIDEs (1.1), there are two cases that will be discussed as follows:
Case 1: If y ( x ) is Caputo [(1)-   β ]-differentiable, we convert the FFVIDEs (1) and (2) to the following OFVIDEs system:
D a + β y 1 r ( x ) + λ   a x φ ( x , t )   y 1 r ( t ) d t   = g 1 r ( x ) D a + β y 2 r ( x ) + λ a x φ ( x , t )   y 2 r ( t ) d t = g 2 r ( x ) ,
with initial conditions
y 1 r ( a ) = y 0 , 1 r   and   y 2 r ( a ) = y 0 , 2 r
Then, do the following steps:
  • Step 1: Solve the system (4) and (5) for y 1 r ( x ) and y 2 r ( x ) .
  • Step 2: Ensure that [ y 1 r ( x ) , y 2 r ( x ) ] and [ D a + β y 1 r ( x ) ,   D a + β y 2 r ( x ) ] are valid level sets on [ a , b ] or on a partial interval in   [ a , b ] .
  • Step 3: Construct a ( 1 ) -differentiable solution y ( x ) whose r -cut representation is [ y 1 r ( x ) , y 2 r ( x ) ] .
Case 2: If y ( x ) is Caputo [(2)-   β ]-differentiable, we convert the FFVIDEs (1) and (2) to the following OFVIDEs system:
D a + β y 1 r ( x ) + λ   a x φ ( x , t )   y 2 r ( t ) d t = g 2 r ( x ) , D a + β y 2 r ( x ) + λ   a x φ ( x , t )   y 1 r ( t ) d t = g 1 r ( x ) ,
with initial conditions
y 1 r ( a ) = y 0 , 1 r   and   y 2 r ( a ) = y 0 , 2 r .
Then, do the following steps:
  • Step 1: Solve the system (6) and (7) for y 1 r ( x )   and y 2 r ( x ) .
  • Step 2: Ensure that [ y 1 r ( x ) , y 2 r ( x ) ] and [ D a + β y 2 r ( x ) ,   D a + β y 1 r ( x ) ] are valid level sets on [ a , b ] or on a partial interval in   [ a , b ] .
  • Step 3: Construct a ( 2 ) -differentiable solution y ( x ) whose r -cut representation is [ y 1 r ( x ) , y 2 r ( x ) ] .
The previous formulation of FFVIDEs (1) and (2) along with Theorem 1 shows us how to deal with the numerical solutions of FFVIDEs. The original FFVIDEs can be converted into a system of OFVIDEs equivalently. As a consequence of this, we can employ the numerical approaches directly to solve the system OFVIDEs obtained without the need for reformulation in the fuzzy setting.

4. Description of the FRPS Technique

In the present section, we show the procedure of FRPS algorithm in order to study and construct analytic-numeric approximated solutions for FFVIDEs (1) and (2) through substituting the expansion of its fractional power series (FPS) among its truncated residual functions. In view of that, the resultant equation helps us derive a recursion formula for the coefficients’ computation, where the FPS expansion’s coefficients can be computed recursively through recurrent fractional differentiating of the truncated residual function.
Definition 7.
[22] A fractional power series (FPS) representation at x 0 has the following form
n = 0 c n ( x x 0 ) n β = c 0 + c 1 ( x x 0 ) β + c 2 ( x x 0 ) 2 β + ,
where 0 m 1 < β m , x x 0 and c n ’s are the coefficients of the series.
Theorem 4.
[22] Suppose that f has the following FPS representation at x 0
f ( x ) = n = 0 c n ( x x 0 ) n β .
where f ( x ) C [ x 0 , x 0 + R ) and D x 0 n β f ( x ) C ( x 0 , x 0 + R ) for n = 0 ,   1 ,   2 , , then the coefficients c n will be in the form c n = D x 0 n β f ( x 0 ) Γ ( 1 + n β )   such that D x 0 n β = D x 0 β · D x 0 β ·     · D x 0 β ( n -times).
Conveniently, to obtain the ( n ) -solution of FFVIDEs (1) and (2), we will explain the fashion to determine (1)-solution equivalent to the solution for the system of OFVIDEs (4) and (5). Further, the same manner can be applied to construct (2)-solution. To reach our purpose, assume that the approximated solutions of OFVIDEs (4) and (5) at x 0 = 0 have the following form:
y 1 r ( x ) = n = 0 c n x n β Γ ( n β + 1 ) , y 2 r ( x ) = n = 0 d n x n β Γ ( n β + 1 ) .
Since y 1 r ( x ) and y 2 r ( x ) satisfy the initial condition (5), then y 1 r ( 0 ) = y 0 , 1 r = c 0 and y 2 r ( 0 ) = y 0 , 2 r = d 0 will be the initial guesses for (5), so the series solutions can be written as
y 1 r ( x ) = c 0 + n = 1 c n x n β Γ ( n β + 1 ) , y 2 r ( x ) = d 0 + n = 1 d n x n β Γ ( n β + 1 ) .
Further, we can approximate y 1 r ( x ) and y 2 r ( x ) by the following k th-truncated series solutions:
y k , 1 r ( x ) = c 0 + n = 1 k c n x n β Γ ( n β + 1 ) , y k , 2 r ( x ) = d 0 + n = 1 k d n x n β Γ ( n β + 1 ) .
Define the so-called k th-residual functions R e s k , 1 r and R e s k , 2 r , for k = 1 ,   2 ,   3 , ,
R e s k , 1 r ( x ) = D 0 + β y k , 1 r ( x ) + λ 0 x φ ( x , t ) y k , 1 r ( t ) d t g 1 r ( x ) , R e s k , 2 r ( x ) = D 0 + β y k , 2 r ( x ) + λ 0 x φ ( x , t ) y k , 2 r ( t ) d t g 2 r ( x ) .
and the following residual functions
R e s 1 r ( x ) = lim k R e s k , 1 r ( x ) = D 0 + β y 1 r ( x ) + λ 0 x φ ( x , t ) y 1 r ( t ) d t g 1 r ( x ) , R e s 2 r ( x ) = lim k R e s k , 2 r ( x ) = D 0 + β y 2 r ( x ) + λ 0 x φ ( x , t ) y 2 r ( t ) d t g 2 r ( x ) .
Clear that, R e s n r ( x ) = 0 , for n = 1 , 2 and for each x 0 ,which leads to D 0 + ( m 1 ) β R e s n r ( x ) = 0 . Also, we have D 0 + ( m 1 ) β R e s n r ( 0 ) = D 0 + ( m 1 ) β R e s k , n r ( 0 ) = 0 , for each m = 0 ,   1 ,   2 , , k . However, D 0 + ( k 1 ) β R e s k , n r ( 0 ) = 0 holds for n = 1 ,   2 .
According to applying FRPS technique to find the 1st-unknown coefficients c 1 and d 1 , substitute the 1st-approximated solutions y 1 , 1 r ( x ) = c 0 + c 1 x β Γ ( β + 1 ) and y 1 , 2 r ( x ) = d 0 + d 1 x β Γ ( β + 1 ) into the 1st residual functions R e s 1 , 1 r ( x ) and R e s 1 , 2 r ( x ) of (12) such that
R e s 1 , 1 r ( x ) = D 0 +   β y 1 , 1 r ( x ) + λ 0 x φ ( x , t ) y 1 , 1 r ( t ) d t g 1 r ( x ) ,   = c 1 + λ a x φ ( x , t ) ( c 0 + c 1 t β Γ ( β + 1 ) ) d t g 1 r ( x ) , R e s 1 , 2 r ( x ) = D 0 +   β y 1 , 2 r ( x ) + λ 0 x φ ( x , t ) y 1 , 2 r ( t ) d t g 2 r ( x ) , = d 1 + λ a x φ ( x , t ) ( d 0 + d 1 t β Γ ( β + 1 ) ) d t g 2 r ( x ) .
Using the facts R e s 1 , 1 r ( 0 ) = R e s 1 , 2 r ( 0 ) = 0 , yields
c 1 = ( g 1 r ( x ) λ c 0 0 x φ ( x , t ) d t 1 + λ Γ ( β + 1 ) 0 x t β φ ( x , t ) d t ) x = 0 and   d 1 = ( g 2 r ( x ) λ d 0 0 x φ ( x , t ) d t 1 + λ Γ ( β + 1 ) 0 x t β φ ( x , t ) d t ) x = 0 .
Again, to determine the 2nd-unknown coefficients c 2 and d 2 , substitute the 2nd-approximated solutions y 2 , 1 r ( x ) = c 0 + c 1 x β Γ ( β + 1 ) + c 2 x 2 β Γ ( 2 β + 1 ) and y 2 , 2 r ( x ) = d 0 + d 1 x β Γ ( β + 1 ) + d 2 x 2 β Γ ( 2 β + 1 ) into the 2nd-residual functions R e s 2 , 1 r ( x ) and R e s 2 , 2 r ( x ) of (11) such that
R e s 2 , 1 r ( x ) = D 0 +   β y 2 , 1 r ( x ) + λ 0 x φ ( x , t ) y 2 , 1 r ( t ) d t g 1 r ( x ) , R e s 2 , 2 r ( x ) = D 0 + β y 2 , 2 r ( x ) + λ 0 x φ ( x , t ) y 2 , 2 r ( t ) d t g 2 r ( x ) .
Then, by applying the fractional derivative D 0 +   β on both sides of R e s 2 , 1 r ( x ) and R e s 2 , 2 r ( x ) , as well using the facts D 0 +   β R e s 2 , 1 r ( 0 ) = D 0 +   β R e s 2 , 2 r ( 0 ) = 0 , the values of   c 2   and d 2   will be given by
c 2 = ( D 0 +   β g 1 r ( x ) λ 0 x ( c 0 + c 1 t β Γ ( β + 1 ) )   D 0 +   β φ ( x , t ) d t 1 + λ Γ ( 2 β + 1 ) 0 x t 2 β   D 0 +   β φ ( x , t ) d t ) x = 0 , d 2 = ( D 0 +   β g 2 r ( x ) λ 0 x ( d 0 + d 1 t β Γ ( β + 1 ) )   D 0 +   β φ ( x , t ) d t 1 + λ Γ ( 2 β + 1 ) 0 x t 2 β   D 0 +   β φ ( x , t ) d t ) x = 0 .
In the same way for the 3rd-unknown coefficients, c 3 and d 3 , substitute the 3rd-approximated solutions y 3 , 1 r ( x ) = c 0 + n = 1 3 c n x n β Γ ( n β + 1 ) and y 3 , 2 r ( x ) = d 0 + n = 1 3 d n x n β Γ ( n β + 1 ) into the 3rd-residual functions R e s 3 , 1 r ( x ) and R e s 3 , 2 r ( x ) of (11) and then by computing D 0 +   2 β R e s 3 , 1 r ( x ) and D 0 +   2 β R e s 3 , 2 r ( x ) and using the facts D 0 +   2 β R e s 3 , 1 r ( 0 ) = D 0 +   2 β R e s 3 , 2 r ( 0 ) = 0 , the coefficients, c 3 and d 3 will be given such that
c 3 = ( D 0 +   2 β g 1 r ( x ) λ 0 x ( c 0 + c 1 t β Γ ( β + 1 ) + c 2 t 2 β Γ ( 2 β + 1 ) )   D 0 +   2 β φ ( x , t ) d t 1 + λ Γ ( 3 β + 1 ) 0 x t 3 β   D 0 +   2 β φ ( x , t ) d t ) x = 0 , d 3 = ( D 0 +   2 β g 2 r ( x ) λ 0 x   ( d 0 + d 1 t β Γ ( β + 1 ) + d 2 t 2 β Γ ( 2 β + 1 ) )   D 0 +   2 β φ ( x , t ) d t 1 + λ Γ ( 3 β + 1 ) 0 x t 3 β   D 0 +   2 β φ ( x , t ) d t ) x = 0 .
The process can be repeated, as can the coefficients, until an arbitrary order of the FPS solution for OFVIDEs (4) and (5) is obtained.

5. Applications and Simulations

The dynamic behavior for fractional systems has gained a lot of prominence over the past few years, mainly because fractional operators have become excellent mathematical tools in modeling many real-life problems, helping us understand the mathematical structure and memory effect. The solution of fractional integro-differential equations, in the Volterra sense, is very important to describe the behavior of linear and non-linear physical systems, in particular, the dynamics of nuclear reactors, and systems which are harmonically excited, or to calculate the probabilistic response of randomly-excited analytical models and so forth. As a matter of terminology, the Volterra series refers to a functional expansion of a dynamic, nonlinear, and time-invariant functional. This section tests two FFIDEs of Volterra type in order to demonstrate the efficiency, accuracy, and applicability of the present novel approach. Here, all necessary calculations and analyses are done using Mathematica 10.
Example 1.
Consider the following FFIDE of Volterra type
D 0 + β y ( x ) = [ r + 1 ,   2 r ] ( 1 + x β ) + 1 Γ ( β ) 0 x ( x t ) β 1 y ( t ) d t   ,   0 < β 1 ,   x [ 0 , 1 ] ,
with the fuzzy initial condition
y ( 0 ) = 0 .
Based on the type of differentiability, fuzzy FVIDEs (18) and (19) can be converted to one of the following systems:
Case1: Under Caputo [(1)-   β ]-differentiability, the system of OFVIDEs corresponding to Caputo [(1)-   β ]-differentiable is
D 0 + β y 1 r ( x ) = ( r + 1 ) ( 1 + x β ) + 1 Γ ( β ) 0 x ( x t ) β 1 y 1 r ( t ) d t , D 0 + β y 2 r ( x ) = ( 2 r ) ( 1 + x β ) + 1 Γ ( β ) 0 x ( x t ) β 1 y 2 r ( t ) d t ,
with the initial conditions
y 1 r ( 0 ) = 0 ,   y 2 r ( 0 ) = 0 .
If β = 1 ,   then the exact solutions of OFVIDEs (20) and (21) are given by
[ y ( x ) ] r = [ r + 1 ,   2 r ] ( e x 1 ) .
Accordingly, the last description of FRPS algorithm, starting with y 0 , 1 r ( 0 ) = y 0 , 2 r ( 0 ) = 0 , then the k th-approximated solutions for the system (20) and (21) are given by
y k , 1 r ( x ) = n = 1 k c n x n β Γ ( n β + 1 ) , y k , 2 r ( x ) = n = 1 k d n x n β Γ ( n β + 1 ) .
Consequently, the k th-residual functions R e s k , 1 r , and R e s k , 2 r for (20) will be
R e s k , 1 r ( x ) = D 0 +   β y k , 1 r ( x ) 1 Γ ( β ) 0 x ( x t ) β 1 y k , 1 r ( t ) d t ( r + 1 ) ( 1 + x β ) , R e s k , 2 r ( x ) = D 0 +   β y k , 2 r ( x ) 1 Γ ( β ) 0 x ( x t ) β 1 y k , 2 r ( t ) d t ( 2 r ) ( 1 + x β ) .
The coefficients of the FRPS approximate solutions for the system (20) and (21) can be found in Appendix A.
To demonstrate the agreement between the lower and upper bound of the (1)-approximated solutions, the 7th-FRPS approximated solutions of Example 1, case 1 with different values of β and r = 0.5 , are shown in Table 1.
Case2: Under Caputo [(2)-   β ]-differentiability, the system of OFVIDEs corresponding to Caputo [(2)-   β ]-differentiable is
D 0 + β y 1 r ( x ) = ( r + 1 ) ( 1 + x β ) + 1 Γ ( β ) 0 x ( x t ) β 1 y 2 r ( t ) d t , D 0 + β y 2 r ( x ) = ( 2 r ) ( 1 + x β ) + 1 Γ ( β ) 0 x ( x t ) β 1 y 1 r ( t ) d t ,
with the initial conditions
y 1 r ( 0 ) = 0 ,   y 2 r ( 0 ) = 0 .
If β = 1 ,   then the exact solutions of OFVIDEs (25) and (26) are given by
[ y ( x ) ] r = 3 2 e x [ 2 r , r + 1 ] + 1 2 [ 1 2 r , 2 r 1 ] ( cos x sin x ) .
To apply the FRPS technique, suppose that the k th-truncated series solutions of y 1 r ( x ) and y 2 r ( x )   for the system (25) and (26) about x 0 = 0 have the following form:
y k , 1 r ( x ) = n = 1 k c n x n β Γ ( n β + 1 ) , y k , 2 r ( x ) = n = 1 k d n x n β Γ ( n β + 1 ) .
Apparently, the k th-residual functions R e s k , 1 r , and R e s k , 2 r for (25) will be
R e s k , 1 r ( x ) = D 0 +   β y k , 1 r ( x ) 1 Γ ( β ) 0 x ( x t ) β 1 y k , 2 r ( t ) d t ( r + 1 ) ( 1 + x β ) , R e s k , 2 r ( x ) = D 0 +   β y k , 2 r ( x ) 1 Γ ( β ) 0 x ( x t ) β 1 y k , 1 r ( t ) d t ( 2 r ) ( 1 + x β ) .
The coefficients of the FRPS approximate solutions for the system (25) and (26) can be found in Appendix A.
In Table 2, the 7th-FRPS approximated solutions of Example 1, case 2 with different values of β and r = 0.5 are computed.
Next, the absolute errors for Example 1, at β = 1 , and n = 10 , have been obtained in Table 3 and Table 4 with different values of r and step size 0.2 on [ 0 , 1 ] . Obviously, from these tables, it can be observed that the FRPS-approximated solutions are in high agreement with the exact solutions.
To show the fuzzy behavior of the proposed algorithm, the core and support of fuzzy ( m )-approximated solutions, m = 1 , 2 , and its Caputo derivative at β = 1 of Example 1, are plotted in Figure 1 and Figure 2, where c o r e ( y ( x ) ) = { τ : y ( x ) ( τ ) = 1 } , and s u p p ( y ( x ) ) = { τ : y ( x ) ( τ ) 0 } .
Example 2.
Consider the following FFIDE of Volterra type
D 0 + β y ( x ) = [ r 1 ,   1 r ] + 0 x y ( t ) d t   ,   0 < β 1 ,   x [ 0 , 1 ] ,
with the fuzzy initial condition
y ( 0 ) = 0.
Based on the type of differentiability, fuzzy FVIDEs (30) and (31) can be converted to one of the following systems:
Case 1: Under Caputo [(1)-   β ]-differentiability, the system of OFVIDEs corresponding to Caputo [(1)-   β ]-differentiable is
D 0 + β y 1 r ( x ) = ( r 1 ) + 0 x y 1 r ( t ) d t , D 0 + β y 2 r ( x ) = ( 1 r ) + 0 x y 2 r ( t ) d t ,
with the fuzzy initial condition
y 1 r ( 0 ) = 0 , y 2 r ( 0 ) = 0 .
If β = 1 ,   then the exact solution of OFVIDEs (32) and (33) is given by
[ y ( x ) ] r = [ r 1 ,   1 r ] s i n h ( x ) .
By utilizing the FRPS algorithm as presented previously, beginning with y 0 , 1 r = y 1 r ( 0 ) = 0 ,   y 0 , 2 r = y 2 r ( 0 ) = 0 , and depending on Equation (10), the k th-FRPS approximated solutions y k , 1 r ( x ) and y k , 2 r ( x ) for the system (32) and (33) are given by
y k , 1 r ( x ) = n = 1 k c n x n β Γ ( n β + 1 ) , y k , 2 r ( x ) = n = 1 k d n x n β Γ ( n β + 1 ) .
Consequently, the k th-residual functions R e s k , 1 r , and R e s k , 2 r for (32) will be
R e s k , 1 r ( x ) = D 0 +   β y k , 1 r ( x ) 0 x y k , 1 r ( t ) d t ( r 1 ) , R e s k , 2 r ( x ) = D 0 +   β y k , 2 r ( x ) 0 x y k , 2 r ( t ) d t ( 1 r ) .
The coefficients of the FRPS approximate solutions for the system (32) and (33) can be found in Appendix B.
Case 2: Under Caputo [(2)-   β ]-differentiability, the system of OFVIDEs corresponding to Caputo [(2)-   β ]-differentiable is
D 0 + β y 1 r ( x ) = ( 1 r ) + 0 x y 2 r ( t ) d t , D 0 + β y 2 r ( x ) = ( r 1 ) + 0 x y 1 r ( t ) d t ,
with the initial conditions
y 1 r ( 0 ) = 0 ,   y 2 r ( 0 ) = 0.
If β = 1 ,   then the exact solutions of OFVIDEs (38) and (39) are given by
[ y ( x ) ] r = [ 1 r ,   r 1 ] sin ( x ) .
In view of the description for FRPS technique and starting with y 0 , 1 r = y 1 r ( 0 ) = 0 , and y 0 , 2 r = y 2 r ( 0 ) = 0 , then the k th-residual functions of Equation (37) will be as follows
R e s k , 1 r ( x ) = D 0 +   β y k , 1 r ( x ) 0 x y k , 2 r ( t ) d t ( 1 r ) , R e s k , 2 r ( x ) = D 0 +   β y k , 2 r ( x ) 0 x y k , 1 r ( t ) d t ( r 1 ) .
where y k , 1 r ( x ) and y k , 2 r ( x ) are given by
y k , 1 r ( x ) = n = 1 k c n x n β Γ ( n β + 1 ) , y k , 2 r ( x ) = n = 1 k d n x n β Γ ( n β + 1 ) .
The coefficients of the FRPS approximate solutions for the system (37) and (38) can be found in Appendix B.
To show the accuracy of the FRPS technique, numerical results for Example 2 at n = 10 , r { 0.5 , 0.75 } with some selected grid points in [ 0 , 1 ] are given in Table 5 and Table 6 with different values of β .
For the purpose of comparison, the absolute errors of y 1 r ( x ) and y 2 r ( x ) of Example 2, case 1, by using the FRPS method and Haar Wavelet (HW) method [35], for n = 8 , with fixed value of β = 1 are given in Table 7. From this table, it can be observed that our method provides us with accurate approximated solutions to OFVIDEs (32) and (33).
To show the fuzzy behavior of the proposed algorithm, the exact and 10th-FRPS approximated solutions of Example 2, case1 and case 2, are plotted in Figure 3. The plots of r -cut representation of exact and FRPS approximated solutions for Example 2, case 1 and case 2, are plotted in Figure 4 at β = 1 , n = 10 and r { 0 , 0.25 , 0.5 , 0.75 , 1 } .

6. Conclusions

In this paper, we have successfully used the FRPS technique to construct and investigate fuzzy approximated solutions for a class of fuzzy fractional Volterra integro-differential equations of order 0 < β 1 , with appropriate fuzzy initial conditions under strongly generalized differentiability. This method is effectively applied without being exposed to any restriction such as conversion, discretization or perturbation. The methodology of the solution basically depends on constructing the residual function and employing the generalized Taylor formula under Caputo sense. Two illustrative examples are given to validate the capability and performance of our algorithm. The numerical results obtained via FRPS algorithm have shown that this technique is an easy, powerful tool which can effectively solve different types of fuzzy fractional integro-differential equations.

Author Contributions

Conceptualization, M.A.; Investigation, M.A.-S.; Methodology, M.A.; Software, M.A.-S.; Supervision, R.R.A. and U.K.S.D.; Writing-original draft, M.A.; Writing-review & editing, R.R.A. and U.K.S.D.

Acknowledgments

Authors gratefully acknowledge to Universiti Kebangsaan Malaysia for providing facilities for our research under the grant no. GP-K007788 and GP-K006926.

Conflicts of Interest

The authors declare no conflict of interest regarding the publication of this paper.

Appendix A

The coefficients of the FRPS approximate solutions for the system (20) and (21): to obtain the value of the coefficients c 1 and d 1 , substitute 1st approximated solutions y 1 , 1 r ( x ) and y 2 , 1 r ( x ) of (23) into R e s 1 , 1 r ( x ) and R e s 1 , 2 r ( x ) , to get, R e s 1 , 1 r ( x ) = ( r + 1 ) ( 1 + x β ) + c 1 ( 1 x 2 β Γ ( 2 β + 1 ) ) , and R e s 1 , 2 r ( x ) = ( 2 r ) ( 1 + x β ) + d 1 ( 1 x 2 β Γ ( 2 β + 1 ) ) . Then, by using the facts R e s 1 , 1 r ( 0 ) = R e s 1 , 2 r ( 0 ) = 0 , yields c 1 = r + 1 , d 1 = 2 r . Therefore, the 1st-FRPS approximated solutions of system (20) and (21) are given by y 1 , 1 r ( x ) = ( r + 1 ) x β Γ ( β + 1 ) and y 1 , 2 r ( x ) = ( 2 r ) x β Γ ( β + 1 ) .
Again, to determine the 2nd unknown coefficients c 2 and d 2 , we have R e s 2 , 1 r ( x ) = ( r + 1 ) ( 1 + x β ) + c 1 ( 1 x 2 β Γ ( 2 β + 1 ) ) + c 2 x β ( 1 Γ ( β + 1 ) x 2 β Γ ( 3 β + 1 ) ) , and R e s 2 , 2 r ( x ) = ( 2 r ) ( 1 + x β ) + d 1 ( 1 x 2 β Γ ( 2 β + 1 ) ) + d 2 x β ( 1 Γ ( β + 1 ) x 2 β Γ ( 3 β + 1 ) ) . Then, through applying the operator D 0 + β on the both sides of R e s 2 , 1 r ( x ) , and R e s 2 , 2 r ( x ) , we have D 0 + β R e s 2 , 1 r ( x ) = ( r + 1 ) β Γ ( β ) c 1 x β Γ ( β + 1 ) + c 2 ( 1 x 2 β Γ ( 2 β + 1 ) ) , and D 0 + β R e s 2 , 2 r ) ( x ) = ( 2 + r ) β Γ ( β ) d 1 x β Γ ( 1 + β ) + d 2 ( 1 x 2 β Γ ( 2 β + 1 ) ) . Finally, by solving the resulting fractional equations at x = 0 , one can get c 2 = ( r + 1 ) β Γ ( β ) ,   d 2 = ( 2 r ) β Γ ( β ) . Hence, y 2 , 1 r ( x ) = ( r + 1 ) ( x β Γ ( β + 1 ) + β Γ ( β ) x 2 β Γ ( 2 β + 1 ) ) and y 2 , 2 r ( x ) = ( 2 r ) ( x β Γ ( β + 1 ) + β Γ ( β ) x 2 β Γ ( 2 β + 1 ) ) .
Likewise, for the 3rd unknown coefficients, c 3 and d 3 , write y 3 , 1 r ( x ) and y 3 , 2 r ( x ) into 3rd residual functions of (24), to conclude that R e s 3 , 1 r ( x ) = ( r + 1 ) ( 1 + x β ) + c 1 + c 2 x β Γ ( β + 1 ) c 2 x 3 β Γ ( 3 β + 1 ) c 3 x 4 β Γ ( 4 β + 1 ) + ( c 1 + c 3 ) x 2 β Γ ( 2 β + 1 ) , and R e s 3 , 2 r ( x ) = ( 2 r ) ( 1 + x β ) + d 1 + d 2 x β Γ ( β + 1 ) d 2 x 3 β Γ ( 3 β + 1 ) d 3 x 4 β Γ ( 4 β + 1 ) + ( d 1 + d 3 ) x 2 β Γ ( 2 β + 1 ) . Thereafter, evaluate D 0 + 2 β of R e s 3 , 1 r ( x ) , and R e s 3 , 2 r ( x ) , to get D 0 + 2 β R e s 3 , 1 r ( x ) = c 1 c 2 x β Γ ( β + 1 ) + c 3 ( 1 x 2 β Γ ( 2 β + 1 ) ) , and D 0 + 2 β R e s 3 , 2 r ( x ) = d 1 d 2 x β Γ ( β + 1 ) + d 3 ( 1 x 2 β Γ ( 2 β + 1 ) ) . Then, based on the fact D 0 + 2 β R e s 3 , 1 r ( 0 ) = D 0 + 2 β R e s 3 , 2 r ( 0 ) = 0 , will yield c 3 = r + 1 , d 3 = 2 r . Thus, y 3 , 1 r ( x ) = ( r + 1 ) ( x β Γ ( β + 1 ) + β Γ ( β ) x 2 β Γ ( 2 β + 1 ) + x 3 β Γ ( 3 β + 1 ) ) and y 3 , 2 r ( x ) = ( 2 r ) ( x β Γ ( β + 1 ) + β Γ ( β ) x 2 β Γ ( 2 β + 1 ) + x 3 β Γ ( 3 β + 1 ) ) .
Using the same argument and the fact D 0 + ( k 1 ) β R e s k , 1 r ( 0 ) = D 0 + ( k 1 ) β R e s k , 2 r ( 0 ) = 0 , for k = 4 , one can get c 4 = ( r + 1 ) β Γ ( β ) , d 4 = ( 2 r ) β Γ ( β ) and hence the 4th-FRPS approximated solutions of system (20) and (21) are given by y 4 , 1 r ( x ) = ( r + 1 ) ( x β Γ ( β + 1 ) + β Γ ( β ) x 2 β Γ ( 2 β + 1 ) + x 3 β Γ ( 3 β + 1 ) + β Γ ( β ) x 4 β Γ ( 4 β + 1 ) ) , and y 4 , 2 r ( x ) = ( 2 r ) ( x β Γ ( β + 1 ) + β Γ ( β ) x 2 β Γ ( 2 β + 1 ) + x 3 β Γ ( 3 β + 1 ) β Γ ( β ) x 4 β Γ ( 4 β + 1 ) ) .
Moreover, the 6th-FRPS approximated solutions of system (20) and (21) are given by
y 6 , 1 r ( x ) = ( r + 1 ) ( 1 Γ ( β + 1 ) x β + β Γ ( β ) Γ ( 2 β + 1 ) x 2 β + 1 Γ ( 3 β + 1 ) x 3 β + β Γ ( β ) Γ ( 4 β + 1 ) x 4 β + 1 Γ ( 5 β + 1 ) x 5 β + β Γ ( β ) Γ ( 63 β + 1 ) x 6 β ) ,
y 6 , 2 r ( x ) = ( 2 r ) ( 1 Γ ( β + 1 ) x β + β Γ ( β ) Γ ( 2 β + 1 ) x 2 β + 1 Γ ( 3 β + 1 ) x 3 β + β Γ ( β ) Γ ( 4 β + 1 ) x 4 β + 1 Γ ( 5 β + 1 ) x 5 β + β Γ ( β ) Γ ( 6 β + 1 ) x 6 β ) .
Therefore, at β = 1 , the approximated solutions of OFVIDEs (20) and (21) can be written as:
y 1 r ( x ) = lim k y k , 1 r ( x ) = ( r + 1 ) ( x + x 2 2 + x 3 6 + x 4 24 + x 5 120 + x 6 720 + ) ,
y 2 r ( x ) = lim k y k , 2 r ( x ) = ( 2 r ) ( x + x 2 2 + x 3 6 + x 4 24 + x 5 120 + x 6 720 + ) .
which coincides precisely with the Taylor series expansion of the exact solutions [ y ( x ) ] r = [ r + 1 ,   2 r ] ( e x 1 ) .
The coefficients of the FRPS approximate solutions for the system (25) and (26): To obtain the value of the coefficients c n and d n , n = 1 , 2 , k , in expansion (28), solve the algebraic system in c n and d n that obtained by the fact D 0 + ( k 1 ) β R e s k , 1 r ( 0 ) = D 0 + ( k 1 ) β R e s k , 2 r ( 0 ) = 0 , 0 < β 1 , k = 1 , 2 , 3 , .
Following the procedure of FRPS-algorithm, the first few coefficients c n and d n are
c 1 = r + 1 ,   d 1 = 2 r , c 2 = ( r + 1 ) β Γ ( β ) ,   d 2 = ( 2 r ) β Γ ( β ) , c 3 = r + 1 , d 3 = 2 r , c 4 = ( 2 r ) β Γ ( β ) , d 4 = ( r + 1 ) β Γ ( β ) , c 5 = r + 1 , d 5 = 2 r ,     and   so   on .
Continue with this procedure to get the 6th-FRPS approximated solutions of (25) and (26) as
y 6 , 1 r ( x ) = ( r + 1 ) Γ ( β + 1 ) x β + ( r + 1 ) β Γ ( β ) Γ ( 2 β + 1 ) x 2 β + ( 2 r ) Γ ( 3 β + 1 ) x 3 β + ( 2 r ) β Γ ( β ) Γ ( 4 β + 1 ) x 4 β + ( r + 1 ) Γ ( 5 β + 1 ) x 5 β + ( r + 1 ) β Γ ( β ) Γ ( 6 β + 1 ) x 6 β ,
y 6 , 2 r ( x ) = ( 2 r ) Γ ( β + 1 ) x β + ( 2 r ) β Γ ( β ) Γ ( 2 β + 1 ) x 2 β + ( r + 1 ) Γ ( 3 β + 1 ) x 3 β + ( r + 1 ) β Γ ( β ) Γ ( 4 β + 1 ) x 4 β + ( 2 r ) Γ ( 5 β + 1 ) x 5 β + ( 2 r ) β Γ ( β ) Γ ( 6 β + 1 ) x 6 β .
Consequently, the approximated solutions of (25) and (26) at β = 1 , can be written as
y 1 r ( x ) = lim k y k , 1 r ( x ) = ( ( r + 1 ) x + ( r + 1 ) 2 x 2 + ( 2 r ) 6 x 3 + ( 2 r ) 24 x 4 + ( r + 1 ) 120 x 5 + ( r + 1 ) 720 x 6 + ) ,
y 2 r ( x ) = lim k y k , 2 r ( x ) = ( ( 2 r ) x + ( 2 r ) 2 x 2 + ( r + 1 ) 6 x 3 + ( r + 1 ) 24 x 4 + ( 2 r ) 120 x 5 + ( 2 r ) 720 x 6 + ) .
which agreement well with the Taylor series expansion of the exact solutions [ y ( x ) ] r = 3 2 e x [ 2 r , r + 1 ] + 1 2 [ 1 2 r , 2 r 1 ] ( cos x sin x ) .

Appendix B

The FRPS approximate solutions of system (32) and (33) can be obtained by determining the unknown coefficients c n and d n of expansion (35).
Anyway, to determine c 1 and d 1 , let k = 1 , in Equation (35), then substitute y 1 , 1 r ( x ) = c 1 x β Γ ( β + 1 ) and y 2 , 1 r ( x ) = d 1 x β Γ ( β + 1 ) , into 1st-residual functions of Equation (36), that is, R e s 1 , 1 r ( x ) = D 0 +   β ( c 1 x β Γ ( β + 1 ) ) 0 x c 1 t β Γ ( β + 1 ) d t ( r 1 ) = 1 r + c 1 ( 1 x 2 β Γ ( 2 β + 1 ) ) , and R e s 1 , 2 r ( x ) = D 0 +   β ( d 1 x β Γ ( β + 1 ) ) 0 x d 1 t β Γ ( β + 1 ) d t ( 1 r ) = r 1 + d 1 ( 1 x 2 β Γ ( 2 β + 1 ) ) . Then, by using the facts R e s 1 , 1 r ( 0 ) = R e s 1 , 2 r ( 0 ) = 0 , yields c 1 = r 1 , d 1 = 1 r . Therefore, the 1st-FRPS approximated solutions of system (33) and (34) is given by y 1 , 1 r ( x ) = ( r 1 ) x β Γ ( β + 1 ) , and y 1 , 2 r ( x ) = ( 1 r ) x β Γ ( β + 1 ) .
To obtain the 2nd unknown coefficients c 2 and d 2 , we have R e s 2 , 1 r ( x ) = r 1 + c 1 ( 1 x 2 β Γ ( 2 β + 1 ) ) + c 2 x β ( 1 Γ ( β + 1 ) x 2 β Γ ( 3 β + 1 ) ) , and R e s 2 , 2 r ( x ) = 1 r + d 1 ( 1 x 2 β Γ ( 2 β + 1 ) ) + d 2 x β ( 1 Γ ( β + 1 ) x 2 β Γ ( 3 β + 1 ) ) . Through utilizing the fact D 0 + β R e s 2 , 1 r ( 0 ) = D 0 + β R e s 2 , 2 r ( 0 ) = 0 , one can get c 2 = 0 ,   d 2 = 0 . Hence, y 2 , 1 r ( x ) = ( r 1 ) x β Γ ( β + 1 ) , and y 2 , 2 r ( x ) = ( 1 r ) x β Γ ( β + 1 ) .
Likewise, for the 3rd unknown coefficients, c 3 and d 3 , write y 3 , 1 r ( x ) and y 3 , 2 r ( x ) into 3rd residual functions of Equation (36), to get R e s 3 , 1 r ( x ) = 1 r + c 1 + c 2 x β Γ ( β + 1 ) c 2 x 3 β Γ ( 3 β + 1 ) c 3 x 4 β Γ ( 4 β + 1 ) + ( c 1 + c 3 ) x 2 β Γ ( 2 β + 1 ) , and R e s 3 , 2 r ( x ) = r 1 + d 1 + d 2 x β Γ ( β + 1 ) d 2 x 3 β Γ ( 3 β + 1 ) d 3 x 4 β Γ ( 4 β + 1 ) + ( d 1 + d 3 ) x 2 β Γ ( 2 β + 1 ) . According to the fact D 0 + 2 β R e s 3 , 1 r ( 0 ) = D 0 + 2 β R e s 3 , 2 r ( 0 ) = 0 , gives c 3 = r 1 , d 3 = 1 r . Thus, y 3 , 1 r ( x ) = ( r 1 ) ( x β Γ ( β + 1 ) + x 3 β Γ ( 3 β + 1 ) ) , and y 3 , 2 r ( x ) = ( 1 r ) ( x β Γ ( β + 1 ) + x 3 β Γ ( 3 β + 1 ) ) .
Using the same argument and the fact D 0 + ( k 1 ) β R e s k , 1 r ( 0 ) = D 0 + ( k 1 ) β R e s k , 2 r ( 0 ) = 0 , for k = 4 , one can get c 4 = 0 , d 4 = 0 . For k = 5 , the 5th-FRPS approximated solutions can be obtained by utilizing the 5th residual functions R e s 5 , 1 r ( x ) , R e s 5 , 2 r ( x ) and the fact D 0 + 4 β R e s 5 , 1 r ( 0 ) = D 0 + 4 β R e s 5 , 2 r ( 0 ) = 0 , such that D 0 + 4 β R e s 5 , 1 r ( x ) = c 3 c 4 x β Γ ( β + 1 ) + c 5 ( 1 x 2 β Γ ( 2 β + 1 ) ) and D 0 + 4 β R e s 5 , 2 r ( x ) = d 3 d 4 x β Γ ( β + 1 ) + d 5 ( 1 x 2 β Γ ( 2 β + 1 ) ) . So, c 5 = r 1 , d 5 = 1 r , and further y 5 , 1 r ( x ) = ( r 1 ) ( x β Γ ( β + 1 ) + x 3 β Γ ( 3 β + 1 ) + x 5 β Γ ( 5 β + 1 ) ) , and y 5 , 2 r ( x ) = ( 1 r ) ( x β Γ ( β + 1 ) + x 3 β Γ ( 3 β + 1 ) + x 5 β Γ ( 5 β + 1 ) ) .
Moreover, based upon the fact D 0 + ( k 1 ) β R e s k , 1 r ( 0 ) = D 0 + ( k 1 ) β R e s k , 2 r ( 0 ) = 0 , k = 5 , 6 , , 10 , the 10th-FRPS approximated solutions of system (32) and (33) are given by
y 10 , 1 r ( x ) = ( r 1 ) ( x β Γ ( β + 1 ) + x 3 β Γ ( 3 β + 1 ) + x 5 β Γ ( 5 β + 1 ) + x 7 β Γ ( 7 β + 1 ) + x 9 β Γ ( 9 β + 1 ) ) ,
y 10 , 2 r ( x ) = ( 1 r ) ( x β Γ ( β + 1 ) + x 3 β Γ ( 3 β + 1 ) + x 5 β Γ ( 5 β + 1 ) + x 7 β Γ ( 7 β + 1 ) + x 9 β Γ ( 9 β + 1 ) ) .
For β = 1 , the approximated solutions of OFVIDEs (32) and (33) can be written in the form
y 1 r ( x ) = lim k y k , 1 r ( x ) = ( r 1 ) ( x + x 3 6 + x 5 120 + x 7 5040 + x 9 362880 + ) ,
y 2 r ( x ) = lim k y k , 2 r ( x ) = ( 1 r ) ( x + x 3 6 + x 5 120 + x 7 5040 + x 9 362880 + ) .
which coincides precisely with Taylor expansion for the exact solutions y 1 r ( x ) = ( r 1 )   sin h ( x ) , and y 2 r ( x ) = ( r 1 )   sin h ( x ) .
The coefficients of the FRPS approximate solutions for the system (37) and (38): To obtain the value of the coefficients c n and d n , n = 1 , 2 , k , in expansion (41), do the following manner:
Through Equations (40) and (41), put k = 1 , we have R e s 1 , 1 r ( x ) = D 0 +   β ( c 1 x β Γ ( β + 1 ) ) 0 x d 1 x β Γ ( β + 1 ) d t ( 1 r ) = r 1 + c 1 d 1 ( 1 x 2 β Γ ( 2 β + 1 ) ) , and R e s 1 , 2 r ( x ) = D 0 +   β ( d 1 x β Γ ( β + 1 ) ) 0 x c 1 x β Γ ( β + 1 ) d t ( r 1 ) = 1 r c 1 ( 1 x 2 β Γ ( 2 β + 1 ) ) + d 1 . Then, we solve R e s 1 , 1 r ( 0 ) = R e s 1 , 2 r ( 0 ) = 0 , it gives c 1 = 1 r , d 1 = r 1 . So, the 1st-FRPS approximated solutions of system (38) and (39) are given by y 1 , 1 r ( x ) = ( 1 r ) x β Γ ( β + 1 ) , and y 1 , 2 r ( x ) = ( r 1 ) x β Γ ( β + 1 ) .
Similarly, for k = 2 , and based on the fact D 0 + β R e s 2 , 1 r ( 0 ) = D 0 + β R e s 2 , 2 r ( 0 ) = 0 , to get c 2 = 0 ,   d 2 = 0 . Thus, y 2 , 1 r ( x ) = ( 1 r ) x β Γ ( β + 1 ) , and y 2 , 2 r ( x ) = ( r 1 ) x β Γ ( β + 1 ) .
To get the 3rd-FRPS approximated solutions, consider k = 3 in Equations (40) and (41) to get R e s 3 , 1 r ( x ) = r 1 + c 1 + c 2 x β Γ ( β + 1 ) + ( c 3 d 1 ) x 2 β Γ ( 2 β + 1 ) d 2 x 3 β Γ ( 3 β + 1 ) d 3 x 4 β Γ ( 4 β + 1 ) , and R e s 3 , 2 r ( x ) = 1 r + d 1 + d 2 x β Γ ( β + 1 ) c 2 x 3 β Γ ( 3 β + 1 ) c 3 x 4 β Γ ( 4 β + 1 ) + ( c 1 + d 3 ) x 2 β Γ ( 2 β + 1 ) . Thereafter, computing 2 β th-Caputo fractional derivative of R e s 3 , 1 r and R e s 3 , 2 r , i.e. D 0 + 2 β R e s 3 , 1 r ( x ) = c 3 d 1 + x β ( d 2 Γ ( β + 1 ) x β d 3 Γ ( 2 β + 1 ) ) , and D 0 + 2 β R e s 3 , 2 r ( x ) = c 1 + d 3 c 2 x β Γ ( β + 1 ) c 3 x β Γ ( 2 β + 1 ) . Lastly, by solving the obtained fractional equations at x = 0 , and considering the values of c 1 ,   d 1 , c 2 and d 2 from previous steps, we get c 3 = r 1 , d 3 = 1 r . Hence, y 3 , 1 r ( x ) = ( 1 r ) x β Γ ( β + 1 ) + ( r 1 ) x 3 β Γ ( 3 β + 1 ) , and y 3 , 2 r ( x ) = ( r 1 ) x β Γ ( β + 1 ) + ( 1 r ) x 3 β Γ ( 3 β + 1 ) .
Continue with this fashion for k = 4 , and depending on the fact D 0 + 3 β R e s 4 , 1 r ( 0 ) = D 0 + 3 β R e s 4 , 2 r ( 0 ) = 0 , will yield c 4 = 0 ,   d 4 = 0 , and hence y 4 , 1 r ( x ) = ( 1 r ) x β Γ ( β + 1 ) + ( r 1 ) x 3 β Γ ( 3 β + 1 ) , and y 4 , 2 r ( x ) = ( r 1 ) x β Γ ( β + 1 ) + ( 1 r ) x 3 β Γ ( 3 β + 1 ) . Next, for k = 5 , the 5th-FRPS approximated solutions can be constructed by considering the 5th-residual functions of Equation (40) and based on the fact D 0 + 4 β R e s 5 , 1 r ( 0 ) = D 0 + 4 β R e s 5 , 2 r ( 0 ) = 0 , such that D 0 + 4 β R e s 5 , 1 r ( x ) = c 5 d 3 d 4 x β Γ ( β + 1 ) d 5 x 2 β Γ ( 2 β + 1 ) and D 0 + 4 β R e s 5 , 2 r ( x ) = d 5 c 3 c 4 x β Γ ( β + 1 ) c 5 x 2 β Γ ( 2 β + 1 ) . Then, by solving the resultant fractional equations at x = 0 , and taking into account the values of c 3 ,   d 3 , c 4 and d 4 , from previous steps, gives c 5 = 1 r , d 5 = r 1 , as well y 5 , 1 r ( x ) = ( 1 r ) x β Γ ( β + 1 ) + ( r 1 ) x 3 β Γ ( 3 β + 1 ) + ( 1 r ) x 5 β Γ ( 5 β + 1 ) , and y 5 , 2 r ( x ) = ( r 1 ) x β Γ ( β + 1 ) + ( 1 r ) x 3 β Γ ( 3 β + 1 ) + ( r 1 ) x 5 β Γ ( 5 β + 1 ) .
Furthermore, through utilizing the fact D 0 + ( k 1 ) β R e s k , 1 r ( 0 ) = D 0 + ( k 1 ) β R e s k , 2 r ( 0 ) = 0 , k = 6 , 7 , , 10 , the 10th-FRPS approximated solutions of system (38) and (39) can be written as
y 10 , 1 r ( x ) = ( 1 r ) x β Γ ( β + 1 ) + ( r 1 ) x 3 β Γ ( 3 β + 1 ) + ( 1 r ) x 5 β Γ ( 5 β + 1 ) + ( r 1 ) x 7 β Γ ( 7 β + 1 ) + ( 1 r ) x 9 β Γ ( 9 β + 1 ) .
y 10 , 2 r ( x ) = ( r 1 ) x β Γ ( β + 1 ) + ( 1 r ) x 3 β Γ ( 3 β + 1 ) + ( r 1 ) x 5 β Γ ( 5 β + 1 ) + ( 1 r ) x 7 β Γ ( 7 β + 1 ) + ( r 1 ) x 9 β Γ ( 9 β + 1 ) .
Consequently, the approximated solutions of system (37) and (38) have general form which are coinciding well with the exact solutions for β = 1 , such that
y 1 r ( x ) = lim k y k , 1 r ( x ) = ( 1 r ) x + ( r 1 ) x 3 6 + ( 1 r ) x 5 120 + ( r 1 ) x 7 5040 + ( 1 r ) x 9 362880 + = ( 1 r ) sin x ,
y 2 r ( x ) = lim k y k , 2 r ( x ) = ( r 1 ) x + ( 1 r ) x 3 6 + ( r 1 ) x 5 120 + ( 1 r ) x 7 5040 + ( r 1 ) x 9 362880 + = ( r 1 ) sin x .

References

  1. Zhang, H.; Liao, X.; Yu, J. Fuzzy modeling and synchronization of hyperchaotic systems. Chaos Solitons Fractals 2005, 26, 835–843. [Google Scholar] [CrossRef]
  2. El Naschie, M.S. From experimental quantum optics to quantum gravity via a fuzzy Kähler manifold. Chaos Solitons Fractals 2005, 25, 969–977. [Google Scholar] [CrossRef]
  3. Diamond, P. Time-dependent differential inclusions, cocycle attractors and fuzzy differential equations. IEEE Trans. Fuzzy Syst. 1999, 7, 734–740. [Google Scholar] [CrossRef]
  4. Feng, G.; Chen, G. Adaptive control of discrete time chaotic systems: a fuzzy control approach. Chaos Solitons Fractals 2005, 23, 459–467. [Google Scholar] [CrossRef]
  5. Priyadharsini, S.; Parthiban, V.; Manivannan, A. Solution of fractional integro-differential system with fuzzy initial condition. Int. J. Pure Appl. Math. 2016, 106, 107–112. [Google Scholar]
  6. Shabestari, M.R.; Ezzati, R.; Allahviranloo, T. Numerical solution of fuzzy fractional integro-differential equation via two-dimensional Legendre wavelet method. J. Intell. Fuzzy Syst. 2018, 34, 2453–2465. [Google Scholar] [CrossRef]
  7. Padmapriya, V.; Kaliyappan, M.; Parthiban, V. Solution of Fuzzy Fractional Integro-Differential Equations Using A Domian Decomposition Method. J. Inform. Math. Sci. 2017, 9, 501–507. [Google Scholar]
  8. Matinfar, M.; Ghanbari, M.; Nuraei, R. Numerical solution of linear fuzzy Volterra integro-differential equations by variational iteration method. J. Intell. Fuzzy Syst. 2013, 24, 575–586. [Google Scholar]
  9. Gumah, G.; Moaddy, K.; Al-Smadi, M.; Hashim, I. Solutions to uncertain Volterra integral equations by fitted reproducing kernel Hilbert space method. J. Funct. Spaces 2016, 2016, 2920463. [Google Scholar] [CrossRef]
  10. Alikhani, R.; Bahrami, F. Global solutions for nonlinear fuzzy fractional integral and integrodifferential equations. Commun. Nonlinear Sci. Numer. Siml. 2013, 18, 2007–2017. [Google Scholar] [CrossRef]
  11. Saadeh, R.; Al-Smadi, M.; Gumah, G.; Khalil, H.; Khan, R.A. Numerical Investigation for Solving Two-Point Fuzzy Boundary Value Problems by Reproducing Kernel Approach. Appl. Math. Inf. Sci. 2016, 10, 1–13. [Google Scholar] [CrossRef]
  12. Gumah, G.N.; Naser, M.F.M.; Al-Smadi, M.; Al-Omari, S.K. Application of reproducing kernel Hilbert space method for solving second-order fuzzy Volterra integro-differential equations. Adv. Diff. Equ. 2018, 2018, 475. [Google Scholar] [CrossRef]
  13. Salahshour, S.; Allahviranloo, T.; Abbasbandy, S. Solving fuzzy fractional differential equations by fuzzy Laplace transforms. Commun. Nonlinear Sci. Numer. Siml. 2012, 17, 1372–1381. [Google Scholar] [CrossRef]
  14. Abu Arqub, O.; Momani, S.; Al-Mezel, S.; Kutbi, M. Existence, uniqueness, and characterization theorems for nonlinear fuzzy integrodifferential equations of Volterra type. Math. Probl. Eng. 2015, 2015, 835891. [Google Scholar] [CrossRef]
  15. Al-Smadi, M.; Abu Arqub, O. Computational algorithm for solving fredholm time-fractional partial integrodifferential equations of dirichlet functions type with error estimates. Appl. Math. Comput. 2019, 342, 280–294. [Google Scholar] [CrossRef]
  16. Abu Arqub, O.; Odibat, Z.; Al-Smadi, M. Numerical solutions of time-fractional partial integrodifferential equations of Robin functions types in Hilbert space with error bounds and error estimates. Nonlinear Dyn. 2018, 94, 1819–1834. [Google Scholar] [CrossRef]
  17. Abu Arqub, O.; Al-Smadi, M. Numerical algorithm for solving time-fractional partial integrodifferential equations subject to initial and Dirichlet boundary conditions. Num. Methods Partial Diff. Equ. 2018, 34, 1577–1597. [Google Scholar] [CrossRef]
  18. Abu Arqub, O.; Al-Smadi, M. Numerical algorithm for solving two-point, second-order periodic boundary value problems for mixed integro-differential equations. Appl. Math. Comput. 2014, 243, 911–922. [Google Scholar] [CrossRef]
  19. Abu Arqub, O. Series solution of fuzzy differential equations under strongly generalized differentiability. J. Adv. Res. Appl. Math. 2013, 5, 31–52. [Google Scholar] [CrossRef]
  20. Moaddy, K.; Al-Smadi, M.; Hashim, I. A novel representation of the exact solution for differential algebraic equations system using residual power-series method. Discrete Dyn. Nat. Soc. 2015, 2015, 205207. [Google Scholar] [CrossRef]
  21. El-Ajou, A.; Abu Arqub, O.; Momani, S.; Baleanu, D.; Alsaedi, A. A novel expansion iterative method for solving linear partial differential equations of fractional order. Int. J. Appl. Math. Comput. 2015, 257, 119–133. [Google Scholar] [CrossRef]
  22. Elajou, A.; Abu Arqub, O.; Zhour, Z.A. New Results on Fractional Power Series: Theory and Applications. Entropy 2013, 12, 5305–5323. [Google Scholar] [CrossRef]
  23. Komashynska, I.; Al-Smadi, M.; Abu Arqub, O.; Momani, S. An efficient analytical method for solving singular initial value problems of nonlinear systems. Appl. Math. Inf. Sci. 2016, 10, 647–656. [Google Scholar] [CrossRef]
  24. Komashynska, I.; Al-Smadi, M.; Ateiwi, A.; Al-Obaidy, S. Approximate analytical solution by residual power series method for system of Fredholm integral equations. Appl. Math. Inf. Sci. 2016, 10, 1–11. [Google Scholar] [CrossRef]
  25. El-Ajou, A.; Abu Arqub, O.; Al-Smadi, M. A general form of the generalized Taylor’s formula with some applications. Appl. Math. Comput. 2015, 256, 851–859. [Google Scholar] [CrossRef]
  26. Alaroud, M.; Al-Smadi, M.; Ahmad, R.R.; Din, S.K.U. Computational optimization of residual power series algorithm for certain classes of fuzzy fractional differential equations. Int. J. Differ. Equ. 2018, 2018, 8686502. [Google Scholar] [CrossRef]
  27. Moaddy, K.; Al-Smadi, M.; Abu Arqub, O.; Hashim, I. Analytic-numeric treatment for handling system of second-order, three-point BVPs. In Proceedings of the AIP Conference, Seoul, Korea, 23–27 May 2017; Volume 1830, p. 020025. [Google Scholar]
  28. Kaleva, O. Fuzzy differential equations. Fuzzy Sets Syst. 1987, 24, 301–317. [Google Scholar] [CrossRef]
  29. Goetschel, J.R.; Voxman, W. Elementary fuzzy calculus. Fuzzy Sets Syst. 1986, 18, 31–43. [Google Scholar] [CrossRef]
  30. Puri, M.L.; Ralescu, D.A. Differentials of fuzzy functions. J. Math. Anal. Appl. 1983, 91, 552–558. [Google Scholar] [CrossRef]
  31. Bede, B.; Gal, S.G. Generalizations of the differentiability of fuzzy-number-valued functions with applications to fuzzy differential equations. Fuzzy Sets Syst. 2005, 151, 581–599. [Google Scholar] [CrossRef]
  32. Moaddy, K.; Freihat, A.; Al-Smadi, M.; Abuteen, E.; Hashim, I. Numerical investigation for handling fractional-order Rabinovich–Fabrikant model using the multistep approach. Soft Comput. 2018, 22, 773–782. [Google Scholar] [CrossRef]
  33. Chalco-Cano, Y.; Roman-Flores, H. On new solutions of fuzzy differential equations. Chaos Solitons Fractals 2008, 38, 112–119. [Google Scholar] [CrossRef]
  34. Salahshour, S.; Allahviranloo, T.; Abbasbandy, S.; Baleanu, D. Existence and uniqueness results for fractional differential equations with uncertainty. Adv. Differ. Equ. 2012, 1, 112. [Google Scholar] [CrossRef]
  35. Ali, M.R. Approximate solutions for fuzzy Volterra integrodifferential equations. J. Abst. Comput. Math. 2018, 3, 11–22. [Google Scholar]
Figure 1. (a) The core and the support of fuzzy (1)-approximated solutions under Caputo [(1)-   β ]-differentiable; (b) The core and the support of derivative of fuzzy (1)-approximated solution under Caputo [(1)-   β ]-differentiable, at β = 1 of Example 1: gray support and red core.
Figure 1. (a) The core and the support of fuzzy (1)-approximated solutions under Caputo [(1)-   β ]-differentiable; (b) The core and the support of derivative of fuzzy (1)-approximated solution under Caputo [(1)-   β ]-differentiable, at β = 1 of Example 1: gray support and red core.
Symmetry 11 00205 g001
Figure 2. (a) The core and the support of fuzzy (2)-approximated solutions under Caputo [(2)-   β ]-differentiable; (b) The core and the support of fuzzy (2)-approximated Caputo derivatives solution under Caputo [(2)-   β ]-differentiable, at β = 1 of Example 1: gray support and red core.
Figure 2. (a) The core and the support of fuzzy (2)-approximated solutions under Caputo [(2)-   β ]-differentiable; (b) The core and the support of fuzzy (2)-approximated Caputo derivatives solution under Caputo [(2)-   β ]-differentiable, at β = 1 of Example 1: gray support and red core.
Symmetry 11 00205 g002
Figure 3. (a) Surface plot exact solution, case 1; (b) Surface plot of 10th-FRPS approximated solutions, case 1; (c) Surface plot exact solution, case 2; (d) Surface plot of 10th-FRPS approximated solutions, case 2, of Example 2 at β = 1 , for all t [ 0 , 1 ] and r [ 0 , 1 ] : (yellow and blue are the upper and lower solution, respectively).
Figure 3. (a) Surface plot exact solution, case 1; (b) Surface plot of 10th-FRPS approximated solutions, case 1; (c) Surface plot exact solution, case 2; (d) Surface plot of 10th-FRPS approximated solutions, case 2, of Example 2 at β = 1 , for all t [ 0 , 1 ] and r [ 0 , 1 ] : (yellow and blue are the upper and lower solution, respectively).
Symmetry 11 00205 g003
Figure 4. (a) Plots of r -cut representations of exact and φ 10 , 1 r ( t ) , of Example 2, case 1; (b) Plots of r -cut representations of exact and φ 10 , 2 r ( t ) , of Example 2, case 2, for β = 1 , t [ 0 , 1 ] with different values of r . in parametric form: red r = 0 , dashed blue r = 0.25 , dashed green r = 0.5 , darker red-dashed r = 0.75 , blue r = 1 .
Figure 4. (a) Plots of r -cut representations of exact and φ 10 , 1 r ( t ) , of Example 2, case 1; (b) Plots of r -cut representations of exact and φ 10 , 2 r ( t ) , of Example 2, case 2, for β = 1 , t [ 0 , 1 ] with different values of r . in parametric form: red r = 0 , dashed blue r = 0.25 , dashed green r = 0.5 , darker red-dashed r = 0.75 , blue r = 1 .
Symmetry 11 00205 g004
Table 1. The (1)-approximated solution of Example 1, case 1 for different values of β with r = 0.5 .
Table 1. The (1)-approximated solution of Example 1, case 1 for different values of β with r = 0.5 .
β i 7th-FRPS Approximated Solutions
1 4 y 7 ( x ) = 3 Γ ( 1 / 4 ) x 4 π + 3 Γ ( 1 / 4 ) x 8 + Γ ( 1 / 4 ) x 3 / 2 2 π + 3 x 1 / 4 2 Γ ( 5 / 4 ) + 3 x 3 / 4 2 Γ ( 7 / 4 ) + 3 x 5 / 4 2 Γ ( 9 / 4 ) + 3 x 7 / 4 2 Γ ( 11 / 4 )
1 2 y 7 ( x ) = 3 x π + 3 π x 4 + 2 x 3 / 2 π + 3 π x 2 8 + 4 x 5 / 2 5 π + π x 3 8 + 8 x 7 / 2 35 π  
3 4 y 7 ( x ) = 3 Γ ( 3 / 4 ) x 3 / 2 2 π + 3 Γ ( 3 / 4 ) x 3 16 + 4 Γ ( 3 / 4 ) x 9 / 2 105 π + 3 x 3 / 4 2 Γ ( 7 / 4 ) + 3 x 9 / 4 2 Γ ( 13 / 4 ) + 3 x 15 / 4 2 Γ ( 19 / 4 ) + 3 x 21 / 4 2 Γ ( 25 / 4 )  
1 y 7 ( x ) = 3 2 x + 3 4 x 2 + 1 4 x 3 + 1 16 x 4 + 1 80 x 5 + 1 480 x 6 + 1 3360 x 7
Table 2. The (2)-approximated solution of Example 5.1, case 2 for different values of β with r = 0.5 .
Table 2. The (2)-approximated solution of Example 5.1, case 2 for different values of β with r = 0.5 .
β i 7th-FRPS Approximated Solutions
1 4 y 7 ( x ) = 3 Γ ( 1 / 4 ) x 4 π + 3 Γ ( 1 / 4 ) x 8 + Γ ( 1 / 4 ) x 3 / 2 2 π + 3 x 1 / 4 2 Γ ( 5 / 4 ) + 3 x 3 / 4 2 Γ ( 7 / 4 ) + 3 x 5 / 4 2 Γ ( 9 / 4 ) + 3 x 7 / 4 2 Γ ( 11 / 4 )
1 2 y 7 ( x ) = 3 x π + 3 π x 4 + 2 x 3 / 2 π + 3 π x 2 8 + 4 x 5 / 2 5 π + π x 3 8 + 8 x 7 / 2 35 π  
3 4 y 7 ( x ) = 3 Γ ( 3 / 4 ) x 3 / 2 2 π + 3 Γ ( 3 / 4 ) x 3 16 + 4 Γ ( 3 / 4 ) x 9 / 2 105 π + 3 x 3 / 4 2 Γ ( 7 / 4 ) + 3 x 9 / 4 2 Γ ( 13 / 4 ) + 3 x 15 / 4 2 Γ ( 19 / 4 ) + 3 x 21 / 4 2 Γ ( 25 / 4 )  
1 y 7 ( x ) = 3 2 x + 3 4 x 2 + 1 4 x 3 + 1 16 x 4 + 1 80 x 5 + 1 480 x 6 + 1 3360 x 7
Table 3. Absolute errors for Example 1, case 1, at β = 1 and n = 10 .
Table 3. Absolute errors for Example 1, case 1, at β = 1 and n = 10 .
y 1 r ( x )
x i r = 0 r = 0.5 r = 1
0.2 5.273559366 × 10 16 7.216449660 × 10 16 1.054711873 × 10 15
0.4 1.086908341 × 10 12 1.630362511 × 10 12 2.173816682 × 10 12
0.6 9.565170878 × 10 11 1.434776741 × 10 10 1.913034175 × 10 10
0.8 2.304785251 × 10 9 3.457177877 × 10 9 4.609570503 × 10 9
y 2 r ( x )
x i r = 0 r = 0.5 r = 1
0.2 1.054711873 × 10 15 7.216449660 × 10 16 5.273559366 × 10 16
0.4 2.173816682 × 10 12 1.630362511 × 10 12 1.086908341 × 10 12
0.6 1.913034175 × 10 10 1.434776741 × 10 10 9.565170877 × 10 11
0.8 4.609570503 × 10 9 3.457177877 × 10 9 2.304785251 × 10 9
Table 4. Absolute errors for Example 1, case 2, at β = 1 and n = 10 .
Table 4. Absolute errors for Example 1, case 2, at β = 1 and n = 10 .
y 1 r ( x )
x i r = 0 r = 0.5 r = 1
0.2 1.137978600 × 10 15 8.326672685 × 10 16 4.996003611 × 10 16
0.4 2.172595439 × 10 12 1.630362512 × 10 12 1.088018564 × 10 12
0.6 1.910848146 × 10 10 1.434776742 × 10 10 9.587064476 × 10 11
0.8 4.600237524 × 10 9 3.457178099 × 10 9 2.314118674 × 10 9
y 2 r ( x )
x i r = 0 r = 0.5 r = 1
0.2 5.551115123 × 10 16 8.326672685 × 10 16 1.193489751 × 10 15
0.4 1.088018564 × 10 12 1.630362512 × 10 12 2.172595437 × 10 12
0.6 9.587064476 × 10 11 1.434776742 × 10 10 1.910848146 × 10 10
0.8 2.314118674 × 10 9 3.457178099 × 10 9 4.600237524 × 10 9
Table 5. Numerical results for Example 2, case 1, with different values of β .
Table 5. Numerical results for Example 2, case 1, with different values of β .
y 1 r ( x )
r i x i β = 1 β = 0.9 β = 0.8 β = 0.7
0.5 0.2 0.100668001 0.123692685 0.151695235 0.186264438
0.4 0.205376163 0.238161466 0.277064940 0.324760477
0.6 0.318326791 0.359432909 0.408753265 0.470283639
0.8 0.444052990 0.494515208 0.556205023 0.634532156
0.75 0.2 0.050334001 0.061846342 0.075847617 0.093132219
0.4 0.102688081 0.119080733 0.138532470 0.162380238
0.6 0.159163396 0.179716454 0.204376633 0.235141819
0.8 0.222026495 0.247257604 0.278102512 0.317266078
y 2 r ( x )
r i x i β = 1 β = 0.9 β = 0.8 β = 0.7
0.5 0.2 0.100668001 0.123692685 0.151695235 0.186264438
0.4 0.205376163 0.238161466 0.277064940 0.324760477
0.6 0.318326791 0.359432909 0.408753265 0.470283639
0.8 0.444052990 0.494515208 0.556205023 0.634532156
0.75 0.2 0.050334001 0.061846342 0.075847617 0.093132219
0.4 0.102688081 0.119080733 0.138532470 0.162380238
0.6 0.159163396 0.179716454 0.204376633 0.235141819
0.8 0.222026495 0.247257604 0.278102512 0.317266078
Table 6. Numerical results for Example 2, case 2, with different values of β .
Table 6. Numerical results for Example 2, case 2, with different values of β .
y 1 r ( x )
r i x i β = 1 β = 0.9 β = 0.8 β = 0.7
0.5 0.2 0.099334665 0.120583966 0.144646654 0.170765057
0.4 0.194709171 0.217958708 0.239847027 0.258218270
0.6 0.282321238 0.299033752 0.310147377 0.313819870
0.8 0.358678047 0.363060372 0.359026685 0.346425845
0.75 0.2 0.049667333 0.060291983 0.072323327 0.085382529
0.4 0.097354586 0.108979354 0.119923514 0.129109135
0.6 0.141160618 0.149516876 0.155073688 0.156909935
0.8 0.179339023 0.181530186 0.179513342 0.173212923
y 2 r ( x )
r i x i β = 1 β = 0.9 β = 0.8 β = 0.7
0.5 0.2 0.099334665 0.120583966 0.144646654 0.170765057
0.4 0.194709171 0.217958708 0.239847027 0.258218270
0.6 0.282321238 0.299033752 0.310147377 0.313819870
0.8 0.358678047 0.363060372 0.359026685 0.346425845
0.75 0.2 0.049667333 0.060291983 0.072323327 0.085382529
0.4 0.097354586 0.108979354 0.119923514 0.129109135
0.6 0.141160618 0.149516876 0.155073688 0.156909935
0.8 0.179339023 0.181530186 0.179513342 0.173212923
Table 7. Comparisons of absolute errors of Example 2, case 1.
Table 7. Comparisons of absolute errors of Example 2, case 1.
y 1 r ( x )
x i FRPS MethodHW Method
0.1 2.76167 × 10 15 4.205 × 10 6
0.2 1.41145 × 10 12 5.305 × 10 6
0.3 5.42854 × 10 11 2.07 × 10 6
y 2 r ( x )
x i FRPS MethodHW Method
0.1 2.76167 × 10 15 2.876 × 10 6
0.2 1.41145 × 10 12 1.836 × 10 6
0.3 5.42854 × 10 11 1.295 × 10 6
1 Results of HW method referred in [35].

Share and Cite

MDPI and ACS Style

Alaroud, M.; Al-Smadi, M.; Rozita Ahmad, R.; Salma Din, U.K. An Analytical Numerical Method for Solving Fuzzy Fractional Volterra Integro-Differential Equations. Symmetry 2019, 11, 205. https://0-doi-org.brum.beds.ac.uk/10.3390/sym11020205

AMA Style

Alaroud M, Al-Smadi M, Rozita Ahmad R, Salma Din UK. An Analytical Numerical Method for Solving Fuzzy Fractional Volterra Integro-Differential Equations. Symmetry. 2019; 11(2):205. https://0-doi-org.brum.beds.ac.uk/10.3390/sym11020205

Chicago/Turabian Style

Alaroud, Mohammad, Mohammed Al-Smadi, Rokiah Rozita Ahmad, and Ummul Khair Salma Din. 2019. "An Analytical Numerical Method for Solving Fuzzy Fractional Volterra Integro-Differential Equations" Symmetry 11, no. 2: 205. https://0-doi-org.brum.beds.ac.uk/10.3390/sym11020205

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop