AFS Advances in Fuzzy Systems 1687-711X 1687-7101 Hindawi 10.1155/2019/6457548 6457548 Research Article Existence and Continuous Dependence on Initial Data of Solution for Initial Value Problem of Fuzzy Multiterm Fractional Differential Equation Choi Huichol 1 https://orcid.org/0000-0003-0436-180X Sin Kinam 1 Pak Sunae 1 So Sungryol 2 Mamlook Rustom M. 1 Faculty of Mathematics Kim Il Sung University Pyongyang Democratic People’s Republic of Korea ryongnamsan.edu.kp 2 Institute of Advanced Science Kim Il Sung University Pyongyang Democratic People’s Republic of Korea ryongnamsan.edu.kp 2019 1962019 2019 17 03 2019 13 05 2019 15 05 2019 1962019 2019 Copyright © 2019 Huichol Choi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

In this paper, the fuzzy multiterm fractional differential equation involving Caputo-type fuzzy fractional derivative of order 0<α<1 is considered. The uniqueness of solution is established by using the contraction mapping principle and the existence of solution is obtained by Schauder fixed point theorem.

1. Introduction

Nowadays the fractional differential equations (FDEs) are powerful tools representing many problems in various areas such as control engineering, diffusion processes, signal processing, and electromagnetism.

Recently the fuzzy fractional differential equations (FFDEs) have been studied by many researchers in order to analyze some systems with fuzzy initial conditions. But, it is too difficult to find the exact solutions of most FFDEs representing real-world phenomena. Therefore, research about FFDEs can be classified into two classes, namely, existence of solution and numerical methods. Many theoretical researches have been advanced on the existence, uniqueness, and stability of solution of FFDEs .

Also the analytical method and the numerical method are typical methods for solving FFDEs. The analytical method includes the Laplace transform method, monotone iterative method, variation of constant formula, and so on . Typical numerical methods are the operational matrix method, fractional Euler method, predictor-corrector method, and so on .

In , the existence and uniqueness of the solutions of fuzzy initial value problems of fractional differential equations with the Caputo-type fuzzy fractional derivative have been proved. Under the conditions which the right sides of equations satisfy Hölder continuity or Lipschitz continuity in its all variables, the existence of a solution to the Cauchy problem for fuzzy fractional differential equations was discussed in . In , by employing the contraction mapping principle on the complete metric space, the existence and uniqueness result for fuzzy fractional functional integral equation has been proved. The existence results of solutions for fuzzy fractional initial value problem under generalized differentiability conditions are obtained by Banach fixed point theorem in . In , researchers discussed the uniqueness and existence of the solutions for FFDEs with Riemann-Liouville H-differentiability of arbitrary order by using Krasnoselskii-Krein type conditions, Kooi type conditions, and Rogers conditions. But, the considerations of researchers in  were restricted to the case of FFDEs with single derivative term. Ngo et al.  presented that the existence and uniqueness results of the solution for fuzzy Caputo-Katugampola (CK) fractional differential equations with initial value and in  proved that the fractional fuzzy differential equation is not equal to the fractional fuzzy integral equation in general.

Based on the above facts, in this paper, we study the existence and uniqueness of solutions for fuzzy multiterm fractional differential equations of order 0<α<1 with fuzzy initial value under Caputo-type H-differentiability.

The paper is organized as follows. In Section 2, we introduced some definitions and properties of fuzzy fractional calculus. The existence result of solution for proposed problem is described in Section 3. Section 4 presented the continuous dependence on initial data of solution. Finally, the conclusion is summarized in Section 5.

2. Preliminaries and Basic Results

We introduce some definitions and notations which will be used throughout our paper.

Definition 1 (see [<xref ref-type="bibr" rid="B27">29</xref>]).

Let us denote by RF the class of fuzzy subsets u:R[0,1] satisfying the following properties:

u is normal, i.e., x0R for which u(x0)=1,

u is fuzzy convex, i.e., u(λx+(1-λ)y)min{u(x),u(y)} for all x,yR,   λ[0,1],

u is upper semicontinuous on R,

suppu=xRux>0 is the support of the u, and its closure cl(suppu) is compact.

Then RF is called the space of fuzzy number and any uRF is called fuzzy number.

We denote the r-cut form of fuzzy number uRF, 0r1, by [u]r[u1(r),u2(r)].

Also let us u,vRF. The metric d:RF×RF[0,+) on RF is defined as follows:(1)du,vsupr0,1maxu1r-v1r,u2r-v2r.

Definition 2 (see [<xref ref-type="bibr" rid="B19">19</xref>]).

Let u,vRF. If there exists wRF such that u=vw, then w is called the H-difference of uu and v, and it is denoted by w=uHv. Note that uHvu+(-1)v.

Definition 3 (see [<xref ref-type="bibr" rid="B19">19</xref>]).

Let f:(a,b)RF and x0(a,b). We say that f is H-differentiable at x0, if for h>0 sufficiently near to 0, there exist the H-differences f(x0+h)Hf(x0), f(x0)Hf(x0-h), and the limits(2)limh0+fx0+hHfx0h=limh0+fx0Hfx0-hh.

Then the limit is denoted by D(1)f(x0).

Theorem 4 (see [<xref ref-type="bibr" rid="B14">14</xref>]).

Let f:[a,b]RF be H-differentiable and [f(x)]r[f1(x,r),f2(x,r)]. Then f1(x,r),f2(x,r) are all differentiable and (3)D1fxr=f1x,r,f2x,r.

Definition 5 (see [<xref ref-type="bibr" rid="B14">14</xref>]).

A function f:[a,b]RF is said to be Riemann integrable on [a,b], if IRRF, ε>0, δ>0, for any division of [a,b], Δ:a=t0<t1<<tn=b with norm λ(Δ)<δ, and for any points ξi[ti,ti+1], i=0,1,,n-1,(4)di=0n-1fξi·ti+1-ti,IR<ε.

We denote the fuzzy Riemann integral of f from a to b by IR=abf(t)dt.

Lemma 6 (see [<xref ref-type="bibr" rid="B14">14</xref>]).

Suppose that f,g:[a,b]RF are continuous, then

(i) f is the fuzzy Riemann integrable on [a,b] and F(t)=atf(s)ds is differentiable as in Definition 3, namely, D(1)F(t)=f(t).

(ii) dabf(s)ds,abg(s)dsabdf(s),g(s)ds.

We introduce following notations:

C F ( I ) is the set of all continuous fuzzy-valued functions on I.

A C F ( I ) is the set of all absolutely continuous fuzzy-valued functions on I.

L F ( I ) is the space of all Lebesque integrable fuzzy-valued functions on I, where I=[0,L] and without losing generality, we promise that L=1.

Definition 7 (see [<xref ref-type="bibr" rid="B19">19</xref>]).

Let fCF(I)LF(I). The fuzzy Riemann-Liouville fractional integral of the fuzzy-valued function f is defined as follows:(5)I0+βfx=1Γβ0xx-sβ-1fsds,x>0.where I0+β is the Riemann-Liouville integral operator of β and Γ(β) is the Gamma function.

Lemma 8 (see [<xref ref-type="bibr" rid="B19">19</xref>]).

Let fCF(I)LF(I). Then Riemann-Liouville integral of the fuzzy-valued function f, based on its r-cut form, can be expressed as follows:(6)I0+βfxr=I0+βf1x,r,I0+βf2x,r,0r1,where (7)I0+βfix,r1Γβ0xfis,rx-s1-βds,i=1,2.

Lemma 9.

Let fCF(I), α,β>0. Then the following relations are satisfied:(8)I0+αI0+βfx=I0+βI0+αfx=I0+α+βfx.

Proof.

Let us denote the r-cut form of f by [f(x)]r[f1(x,r),f2(x,r)].

Then we have(9)I0+αI0+βfxr=I0+αI0+βf1x,r,I0+αI0+βf2x,r=I0+βI0+αf1x,r,I0+βI0+αf2x,r=I0+βI0+αfxr.

Moreover, since fCF(I), we get(10)I0+αI0+βfxr=I0+α+βf1x,r,I0+α+βf2x,r=I0+α+βfxr.

Definition 10 (see [<xref ref-type="bibr" rid="B19">19</xref>]).

Let 0<β<1. We say that f is fuzzy Riemann-Liouville H-differentiable of order β if (11)I0+1-βfx=1Γ1-β0xfsx-sβds,x>0is H-differentiable. Then fuzzy Riemann-Liouville H-derivative of order β of function f is denoted by D0+βRLf(x)D(1)I0+1-βf(x).

Definition 11 (see [<xref ref-type="bibr" rid="B19">19</xref>]).

Let 0<β<1. We say that f is a fuzzy Caputo-type differentiable function if H-difference f(x)Hf(0) exists and I0+1-β(f(x)Hf(0))ACF(I) satisfies. Then fuzzy Caputo-type derivative of order β of function f is denoted by(12)D0+βcfxD0+βRLfxHf0,x>0.

Lemma 12 (see [<xref ref-type="bibr" rid="B19">19</xref>]).

Let fCF(I)LF(I) and [f(x)]r[f1(x,r),f2(x,r)] for r[0,1]. If f is a fuzzy Caputo-type fractional differentiable function, then(13)D0+βcfxr=I0+1-βDfxr=I0+1-βDf1x,r,I0+1-βDf2x,r.

Lemma 13.

Let f:[a,b]RF be H-differentiable. Then the following relations hold:

fCF(0,1]D0+βcI0+βf(x)=f(x).

fCF(0,1]D0+βRLI0+βf(x)=f(x).

0<λ<β, fCF(0,1]D0+λcI0+βf(x)=I0+β-λf(x).

Proof.

First we prove (i). From the assumption of Lemma 13, we have (14)D0+βcI0+βfx=D0+βRLI0+βfxHI0+βf0=D1I0+1-βI0+βfxHI0+βf0,=D1I0+1-βI0+βfx=D1I0+1fx.

By Lemma 6 (i), we get (15)D1I0+1fx=fx.

From the result (i) of lemma, it is obvious that (ii) holds.

Next let prove (iii). (16)D0+λcI0+βfx=D0+λRLI0+βfxHI0+βf0=D1I0+1-λI0+βfxHI0+βf0=D1I0+1-λI0+βfx=D1I0+1+β-λfx.

By Lemma 6 (i), we obtain (17)D1I0+1+β-λfx=I0+β-λfx.

Lemma 14.

The following facts are true:

(i) Let zCF(I). For any positive number σ, the fractional integral I0+σz(x) is continuous in x.

(ii) For any positive numbers σ and k, it holds that I0+σekxekx/kα.

(iii) Let u,vCF(I). For any positive number σ, it holds that(18)dI0+σux,I0+σvxI0+σdux,vx.

Proof.

Let us consider the assertion (i). For any x0(0,1], it is enough to prove that (19)limxx0dI0+σzx,I0+σzx0=0.

We use the notation [z(x)]r[z1(x,r),z2(x,r)] for r-cut representation of z(x).

Since zCF(I), z1(r,x), and z2(r,x) are continuous in x. And by Lemma 8, the following expression holds:(20)dI0+σzx,I0+σzx0=supr0,1maxI0+σz1x,r-I0+σz1x0,r,I0+σz2x,r-I0+σz2x0,r.

It is well known that I0+σzi(x,r), i=1,2, are continuous with respective to x. Thus we can see that limxx0d(I0+σz(x),I0+σz(x0))=0.

Now we prove (ii).

From the definition of fractional integral, the following is true:(21)I0+αekt=1Γα0tt-sα-1eksds=1Γα0tzα-1ekt-zdz=ektΓα0tzα-1e-ktdz=ektΓα0ktskα-1e-s1kds=ektΓα·kα·0ktsα-1e-sds

By the definition of Gamma function, we have(22)ektΓα·kα·0sα-1e-sds=ektΓα·kαΓα=ektkα.

Let us consider assertion (iii).

We use the notations [u(x)]r[u1(x,r),u2(x,r)],[v(x)]r[v1(x,r),v2(x,r)] for r-cut representations of u(x),v(x), respectively.

Since u,vCF(I), u1(x,r),u2(x,r),v1(x,r),v2(x,r) are continuous in x. Therefore the following evaluations are true:(23)dI0+σux,I0+σvx=supr0,1maxI0+σu1x,r-I0+σv1x,r,I0+σu2x,r-I0+σv2x,r=supr0,1maxI0+σu1x,r-v1x,r,I0+σu2x,r-v2x,rsupr0,1maxI0+σu1x,r-v1x,r,I0+σu2x,r-v2x,rI0+σsupr0,1maxu1x,r-v1x,r,u2x,r-v2x,r=I0+σdz1x,z2x.

The proof is completed.

3. Existence of Solution for Fuzzy Multiterm Fractional Differential Equation

Let us consider the existence of solution for initial value problem of following fuzzy multiterm fractional differential equation: (24)D0+αcyx=fx,yx,D0+βcyx,x0,1,y0=y0,0<β<α<1,where f:I×RF×RFRF, y0RF, and D0+αc is fuzzy Caputo-type derivative.

Definition 15.

Let y:IRF. We say that y is the solution of initial value problem (24) if D0+β  cy(x)CF(I) holds and y satisfies (24).

Now let us consider the following:(25)D0+αcyx=zx,x0,1,y0=y0,y0RF,where z(x)CF(I).

Lemma 16.

The solution of initial value problem of fuzzy fractional differential equation (25) is represented as (26)yx=y01Γα0xzsx-s1-αds.

Proof.

Let y be the solution of initial value problem (25). Then we have (27)D0+αcyxzx,xI.

Also since z is the fuzzy continuous, it is the fuzzy integrable and I0+αz(x) exists for xI.

Therefore the following relations hold:(28)I0+αD0+αcyxI0+αzx,xI,I0+αD0+αRLyxHy0I0+αzx,I0+αD1I0+1-αyxHy0I0+αzx.

From the Caputo-type differentiability of fuzzy-valued function y, we get (29)I0+1-αyxHy0ACFIand since the space of the absolutely continuous functions coincides with the space of primitive functions of Lebesque integrable functions, the following relation is satisfied: (30)φLFI;I0+1-αyxHy0=I0+1φx.

Therefore we obtain(31)yxHy0=I0+αφx.

By (30) and Lemma 6 (i), the left side of (28) exchanges as (32)I0+αD1I0+1-αyxHy0=I0+αD1I0+1φx=I0+αφx.

From (28) and (31), y(x)Hy0=I0+αz(x) is satisfied. Namely, we have (33)yx=y0I0+αzx.

Conversely, we prove that y denoted by (26) is the solution of fuzzy initial value problem (25). Since (34)yx=y01Γα0xzsx-s1-αds=y0I0+αzx,y(x)Hy0=I0+αz(x) holds. Also as zCF(I), we get(35)I0+1-αyxHy0=I0+1zx.

On the other hand, when [z(x)]r[z1(x,r),z2(x,r)], by Lemma 8(36)I0+1zxr=0xz1s,rds,0xz2s,rdsholds r[0,1]. Therefore (37)I0+1zx+hr=0x+hz1s,rds,0x+hz2s,rds=0xz1s,rds+xx+hz1s,rds,0xz2s,rds+xx+hz2s,rds=0xz1s,rds,0xz2s,rds+xx+hz1s,rds,xx+hz2s,rds.

Since the interval family xx+hz1(s,r)ds,xx+hz2(s,r)dsr generates obviously a fuzzy number, we can see that there exist the H-differences as (38)I0+1zx+hHI0+1zx,I0+1zxHI0+1zx-h.

Consequently the H-differentiability of y is leaded.

From (35), we obtain (39)D1I0+1-αyxHy0=D1I0+1zx=zx.

Also it is obvious that y satisfies the initial condition.

Theorem 17.

Let f of (24) be a fuzzy continuous with respect to every variable. If y(x) is the solution of initial value problem (24), the fuzzy-valued function z(x) which is constructed by z(x)D0+αcy(x) is the solution in CF(I) of fuzzy integral equation as(40)zx=fx,y01Γα0xzsx-s1-αds,I0+α-βzx,x0,1.Conversely if z(x) is the solution in CF(I) of fuzzy integral equation (40), y(x) which is constructed by (26) is the solution of initial value problem (24).

Proof.

Let y(x) be the solution of initial value problem (24). Namely, let assume that y(x) satisfies (41)D0+αcyx=fx,yx,D0+βcyx,x0,1,y0=y0,y0RF.

Then if z(x) is denoted by z(x)D1αcy(x), by Lemma 16, the following relation is leaded:(42)yx=y01Γα0xzsx-s1-αds=y0I0+αzx.

Applying the operator D0+βc to the both side of above equation, by Lemma 13 (iii), we get (43)D0+βcyx=I0+α-βzx.

Therefore substituting the above results to D0+αcy(x), D0+βcy(x),y(x) of (41), we obtain (40).

Next let z(x) be the solution of fuzzy integral equation (40). Then it is obvious that (44)yx=y0I0+αzxsatisfies the initial condition of problem (24) from the continuousness of z(x).

Namely (45)y0=y0.

Consequently from (44)(46)yxHy0=yxHy0=I0+αzxis leaded and regarding the above equation, we get(47)I0+1-αyxHy0=I0+1-αI0+αzx=I0+1zx,D1I0+1-αyxHy0=D1I0+1zx=zx,D0+αcyx=zx=fx,ys,I0+α-βzx.

Moreover since (48)D0+βcyx=D0+βcI0+αzxy0=D0+βcI0+αzx=I0+α-βzx,we obtain (49)D0+αcyx=fx,ys,D0+βcyx.

Now we employ the following metric structure in CF(I):(50)u,vCFI,du,vmaxtIdut,vt.

Obviously we can see that (CF(I),d) is a complete metric space (see ).

For any positive number k, we can consider the metric structure as(51)u,vCFI,dku,vmaxtIe-ktdut,vt.

Then the metric dk is equivalent to the metric d. Namely, (52)M,m>0;u,vCFI,mdku,vdu,vMdku,v.

Theorem 18.

Assume that the function f in (40) is continuous in its all variables and especially, for any y1,y2,z1,z2RF, f satisfies the following condition:(53)dfx,y1,z1,fx,y2,z2L1·dy1,y2+L2·dz1,z2.Then the fuzzy integral equation (40) has a unique solution.

Proof.

Since α,α-β>0, there exists k>0 that inequality L11/kα+L21/kα-β<1 is true.

Therefore for any k satisfying this inequality, we put as follows:(54)qL11kα+L21kα-β.

Also we define the operator T by (55)Tzxfx,y01Γα0xzsx-s1-αds,I0+α-βzx.

For any zCF(I), by (i) of Lemma 14, I0+αz(x),I0+α-βz(x) are continuous and f is continuous by assumptions of theorem. Thus the operator T is a map from CF(I) to CF(I).

Thus the fuzzy integral equation (40) is represented as (56)z=Tz,zCFI.

The existence of solution for the fuzzy integral equation (40) is equivalent to the existence of the fixed point of the operator T in CF(I).

For any z1,z2CF(0,1], let us evaluate d(Tz1(x),Tz2(x)).(57)dTz1x,Tz2x=dfx,y01Γα0xz1sx-s1-αds,I0+α-βz1x,fx,y01Γα0xz2sx-s1-αds,I0+α-βz2xL1dy01Γα0xz1sx-s1-αds,y01Γα0xz2sx-s1-αds+L2·dI0+α-βz1x,I0+α-βz2xL1d1Γα0xz1sx-s1-αds,1Γα0xz2sx-s1-αds+L2·dI0+α-βz1x,I0+α-βz2xBy (iii) of Lemma 14, we have(58)L11Γα0xdz1s,z2sx-s1-αds+L2·I0+α-βdz1x,z2x=L11Γα0xe-kseksdz1s,z2sx-s1-αds+L2·I0+α-βe-kxekxdz1x,z2xL11Γα0xeksx-s1-αds·dkz1,z2+L2·I0+α-βeks·dkz1,z2By (ii) of Lemma 14, we have(59)L1ekxkα+L2·ekxkα-β·dkz1,z2.

Consequently the following inequality is obtained:(60)dTz1x,Tz2xL1ekxkα+L2·ekxkα-β·dkz1,z2.

Multiplying the both sides of the above equation by e-kx, we obtain that(61)e-kxdTz1x,Tz2xL11kα+L2·1kα-β·dkz1,z2and(62)dkTz1,Tz2L11kα+L2·1kα-β·dkz1,z2=q·dkz1,z2.

Thus the operator T is contractive on CF(I) with respect to the distance dk and we obtain the unique fixed point CF(I) of the operator T by contraction mapping principle. By the way, since the distance dk is equivalent to d in CF(I), z is also the unique fixed point in sense of the distance d. This completes the proof of theorem.

Next let us consider the existence of solution in case which have not satisfied the Lipschitz condition.

Lemma 19 (Schauder fixed point theorem).

Assume that (E,d) is the complete metric space, U is a nonempty convex closed subset of E, and A is a continuous mapping of U into itself such that A(U) is contained in a compact subset of U, then A has a fixed point in U.

Theorem 20.

Suppose that the followings conditions are satisfied:

(i) r>d(f(·,y0,0^),y0),L1,L2,L3>0,x1,x2(0,1],y1,y2Ur(y0),z1,z2Ur(0^), (63)dfx1,y1,z1,fx2,y2,z2L1·dy1,y2+L2·dz1,z2+L3·x1-x2,where Ur(y0)ud(u,y0)r,uCF(I) and 0^ is zero fuzzy number.

(ii) qL11/Γ(α+1)+L2·1/Γ(α-β+1)<1.

Let (64)xmaxxL1·xα+L2·xα-βminΓα+1,Γα-β+1r-dfx,y0,0^,y0r,xαΓα+11,xαΓα-β+11and(65)U-ry0udu,y0r,uCF0,x,U-r0^udu,0^r,uCF0,x.Then the following results hold:

T:U-r(0^)U-r(0^),

TuuU-r0^ is relatively compact.

Proof.

Firstly, let us prove that T:U-r(0^)U-r(0^).

We have that, for any uU-r(0^),(66)dy01Γα0xusx-s1-αds,y0d1Γα0xusx-s1-αds,0^1Γα0xdus,0^x-s1-αdsxαdu,0^Γα+1r.

Thus we can get that y01/Γ(α)0xu(s)/(x-s)1-αdsU-r(y0).

Also by (iii) of Lemma 14, we have that, for any uU-r(0^),(67)dI0+α-βux,0^I0+α-βdux,0^xα-βΓα-β+1dux,0^r.

Thus we can see that I0+α-βu(x)U-r(0^).

Therefore we can get that for any uU-r(0^),(68)dTux,y0=dfx,y01Γα0xusx-s1-αds,I0+α-βux,y0dfx,y01Γα0xusx-s1-αds,I0+α-βux,fx,y0,0^+dfx,y0,0^,y0L1·dy01Γα0xusx-s1-αds,y0+L2·dI0+α-βux,0^+dfx,y0,0^,y0L1·d1Γα0xusx-s1-αds,0^+L2·dI0+α-βux,0^+dfx,y0,0^,y0 By (iii) of Lemma 14, we have(69)L11Γα0xdus,0^x-s1-αds+L2·I0+α-βdux,0^+dfx,y0,0^,y0L1·xα·rΓα+1+L2·xα-β·rΓα-β+1+dfx,y0,0^,y0L1·xα·rminΓα+1,Γα-β+1+L2·xα-β·rminΓα+1,Γα-β+1+dfx,y0,0^,y0r.

Therefore the operator T is a continuous mapping from convex closed subset U-r(0^) into itself; namely,(70)T:U-r0^U-r0^.

Next let us prove (ii). The uniformly boundness of TuuU-r0^ is obvious from (70).

Let us consider the equicontinuity of TuuU-r0^. For any yTuuU-r0^, there exists zU-r(0^) which satisfies y=Tz.

For x1,x2(0,x], we estimate d(y(x1),y(x2))=d(Tz(x1),Tz(x2)).

Without losing generality, let x1<x2x.

Then we have (71) d T z x 1 , T z x 2 = d f x 1 , y 0 1 Γ α 0 x 1 z s x 1 - s 1 - α d s , I 0 + α - β z x x = x 1 , (72) f x 2 , y 0 1 Γ α 0 x 2 z s x 2 - s 1 - α d s , I 0 + α - β z x x = x 2 = d f x 1 , y 0 1 Γ α 0 x 1 z s x 1 - s 1 - α d s , 1 Γ α - β 0 x 1 z s x 1 - s 1 - α - β d s , f x 2 , y 0 1 Γ α 0 x 2 z s x 2 - s 1 - α d s , 1 Γ α - β 0 x 2 z s x 2 - s 1 - α - β d s d f x 1 , y 0 1 Γ α 0 x 1 z s x 1 - s 1 - α d s , 1 Γ α - β 0 x 1 z s x 1 - s 1 - α - β d s , f x 2 , y 0 1 Γ α 0 x 1 z s x 1 - s 1 - α d s , 1 Γ α - β 0 x 1 z s x 1 - s 1 - α - β d s + d f x 2 , y 0 1 Γ α 0 x 1 z s x 1 - s 1 - α d s , 1 Γ α - β 0 x 1 z s x 1 - s 1 - α - β d s , f x 2 , y 0 1 Γ α 0 x 2 z s x 2 - s 1 - α d s , 1 Γ α - β 0 x 2 z s x 2 - s 1 - α - β d s L 3 x 1 - x 2 + L 1 d y 0 1 Γ α 0 x 1 z s x 1 - s 1 - α d s , y 0 1 Γ α 0 x 2 z s x 2 - s 1 - α d s + L 2 d 1 Γ α - β 0 x 1 z s x 1 - s 1 - α - β d s , 1 Γ α - β 0 x 2 z s x 2 - s 1 - α - β d s = L 3 x 1 - x 2 + L 1 d 1 Γ α 0 x 1 z s x 1 - s 1 - α d s , 1 Γ α 0 x 2 z s x 2 - s 1 - α d s + L 2 d 1 Γ α - β 0 x 1 z s x 1 - s 1 - α - β d s , 1 Γ α - β 0 x 2 z s x 2 - s 1 - α - β d s .

Now we estimate the second term of the above expression.(73)d1Γα0x1zsx1-s1-αds,1Γα0x2zsx2-s1-αds=d1Γα0x1zsx1-s1-αds,1Γα0x1zsx2-s1-αds1Γαx1x2zsx2-s1-αdsd1Γα0x1zsx1-s1-αds,1Γα0x1zsx2-s1-αds+d1Γαx1x2zsx2-s1-αds,0^rΓα0x11x1-s1-α-1x2-s1-αds+rΓαx1x21x2-s1-αds2rΓ1+αx2-x1α+rΓ1+αx2α-x1α.

Estimating the third term of above equation similarly to above, we have(74)d1Γα-β0x1zsx1-s1-α-βds,1Γα-β0x2zsx2-s1-α-βds2rΓ1+α-βx2-x1α-β+rΓ1+α-βx2α-β-x1α-β.

So we obtain as follows: (75)dTzx1,Tzx2L3x1-x2+2rL1Γ1+αx2-x1α+rL1Γ1+αx2α-x1α+2rL2Γ1+α-βx2-x1α-β+rL2Γ1+α-βx2α-β-x1α-β.

From right side of the above expression, we can see that TuuU-r0^ is equicontinuous.

By Arzelà-Ascoli theorem, TuuU-r0^ is relatively compact.

By Theorem 20 and Schauder fixed point theorem, we can see the operator T has a fixed point in U-r(0^).

4. Continuous Dependence on Initial Condition of Solution

We consider the continuous dependence on initial condition of solution for proposed problem (24). Let y1(t,y1,0), y2(t,y2,0) are solutions with initial values y1,0,y2,0RF and z1(t,y1,0), z2(t,y2,0) are solutions of corresponding integral equation.

The following theorem shows the continuous dependence on initial condition of solution.

Theorem 21.

Suppose that the following conditions are satisfied:

L1,L2>0;x[0,1],y1,y2,z1,z2CF(I),(76)dfx,y1,z1,fx,y2,z2L1·dy1,y2+L2·dz1,z2.

qL11/Γ(α+1)+L2·1/Γ(α-β+1)<1.

Then following relation holds:(77)dy1t,y1,0,y2t,y2,01+L11-qΓα+1dy1,0,y2,0.

Proof.

Since z1(t,y1,0),z2(t,y2,0) are solutions of the integral equation (40), we have (78)z1x,y1,0fx,y1,0I0+αz1x,y1,0,I0+α-βz1x,y1,0,z2x,y2,0fx,y2,0I0+αz2x,y2,0,I0+α-βz2x,y2,0.

Now we estimate d(z1(x,y1,0),z2(x,y2,0)).(79)dz1x,y1,0,z2x,y2,0=dfx,y1,0I0+αz1x,y1,0,I0+α-βz1x,y1,0,fx,y2,0I0+αz2x,y2,0,I0+α-βz2x,y2,0L1dy1,0I0+αz1x,y1,0,y2,0I0+αz2x,y2,0+L2dI0+α-βz1x,y1,0,I0+α-βz2x,y2,0L1dy1,0,y2,0+dI0+αz1x,y1,0,I0+αz2x,y2,0+L2dI0+α-βz1x,y1,0,I0+α-βz2x,y2,0 By (iii) of Lemma 14, we have(80)L1dy1,0,y2,0+I0+αdz1x,y1,0,z2x,y2,0+L2I0+α-βdz1x,y1,0,z2x,y2,0L1dy1,0,y2,0+L1Γα+1dz1,z2+L2Γα-β+1dz1,z2=L1dy1,0,y2,0+qdz1,z2.

Consequently, we get(81)maxx0,1dz1x,y1,0,z2x,y2,0=dz1,z2L1dy1,0,y2,0+qdz1,z2,dz1,z2L11-qdy1,0,y2,0.

Next we estimate d(y1(t,y1,0),y2(t,y2,0)).(82)dy1t,y1,0,y2t,y2,0=dy1,01Γα0xz1s,y1,0x-s1-αds,y2,01Γα0xz2s,y1,0x-s1-αdsdy1,0,y2,0+d1Γα0xz1s,y1,0x-s1-αds,1Γα0xz2s,y1,0x-s1-αdsBy (iii) of Lemma 14, we have(83)dy1,0,y2,0+1Γα0xdz1s,y1,0,z2s,y1,0x-s1-αdsdy1,0,y2,0+1Γα+1dz1,z2.

By (81), we obtain(84)dy1t,y1,0,y2t,y2,01+L11-qΓα+1dy1,0,y2,0.

5. Conclusions

We obtained the uniqueness results by employing contraction mapping principle and the existence results by using Schauder fixed point theorem for the solution of the fuzzy multiterm fractional differential equation involving Caputo-type fuzzy fractional derivative of order 0<α<1 on complete metric space of continuous fuzzy number value functions. We also established the continuous dependence of solution on its initial condition. Our results can be expanded to the case in which the right side of equation involves more than two derivative terms.

Data Availability

No data were used to support this study.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.