Next Article in Journal
Quasi Association Models for Square Contingency Tables with Ordinal Categories
Previous Article in Journal
Applications of Symmetric Conic Domains to a Subclass of q-Starlike Functions
Previous Article in Special Issue
On Conditions for L2-Dissipativity of an Explicit Finite-Difference Scheme for Linearized 2D and 3D Barotropic Gas Dynamics System of Equations with Regularizations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Approximate Analytic–Numeric Fuzzy Solutions of Fuzzy Fractional Equations Using a Residual Power Series Approach

1
Department of Mathematics, Faculty of Arts and Science, Amman Arab University, Amman 11953, Jordan
2
Mathematics Department, Al-Qunfudah University College, Umm Al-Qura University, Mecca 24382, Saudi Arabia
3
Department of Scientific Basic Sciences, Faculty of Engineering Technology, Al-Balqa Applied University, Amman 11134, Jordan
*
Author to whom correspondence should be addressed.
Symmetry 2022, 14(4), 804; https://doi.org/10.3390/sym14040804
Submission received: 3 March 2022 / Revised: 4 April 2022 / Accepted: 7 April 2022 / Published: 12 April 2022

Abstract

:
In this article, we consider a reliable analytical and numerical approach to create fuzzy approximated solutions for differential equations of fractional order with appropriate uncertain initial data by the means of a residual error function. The concept of strongly generalized differentiability is utilized to introduce the fuzzy fractional derivatives. The proposed method provides a systematic scheme based on generalized Taylor expansion and minimization of the residual error function, so as to obtain the coefficients values of a fractional series based on the given initial data of triangular fuzzy numbers in the parametric form. The obtained approximated solutions are provided within an appropriate radius to the requisite domain in the form of rapidly convergent fractional series according to their parametric form. The method’s performance and applicability are verified by applying it on some numerical examples. The impact of r -levels and fractional order γ is presented quantitatively and graphically, showing the coincidence between the exact and the fuzzy approximated solutions. Moreover, for reliability and accuracy, our obtained results are numerically compared with the exact solutions and with results obtained using other methods described in the literature. This indicates that the proposed approach overcomes the difficulties that appear in other approaches to create fractional series solutions for varied uncertain natural problems arising within the fields of applied physics and engineering.

1. Introduction

Uncertain models are one of the most significant parts of the fuzzy analysis theory and have rapidly developed in the last decades. With this, recent theoretical and applied aspects have been discussed by many mathematicians, including measure, symmetric and control theories, radiation transfer in a semi-infinite atmosphere, and so forth. They are considered influential tools for modeling several real-life situations and phenomena in which uncertainty results from several factors such as measurement errors, deficient data, and initial guesses. Recently, numerous publications have shown that fractional differential equations (FDEs) are a powerful and applicable instrument to describe the exact results of physical, applied mathematics, and engineering phenomena such as control systems, aerodynamics, signal processing, bio-mathematical problems, and others [1,2,3,4,5,6]. However, in some cases, the raw initial data are imprecise and could be replaced by uncertain initial data to obtain fuzzy FDEs. Therefore, fuzzy FDEs are crucial in fuzzy calculus and are a widespread model in various natural scientific areas, including population analysis, evaluation of weapon systems, civil engineering, and modeling in electro-hydraulics.
The concept of the fuzzy FDEs was first introduced by Agarwal et al. [7], who investigated fuzzy solutions for a certain class of fuzzy FDEs under the Hukuhara differentiability in the sense of Riemann–Liouville differentiability. Thereafter, many researchers have investigated solutions of ordinary fuzzy DEs and fuzzy FDEs (for more details, we refer to [8,9,10,11,12,13]). In addition, some mathematicians showed an interest in the existence and uniqueness of solutions for fuzzy FDEs. The authors of [14] proposed the existence and uniqueness of the solution to a fuzzy FDE under Hukuhara fractional Riemann–Liouville differentiability. Later, new and different techniques and methods were presented so that the existence and uniqueness of solutions for fuzzy FDEs were proved. Alikhani et al. [15] also confirmed the results of the existence and uniqueness of nonlinear fuzzy fractional integral and integro-differential equations by using the technique of upper and lower solutions. Additionally, these authors examined some related results about the existence and uniqueness of solutions to fuzzy FDEs under Caputo type-2 fuzzy fractional derivative and the definition of Laplace transform of type-2 fuzzy number-valued functions [16]. For instance, in [17], Salahshour et al. recommended some novel and different results for the existence and uniqueness of solutions of fuzzy FDEs.
Providing exact solutions to fuzzy FDEs is a difficult task. As a result, it is necessary to develop a robust numeric–analytic approach to deal with the complications of uncertain models and attain a precise mathematical framework for processing fuzzy initial value problems (IVPs) [18,19,20,21]. This analysis aimed to apply a recent treatment method, called the residual power series (RPS) method, to provide fuzzy approximated analytical solutions for a class of fuzzy FDEs under the concept of strongly generalized differentiability. This concept was introduced and discussed by Bede and Gal [22]. Later, it was developed and investigated (for more details see [23,24]). In fact, by utilizing strongly generalized differentiability, it is possible to find solutions for larger classes of fuzzy FDEs than by using other types of differentiability. More specifically, we here provide fuzzy approximated analytical solutions for the fuzzy fractional initial value problem (FFIVP) of the general form:
{ D a + γ φ ( t ) = F ( t , φ ( t ) ) , a   t b , φ ( a ) = μ ,
where D a + γ is the fuzzy Caputo fractional derivative of order γ : 0 < γ 1 , F : [ a , b ] × R F R F is a continuous fuzzy-valued function, φ ( t ) is an unknown fuzzy analytical function to be determined, and μ R F , where R F stands for the set of fuzzy numbers on a real line.
In 2013, the scholar Abu Arqub [25] proposed the RPS as an effectively numeric–analytic approach and easily applied it to define the components of the suggested series solutions to a certain class of classical fuzzy DEs. Later, the RPS approach was developed for handling various kinds of FDEs [26,27,28,29]. This approach produces solutions to the given problem in the convergent generalized Taylor’s series formula without involving discretization, linearization, or perturbation [30,31,32]. It may be applied directly to given problems by selecting an appropriate value for the initial guess approximation. Recently, applications of the RPS approach for the simulation and creation of analytical solutions of FDEs, partial FDEs, and fuzzy FDEs have become popular and diverse, and numerous real-world problems have been studied and analyzed using the RPS approach, such as fractional stiff systems [33], time-fractional Fokker–Planck equations [34], time-fractional Whitham–Broer–Kaup equations [35], time-fractional Sharma–Tasso–Olever equation [34], fractional Newell–Whitehead–Segel equation [36], coupled fractional resonant Schrödinger equation [37], fractional foam drainage equation [38], and certain class of fractional systems of partial differential equations [39].
Approximate analytic–numeric techniques are considered to deal with fuzzy models of fractional PDEs, systems of fractional ODEs, and delay differential models. However, fuzzy fractional differential equations have not been investigated using the fractional power series method. Motivated by this, the primary objective of this work was to provide approximate analytic numerical solutions to fuzzy fractional initial value problems (IVPs) utilizing the RPS. The current article is organized as follows: Section 2 presents some of the well-known concepts and primary results of the fuzzy set theory and fuzzy fractional calculus theory. Section 3 discusses the formulation of the FFIVP (1) in the parametric form. Some FFIVPs are considered to demonstrate the efficiency and applicability of the RPS scheme presented in Section 4. Concluding remarks are outlined in Section 5.

2. Overview of the Fuzzy Fractional Calculus

In the subsequent section, we revise the most significant definitions and preliminary results related to the fuzzy fractional calculus. Assume that R F = { σ : R [ 0 , 1 ] } , where σ is a fuzzy number satisfying the following conditions:
  • σ is normal; that is, there is an element ζ R m such that σ ( ζ ) = 1 .
  • σ is convex; that is, for each ζ 1 , ζ 2 R m and 0 θ 1 , we have σ ( θ ζ 1 + ( 1 θ ) ζ 2 ) min ( σ ( ζ 1 ) , σ ( ζ 2 ) ) .
  • σ is upper semi-continuous.
  • The closure of supp ( σ ) is a compact subset, where supp ( σ ) = { ζ R m : σ ( ζ ) > 0 } .
  • R F is the fuzzy numbers set.
For r ( 0 , 1 ] , the r -level representation of the fuzzy number σ is defined by [ σ ] r = { ζ R m : σ ( ζ ) > r } . Then, σ R F if the r -level representation is a compact convex subset of R . So, if σ R F , then [ σ ] r = [ σ 1 r , σ 2 r ] , such that σ 1 r = min { ζ : ζ [ σ ] r } , and σ 2 r = max { ζ : ζ [ σ ] r } .
Definition 1.
Ref. [40] A triangular fuzzy number σ is defined as a fuzzy set in R F , which is given by σ = ( s , t , u ) R 3 , with s t u , where the lower bound σ 1 r = s + ( t s ) , and the upper bound σ 2 r = u ( u t ) r are the endpoints of the r -level representation for each r [ 0 , 1 ] .
The Hausdorff distance between two arbitrary fuzzy numbers σ and ϑ is defined as a mapping d H :   R F × R F R + { 0 } and given by d h ( σ , ϑ ) = sup 0   r 1   max { σ 1 r ϑ 1 r , σ 2 r ϑ 2 r } , where the r -level representations of σ and ϑ are [ σ ] r = [ σ 1 r , σ 2 r ] and [ ϑ ] r = [ ϑ 1 r , ϑ 2 r ] , respectively.
Definition 2.
Ref. [22] Let φ : ( a , b ) R F and fix t 0 [ a , b ] . One can say φ is strongly generalized differentiable at t 0 , if there is an element φ ( t 0 ) R F such that either:
  • The -differences φ ( t 0   + η ) φ ( t 0 )   ,   φ ( t 0 ) φ ( t 0   η ) exist,   η > 0 , sufficiently approach to 0, and lim η 0 + φ ( t 0 + η ) φ ( t 0 ) η = φ ( t 0 ) = lim η 0 + φ ( t 0 ) φ ( t 0 η ) η ,
    or
  • The -differences φ ( t 0 ) φ ( t 0   + η ) , φ ( t 0   η ) φ ( t 0 ) exist   η > 0 , sufficiently approach to 0, and lim η 0 + φ ( t 0 )   φ ( t 0   + η ) η = φ ( t 0 ) = lim η 0 + φ ( t 0 η ) φ ( t 0 ) η .
    where the limits here are taken in the complete metric space ( R F ,   d H ) .
Remark 1.
One can say φ is differentiable on ( a , b ) , when φ is differentiable for any point t ( a , b ) . Furthermore:
  • φ is (1)-differentiable on ( a , b ) , and its derivative of φ at t = t 0   is given by φ ( t 0   ) = D 1 1 φ ( t 0   ) , when φ is differentiable in terms of the first condition of Definition 2.
  • φ is (2)-differentiable on ( a , b ) , and its derivative of φ at t = t 0   is given by φ ( t 0   ) = D 2 1 φ ( t 0   ) , when φ is differentiable in terms of the second condition of Definition 2.
Definition 3.
Ref. [17] Let φ : [ a , b ]   R F and φ C F [ a , b ] L F [ a , b ] . One can say φ is Caputo fuzzy H -differentiable at t when D a + γ φ ( t ) = 1 Γ ( 1 γ ) a t φ ( τ ) ( t τ ) γ   d τ exists, where 0 < γ 1 .
Remark 2.
We say that φ is Caputo [ ( 1 ) - γ ] - differentiable if φ is (1)-differentiable, and φ is Caputo [ ( 2 ) - γ ] - differentiable if φ is (2)-differentiable.
Theorem 1.
Ref. [17] Let φ : [ a , b ]   R F and φ C F [ a , b ] , for γ ( 0 , 1 ]   . Then, the fuzzy Caputo fractional derivative exists on ( a , b ) for each r [ 0 , 1 ] , such that
  • [ D a + γ φ ( t ) ] r = [ D a + γ φ 1 r ( t ) ,   D a + γ φ 2 r ( t ) ] , when φ is (1)-differentiable.
  • [ D a + γ φ ( t ) ] r = [ D a + γ φ 2 r ( t ) ,   D a + γ φ 1 r ( t ) ] , when φ is (2)-differentiable.
Definition 4.
Ref. [29] A power series representation at t = a has the following form
n = 0 c n ( t a )   n γ , 0 n 1 < γ n ,
It is called a fractional series (FS), where t a and c n ’s are called the coefficients of the series.

3. Fuzzy Fractional Initial Value Problems

In this section, we study a certain class of FFIVPs in the meaning of Caputo’s fuzzy H-differentiability throughout converting the main problem from the fuzzy environment into a crisp environment based on the differentiability type. Furthermore, we present an algorithm to solve the new system which consists of two fractional initial value problems (FIVPs).
The formulation of the target problem is the significant part of the procedure. Anyhow, to create the fuzzy solution of the FFIVPs, we reformulate (1) based on the type of differentiability in the r -level representation as follows:
{ [ D a + γ φ ( t ) ] r = [ F ( t , φ ( t ) ) ] r ,   a t b , [ φ ( a ) ] r = [ μ ] r ,
where
[ D a + γ φ ( t ) ] r = [ D a + γ φ 1 r ( t ) , D a + γ φ 2 r ( t ) ] , [ μ ] r = [ μ 1 r , μ 2 r ] , and
[ F ( t , φ ( t ) ) ] r = [ F 1 r ( t , φ 1 r ( t ) , φ 2 r ( t ) ) ,   F 2 r ( t , φ 1 r ( t ) , φ 2 r ( t ) ] , for r [ 0 , 1 ] , and γ ( 0 , 1 ] .
For n { 1 , 2 } , the ( n ) -solution of FFIVPs (1) is a fuzzy function φ : [ a , b ]   R F that has Caputo [ ( n ) - γ ]-differentiable and satisfies the fuzzy FIVPs (1). The next algorithm (Algorithm 1) along with Theorem 1 assisted us to find these solutions, ignoring the fuzzy settings approach:
Algorithm   1 :   To   determine   the   ( n ) -solutions of FFIVPs (1), we considered the following cases:
Case 1: Under Caputo [(1)- γ ]-differentiable, the FFIVPs (1) converts to the following FIVPs system
{ D a + γ φ 1 r ( t ) = F 1 r ( t , φ 1 r ( t ) , φ 2 r ( t ) ) , D a + γ φ 2 r ( t ) = F 2 r ( t , φ 1 r ( t ) , φ 2 r ( t ) ) ,   φ 1 r ( a ) = μ 1 r ,   and   φ 2 r ( a ) = μ 2 r .
 Then, we used the following procedure:
  •    First: Solve the system (3) for φ 1 r ( t ) and φ 2 r ( t ) .
  •    Second: Ensure that [ φ 1 r ( t ) , φ 2 r ( t ) ] and [ D a + γ φ 1 r ( t ) ,   D a + γ φ 2 r ( t ) ] are valid level sets for each r [ 0 , 1 ] .
  •    Third: Construct the (1)-solution, φ ( t ) ,   whose   r - level   representation   is   [ φ 1 r ( t ) , φ 2 r ( t ) ] .
Case 2: Under Caputo [(2)- γ ]-differentiable, the FFIVPs (1) converts to the following FIVPs system:
{ D a + γ φ 1 r ( t ) = F 2 r ( t , φ 1 r ( t ) , φ 2 r ( t ) ) , D a + γ φ 2 r ( t ) = F 1 r ( t , φ 1 r ( t ) , φ 2 r ( t ) ) ,   φ 1 r ( a ) = μ 1 r ,   and   φ 2 r ( a ) = μ 2 r .
 Then, we performed the following procedure:
  •    First: Solve the system (4) for φ 1 r ( t ) and φ 2 r ( t ) .
  •    Second: Ensure that [ φ 1 r ( t ) , φ 2 r ( t ) ]   and   [ D a + γ φ 2 r ( t ) ,   D a + γ φ 1 r ( t ) ]   are   valid   level   sets   for   each   r [ 0 , 1 ] .
  •    Third: Construct the (2)-solution, φ ( t ) ,   whose   r - level   representation   is   [ φ 1 r ( t ) , φ 2 r ( t ) ] .
Remark 3.
Let n { 1 , 2 } and let [ φ ( t ) ] r = [ φ 1 r ( t ) ,   φ 2 r ( t ) ] be an ( n ) -solution of FFIVPs (1) on [ a , b ] . Then, φ 1 r ( t ) and φ 2 r ( t ) will be the solutions to the ( n ) -corresponding FIVPs systems.
Remark 4.
Let n { 1 , 2 } and let φ 1 r ( t )   , φ 2 r ( t ) represent the solutions of ( n ) -corresponding FIVPs systems for each r [ 0 , 1 ] . If [ φ ( t ) ] r = [ φ 1 r ( t ) ,   φ 2 r ( t ) ] has valid level sets and φ ( t ) is Caputo [ ( n ) - γ ]-differentiable, then φ ( t ) is an ( n ) -solution of FFIVPs (1) on [ a , b ] .

4. Application of the RPS Method to Solve FFIVPs

In this section, the fundamental principle of the proposed technique is introduced to predict and obtain analytical solutions for FFIVPs (1). The RPS approach provides an approximate solution by substituting the FPS expansion in its fractional truncated residual function.
Theorem 2.
Suppose that φ 1 r ( t ) and φ 2 r ( t ) have the following FS expansions about t = a :
φ 1 r ( t ) = n = 0 ω n ( t a ) n γ Γ ( n γ + 1 ) , φ 2 r ( t ) = n = 0 ρ n ( t a ) n γ Γ ( n γ + 1 ) .
where γ ( 0 , 1 ] and t [ a , a + R ) . If D a + γ φ m r ( t ) C [ a , a + R ) , for m { 1 , 2 } , then the coefficients ω n and ρ n will be written as ω n = D a + n γ φ 1 r ( t ) , and ρ n = D a + n γ φ 2 r ( t ) , so that D a + n γ = D a + γ · D a + γ · · D a + γ ( n -times).
Proof: 
Let φ 1 r ( t ) and φ 2 r ( t ) be two arbitrary functions that could be expressed by an FS expansion (4). If we substitute = , into (5), one can notice that ω 0 = φ 1 r ( a ) , ρ 0 = φ 2 r ( a ) , and ω n = ρ n = 0 , for n = 1 , 2 , .
On other hand, operating D a + γ on both sides of (5) gives
D a + γ φ 1 r ( t ) = ω 1 + ω 2 ( t a ) γ Γ ( γ + 1 ) + ω 3 ( t a ) 2 γ Γ ( 2 γ + 1 ) + , D a + γ φ 2 r ( t ) = ρ 1 + ρ 2 ( t a ) γ Γ ( γ + 1 ) + ρ 3 ( t a ) 2 γ Γ ( 2 γ + 1 ) + .
Then, by substituting t = a into (6), we obtain ω 1 = D a + γ φ 1 r ( a ) and ρ 1 = D a + γ φ 2 r ( a ) .
Next, by applying D a + γ once on the resulting Equation (6):
D a + 2 γ φ 1 r ( t ) = D a + γ ( D a + γ φ 1 r ( t ) ) = ω 2 + ω 3 ( t a ) γ Γ ( γ + 1 ) + ω 4 ( t a ) 3 γ Γ ( 3 γ + 1 ) + , D a + 2 γ φ 2 r ( t ) = D a + γ ( D a + γ φ 2 r ( t ) ) = ρ 2 + ρ 3 ( t a ) γ Γ ( γ + 1 ) + ρ 4 ( t a ) 3 γ Γ ( 3 γ + 1 ) + .
Here, if t = a in (7), then the second coefficients of (5) will be ω 2 = D a + 2 γ φ 1 r ( a ) and ρ 2 = D a + 2 γ φ 2 r ( a ) . Likewise, by operating D a + γ on both sides of (6) and substituting t = a into the resultant fractional equation the result is ω 3 = D a + 3 γ φ 1 r ( a ) , and ρ 3 = D a + 3 γ φ 2 r ( a ) .
In the same way, we applied D a + γ , n -times, and then considered t = a in the resultant fractional equation, then the pattern of the unknown coefficients were obtained, and hence ω n and ρ n , in the FS expansions (5) had the general forms ω n = D a + n γ φ 1 r ( t ) and ρ n = D a + n γ φ 2 r ( t ) , for n = 0 , 1 , 2 , 3 , . □
To reach our purpose, the following approach was used under Caputo [(1)- γ ]-differentiable. Likewise, it can be applied to solve FFIVPs (1) under Caputo [(2)- γ ]-differentiable.
Step A: According to Theorem 2, the RPS solutions of FIVPs system (3) at a = 0 have the following FS forms:
φ 1 r ( t ) = n = 0 ω n t n γ Γ ( n γ + 1 ) , φ 2 r ( t ) = n = 0 ρ n t n γ Γ ( n γ + 1 ) .
It is clear that φ 1 r ( t ) and φ 2 r ( t ) satisfy the initial condition of (3), then φ 1 r ( 0 ) = ω 0 = μ 1 r and φ 2 r ( 0 ) = ρ 0 = μ 2 r will be the initial guess approximations for (3). So, the series solutions can be written as:
φ 1 r ( t ) = μ 1 r + n = 1 ω n t n γ Γ ( n γ + 1 ) , φ 2 r ( t ) = μ 2 r + n = 1 ρ n t n γ Γ ( n γ + 1 ) .
Then, the k t h -truncated series of the solutions φ 1 r ( t ) and φ 2 r ( t ) can be given as:
φ k , 1 r ( t ) = μ 1 r + n = 1 k ω n t n γ Γ ( n γ + 1 ) , φ k , 2 r ( t ) = μ 2 r + n = 1 k ρ n t n γ Γ ( n γ + 1 ) .
Step B: Identify the so-called k -th residual functions of (3) as follows:
r e s k , 1 r ( t ) = D 0 + γ φ k , 1 r ( t ) F 1 r ( t , φ k , 1 r ( t ) , φ k , 2 r ( t ) ) , r e s k , 2 r ( t ) = D 0 + γ φ k , 2 r ( t ) F 2 r ( t , φ k , 1 r ( t ) , φ k , 2 r ( t ) ) .
As in [18], we note that lim k r e s k , m r ( t ) = r e s m r = 0 , for m = { 1 , 2 } , and each t 0 . In fact, this leads to D 0 + n γ r e s m r ( t ) = 0 , because of D 0 + γ C = 0 , for any constant C . Further, D 0 + n γ r e s m r ( t ) and D 0 + n γ r e s k , m r ( t ) are equivalent at t = 0 , for each n = 0 , 1 , 2 , , k .
Step C: Substitute the k t h -truncated series of the solutions φ 1 r ( t ) and φ 2 r ( t ) of (7) into the k -th residual functions r e s k , 1 r and r e s k , 2 r ( t ) .
Step D: Apply the fractional operator D 0 + ( k 1 ) γ , for k = 1 , 2 , to both sides of the obtained fractional equations in Step C and then solve the following fractional systems for the target unknown coefficients:
  D 0 + ( k 1 ) γ r e s k , 1 r ( 0 ) = 0 , D 0 + ( k 1 ) γ r e s k , 2 r ( 0 ) = 0 .
Step E: After solving (12), we obtained the forms of ω n and ρ n in the expansions (10), and hence the k t h -truncated series solutions were found.
Now, to find ω 1 and ρ 1 , we considered k = 1 , in (10), then substituted φ 1 , 1 r ( t ) = μ 1 r + ω 1 t γ Γ ( γ + 1 ) and φ 1 , 2 r ( t ) = μ 2 r + ρ 1 t γ Γ ( γ + 1 ) into r e s 1 , 1 r ( t ) and r e s 1 , 2 r ( t ) of (11), that is,
r e s 1 , 1 r ( t ) = D 0 + γ φ 1 , 1 r ( t ) F 1 r ( t , φ 1 , 1 r ( t ) , φ 1 , 2 r ( t ) ) = ω 1 F 1 r ( t , μ 1 r + ω 1 t γ Γ ( γ + 1 ) ,   μ 2 r + ρ 1 t γ Γ ( γ + 1 ) ) , r e s 1 , 2 r ( t ) = D 0 + γ φ 1 , 2 r ( t ) F 2 r ( t , φ 1 , 1 r ( t ) , φ 1 , 2 r ( t ) ) = ρ 1 F 2 r ( t , μ 1 r + ω 1 t γ Γ ( γ + 1 ) ,   μ 2 r + ρ 1 t γ Γ ( γ + 1 ) ) .
Then, by solving the system r e s 1 , 1 r ( 0 ) = 0 and r e s 1 , 2 r ( 0 ) = 0 , we obtained ω 1 = F 1 r ( 0 ,   μ 1 r ,   μ 2 r ) and ρ 1 = F 2 r ( 0 ,   μ 1 r ,   μ 2 r ) . Thus, the 1st-FS approximated solutions for the system of FIVPs (3) can be written as:
φ 1 , 1 r ( t ) = μ 1 r + F 1 r ( 0 , μ 1 r ,   μ 2 r ) t γ Γ ( γ + 1 ) , φ 1 , 2 r ( t ) = μ 2 r + F 2 r ( 0 , μ 1 r ,   μ 2 r ) t γ Γ ( γ + 1 ) .
Similarly, to determine   ω 2 and ρ 2 , we set k = 2 in (10), then substituted φ 2 , 1 r ( t ) = μ 1 r + ω 1 t γ Γ ( γ + 1 ) + ω 2 t 2 γ Γ ( 2 γ + 1 ) , and φ 2 , 2 r ( t ) = μ 2 r + ρ 1 t γ Γ ( γ + 1 ) + ρ 2 t 2 γ Γ ( 2 γ + 1 ) into r e s 2 , 1 r ( t )   and   r e s 2 , 2 r ( t ) of Equation (11), as follows:
r e s 2 , 1 r ( t ) = D 0 + γ ( μ 1 r + ω 1 t γ Γ ( γ + 1 ) + ω 2 t 2 γ Γ ( 2 γ + 1 ) ) F 1 r ( t , μ 1 r + ω 1 t γ Γ ( γ + 1 ) + ω 2 t 2 γ Γ ( 2 γ + 1 ) , μ 2 r + ρ 1 t γ Γ ( γ + 1 ) + ρ 2 t 2 γ Γ ( 2 γ + 1 ) ) = ω 1 + ω 2 t γ Γ ( γ + 1 ) F 1 r ( t , μ 1 r + ω 1 t γ Γ ( γ + 1 ) + ω 2 t 2 γ Γ ( 2 γ + 1 ) , μ 2 r + ρ 1 t γ Γ ( γ + 1 ) + ρ 2 t 2 γ Γ ( 2 γ + 1 ) ) , r e s 1 , 2 r ( t ) = D 0 + γ ( μ 2 r + ρ 1 t γ Γ ( γ + 1 ) + ρ 2 t 2 γ Γ ( 2 γ + 1 ) ) F 2 r ( t , μ 1 r + ω 1 t γ Γ ( γ + 1 ) + ω 2 t 2 γ Γ ( 2 γ + 1 ) , μ 2 r + ρ 1 t γ Γ ( γ + 1 ) + ρ 2 t 2 γ Γ ( 2 γ + 1 ) ) = ρ 1 + ρ 2 t γ Γ ( γ + 1 ) F 2 r ( t , μ 1 r + ω 1 t γ Γ ( γ + 1 ) + ω 2 t 2 γ Γ ( 2 γ + 1 ) , μ 2 r + ρ 1 t γ Γ ( γ + 1 ) + ρ 2 t 2 γ Γ ( 2 γ + 1 ) ) .
By applying the operator D 0 + γ on both sides of (15), we obtained the γ -th Caputo fractional derivative of r e s 2 , 1 r ( t ) and r e s 2 , 2 r ( t ) and then we solved the obtained algebraic equations D 0 + γ r e s 2 , 1 r ( 0 ) = 0 and D 0 + γ r e s 2 , 2 r ( 0 ) = 0 , obtaining ω 2 = F 1 r ( 0 ,   ω 1 ,   ρ 1 ) and ρ 2 = F 2 r ( 0 ,   ω 1 ,   ρ 1 ) . Therefore, the 2nd-FS approximated solutions for the system of FIVPs (3) can be written as:
φ 2 , 1 r ( t ) = μ 1 r + F 1 r ( 0 , μ 1 r ,   μ 2 r ) t γ Γ ( γ + 1 ) + F 1 r ( 0 , ω 1 ,   ρ 1 ) t 2 γ Γ ( 2 γ + 1 ) , φ 2 , 2 r ( t ) = μ 2 r + F 2 r ( 0 , μ 1 r ,   μ 2 r ) t γ Γ ( γ + 1 ) + F 2 r ( 0 , ω 1 ,   ρ 1 ) t 2 γ Γ ( 2 γ + 1 ) .
Thirdly, to obtain the coefficients c 3 and d 3 , we considered k = 3 , in (10), then substituted φ 3 , 1 r ( t ) = μ 1 r + n = 1 k ω n t n γ Γ ( n γ + 1 ) , and φ 3 , 2 r ( t ) = μ 2 r + n = 1 k ρ n t n γ Γ ( n γ + 1 ) in r e s 3 , 1 r   and   r e s 3 , 2 r of (8); then, by computing D 0 +   2 γ r e s 3 , 1 r ( t ) and D 0 +   2 γ r e s 3 , 2 r ( t ) 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 ω 3 , and ρ 3 were obtained such that ω 3 = F 1 r ( 0 ,   ω 2 ,   ρ 2 ) and ρ 3 = F 2 r ( 0 ,   ω 2 ,   ρ 2 ) . Hence, the 3rd-FS approximated solutions for the system of FIVPs (3) can be summarized in the following expansions:
φ 3 , 1 r ( t ) = μ 1 r + F 1 r ( 0 ,   μ 1 r ,   μ 2 r ) t γ Γ ( γ + 1 ) + F 1 r ( 0 ,   ω 1 ,   ρ 1 ) t 2 γ Γ ( 2 γ + 1 ) + F 1 r ( 0 ,   ω 2 ,   ρ 2 ) t 3 γ Γ ( 3 γ + 1 )  
Using the same argument, the process can be repeated till the arbitrary order coefficients of the FS solutions for the system of FIVPs (3) are obtained. Hence, a higher degree of approximated solutions was achieved.

5. Numerical Experiments

In this section, we considered two FFIVPs of order γ to demonstrate the efficiency and applicability of our algorithm. Here, all the symbolic and numerical computations were performed by using Mathematica 12.
Example 1.
Consider the following FFIVPs:
{ D a + γ φ ( t ) = φ ( t ) , 0 < γ 1 , t [ 0 , 1 ] , φ ( 0 ) = μ ,
where μ = ( 2 , 3 , 4 ) is a fuzzy triangular number and has the r -level representations [ 2 + r , 4 r ] for r [ 0 , 1 ] . Based on Algorithm 1, the FFIVPs (18) will be transformed to one of the subsequent FIVPs systems:
Case 1:
The system of FIVPs corresponding to Caputo [(1)- γ ]-differentiable is
{ D 0 + γ φ 1 r ( t ) = φ 1 r ( t ) , D 0 + γ φ 2 r ( t ) = φ 2 r ( t ) ,   φ 1 r ( 0 ) = r + 2 ,   φ 2 r ( 0 ) = 4 r .
If γ = 1 , then the r -level representations of the exact solutions for the FIVPs system (19) are given by:
φ 1 r ( t ) = ( r + 2 ) e t , φ 2 r ( t ) = ( 4 r ) e t .
In light of the previous steps for the RPS algorithm, starting with φ 0 , 1 r ( 0 ) = r + 2 and φ 0 , 2 r ( 0 ) = 4 r , the k t h -residual functions r e s k , 1 r and r e s k , 2 r for (19) will be defined as:
r e s k , 1 r ( t ) = D 0 + γ φ k , 1 r ( t ) + φ k , 1 r ( t ) , r e s k , 2 r ( t ) = D 0 + γ φ k , 2 r ( t ) + φ k , 2 r ( t ) .
where φ k , 1 r and φ k , 2 r , indicating the k th-FS approximated solutions for (19), have the following forms:
φ k , 1 r ( t ) = ( r + 2 ) + n = 1 k ω n t n γ Γ ( n γ + 1 ) , φ k , 2 r ( t ) = ( 4 r ) + n = 1 k ρ n t n γ Γ ( n γ + 1 ) .
Now, to construct the 1st-FRPS-approximated solutions, consider k = 1 in the residual Equation (21) to obtain r e s 1 , 1 r ( t ) = D 0 + γ ( ( r + 2 ) + ω 1 t γ Γ ( γ + 1 ) ) + ( r + 2 ) + ω 1 t γ Γ ( γ + 1 ) = ω 1 + ( r + 2 ) + ω 1 t γ Γ ( γ + 1 ) and R e s 2 , 1 r ( t ) = D 0 + γ ( ( 4 r ) + ρ 1 t γ Γ ( γ + 1 ) ) + ( 4 r ) + ρ 1 t γ Γ ( γ + 1 ) = ρ 1 + ( 4 r ) + ρ 1 t γ Γ ( γ + 1 ) . Then, by using the fact r e s 1 , 1 r ( 0 ) = r e s 1 , 2 r ( 0 ) = 0 , we obtained ω 1 = ( r + 2 ) ,   ρ 1 = ( 4 r ) . So, the 1st-FS approximated solutions for the FIVPs (19) can be expressed as φ 1 , 1 r ( t ) = ( r + 2 ) ( r + 2 ) t γ Γ ( γ + 1 ) and φ 1 , 2 r ( t ) = ( 4 r ) ( 4 r ) t γ Γ ( γ + 1 ) .
Again, to find the 2nd-FS approximated solutions, put k = 2 , in (22), taking into account ω 1 = ( r + 2 ) ,   ρ 1 = ( 4 r ) and applying D 0 + γ in the resulting equations to obtain D 0 + γ r e s 2 , 1 r ( t ) = ω 2 ( r + 2 ) + ω 2 t γ Γ ( γ + 1 ) and D 0 + γ r e s 2 , 2 r ( t ) = ρ 2 ( 4 r ) + ρ 2 t γ Γ ( γ + 1 ) . Then, by considering that D 0 + γ r e s 2 , 1 r ( 0 ) = D 0 + γ r e s 2 , 2 r ( 0 ) = 0 , the coefficients ω 2 and ρ 2 will be obtained, such that ω 2 = r + 2 ,   ρ 2 = 4 r . Hence, the 2nd-FS approximated solutions could be given as φ 2 , 1 r ( t ) = ( r + 2 ) ( 1 t γ Γ ( γ + 1 ) + t 2 γ Γ ( 2 γ + 1 ) ) and φ 2 , 2 r ( t ) = ( 4 r ) ( 1 t γ Γ ( γ + 1 ) + t 2 γ Γ ( 2 γ + 1 ) ) .
Similarly, by computing the operator D 0 + 2 γ of the 3rd-residual functions, one can get D 0 + 2 γ r e s 3 , 1 r ( t ) = ω 3 + ( r + 2 ) + ω 3 t γ Γ ( γ + 1 ) and D 0 + 2 γ r e s 3 , 2 r ( t ) = ρ 3 + ( 4 r ) + ρ 3 t γ Γ ( γ + 1 ) . Then, by solving the resultant fractional equations at t = 0 , we obtained ω 3 = r + 1 ,   ρ 3 = 3 r . Therefore, the 3rd-FS approximated solutions could be given as φ 3 , 1 r ( t ) = ( r + 2 ) ( 1 t γ Γ ( γ + 1 ) + t 2 γ Γ ( 2 γ + 1 ) t 3 γ Γ ( 3 γ + 1 ) ) and φ 3 , 2 r ( t ) = ( 4 r ) ( 1 t γ Γ ( γ + 1 ) + t 2 γ Γ ( 2 γ + 1 ) t 3 γ Γ ( 3 γ + 1 ) ) .
Using the same approach for k = 4 and based on the fact that D 0 + 3 γ r e s 4 , 1 r ( 0 ) = D 0 + 3 γ r e s 4 , 2 r ( 0 ) = 0 , we obtained ω 4 = r + 1 ,   ρ 4 = 4 r . Depending on this, the 4th-FS approximated solutions can be written as φ 4 , 1 r ( t ) = ( r + 2 ) ( 1 t γ Γ ( γ + 1 ) + t 2 γ Γ ( 2 γ + 1 ) t 3 γ Γ ( 3 γ + 1 ) + t 4 γ Γ ( 4 γ + 1 ) ) and φ 4 , 2 r ( t ) = ( 4 r ) ( 1 t γ Γ ( γ + 1 ) + t 2 γ Γ ( 2 γ + 1 ) t 3 γ Γ ( 3 γ + 1 ) + t 4 γ Γ ( 4 γ + 1 ) ) . Moreover, depending on the fact that 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 = 5 , 6 , 7 , , the FS approximated solutions for (19) could be reformulated as:
φ 1 r ( t ) = ( r + 2 ) ( 1 t γ Γ ( γ + 1 ) + t 2 γ Γ ( 2 γ + 1 ) t 3 γ Γ ( 3 γ + 1 ) + t 4 γ Γ ( 4 γ + 1 ) ) = ( ( r + 2 ) n = 0 t n β Γ ( n β + 1 ) ) = ( r + 2 ) E γ ( t ) , φ 2 r ( t ) = ( 4 r ) ( 1 t γ Γ ( γ + 1 ) + t 2 γ Γ ( 2 γ + 1 ) t 3 γ Γ ( 3 γ + 1 ) + t 4 γ Γ ( 4 γ + 1 ) ) = ( ( r + 2 ) n = 0 t n β Γ ( n β + 1 ) ) = ( 4 r ) E γ ( t ) .
where E γ ( t ) is the Mittag–Leffler function.
In the case of γ = 1 , the FS expansions (23) could be reduced to the following forms:
φ 1 r ( t ) = ( r + 2 ) ( 1 t + t 2 2 t 3 3 ! + t 4 4 ! t 5 5 ! + t 6 6 ! ) = ( r + 2 ) e t ,   φ 2 r ( t ) = ( 4 r ) ( 1 t + t 2 2 t 3 3 ! + t 4 4 ! t 5 5 ! + t 6 6 ! ) = ( 4 r ) e t .
which coincide with the Taylor series expansions of the exact solutions φ 1 r ( t ) = ( r + 2 ) e t and φ 2 r ( t ) = ( 4 r ) e t .
Table 1 shows the lower and the upper bound solutions,   φ 1 r ( t ) and φ 2 r ( t ) , of the 7th-FS approximated solutions for FIVPs (19) for different values of γ , when r = 1 .
Case 2:
The system of FIVPs corresponding to Caputo [(2)- γ ]-differentiable is
{ D 0 + γ φ 1 r ( t ) = φ 2 r ( t ) , D 0 + γ φ 2 r ( t ) = φ 1 r ( t ) ,   φ 1 r ( 0 ) = r + 2 ,   φ 2 r ( 0 ) = 4 r .
If γ = 1 , then the r -level representations of the exact solution for the FIVPs system (25) are given by:
φ 1 r ( t ) = 3 e t + ( r 1 ) e t , φ 2 r ( t ) = 3 e t + ( 1 r ) e t   .
According to the RPS approach, starting with the 0th-FS approximated solutions φ 0 , 1 r ( 0 ) = r + 2 ,   φ 0 , 2 r ( 0 ) = 4 r , the k th-FS approximated solutions of (25) take the forms:
φ k , 1 r ( t ) = ( r + 2 ) + n = 1 k ω n t n γ Γ ( n γ + 1 ) , φ k , 2 r ( t ) = ( 4 r ) + n = 1 k ρ n t n γ Γ ( n γ + 1 ) .
Thus, the k th-residual functions of (25) will be
r e s k , 1 r ( t ) = D 0 + γ ( ( r + 2 ) + n = 1 k ω n t n γ Γ ( n γ + 1 ) ) + ( ( 4 r ) + n = 1 k ρ n t n γ Γ ( n γ + 1 ) ) , r e s k , 2 r ( t ) = D 0 + γ ( ( 4 r ) + n = 1 k ρ n t n γ Γ ( n γ + 1 ) ) + ( ( r + 2 ) + n = 1 k ω n t n γ Γ ( n γ + 1 ) )
To obtain the values of the coefficients ω n and ρ n , n = 1 , 2 , 3 , , k , in FS expansions (27), solve the algebraic fractional system in ω n and ρ n that was obtained considering 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 the RPS algorithm, the values of ω n , and ρ n , n = 1 , 2 , 3 , , k in (27) can be obtained as follows:
  • For k = 1 , we had r e s 1 , 1 r ( t ) = D 0 + γ ( ( r + 2 ) + ω 1 t γ Γ ( γ + 1 ) ) + ( 4 r ) + ρ 1 t γ Γ ( γ + 1 ) = ω 1 + ( 4 r ) + ρ 1 t γ Γ ( γ + 1 ) , and r e s 1 , 2 r ( t ) = D 0 + γ ( ( 4 r ) + ρ 1 t γ Γ ( γ + 1 ) ) + ( r + 2 ) + ω 1 t γ Γ ( γ + 1 ) = ρ 1 + ( r + 2 ) + ω 1 t γ Γ ( γ + 1 ) . Then, considering r e s 1 , 1 r ( 0 ) = r e s 1 , 2 r ( 0 ) = 0 , we obtained ω 1 = ( 4 r ) , ρ 1 = ( r + 2 ) .
  • For k = 2 , we had D 0 + γ r e s 2 , 1 r ( t ) = D 0 + γ ( D 0 + γ ( ( r + 2 ) ( 4 r ) t γ Γ ( γ + 1 ) + ω 2 t 2 γ Γ ( 2 γ + 1 ) ) + ( 4 r ) ( r + 2 ) t γ Γ ( γ + 1 ) + ρ 2 t 2 γ Γ ( 2 γ + 1 ) ) = ω 2 ( r + 2 ) + ρ 2 t γ Γ ( γ + 1 ) , and D 0 + γ r e s 2 , 2 r ( t ) = D 0 + γ ( D 0 + γ ( ( 4 r ) ( r + 2 ) t γ Γ ( γ + 1 ) + ρ 2 t 2 γ Γ ( 2 γ + 1 ) ) + ( r + 2 ) ( 4 r ) t γ Γ ( γ + 1 ) + ω 2 t 2 γ Γ ( 2 γ + 1 ) ) = ρ 2 ( 4 r ) + ω 2 t γ Γ ( γ + 1 ) . Lastly, by considering D 0 + γ r e s 2 , 1 r ( 0 ) = D 0 + γ r e s 2 , 2 r ( 0 ) = 0 , we obtained ω 2 = ( r + 2 ) , ρ 2 = ( 4 r ) .
  • For k = 3 , we had D 0 + 2 γ r e s 3 , 1 r ( t ) = D 0 + 2 γ ( D 0 + γ ( ( r + 2 ) ( 4 r ) t γ Γ ( γ + 1 ) + ( r + 2 ) t 2 γ Γ ( 2 γ + 1 ) + ω 3 t 3 γ Γ ( 3 γ + 1 ) ) + ( 4 r ) ( r + 2 ) t γ Γ ( γ + 1 ) + ( 4 r ) t 2 γ Γ ( 2 γ + 1 ) + ρ 3 t 3 γ Γ ( 3 γ + 1 ) ) = ω 3 + ( 4 r ) + ρ 3 t γ Γ ( γ + 1 ) , and D 0 + 2 γ R e s 3 , 2 r ( t ) = D 0 + 2 γ ( D 0 + γ ( ( 4 r ) ( r + 2 ) t γ Γ ( γ + 1 ) + ( 4 r ) t 2 γ Γ ( 2 γ + 1 ) + ρ 3 t 3 γ Γ ( 3 γ + 1 ) ) + ( r + 2 ) ( 4 r ) t γ Γ ( γ + 1 ) + ( r + 2 ) t 2 γ Γ ( 2 γ + 1 ) + ω 3 t 3 γ Γ ( 3 γ + 1 ) ) = ρ 3 + ( r + 2 ) + ω 3 t γ Γ ( γ + 1 ) . Thus, by considering D 0 + 2 γ r e s 2 , 1 r ( 0 ) = D 0 + 2 γ r e s 2 , 2 r ( 0 ) = 0 , we obtained ω 3 = ( 4 r ) , ρ 3 = ( r + 2 ) .
  • For k = 4 , we had D 0 + 3 γ r e s 4 , 1 r ( t ) = D 0 + 3 γ ( D 0 + γ ( ( r + 2 ) ( 4 r ) t γ Γ ( γ + 1 ) + ( r + 2 ) t 2 γ Γ ( 2 γ + 1 ) ( 4 r ) t γ Γ ( γ + 1 ) + ω 4 t 4 γ Γ ( 4 γ + 1 ) ) + ( 4 r ) ( r + 2 ) t γ Γ ( γ + 1 ) + ( 4 r ) t 2 γ Γ ( 2 γ + 1 ) ( r + 2 ) t 3 γ Γ ( 3 γ + 1 ) + ρ 4 t 4 γ Γ ( 4 γ + 1 ) ) = ω 4 ( r + 2 ) + ρ 4 t γ Γ ( γ + 1 ) , and D 0 + 3 γ r e s 4 , 2 r ( t ) = D 0 + 3 γ ( D 0 + γ ( ( 4 r ) ( r + 2 ) t γ Γ ( γ + 1 ) + ( 4 r ) t 2 γ Γ ( 2 γ + 1 ) ( r + 2 ) t 3 γ Γ ( 3 γ + 1 ) + ρ 4 t 4 γ Γ ( 4 γ + 1 ) ) + ( r + 2 ) ( 4 r ) t γ Γ ( γ + 1 ) + ( r + 2 ) t 2 γ Γ ( 2 γ + 1 ) + ( 4 r ) t γ Γ ( γ + 1 ) + ω 4 t 4 γ Γ ( 4 γ + 1 ) ) = ρ 4 ( 4 r ) + ω 4 t γ Γ ( γ + 1 ) . Thus, by considering D 0 + 3 γ r e s 4 , 2 r ( 0 ) = D 0 + 3 γ r e s 4 , 2 r ( 0 ) = 0 , we obtained ω 4 = ( r + 2 ) , ρ 4 = ( 4 r ) .
  • Likewise, for k = 5 and considering D 0 + 4 γ r e s 5 , 1 r ( 0 ) = D 0 + 4 γ r e s 5 , 2 r ( 0 ) = 0 , the coefficients ω 5 and ρ 5 will be obtained such that ω 5 = ( 4 r ) , ρ 5 = ( r + 2 ) .
  • Continuing with this procedure and based upon 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 , 8 , the 8th-FS approximated solutions for IVPs (25) were obtained:
    φ 8 , 1 r ( t ) = ( ( r + 2 ) ( 4 r ) t γ Γ ( γ + 1 ) + ( r + 2 ) t 2 γ Γ ( 2 γ + 1 ) ( 4 r ) t 3 γ Γ ( 3 γ + 1 ) + ( r + 2 ) t 4 γ Γ ( 4 γ + 1 ) ( 4 r ) t 5 γ Γ ( 5 γ + 1 ) + ( r + 2 ) t 6 γ Γ ( 6 γ + 1 ) ( 4 r ) t 7 γ Γ ( 7 γ + 1 ) + ( r + 2 ) t 8 γ Γ ( 8 γ + 1 ) ) , φ 8 , 2 r ( t ) = ( ( 4 r ) ( r + 2 ) t γ Γ ( γ + 1 ) + ( 4 r ) t 2 γ Γ ( 2 γ + 1 ) ( r + 2 ) t 3 γ Γ ( 3 γ + 1 ) + ( 4 r ) t 4 γ Γ ( 4 γ + 1 ) ( r + 2 ) t 5 γ Γ ( 5 γ + 1 ) + ( 4 r ) t 6 γ Γ ( 6 γ + 1 ) ( r + 2 ) t 7 γ Γ ( 7 γ + 1 ) + ( 4 r ) t 8 γ Γ ( 8 γ + 1 ) ) .
In particular, when γ = 1 , the FS approximated solutions for (25) could be expressed as
φ 1 r ( t ) = lim k φ k , 1 r ( t ) = ( ( r + 2 ) ( 4 r ) t + ( r + 2 ) t 2 2 ! ( 4 r ) t 3 3 ! + ( r + 2 ) t 4 4 ! + ) , φ 2 r ( t ) = lim k φ k , 2 r ( t ) = ( ( 4 r ) ( r + 2 ) t + ( 4 r ) t 2 2 ! ( r + 2 ) t 3 3 ! + ( 4 r ) t 4 4 ! + ) .
which agrees with the Maclurain series expansions of the exact solutions φ 1 r ( t ) = 3 e t + ( r 1 ) e t and φ 2 r ( t ) = 3 e t + ( 1 r ) e t .
Utilizing the RPS method, the numerical results of the fuzzy 8th-FS approximated solutions [ φ 8 , 1 r ( t ) , φ 8 , 2 r ( t ) ] are shown in Table 2 of Example 1, case 1, for different values of γ and a fixed value of the r -level. The effectiveness and reliability of the present method were also demonstrated via computing the absolute errors of the lower and upper approximated solutions and are presented in Table 3 of Example 1, case 2. From the table, we note the agreement between the obtained and the exact solutions at standard order γ = 1 .
The behavior of the fuzzy 8th-FS approximated solutions [ φ 8 , 1 r ( t ) , φ 8 , 2 r ( t ) ] of Example 1, case 1, with various values of fractional order γ are shown in Figure 1. The impact of the parameter r -level on the behavior of the lower and upper 8th-FS approximated solutions of Example 1, case 2, are illustrated in Figure 2. Moreover, the effect of the fractional order γ on the behavior of the lower and upper 8th-FS approximated solutions is demonstrated in Figure 3. Notice that, for different values of γ and r , the approximated solutions are continuously approaching to the exact solutions when γ = 1 . Therefore, we expect a veracious solution to such problems with various values of γ .
Example 2.
Consider the following FFIVPs:
{ D 0 + γ φ ( t ) = 2 t γ φ ( t ) + λ t γ   ,   0 < γ 1   ,   t [ 0 , 1 ] , φ ( 0 ) = μ ,
where λ and μ are two fuzzy triangular numbers and have the r -level representations [ r + 1 , 3 r ] .
Indeed, the FFIVPs (29) will be transformed to one of the subsequent systems with respect to type of Caputo differentiability:
Case 1: The system of FIVPs corresponding to Caputo [(1)- γ ]-differentiable is
{ D 0 + γ φ 1 r ( t ) = 2 t γ φ 1 r ( t ) + ( r + 1 ) t γ , D 0 + γ φ 2 r ( t ) = 2 t γ φ 2 r ( t ) + ( 3 r ) t γ   φ 1 r ( 0 ) = r + 1 ,   φ 2 r ( 0 ) = 3 r . ,
The exact solutions of FIVPs system (30) when γ = 1 could be obtained as
φ 1 r ( t ) = 1 2 ( r + 1 ) ( 3 e t 2 1 ) ,   φ 2 r ( t ) = 1 2 ( 3 r ) ( 3 e t 2 1 ) .
In view of the last described FS technique, we took into account φ 1 r ( 0 ) = ω 0 = r + 1 and φ 2 r ( 0 ) = ρ 0 = 3 r . Suppose that the k -th approximated solutions for FIVPs (30) have the following FS expansions form
φ k , 1 r ( t ) = ( r + 1 ) + n = 1 k ω n t n γ Γ ( n γ + 1 ) , φ k , 2 r ( t ) = ( 3 r ) + n = 1 k ρ n t n γ Γ ( n γ + 1 ) .
To determine ω n and ρ n ,   n = 1 , 2 , 3 , , k , we considered the solutions of D 0 + ( k 1 ) γ r e s k , 1 r ( 0 ) = D 0 + ( k 1 ) γ r e s k , 2 r ( 0 ) = 0 , k = 1 , 2 , 3 , , in which r e s k , 1 r and r e s k , 2 r are the k th residual functions of (30), defined as
r e s k , 1 r ( t ) = D 0 + γ φ k , 1 r ( t ) 2 t γ φ k , 1 r ( t ) ( r + 1 ) t γ , r e s k , 2 r ( t ) = D 0 + γ φ , k 2 r ( t ) 2 t γ φ k , 2 r ( t ) ( 3 r ) t γ .
Anyhow, by using the FS algorithm, the first few coefficients ω n and ρ n are:
ω 1 = 0 ,   ρ 1 = 0 , ω 2 = 3 ( r + 1 ) γ Γ ( γ ) ,   ρ 2 = 3 ( 3 r ) γ Γ ( γ ) , ω 3 = 0 ,   ρ 3 = 0 , ω 4 = 6 ( r + 1 ) γ Γ ( γ ) Γ ( 3 γ + 1 ) Γ ( 2 γ + 1 ) , ρ 4 = 6 ( 3 r ) γ Γ ( γ ) Γ ( 3 γ + 1 ) Γ ( 2 γ + 1 ) , ω 5 = 0 ,   ρ 5 = 0 , ω 6 = 12 ( r + 1 ) γ Γ ( γ ) Γ ( 3 γ + 1 ) Γ ( 5 γ + 1 ) Γ ( 2 γ + 1 ) Γ ( 4 γ + 1 ) , ρ 6 = 12 ( 3 r ) γ Γ ( γ ) Γ ( 3 γ + 1 ) Γ ( 5 γ + 1 ) Γ ( 2 γ + 1 ) Γ ( 4 γ + 1 ) , and so on .
Consequently, the 7th-FS approximated solutions of FFIVPs (30) can be represented as
φ 7 , 1 r ( t ) = ( r + 1 ) ( 1 + 3 γ Γ ( γ ) t 2 γ Γ ( 2 γ + 1 ) + 6 γ Γ ( γ ) Γ ( 3 γ + 1 ) t 4 γ Γ ( 2 γ + 1 ) Γ ( 4 γ + 1 ) + 12 γ Γ ( γ ) Γ ( 3 γ + 1 ) Γ ( 5 γ + 1 ) t 6 γ Γ ( 2 γ + 1 ) Γ ( 4 γ + 1 ) Γ ( 6 γ + 1 ) ) , φ 7 , 2 r ( t ) = ( 3 r ) ( 1 + 3 γ Γ ( γ ) t 2 γ Γ ( 2 γ + 1 ) + 6 γ Γ ( γ ) Γ ( 3 γ + 1 ) t 4 γ Γ ( 2 γ + 1 ) Γ ( 4 γ + 1 ) + 12 γ Γ ( γ ) Γ ( 3 γ + 1 ) Γ ( 5 γ + 1 ) t 6 γ Γ ( 2 γ + 1 ) Γ ( 4 γ + 1 ) Γ ( 6 γ + 1 ) ) .
For the particular case of γ = 1 , the FS-approximated solutions for (30) can be written as
φ 1 r ( t ) = lim k φ k , 1 r ( t ) = ( r 1 ) + 3 ( r 1 ) t 2 2 + 3 ( r 1 ) t 4 4 + 3 ( r 1 ) t 6 12 + , φ 2 r ( t ) = lim k φ k , 2 r ( t ) = ( 1 r ) + 3 ( 1 r ) t 2 2 + 3 ( 1 r ) t 4 4 + 3 ( 1 r ) t 6 12 + .
and are in agreement with the Taylor series expansions of the exact solutions
φ 1 r ( t ) = 1 2 ( r + 1 ) ( 3 e t 2 1 ) ,   and   φ 2 r ( t ) = 1 2 ( 3 r ) ( 3 e t 2 1 ) .
Case 2:
The system of FIVPs corresponding to Caputo [(2)- γ ]-differentiable is
{ D 0 + γ φ 1 r ( t ) = 2 t γ φ 2 r ( t ) + ( 3 r ) t γ , D 0 + γ φ 2 r ( t ) = 2 t γ φ 1 r ( t ) + ( r + 1 ) t γ ,   φ 1 r ( 0 ) = r + 1 ,   φ 2 r ( 0 ) = 3 r .
The exact solutions of FIVPs system (34) when γ = 1 could be obtained as
φ 1 r ( t ) = 1 2 ( 1 + r ) + 3 e t 2 + 1 2 ( r 1 ) e t 2 ,   φ 2 r ( t ) = 1 2 ( 3 r ) + 3 e t 2 + 1 2 ( 1 r ) e t 2 .
According to the application of the RPS approach, selecting φ 1 r ( 0 ) = r 1 and φ 2 r ( 0 ) = 1 r , we obtained the 0th-FS approximated solutions; then, the k th-truncated FS approximated solutions for FIVPs (34) are given by the following forms:
φ k , 1 r ( t ) = ( r + 1 ) + n = 1 k ω n t n γ Γ ( n γ + 1 ) , φ k , 2 r ( t ) = ( 3 r ) + n = 1 k ρ n t n γ Γ ( n γ + 1 ) .
Next, we defined the k th-residual functions r e s k , 1 r and r e s k , 2 r for (34) as follows:
r e s k , 1 r ( t ) = D 0 + γ ( ( r + 1 ) + n = 1 k ω n t n γ Γ ( n γ + 1 ) ) 2 t γ ( ( 3 r ) + n = 1 k ρ n t n γ Γ ( n γ + 1 ) ) ( 3 r ) t γ , r e s k , 2 r ( t ) = D 0 + γ ( ( 3 r ) + n = 1 k ρ n t n γ Γ ( n γ + 1 ) ) 2 t γ ( ( r + 1 ) + n = 1 k ω n t n γ Γ ( n γ + 1 ) ) ( r + 1 ) t γ .
Following the procedure of the RPS algorithm, the first few coefficients ω n and ρ n are
ω 1 = 0 ,   ρ 1 = 0 , ω 2 = 3 ( 3 r ) γ Γ ( γ ) ,   ρ 2 = 3 ( 1 + r ) γ Γ ( γ ) , ω 3 = 0 ,   ρ 3 = 0 , ω 4 = 6 ( 1 + r ) γ Γ ( γ ) Γ ( 3 γ + 1 ) Γ ( 2 γ + 1 ) , ρ 4 = 6 ( 3 r ) γ Γ ( γ ) Γ ( 3 γ + 1 ) Γ ( 2 γ + 1 ) , ω 5 = 0 ,   ρ 5 = 0 , ω 6 = 12 ( 3 r ) γ Γ ( γ ) Γ ( 3 γ + 1 ) Γ ( 5 γ + 1 ) Γ ( 2 γ + 1 ) Γ ( 4 γ + 1 ) ,   ρ 6 = 12 ( 1 + r ) γ Γ ( γ ) Γ ( 3 γ + 1 ) Γ ( 5 γ + 1 ) Γ ( 2 γ + 1 ) Γ ( 4 γ + 1 ) , ω 7 = 0 ,   ρ 7 = 0 , ω 8 = 24 ( 1 + r ) γ Γ ( γ ) Γ ( 3 γ + 1 ) Γ ( 5 γ + 1 ) Γ ( 7 γ + 1 ) Γ ( 2 γ + 1 ) Γ ( 4 γ + 1 ) Γ ( 6 γ + 1 ) , ρ 8 = 24 ( 3 r ) γ Γ ( γ ) Γ ( 3 γ + 1 ) Γ ( 5 γ + 1 ) Γ ( 7 γ + 1 ) Γ ( 2 γ + 1 ) Γ ( 4 γ + 1 ) Γ ( 6 γ + 1 ) , and so on .
Therefore, the 8th-FS approximated solutions of FIVPs (34) can be represented as
φ 8 , 1 r ( t ) = ( r + 1 ) + 3 γ Γ ( γ ) t 2 γ Γ ( 2 γ + 1 ) ( 3 r + 2 t 2 γ Γ ( 3 γ + 1 ) Γ ( 4 γ + 1 ) Γ ( 6 γ + 1 ) Γ ( 8 γ + 1 ) ( ( 1 + r ) Γ ( 6 γ + 1 ) Γ ( 8 γ + 1 ) + 2 t 2 γ Γ ( 5 γ + 1 ) ( 2 ( 1 + r ) t 2 γ Γ ( 7 γ + 1 ) + ( 3 r ) Γ ( 8 γ + 1 ) ) ) ) , φ 8 , 2 r ( t ) = ( 3 r ) + 3 γ Γ ( γ ) t 2 γ Γ ( 2 γ + 1 ) ( 1 + r + 2 t 2 γ Γ ( 3 γ + 1 ) Γ ( 4 γ + 1 ) Γ ( 6 γ + 1 ) Γ ( 8 γ + 1 ) ( ( 3 r ) Γ ( 6 γ + 1 ) Γ ( 8 γ + 1 ) + 2 t 2 γ Γ ( 5 γ + 1 ) ( 2 ( 3 r ) t 2 γ Γ ( 7 γ + 1 ) + ( 1 + r ) Γ ( 8 γ + 1 ) ) ) ) .
Correspondingly, for γ = 1 , the 8th-FS approximated solutions (38) can be written as
φ 8 , 1 r ( t ) = ( 1 + r ) + 3 2 t 2 ( 3 r + t 2 ( 29030400 ( 1 + r ) + 240 t 2 ( 40320 ( 3 r ) + 10080 ( 1 + r ) t 2 ) ) 58060800 ) , φ 8 , 2 r ( t ) = ( 3 r ) + 3 2 t 2 ( 1 + r + t 2 ( 29030400 ( 3 r ) + 240 t 2 ( 40320 ( 1 + r ) 10080 ( 3 + r ) t 2 ) ) 58060800 ) .
which agrees with the first eighth terms of the MacLaurin series of the exact forms φ 1 r ( t ) = 1 2 ( 1 + r ) + 3 e t 2 + 1 2 ( r 1 ) e t 2 , and φ 2 r ( t ) = 1 2 ( 3 r ) + 3 e t 2 + 1 2 ( 1 r ) e t 2 .
Table 4 presents the absolute errors of the obtained solutions by the RPS method for Example 2, case 2. The results in Table 4 show that the absolute errors of the proposed method were quite small. Further, numerical simulations of the outcomes for Example 2, case 1, were performed and are presented in Table 5.
Next, the numerical comparisons of the errors for Example 2 under Caputo [(1)-   γ ]-differentiability are discussed using our method and the homotopy analysis (HA) method [41] for different values of r , as shown in Table 6. From this table, one can observe that the RPS solutions were more accurate than the HA solutions.
The surface plots in Figure 4 show the 8th-FS approximated solution behavior at various values of γ for ( t ,   r ) [ 0 , 1 ] 2 , for Example 2, case 2. Figure 5 illustrates the effect of the parameter r on the obtained solutions against the exact solutions for Example 2, case 2, at standard order γ = 1 .

6. Conclusions

The major aim of the current analysis was to propose an efficient approach for obtaining fuzzy approximated solutions of FFIVPs with the assumption of strongly generalized differentiability. The main equations can be reformulated in parametric form and then translated into a fractional IVPs system to be solved by the RPS approach. This approach was applied directly by choosing suitable uncertain initial data to construct approximated solutions in the FS expansion with no need for nonphysical restrictive hypotheses. Numerical examples were provided to clarify the compatibility and reliability of the RPS approach. The graphical and numerical results satisfied the convex symmetric triangular fuzzy number and indicated that the proposed approach is an accurate instrument that can be used suitably as an alternative approach for constructing analytical solutions of diverse kinds of fuzzy fractional problems appearing in the fields of physics and engineering. In future work, it will be possible to apply the proposed approach for solving coupled systems of FFIVPs, fuzzy fractional BVPs with different order of γ , and fuzzy delay models.

Author Contributions

Conceptualization, Y.A.-q.; methodology, S.A.-O.; software, S.E.A.; validation, S.A.-O., S.E.A.; formal analysis, S.A.-O.; investigation, M.A.; writing-original draft preparation, M.A.; writing-review and editing, S.A.-O.; visualization; supervision, S.A.-O.; project administration, A.J.; funding acquisition, H.Q. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to thank the Deanship of Scientific Research at Umm Al-Qura University for supporting this work by Grant Code: (22UQU4282396DSR04).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Caputo, M. Linear models of dissipation whose Q is almost frequency independent: Part II. Geophys. J. Int. 1967, 13, 529–539. [Google Scholar] [CrossRef]
  2. Momani, S.; Freihat, A.; Al-Smadi, M. Analytical study of fractional-order multiple chaotic Fitzhugh-Nagumo neurons model using multistep generalized differential transform method. Abstr. Appl. Anal. 2014, 2014, 276279. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. 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]
  4. Altawallbeh, Z.; Al-Smadi, M.; Komashynska, I.; Ateiwi, A. Numerical solutions of fractional systems of two-point BVPs by using the iterative reproducing kernel algorithm. Ukr. Math. J. 2018, 70, 687–701. [Google Scholar] [CrossRef]
  5. Alabedalhadi, M.; Al-Smadi, M.; Al-Omari, S.; Baleanu, D.; Momani, S. Structure of optical soliton solution for nonliear resonant space-time Schrödinger equation in conformable sense with full nonlinearity term. Phys. Scr. 2020, 95, 105215. [Google Scholar] [CrossRef]
  6. Al-Smadi, M.; Abu Arqub, O.; Gaith, M. Numerical simulation of telegraph and Cattaneo fractional-type models using adaptive reproducing kernel framework. Math. Methods Appl. Sci. 2021, 44, 8472–8489. [Google Scholar] [CrossRef]
  7. Agarwal, R.P.; Lakshmikantham, V.; Nieto, J.J. On the concept of solution for fractional differential equations with uncertainty. Nonlinear Anal. Theory Methods Appl. 2010, 72, 2859–2862. [Google Scholar] [CrossRef]
  8. Al-Smadi, M.; Abu Arqub, O.; Zeidan, D. Fuzzy fractional differential equations under the Mittag-Leffler kernel differential operator of the ABC approach: Theorems and applications. Chaos Solitons Fractals 2021, 146, 110891. [Google Scholar] [CrossRef]
  9. Ahmadian, A.; Suleiman, M.; Salahshour, S.; Baleanu, D. A Jacobi operational matrix for solving a fuzzy linear fractional differential equation. Adv. Differ. Equ. 2013, 2013, 104–115. [Google Scholar] [CrossRef] [Green Version]
  10. Jafarian, A.; Golmankhaneh, A.K.; Baleanu, D. On fuzzy fractional Laplace transformation. Adv. Math. Phys. 2014, 2014, 295432. [Google Scholar] [CrossRef]
  11. Hasan, S.; Al-Smadi, M.; El-Ajou, A.; Momani, S.; Hadid, S.; Al-Zhour, Z. Numerical approach in the Hilbert space to solve a fuzzy Atangana-Baleanu fractional hybrid system. Chaos Solitons Fractals 2021, 143, 110506. [Google Scholar] [CrossRef]
  12. Al-Smadi, M.; Dutta, H.; Hasan, S.; Momani, S. On numerical approximation of Atangana-Baleanu-Caputo fractional integro-differential equations under uncertainty in Hilbert Space. Math. Model. Nat. Phenom. 2021, 16, 41. [Google Scholar] [CrossRef]
  13. Harrouche, N.; Momani, S.; Hasan, S.; Al-Smadi, M. Computational algorithm for solving drug pharmacokinetic model under uncertainty with nonsingular kernel type Caputo-Fabrizio fractional derivative. Alex. Eng. J. 2021, 60, 4347–4362. [Google Scholar] [CrossRef]
  14. Arshad, S.; Lupulescu, V. On the fractional differential equations with uncertainty. Nonlinear Anal. 2011, 74, 3685–3693. [Google Scholar] [CrossRef]
  15. Alikhani, R.; Bahrami, F. Global solutions for nonlinear fuzzy fractional integral and integrodifferential equations. Commun. Nonlinear Sci. Numer. Simul. 2013, 18, 2007–2017. [Google Scholar] [CrossRef]
  16. Mazandarani, M.; Najariyan, M. Type-2 fuzzy fractional derivatives. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 2354–2372. [Google Scholar] [CrossRef]
  17. Salahshour, S.; Allahviranloo, T.; Abbasbandy, S.; Baleanu, D. Existence and uniqueness results for fractional differential equations with uncertainty. Adv. Diff. Equ. 2012, 2012, 112. [Google Scholar] [CrossRef] [Green Version]
  18. Alaroud, M.; Al-Smadi, M.; Ahmad, R.R.; Salma Din, U.K. 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]
  19. 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. [Google Scholar] [CrossRef] [Green Version]
  20. Gumah, G.; Naser, M.F.M.; Al-Smadi, M.; Al-Omari, S.K.Q.; Baleanu, D. Numerical solutions of hybrid fuzzy differential equations in a Hilbert space. Appl. Numer. Math. 2020, 151, 402–412. [Google Scholar] [CrossRef]
  21. Alshammari, M.; Al-Smadi, M.; Abu Arqub, O.; Hashim, I.; Alias, M.A. Residual series representation algorithm for solving fuzzy duffing oscillator equations. Symmetry 2020, 12, 572. [Google Scholar] [CrossRef] [Green Version]
  22. Bede, B.; Gal, S.G. Generalizations of the differentiability of fuzzy number value functions with applications to fuzzy differential equations. Fuzzy Sets Syst. 2005, 151, 581–599. [Google Scholar] [CrossRef]
  23. Chalco-Cano, Y.; Román-Flores, H. On new solutions of fuzzy differential equations. Chaos Solitons Fractals 2008, 38, 112–119. [Google Scholar] [CrossRef]
  24. Kaleva, O. Fuzzy differential equations. Fuzzy Sets Syst. 1987, 24, 301–317. [Google Scholar] [CrossRef]
  25. 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]
  26. Al-Smadi, M.; Abu Arqub, O.; Hadid, S. Approximate solutions of nonlinear fractional Kundu-Eckhaus and coupled fractional massive Thirring equations emerging in quantum field theory using conformable residual power series method. Phys. Scr. 2020, 95, 10520. [Google Scholar] [CrossRef]
  27. Al-Smadi, M.; Abu Arqub, O.; Hadid, S. An attractive analytical technique for coupled system of fractional partial differential equations in shallow water waves with conformable derivative. Commun. Theor. Phys. 2020, 72, 085001. [Google Scholar] [CrossRef]
  28. Al-Smadi, M.; Abu Arqub, O.; Momani, S. Numerical computations of coupled fractional resonant Schrödinger equations arising in quantum mechanics under conformable fractional derivative sense. Phys. Scr. 2020, 95, 075218. [Google Scholar] [CrossRef]
  29. 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]
  30. Hasan, S.; Djeddi, N.; Al-Smadi, M.; Al-Omari, S.; Momani, S.; Pulga, A. Numerical solvability of generalized Bagley-Torvik fractional models under Caputo-Fabrizio derivative. Adv. Differ. Equ. 2021, 2021, 469. [Google Scholar] [CrossRef]
  31. Hasan, S.; Al-Smadi, M.; Freihet, A.; Momani, S. Two computational approaches for solving a fractional obstacle system in Hilbert space. Adv. Differ. Equ. 2019, 2019, 55. [Google Scholar] [CrossRef] [Green Version]
  32. Abu-Gaairi, R.; Hasan, S.; Al-Omari, S.; Al-Smadi, M.; Momani, S. Attractive multistep reproducing kernel approach for solving stiffness differential systems of ordinary differential equations and some error analysis. Compu. Model. Eng. Sci. 2022, 130, 299–313. [Google Scholar]
  33. Freihet, A.; Hasan, S.; Al-Smadi, M.; Gaith, M.; Momani, S. Construction of fractional power series solutions to fractional stiff system using residual functions algorithm. Adv. Differ. Equ. 2019, 2019, 95. [Google Scholar] [CrossRef] [Green Version]
  34. Freihet, A.; Hasan, S.; Alaroud, M.; Al-Smadi, M.; Ahmad, R.R.; Salma Din, U.K. Toward computational algorithm for time-fractional Fokker–Planck models. Adv. Mech. Eng. 2019, 11, 1687814019881039. [Google Scholar] [CrossRef]
  35. Al-Smadi, M.; Djeddi, N.; Momani, S.; Al-Omari, S.; Araci, S. An attractive numerical algorithm for solving nonlinear Caputo–Fabrizio fractional Abel differential equation in a Hilbert space. Adv. Differ. Equ. 2021, 2021, 271. [Google Scholar] [CrossRef]
  36. Moaddy, K.; Al-Smadi, M.; Hashim, I. A novel representation of the exact solution for differential algebraic equations system using residual power-series method. Discret. Dyn. Nat. Soc. 2015, 2015, 205207. [Google Scholar] [CrossRef]
  37. Saadeh, R.; Alaroud, M.; Al-Smadi, M.; Ahmad, R.R.; Salma Din, U.K. Application of fractional residual power series algorithm to solve Newell–Whitehead–Segel equation of fractional order. Symmetry 2019, 11, 1431. [Google Scholar] [CrossRef] [Green Version]
  38. Hasan, S.; Al-Smadi, M.; Dutta, H.; Momani, S.; Hadid, S. Multi-step reproducing kernel algorithm for solving Caputo–Fabrizio fractional stiff models arising in electric circuits. Soft Comput. 2022, 26, 3713–3727. [Google Scholar] [CrossRef]
  39. Aljarrah, H.; Alaroud, M.; Ishak, A.; Darus, M. Adaptation of residual-error series algorithm to handle fractional system of partial differential equations. Mathematics 2021, 9, 2868. [Google Scholar] [CrossRef]
  40. Allahviranloo, T.; Noeiaghdam, Z.; Noeiaghdam, S.; Nieto, J.J. A fuzzy method for solving fuzzy fractional differential equations based on the generalized fuzzy taylor expansion. Mathematics 2020, 8, 2166. [Google Scholar] [CrossRef]
  41. Abu Arqub, O.; Abo-Hammour, Z.; Al-Badarneh, R.; Momani, S. A reliable analytical method for solving higher-order initial value problems. Discret. Dyn. Nat. Soc. 2013, 2013, 673829. [Google Scholar] [CrossRef]
Figure 1. Surface plots of the fuzzy approximated solutions of Example 1, case 1, for all t [ 0 , 1 ] , and r [ 0 , 1 ] , at various values of γ : (a) γ = 1 ; (b) γ = 0.85 ; (c) γ = 0.65 ; (d) γ = 0.45 : (green and yellow are the lower and the upper solutions, respectively).
Figure 1. Surface plots of the fuzzy approximated solutions of Example 1, case 1, for all t [ 0 , 1 ] , and r [ 0 , 1 ] , at various values of γ : (a) γ = 1 ; (b) γ = 0.85 ; (c) γ = 0.65 ; (d) γ = 0.45 : (green and yellow are the lower and the upper solutions, respectively).
Symmetry 14 00804 g001
Figure 2. (a) Fuzzy 8th-FS approximated solutions [ φ 8 , 1 r ( t ) , φ 8 , 2 r ( t ) ] , case 1. (b) Fuzzy 8th-FS approximated solutions [ φ 8 , 1 r ( t ) , φ 8 , 2 r ( t ) ] , case 2, for Example 1 in parametric form, when γ = 1 .
Figure 2. (a) Fuzzy 8th-FS approximated solutions [ φ 8 , 1 r ( t ) , φ 8 , 2 r ( t ) ] , case 1. (b) Fuzzy 8th-FS approximated solutions [ φ 8 , 1 r ( t ) , φ 8 , 2 r ( t ) ] , case 2, for Example 1 in parametric form, when γ = 1 .
Symmetry 14 00804 g002
Figure 3. (a) Solution behavior of the Exact and the 8th-FS approximated solutions φ 8 , 1 r ( t ) , case 1; (a,b) Solution behavior of the Exact and the 8th-FS approximated solutions φ 8 , 2 r ( t ) , case 1; (c) Solution behavior of the Exact and the 8th-FS approximated solutions φ 8 , 1 r ( t ) , case 2; (d) Solution behavior of the Exact and the 8th-FS approximated solutions φ 8 , 2 r ( t ) , case 2, for Example 1, when r = 0.5 .
Figure 3. (a) Solution behavior of the Exact and the 8th-FS approximated solutions φ 8 , 1 r ( t ) , case 1; (a,b) Solution behavior of the Exact and the 8th-FS approximated solutions φ 8 , 2 r ( t ) , case 1; (c) Solution behavior of the Exact and the 8th-FS approximated solutions φ 8 , 1 r ( t ) , case 2; (d) Solution behavior of the Exact and the 8th-FS approximated solutions φ 8 , 2 r ( t ) , case 2, for Example 1, when r = 0.5 .
Symmetry 14 00804 g003
Figure 4. Surface plot of the fuzzy 7th-FS approximated solutions, of Example 2, case 2, for all t [ 0 , 1 ] and r [ 0 , 1 ] at various values of γ : (a) γ = 1 ; (b) γ = 0.9 ; (c) γ = 0.8 ; (d) γ = 0.7 ; (green and yellow are the lower and the upper solutions, respectively).
Figure 4. Surface plot of the fuzzy 7th-FS approximated solutions, of Example 2, case 2, for all t [ 0 , 1 ] and r [ 0 , 1 ] at various values of γ : (a) γ = 1 ; (b) γ = 0.9 ; (c) γ = 0.8 ; (d) γ = 0.7 ; (green and yellow are the lower and the upper solutions, respectively).
Symmetry 14 00804 g004
Figure 5. (a) Fuzzy approximated solutions [ φ 7 , 1 r ( t ) , φ 7 , 2 r ( t ) ] , case 1. (b) Fuzzy approximated solutions [ φ 8 , 1 r ( t ) , φ 8 , 2 r ( t ) ] , case 2, for Example 2, in parametric form, when γ = 1 .
Figure 5. (a) Fuzzy approximated solutions [ φ 7 , 1 r ( t ) , φ 7 , 2 r ( t ) ] , case 1. (b) Fuzzy approximated solutions [ φ 8 , 1 r ( t ) , φ 8 , 2 r ( t ) ] , case 2, for Example 2, in parametric form, when γ = 1 .
Symmetry 14 00804 g005
Table 1. The (1)-approximated solutions of Example 1, case 1, for different values of γ , with r = 1 .
Table 1. The (1)-approximated solutions of Example 1, case 1, for different values of γ , with r = 1 .
γ i   7th-FS Approximated Solutions, Case 1
1 4 φ 7 , 1 r ( t ) = φ 7 , 2 r ( t ) = 3 + 3 t + 2 t ( 3 + 2 t ) π 3 t 1 / 4 ( 5 + 4 t ) 5 Γ ( 5 4 ) 3 t 3 / 4 ( 7 + 4 t ) 7 Γ ( 7 4 ) .
1 2 φ 7 , 1 r ( t ) = φ 7 , 2 r ( t ) = 3 + 1 2 t ( 6 + t ( 3 + t ) ) 2 t ( 105 + 2 t ( 35 + 2 t ( 7 + 2 t ) ) ) 35 π .
3 4 φ 7 , 1 r ( t ) = φ 7 , 2 r ( t ) = 3 + t 3 2 + 4 t 3 / 2 ( 315 + 8 t 3 ) 315 π t 3 / 4 ( 1155 + 64 t 3 ) 385 Γ ( 7 4 ) t 9 / 4 ( 4641 + 64 t 3 ) 1547 Γ ( 13 4 ) .
1 φ 7 , 1 r ( t ) = φ 7 , 2 r ( t ) = 3 3 t + 3 t 2 2 t 3 2 + t 4 4 t 5 20 + t 6 120 t 7 840 .
Table 2. Numerical results of the 8th-FS approximated solutions [ φ 8 , 1 r ( t ) , φ 8 , 2 r ( t ) ] , with various values of γ and r , for Example 1, case 1.
Table 2. Numerical results of the 8th-FS approximated solutions [ φ 8 , 1 r ( t ) , φ 8 , 2 r ( t ) ] , with various values of γ and r , for Example 1, case 1.
r i   γ i   t i = 0.2   t i = 0.4   t i = 0.6  
0.5 1 [ 2.0468269 ,   2.8655576 ] [ 1.6758001 ,   2.3461202 ] [ 1.3720292 ,   1.9208408 ]
0.9 [ 1.9644060 ,   2.7501685 ] [ 1.6027476 ,   2.2438467 ] [ 1.3277858 ,   1.8589001 ]
0.8 [ 1.8765632 ,   2.6271884 ] [ 1.5358712 ,   2.1502197 ] [ 1.2939282 ,   1.8114995 ]
0.7 [ 1.7860931 ,   2.5005303 ] [ 1.4768338 ,   2.0675673 ] [ 1.2698638 ,   1.7778093 ]
1 1 [ 2.4561923 ,   2.4561923 ] [ 2.0109601 ,   2.0109601 ] [ 1.6464349 ,   1.6464349 ]
0.9 [ 2.3572873 ,   2.3572873 ] [ 1.9232971 ,   1.9232971 ] [ 1.5933429 ,   1.5933429 ]
0.8 [ 2.2518758 ,   2.2518758 ] [ 1.8430455 ,   1.8430455 ] [ 1.5527139 ,   1.5527139 ]
0.7 [ 2.1433117 ,   2.1433117 ] [ 1.7722006 ,   1.7722006 ] [ 1.5238366 ,   1.5238366 ]
Table 3. Absolute errors of Example 1, case 2, at n = 8 and various r values.
Table 3. Absolute errors of Example 1, case 2, at n = 8 and various r values.
t i r = 0 r = 0.5 r = 1
φ 1 r ( t ) 0.15 4.205524817 × 10 13 3.668176873 × 10 13 3.126388037 × 10 13
0.30 2.138851318 × 10 10 2.495841311 × 10 10 1.579714137 × 10 10
0.45 8.168253463 × 10 9 1.045280973 × 10 8 5.985222318 × 10 9
0.60 1.081084295 × 10 7 1.445286071 × 10 7 7.857506357 × 10 8
0.75 8.007310266 × 10 7 1.097021400 × 10 6 5.771710252 × 10 7
0.90 4.108761045 × 10 6 5.713634080 × 10 6 2.936552756 × 10 6
φ 2 r ( t ) 0.15 2.056133041 × 10 13 4.627409567 × 10 13 3.126388037 × 10 13
0.30 1.020570295 × 10 10 1.743996059 × 10 10 1.579714137 × 10 10
0.45 4.904169737 × 10 8 5.955703042 × 10 9 5.985222318 × 10 9
0.60 4.904169737 × 10 8 7.393833812 × 10 8 7.857506357 × 10 8
0.75 3.536110236 × 10 7 5.235650713 × 10 7 5.771710252 × 10 7
0.90 1.764344467 × 10 6 2.592500720 × 10 6 2.936552756 × 10 6
Table 4. Absolute errors of Example 2, case 2, at n = 8 and various r values.
Table 4. Absolute errors of Example 2, case 2, at n = 8 and various r values.
t i r = 0.25 r = 0.5 r = 0.75 r = 1
φ 1 r ( t ) 0.16 4.000000 × 10 10 4.000000 × 10 10 4.000000 × 10 10 2.000000 × 10 10
0.36 1.269000 × 10 6 1.158000 × 10 6 1.046000 × 10 6 9.340000 × 10 7
0.56 1.069870 × 10 4 9.798300 × 10 4 8.897900 × 10 5 7.997400 × 10 5
0.76 2.324840 × 10 3 2.141800 × 10 3 1.958750 × 10 3 1.775700 × 10 3
0.96 2.494130 × 10 2 2.314519 × 10 2 2.134908 × 10 2 1.955297 × 10 2
φ 2 r ( t ) 0.16 2.000000 × 10 10 3.000000 × 10 10 2.000000 × 10 10 2.000000 × 10 10
0.36 5.980000 × 10 7 7.100000 × 10 7 8.220000 × 10 7 9.340000 × 10 7
0.56 5.296200 × 10 5 6.196600 × 10 5 7.097000 × 10 5 7.997400 × 10 5
0.76 1.226540 × 10 3 1.409600 × 10 3 1.592650 × 10 3 1.775700 × 10 3
0.96 1.416464 × 10 2 1.596075 × 10 2 1.775686 × 10 2 1.955297 × 10 2
Table 5. Numerical results of the fuzzy 8th-FS approximated solutions [ φ 8 , 1 r ( t ) , φ 8 , 2 r ( t ) ] with various values of γ and r for Example 2, case 1.
Table 5. Numerical results of the fuzzy 8th-FS approximated solutions [ φ 8 , 1 r ( t ) , φ 8 , 2 r ( t ) ] with various values of γ and r for Example 2, case 1.
γ i   t i   r = 0.2   r = 0.6   r = 1.0  
1 0.15 [ 1.2409590 ,   2.8955710 ] [ 1.6546121 ,   2.4819181 ] [ 2.0682651 ,   2.0682651 ]
0.30 [ 1.3695087 ,   3.1955203 ] [ 1.8260116 ,   2.7390174 ] [ 2.2825145 ,   2.2825145 ]
0.45 [ 1.6038968 ,   3.7424258 ] [ 2.1385291 ,   3.2077935 ] [ 2.6731612 ,   2.6731612 ]
0.60 [ 1.9786368 ,   4.6168192 ] [ 2.6381824 ,   3.9572736 ] [ 3.2977280 ,   3.2977280 ]
0.85 0.15 [ 1.2900668 ,   3.0101558 ] [ 1.7200890 ,   2.5801335 ] [ 2.1501113 ,   2.1501113 ]
0.30 [ 1.5117406 ,   3.5273948 ] [ 2.0156542 ,   3.0234813 ] [ 2.5195677 ,   2.519567 ]
0.45 [ 1.8810863 ,   4.3892014 ] [ 2.5081151 ,   3.7621726 ] [ 3.1351438 ,   3.1351438 ]
0.60 [ 2.4495733 ,   5.7156709 ] [ 3.2660977 ,   4.8991466 ] [ 4.0826221 ,   4.0826221 ]
0.65 0.15 [ 1.4577914 ,   3.4015132 ] [ 1.9437218 ,   2.9155827 ] [ 2.4296523 ,   2.4296523 ]
0.30 [ 1.9263358 ,   4.4947835 ] [ 2.5684477 ,   3.8526716 ] [ 3.2105596 ,   3.2105596 ]
0.45 [ 2.6433967 ,   6.1679256 ] [ 3.5245289 ,   5.2867934 ] [ 4.4056612 ,   4.4056612 ]
0.60 [ 3.6985342 ,   8.6299130 ] [ 4.9313789 ,   7.3970683 ] [ 6.1642236 ,   6.1642236 ]
Table 6. Numerical comparison of the approximated solutions of Example 2, case 1, at t = 1 and different values of r .
Table 6. Numerical comparison of the approximated solutions of Example 2, case 1, at t = 1 and different values of r .
r i Absolute   Errors   of   φ 1 r ( t )
RPS methodHA method
0 4.09689 × 10 8 4.09689 × 10 5
0.2 3.20000 × 10 8 3.34322 × 10 5
0.4 2.40000 × 10 8 2.50742 × 10 5
0.6 1.60000 × 10 8 1.67161 × 10 5
0.8 8.0000 × 10 9 8.35806 × 10 6
1 0.0 0.0
r i Absolute   Errors   of   φ 1 r ( t )
RPS methodHA method
0 4.09689 × 10 8 3.34475 × 10 5
0.2 3.20000 × 10 8 2.67580 × 10 5
0.4 2.40000 × 10 8 2.00685 × 10 5
0.6 1.60000 × 10 8 1.33790 × 10 5
0.8 8.00000 × 10 9 6.68949 × 10 6
1 0.0 0.0
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Al-qudah, Y.; Alaroud, M.; Qoqazeh, H.; Jaradat, A.; Alhazmi, S.E.; Al-Omari, S. Approximate Analytic–Numeric Fuzzy Solutions of Fuzzy Fractional Equations Using a Residual Power Series Approach. Symmetry 2022, 14, 804. https://doi.org/10.3390/sym14040804

AMA Style

Al-qudah Y, Alaroud M, Qoqazeh H, Jaradat A, Alhazmi SE, Al-Omari S. Approximate Analytic–Numeric Fuzzy Solutions of Fuzzy Fractional Equations Using a Residual Power Series Approach. Symmetry. 2022; 14(4):804. https://doi.org/10.3390/sym14040804

Chicago/Turabian Style

Al-qudah, Yousef, Mohammed Alaroud, Hamza Qoqazeh, Ali Jaradat, Sharifah E. Alhazmi, and Shrideh Al-Omari. 2022. "Approximate Analytic–Numeric Fuzzy Solutions of Fuzzy Fractional Equations Using a Residual Power Series Approach" Symmetry 14, no. 4: 804. https://doi.org/10.3390/sym14040804

APA Style

Al-qudah, Y., Alaroud, M., Qoqazeh, H., Jaradat, A., Alhazmi, S. E., & Al-Omari, S. (2022). Approximate Analytic–Numeric Fuzzy Solutions of Fuzzy Fractional Equations Using a Residual Power Series Approach. Symmetry, 14(4), 804. https://doi.org/10.3390/sym14040804

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