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

: The modeling of fuzzy fractional integro-differential equations is a very signiﬁcant 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 classiﬁcation 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.


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: with the fuzzy initial conditions y(a) = y 0 (2) where λ is a parameter, g : [a, b] → R 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 ∈ R F and y(x) is unknown analytical fuzzy function to be determined.Throughout this paper, R F stands for the set of all fuzzy numbers on the real line R.
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.

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 R 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].
The parametric form or the r-cut representation of a fuzzy number ρ ∈ R F , is defined as: , and [ρ] r = supp(ρ), if r = 0. Clearly, the parametric form of ρ is closed and bounded interval [ρ 1r , ρ 2r ] in which ρ 1r is the lower r-cut representation of ρ, and ρ 2r is the upper r-cut representation of ρ.An equivalent parametric definition is also given in [13] as follows:
The metric structure of R F is defined by the Hausdorff distance mapping , for arbitrary fuzzy numbers ρ and ω.The metric (R 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) → R 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 ) ∈ R F such that either:
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 : Ω → R F be a continuous fuzzy-valued function.For each partition ℘ = x * 0 , x * 1 , . . ., x * n of Ω and for arbitrary point provided that the limit exists in the metric space (R 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. 1.
If y is (1)-differentiable on [a, b], then y 1r and y 2r are two differentiable functions on [a, b], and [D One can say y is Caputo fuzzy H-differentiable at x when D 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

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 1r (x), g 2r (x)] and [y(x)] r = [y 1r (x), y 2r (x)], respectively, as well as [y 1r (a), y 2r (a)] = [y 0,1r , y 0,2r ].Hence, one can write the FFVIDEs (1) and (2) in the r-cut representation as follows: where the r-cut representation of the function, [φ(x, t, y(t))] r is given by: [φ(x, t, y(t))] r = [φ 1r (x, t, y 1r (t), y 2r (t)), φ 2r (x, t, y 1r (t), y 2r (t))], in which The (n)-solution of fuzzy FVIDEs ( 1) and ( 2) is a function y -differentiable and satisfies the FFVIDEs ( 1) and ( 2).To compute it, we perform the next algorithm.

•
Step 2: Ensure that [y 1r (x), y 2r (x)] and D 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.

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 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 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: Since y 1r (x) and y 2r (x) satisfy the initial condition (5), then y 1r (0) = y 0,1r = c 0 and y 2r (0) = y 0,2r = d 0 will be the initial guesses for (5), so the series solutions can be written as Further, we can approximate y 1r (x) and y 2r (x) by the following kth-truncated series solutions: Define the so-called kth-residual functions Res k,1r and Res k,2r , for k = 1, 2, 3, . .., and the following residual functions Clear that, Res nr (x) = 0, for n = 1, 2 and for each x ≥ 0,which leads to Res k,nr (0) = 0, for each m = 0, 1, 2, . . ., k.However, According to applying FRPS technique to find the 1st-unknown coefficients c 1 and d 1 , substitute the 1st-approximated solutions y 1,1r (x) = c 0 + c 1 Γ(β+1) into the 1st residual functions Res 1,1r (x) and Res 1,2r (x) of ( 12) such that Using the facts Res 1,1r (0) = Res 1,2r (0) = 0, yields Again, to determine the 2nd-unknown coefficients c 2 and d 2 , substitute the 2nd-approximated solutions y 2,1r (x) = c 0 + c 1 Γ(2β+1) into the 2nd-residual functions Res 2,1r (x) and Res 2,2r (x) of ( 11) such that In the same way for the 3rd-unknown coefficients, c 3 and d 3 , substitute the 3rd-approximated solutions y 3,1r (x Γ(nβ+1) into the 3rd-residual functions Res 3,1r (x) and Res 3,2r (x) of ( 11) and then by computing D 2β 0 + Res 3,1r (x) and D 2β 0 + Res 3,2r (x) and using the facts 0) = 0, the coefficients, c 3 and d 3 will be given such that The process can be repeated, as can the coefficients, until an arbitrary order of the FPS solution for OFVIDEs (4) and ( 5) is obtained.

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
with the fuzzy initial condition y(0 Based on the type of differentiability, fuzzy FVIDEs ( 18) and ( 19) can be converted to one of the following systems: with the initial conditions y 1r (0) = 0, y 2r (0) = 0.
If β = 1, then the exact solutions of OFVIDEs ( 20) and ( 21) are given by Accordingly, the last description of FRPS algorithm, starting with y 0,1r (0) = y 0,2r (0) = 0, then the kth-approximated solutions for the system (20) and ( 21) are given by Consequently, the kth-residual functions Res k,1r , and Res k,2r for (20) will be 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.
+ 4x with the initial conditions y 1r (0) = 0, y 2r (0 If β = 1, then the exact solutions of OFVIDEs ( 25) and ( 26) are given by To apply the FRPS technique, suppose that the kth-truncated series solutions of y 1r (x) and y 2r (x) for the system ( 25) and ( 26) about x 0 = 0 have the following form: Apparently, the kth-residual functions Res k,1r , and Res k,2r for (25) will be 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 Tables 3 and 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 Figures 1 and 2, where core(y(x)) = {τ ∈ R : y(x)(τ) = 1}, and supp(y(x)) = {τ ∈ R : y(x)(τ) ≥ 0}.To show the fuzzy behavior of the proposed algorithm, the core and support of fuzzy ()approximated solutions,  = 1,2, and its Caputo derivative at  = 1 of Example 1, are plotted in Figure 1    Example 2. Consider the following FFIDE of Volterra type with the fuzzy initial condition y(0 Based on the type of differentiability, fuzzy FVIDEs ( 30) and ( 31) can be converted to one of the following systems: x 0 y 2r (t)dt, (32) with the fuzzy initial condition y 1r (0) = 0, y 2r (0 If β = 1, then the exact solution of OFVIDEs ( 32) and ( 33) is given by By utilizing the FRPS algorithm as presented previously, beginning with y 0,1r = y 1r (0) = 0, y 0,2r = y 2r (0) = 0, and depending on Equation ( 10), the kth-FRPS approximated solutions y k,1r (x) and y k,2r (x) for the system (32) and ( 33) are given by Consequently, the kth-residual functions Res k,1r , and Res k,2r for (32) will be 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
with the initial conditions y 1r (0) = 0, y 2r (0 If β = 1, then the exact solutions of OFVIDEs ( 38) and ( 39) are given by In view of the description for FRPS technique and starting with y 0,1r = y 1r (0) = 0, and y 0,2r = y 2r (0) = 0, then the kth-residual functions of Equation (37) will be as follows where y k,1r (x) and y k,2r (x) are given by 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 Tables 5 and 6 with different values of β.For the purpose of comparison, the absolute errors of y 1r (x) and y 2r (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). 1 Results of HW method referred in [35].
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}.
To show the fuzzy behavior of the proposed algorithm, the exact and 10 th -FRPS approximated solutions of Example 2, case1 and case 2, are plotted in Figure 3.The plots of -cut representation of exact and FRPS approximated solutions for Example 2, case 1 and case 2, are plotted in Figure 4 at  = 1,  = 10 and  ∈ {0,0.25,0.5,0.75,1}.

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.

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).

βa
+ y 1r (x), D β a + y 2r (x) are valid level sets on [a, b] or on a partial interval in [a, b].

Table 5 .
Numerical results for Example 2, case 1, with different values of β.

Table 6 .
Numerical results for Example 2, case 2, with different values of β.

Table 7 .
Comparisons of absolute errors of Example 2, case 1.