1.
Introduction
Fractional calculus has been concerned with integration and differentiation of fractional (non-integer) order of the function. Riemann and Liouville defined the concept of fractional order intgro-differential equations [1]. Fractional calculus has developed an extensive attraction in current years in applied mathematics such as physics, medical, biology and engineering [2,3,4,5,6,7,8]. Whenever dealing with the fractional integro-differential equation many authors consider the terms Caputo fractional derivative, Riemann-Liouville and Grunwald-Letnikvo [9,10,11,12,13]. The subject fractional calculus has many applications in widespread and diverse field of science and engineering such as fractional dynamics in the trajectory control of redundant manipulators, viscoelasticity, electrochemistry, fluid mechanics, optics and signals processing etc.
Fractional integro-differential equations having some uncertainties in the form of boundary conditions, initial conditions and so on [14,15,16]. To resolve these type of uncertainties mathematicians introduced some concepts fuzzy set theory is one of them.
Zadeh introduced the concept of fuzzy set theory [17,18,19,20]. Later on Prade and Dubois [21,22], Nahmias [23], Tanaka and Mizumoto [24]. All of them experienced that the fuzzy number as a location of r-cut 0⩽r⩽1.
Many authors investigated some numerical techniques related to these problem which include the existence of the solution for discontinuous [25], reproducing kernel algorithm [26], integro-differential under generalized Caputo differentiability [27], A domain decomposition method [28], fractional differential transform method [29], Jacobi polynomial operational matrix [30], global solutions for nonlinear fuzzy equations [31], radioactivity decay model [32], Caputo-Katugampola fractional derivative approach [33], two-dimensional legendre wavelet method [34], fuzzy Laplace transform [35], fuzzy sumudu transform [36]. Further we can see [37,38,39,40]
Optimal Homotopy Asymptotic Method (OHAM) is one of the powerful techniques introduced by Marinca at al. [41,42,43] for approximate solution of differential equations. OHAM attracted an enormous importance in solving various problems in different field of science. Iqbal et al. applied this technique to Klein-Gordon equations and singular Lane-Emden type equation [44]. Sheikholeslami et al. used the proposed method for investigation of the laminar viscous flow and magneto hydrodynamic flow in a permeable channel [45]. Hashmi et al. obtained the solution of nonlinear Fredholm integral equations using OHAM [46]. Nawaz at al. applied the proposed method for solution of fractional order integro-differential equations [47], fractional order partial differential equations [48] and three-dimensional integral equations [49].
Aim of our study is to extend OHAM for solution of system of fuzzy Volterra integro differential equation of fractional order of the following form
with the given initial condition
Where Dαx represents the fuzzy fractional derivative in Caputo sense for fractional order of α with respect to x, h:[a,b]→RF is fuzzy valued function, k(x,t) is arbitrary kernel u0(x)∈RF is an unknown solution. RF represent set of all fuzzy valued function on real line.
The remaining paper is structured as follows: A brief overview on some elementary concept, notations and definitions of fuzzy calculus and fuzzy fractional calculus are discussed in section 2. Analysis of the technique is presented in section 3. Proposed method is applied to solve fuzzy fractional order Volterra integro-differential equations in section 4. Result and discussion of the paper is given in section 5 and section 6 is the conclusion of the paper.
2.
Preliminaries
In literature there exist various definitions of fuzzy calculus and fuzzy fractional calculus [50]. Some elementary concept, notations and definitions of fuzzy calculus and fuzzy fractional calculus related to this study are provided in this section.
Definition 2.1. The Riemann-Liouville fractional integral operator Iαx of order α is [50]:
Definition 2.2. Caputo partial fractional Derivative operator Dαx of order α with respect to x is defined as follow [50]:
which clearly shows that
Definition 2.3. A fuzzy number σ is a mapping σ:R→[0,1], satisfy the following property:
a. σ is normal that is, ∃x0∈R with u(x0)=1 [51,52].
b. σ is a convex fuzzy set that is, u(λx+(1−λ)y)⩾min{u(x),u(y)} for all x,y∈R, λ∈[0,1].
c. σ is upper semi-continuous in R.
d. ¯{x∈R:u(x)>0} is compact.
Definition 2.4. Parametric form of fuzzy number σ represented by an order pair (σ_,ˉσ) of the function (σ_(r),ˉσ(r)), satisfies the following conditions [52,53]:
a. σ_(r) is bounded monotonic increasing left continuous ∀r∈[0,1].
b. ˉσ(r) is bounded monotonic decreasing left continuous ∀r∈[0,1].
c. σ_(r)⩽ˉσ(r)∀r∈[0,1].
Definition 2.5. Addition and scalar multiplication of fuzzy number is given as:
a. (σ1⊕σ2)=(σ_1(r)+σ_2(r),ˉσ1(r)+ˉσ2(r))
b. (k⊗σ)={(σ_(r),ˉσ(r)),k⩾0,(σ_(r),ˉσ(r)),k<0.
Definition 2.6. A fuzzy real valued function σ1,σ2:[a,b]→R, then in [54]:
Definition 2.7. Assume u:[a,b]→RF. For every partition P={σ0,σ1,σ2,σ3,....,σn} and arbitrary ℓi:σi−1⩽ℓi⩽σi, 2⩽i⩽n consider
Rp=nΣi=2u(ℓj)(σi−σi−1). The definite integral of u(x) over [α,β] is
which show existence of limit in metric [55].
Definite integral exist if u(x) is continuous in metric D [51]:
3.
Application of OHAM
By considering definition 2.4. as discussed in section 2, Eq (1.1) becomes:
with the given initial condition
The homotopy of OHAM [41,42,43], constructed as follow:
where ρ∈[0,1], H(ρ)=∑m⩾1cmρm for all ρ≠0 is an auxiliary function, if ρ=0 then H(0)=0 where
and cm represent auxiliary constants. Using Taylor's series to expand υ(x,r;ρ) about ρ we get
Inserting Eq (3.4) into Eq (3.3) we get series of the problems by comparing the like power of ρ given as follow:
For calculating the constants c1,c2,c3..., mth order optimum solution becomes
Putting Eq (3.9) into Eq (3.1), we can found our residual given as follow:
If R(x,r;cl)=0, then u−m(x,r;cl)&ˉum(x,r;cl) will be the exact solutions.
Optimum solution contains some auxiliary constants; the optimal values of these constants are obtained through various techniques. In the present work, we have used the least square method [56,57]. The method of least squares is a powerful technique for obtaining the values of auxiliary constants. By putting the optimal values of these constants in Eq (8), we obtain the OHAM solution.
4.
Application and accuracy
Problem 4.1. Consider system of fuzzy fractional order Volterra integro-differential equation as [58]:
subject to the fuzzy initial condition [u(0)]r=[r−1,1−r], and for α=1 fuzzy fractional order Volterra integro-differential equations the exact solution is [u(x)]r=[r−1,1−r]Sinh(x) and 0⩽r⩽1.
By follow the technique as discussed in section 3, we get series of problems and their solutions as:
Their solutions are
Adding (4.6), (4.7), (4.8) and (4.9), one can construct u_(x,r) & \bar u(x, r) :
Values of {c_1}, \, \, {c_2} and {c_3} contain is in Eq (4.10)
Substituting the values from Table 1 into Eq (4.10), the approximate solutions for \underline u (x, r) & \bar u(x, r) at different values of \alpha taking r = 0.75 respectively is as follow
Substituting the values from Table 2 into Eq (4.10), the approximate solutions for \underline u (x, r) & \bar u(x, r) at different values of \alpha taking r = 0.5 respectively is as follow
Problem 4.2. Consider system of fuzzy fractional order Volterra integro-differential equation as [59]:
subject to the fuzzy initial condition {\left[ {u(0)} \right]^r} = \left[ {r - 1, 1 - r} \right], and the exact solution is \underline u (x, r) = \left( {r - 1} \right){E_{\alpha + 1}}\left( { - {t^{\alpha + 1}}} \right), \bar u(x, r) = \left( {1 - r} \right){E_{\alpha + 1}}\left( { - {t^{\alpha + 1}}} \right),
where {E_{\alpha + 1}} is a Mittag-Leffler function and 0 \leqslant r \leqslant 1.
By follow the technique as discussed in section 3, we get series of problems and their solutions as:
And their solutions are
Adding (4.24), (4.25), (4.26) and (4.27), one can construct \underline u (x, r) & \bar u(x, r) :
Values of {c_1}, \, \, {c_2} and {c_3} contain in Eq (4.28)
Substituting the values from Tables 3 and 4 into Eq (4.28), the approximate solutions for \underline u (x, r) & \bar u(x, r) at different values of \alpha taking r = 0.5 is as follow
Substituting the values from Tables 5 and 6 into Eq (4.28), the approximate solutions for \underline u (x, r) & \bar u(x, r) at different values of r taking \alpha = 0.5 is as follow
5.
Result and discussion
Tables 1–6 show the values of auxiliary constant at different values of r & \alpha for both lower and upper solution of OHAM for the solved problems. Tables 7 and 8 show the comparison of absolute error of 3rd order OHAM with Fractional Residual Power Series (FRPS) Method for 5-approximated solution and k = 5 for both lower and upper solutions of OHAM at different value of \alpha for problem 1. Comparison of absolute error of 3rd orders OHAM for both lower and upper solution of OHAM are shown in Tables 9 and 10. Numerical result show that OHAM provide more accuracy as compared to the other method and as \alpha \to 1 the approximate solution become very close to the exact solution. Graphical representation confirmed the convergence of fractional order solution towards the integer order solution. In Figure 1 graphical representation of OHAM at \alpha = 0.7, \, \, 0.8, \, \, 0.9\, , \, \, 1, \, \, r = 0.75 and \alpha = 0.7, \, \, 0.8, \, \, 0.9\, , \, \, 1, \, \, r = 0.50 are discussed for both \underline u (x, r) & \bar u(x, r) for problem 1. Figures 2 and 3 show the comparison of OHAM with the exact solution at different values of and taking r = 0.75 & r = 0.5 respectively for problem 1. Figure 4 represent the comparison of OHAM at \alpha = 0.2, \, 0.4, \, \, 0.6, \, \, 0.8, \, \, 1, \, \, r = 0.5 and r = 0, \, \, 0.2, \, 0.4, \, \, 0.6, \, \, 0.8, \, \, \alpha \, = 0.5 for both \underline u (x, r) and \bar u(x, r) for problem 2. Figure 5 shows the comparison of OHAM with the exact solution at different values of and r = 0.5 while Figure 6 shows the comparison of OHAM with the exact solution at different values of r and = 0.5 for problem 2.
6.
Conclusions
In the research paper, a powerful technique known as Optimal Homotopy Asymptotic Method (OHAM) has been extended to the solution of system of fuzzy integro differential equations of fractional order. The obtained results are quite interesting and are in good agreement with the exact solution. Two numerical equations are taken as test examples which show the behavior and reliability of the proposed method. The extension of OHAM to system of fuzzy integro differential equations of fractional order is more accurate and as a result this technique will more appealing for the researchers for finding out optimum solutions of system of fuzzy integro differential equations of fractional order.
Conflict of interest
The authors declare no conflict of interest.