Open access peer-reviewed chapter

# Existence Theory of Differential Equations of Arbitrary Order

Written By

Kamal Shah and Yongjin Li

Submitted: 01 November 2017 Reviewed: 16 February 2018 Published: 23 May 2018

DOI: 10.5772/intechopen.75523

From the Edited Volume

## Differential Equations - Theory and Current Research

Edited by Terry E. Moschandreou

Chapter metrics overview

View Full Metrics

## Abstract

The aims of this chapter are devoted to investigate a system of fractional-order differential equations (FDEs) with multipoint boundary conditions. Necessary and sufficient conditions are investigated for at most one solution to the proposed problem. Also, results for the existence of at least one or two positive solutions are developed by using a fixed-point theorem of concave-type operator for the considered problem. Further, we extend the conditions for more than two solutions and established some adequate conditions for multiplicity results to the proposed problem. Also, a result devoted to Hyers-Ulam stability is discussed. Suitable examples are provided to verify the established results.

### Keywords

• fractional differential equations
• coupled system
• boundary condition
• concave operator
• Mathematics subject classification: 26A33
• 34A08
• 35B40

## 1. Introduction

Arbitrary-order differential equations are the excellent tools in the description of many phenomena and process in different fields of science, technology, and engineering (see [1, 2]). Therefore, considerable attention has been paid to the subject of differential equations of arbitrary order (see [3, 4, 5] and the references therein). The area devoted to the existence of positive solutions to fractional differential equations and their system especially coupled systems was greatly studied by many authors (for details see [6, 7, 8, 9]). In all these articles, the concerned results were obtained by using classical fixed point theorems like Banach contraction principle, Leray-Schauder fixed point theorem, and fixed point theorems of cone type. The aforesaid area has been very well explored for both ordinary- and arbitrary-order differential equations. Existence and uniqueness results for nonlinear and linear, classical, as well as arbitrary-order differential equations have been investigated in many papers (see few of them as [10, 11, 12, 13]).

Another warm area of research in the theory of fractional-order differential equations (FDEs) is devoted to the multiplicity of solutions. Plenty of research articles are available on this topic in literature. In [14], the author studied the given boundary value problem (BVP) for existence of multiple solutions:

Dθ1pt+Htpt=0,tI,θ112,ptt=0=ptt=1=0.

where D is the Riemann-Liouville derivative of non-integer order and I=01. In same line, Kaufmann and Mboumi [15] studied the given boundary value problem of fractional differential equations for multiplicity of positive solutions:

Dθ1pt+ϕtHtpt=0,tI,θ112,ptt=0=p'tt=1=0,

where D is the Riemann-Liouville derivative and ϕCIR,HCI×RR.

In the last few decades, the theory devoted to the multiplicity of solutions is very well extended to coupled systems of nonlinear FDEs, and we refer to few papers in [16, 17, 18]. Wang et al. [19] established some conditions under which the given system of three point BVP

Dθ1pt=H1tqt;tI,Dθ2qt=H2tpt;tI,ptt=0=0ptt=1=μptt=ξ,qtt=0=0qtt=1=νqtt=ξ,

has a solution, where θ1,θ212 and μ,νI,ξ01,Hi:01×RR for i=1,2 are nonlinear functions.

In the last few decades, another important aspect devoted to stability analysis of FDEs with initial/boundary conditions has been given much attention. This is because stability is very important from the numerical and optimization point of view. Various forms of stabilities were studied for the aforesaid FDEs including exponential, Mittag-Leffler, and Lyapunov stability. Recently, Hyers-Ulam stability has given more attention. This concept was initially introduced by Ulam and then by Hyers (for details see [20, 21, 22]). Now, many articles have been written on this concept (see [23, 24, 25, 26, 27]). So far, the aforementioned stability has not yet well studied for multipoint BVPs of FDEs. Motivated by the aforesaid discussion, we propose the following coupled system of four-point BVP provided as

Dθ1pt=H1tptqt;tI;θ1m1m,Dθ2qt=H2tptqt;tI;θ2m1m,pjtt=0=qjtt=0=0,ptt=1=ptt=η,qtt=1=qtt=ξ.E1

where j=0,1,2,m2,m3,I=01,η,ξ01, H1,H2:01×0R+×0R+0R+ are continuous functions, and Dθ1,Dθ2 stand for Riemann-Liouville fractional derivative of order θ1,θ2 in sequel. We obtain necessary and sufficient conditions for the existence of solution to system (1) by using another type of fixed point result based on a concave-type operator with increasing or decreasing property. The idea then extends to form some conditions which ensure multiplicity of solutions to the considered problem. Also, we discuss some results about the Hyers-Ulam stability for the considered problem. Further by providing examples, we illustrate the established results.

## 2. Preliminaries

In the current section, we review few fundamental lemmas and results found in [2, 4, 6, 28, 29].

Arbitrary-order integral of function ψ:0R is recalled as

Iθ1ψt=1Γθ10ttsθ11ψsds,
where θ1>0 is a real number and also the integral is pointwise defined on R+

Arbitrary-order derivative in Riemann-Liouville sense for a function ψ0R is given by

Dθ1ψt=ddtm0ttsmθ11Γmθ1ψsds,θ1>0,where m=θ1+1.

[16] Let θ1>0, then for arbitrary CjR,j=1,2,,m,m=θ1+1, and the solution of

Dθ1ψt=ft
is provided by
ψt=Iθ1ft+C1tθ11+C2tθ12++Cmtθ1m.

[17, 28] Consider a Banach space E with a closed set CE. Then, C is said to be partially ordered if pq such that qpC. Further, C is said to be a cone if it holds the given conditions:

1. pC and for a real constant κ0 the relation κpC holds.

2. p and pC yield that 0C, where 0 is zero element of Banach space E

[17, 28] A closed and convex set C of E is said to be a normal cone if it obeys the given properties:

1. For 0pqE, there exists β>0, such that pEβqE;

2. pq, for all p,qE yields that there exist constants a,b>0 such that apqbq.

As is an equivalence relation, therefore defines a set Cf=pE:pf for fC. Obviously, one can derive that CfC for f0.

The operator S:CC is said to be λ concave for every θ,λ01,pC, if and only if SλpθλSp.

The operator S:CC is said to be to be increasing if p,qC,pq gives that SpSq.

[17, 28] Assume that S:CC is increasing λconcave operator for a normal cone C produced by Banach space E, such that there exists p0 with SfCf. Then, S has a unique fixed point pCf

[30] Let E be a Banach space with CB, which is closed and convex. Let E be a relatively open subset of C with 0E and S:E¯C be a continuous and compact operator. Then.

1. The operator S has a fixed point in E¯,

2. There exist w∂ℰ and λ01 with w=λSw.

[30] For a Banach space E together with a cone C, there exist two relatively open subsets A1 and A2 of E such that 0A1A¯1A2. Moreover, for a completely continuous operator S:CA¯2\A1C, one of the given conditions holds:

1. Spp for all pCA1;Spp, for all pCA2;

2. Spp for all pCA1;Spp, for all pCA2

Then, S has at least one fixed point in CA¯2\A1.

## 3. Main results

Let φC01R,η01 and λ1=1ηθ11<1, and then the unique solution to BVP of linear FDE

Dθ1pt=φt,tI,θ1m1m,pjtt=0=0,ptt=1=ptt=η,j=0,1,2,m2,m3,E2
is given by
pt=01Gtsφsds,E3
where Gts is the Green’s function defined by
Gts=1Γθ11λ1t1sθ11+tηsθ11+tsθ11,0stη1,1λ1t1sθ11+tηsθ11,0tsη1,1λ1t1sθ11+tsθ11,0ηst1,1λ1t1sθ11,0ηts1.E4

Proof. In view of Lemma 2.3, we may write Eq. (2) as

pt=Iθ1φt+C1tθ11+C2tθ12++Cmtθ1m.E5

In view of conditions pjtt=0=0,j=0,1,m2,m3,, Eq. (5) suffers from singularity; therefore, we have C2=C3==Cn=0. Hence, Eq. (5) becomes

pt=Iθ1φt+C1tθ11.E6

Applying boundary condition ptt=1=ptt=η and d=1η1θ in Eq. (6), one has

pt=Iθ1φt+tθ11λ1Iθ1φηIθ1φ1pt=01Gtsφsds.E7

where Gts is Green’s function given in Eq. (4).

In view of Theorem 3.1 and using λ1=1ηθ11,λ2=1ξθ21, the corresponding coupled system of integral equations to the proposed system (1) is given as

pt=01G1tsH1spsqsds,qt=01G2tsH2spsqsds,E8

where G1ts,G2ts are Green’s functions, which can be similarly computed like in Theorem 3.1. Further, they are continuous on I×I and satisfy the following properties:

1. maxtIG1tsλ1+11sθ11λ1=G11s,forallsI,maxtIG1tsλ2+11sθ21λ2=G21s,forallsI;

2. mintθ1θG1tsγ1s2G1s for every θs01;

mintθ1θG2tsγ2s2G1s for every θs01;

Further, taking that γ=infγ1=θθ11γ2=θθ21.

Let us define a Banach space by E=ptpCI endowed with a norm pE=maxtIpt. Further, in the norm for the product space, we define it as pqE×E=pE+qE. Clearly, E×EE×E is a Banach space. Onward, we define the cone CE×E by

C=pqE×E:mintIpt+qtγ(pq)E×E.

Consider an operator S:E×EE×E defined by

Spqt=01G1tsH1(spsqs)ds01G2(ts)H2(spsqs)ds.=S1ptS2qt.E9

It is to be noted that the fixed points of the operator S correspond with the solution of the system (1) under consideration.

Under the continuity of H1,H2:I×R+0×R+0R+0, the operator S satisfies that SCCandS:CC is completely continuous.

Proof. To derive SCC, let pqC, and then we have

S1ptqt=01G1tsH1spsqsdsγ101G11sH1sps,qsds.E10

Also, we get

S1ptqt=01G1tsH1spsqsds01G11sH1sps,qsds.E11

Thus, from Eqs. (10) and (11), we have

S1ptqtγS1pqE,for everytI.

Similarly, we can obtain

S2ptqtγS2pqE,for everytI.
ThusS1ptqt+S2ptqtγpqE×E,foralltI,
mintIS1ptqt+S2ptqtγpqE×E.

Hence, we have SpqCSCC.

Let us consider

maxtIH1tptqtM1,maxtIH2tptqtM2.

Then, we consider t1<t2I, such that

S1pqt2S1pqt1=01G(t2s)G1t1sH1(spsqs)dsM1Γθ1t2θ11t1θ11λ10ηηsθ11ds011sθ11ds+M1Γθ10t2t2sθ11ds0t1t1sθ11dsM1λ1Γθ1+1t2θ11t1θ11ηθ1λ1+λ1t2θ1t1θ1.E12

By the same fashion, we obtain for S2 as

S2pqt2S2pqt1M2λ2Γθ2+1t2θ21t1θ21ξθ2λ2+λ2t2θ2t1θ2.E13

The right hand sides of Eqs. (12) and (13) are approaching to zero at t1t2. Thus, the operator S is equi-continuous. Therefore, thanks to the Arzelá-Ascoli theorem, we receive that S=S1S2:CC is completely continuous.

Due to continuity of H1 and H2 on I×R+0×R+0R+, there exist φj,ψj,σjj=12:01R+0 for t01,p,q0 such that

H1tptqtφ1t+ψ1tpt+σ1tqt;H2tptqtφ2t+ψ2tpt+σ2tqt,
along with the following conditions:
1. Δ1=01G11sφ1sds<,Λ1=01G11sψ1s+σ1sds<1;

2. Δ2=01G21sφ2sds<,Λ2=01G21sψ2s+σ2sds<1

are satisfied. Then, the system (1) has at least one solution pq which lies in
E=pqC:(pq)E×E<min2Δ112Λ12Δ212Λ2.

Proof. Let E=pqC:pqE×E<r with min2Δ112Λ12Δ212Λ2<r.

Define the operator S:E¯C as in Eq. (9).

Let pqE that is pqE×E<r. Then, we have

S1pqt=maxtI01G1tsH1(spsqs)ds01G11sφ1sds+01G1(1s)ψ1spsds+01G1(1s)σ1sqsds01G11sφ1sds+r01G11sψ1s+σ1sds=Δ1+rΛ1r2.E14

Thus, from Eq. (14), we have

S1pqEr2.E15

Similarly, one can derive that

S2pqEr2.E16

Thus, from Eqs. (15) and (16), we get

SpqE×Er.E17

Therefore, SpqE¯. Hence, by Theorem 3.2 the operator S:E¯E is completely continuous.

Consider the eigenvalue problem:

pq=ρSpq,withρ01.E18

Under the assumption that pq is a solution of Eq. (18) for ρ01, we have

ptρmaxtI01G1tsH1(spsqs)dsρ01G11sφ1sds+01G1(1s)(ψ1sps+σ1sqs)dsρΔ1+rΛ1which implies thatpE<r2.

Similarly, we can obtain that qE<r2, so pqE×E<r, which implies that pq does not belong to E for all ρ01. Therefore, due to Theorem 2.10, S has a fixed point in E¯

Assume that the given hypothesis holds:

(H1) The nonlinear functions H1 and H2 are continuous on I×R+0×R+0R+0

(H2) For all tI, we have

H1tpq=0,H2tpq0,atpq=00

and

H1tpq1,H2t,1,11;atpq=11;

(H3) For all tI such that

0pp1,0qq1H1tpqH1tp1q1,H2tpqH1tp1q1;

(H4) For p,q0, there exist real numbers 0<λ,μ<1, such that for each tI,τ01, we have

H1tτpτqτλH1tpq,H2tτpτqτμH2tpq.

Under the assumptions H1H4, the BVP (1) has a unique solution in Cf where ft=tθ11tθ21.

Proof. Let maxλμ=κ and pqC. For each tI, using H4, we have

S1τpτqt=01G1tsH1sτpsτqsdsτλ01G1tsH1spsqsds=τλS1pqtτκS1pqt,

Analogously, we also get

S2τpτqtτκS2pqt.

In view of partial order on E×E induced by the cone C, we get SτpτqτκSptqt,τ01,pqC. Which yields that S is τ concave and nondecreasing operator with respect to the partial order by using hypothesis H4. Hence, taking fC for each tI defined by

ft=tθ11tθ21=f1tf2t.

Suppose that

w1=max1Γθ101H1s,1,1ds1Γθ201H2s,1,1ds

and

w2=max1Γθ101LsH1s,0,0ds1Γθ201KsH2s,0,0ds.

Also, from Green’s functions, we can obtain that

Ls=1sθ111+λ1λ1,Ks=1sθ211+λ2λ2.E19

Due to nondecreasing property of H1,H2 in view of H3, we get μ>0,ν>0. Therefore, applying (19) together with H4, one has

S1ht=01G1tsH1sf1sf2sds=01G1tsH1ssθ11sθ21ds01G1tsH1s,1,1ds1Γθ1011sθ11H1s,1,1dstθ11μf1t.

Similarly, we can get

S2ftμf2t.

Then, we obtain

Sfμf.E20

Like the aforesaid process, applying Eq. (19) together with H4, for each tI, one has

S1ft=01G1tsH1ssθ11sθ21ds01G1tsH1s,0,0ds1Γθ101LsH1s,0,0dstθ11νh1t,

With same fashion, we can obtain

S2ftνf2t.

Thus, we have

Sftνf.E21

From Eqs. (20) and (21), we produce

νfSfμf,

which implies that SfCf. So, thanks to Lemma 2.9, we see that the operator S is concave; hence, it has at most one fixed point pqCf which is the corresponding solution of BVPs (1).

Now, we define the following:

(C1) Hjj=12:I×R+0×R+0R+0 is uniformly bounded and continuous on I with respect to t.

(C2) Green’s functions G11s,G21s satisfy

0<01G11sds<,0<01G21sds<;

(C3) Let these limits hold:

H1ϱ=limp+qϱmaxtIH1tpqp+q,H2ϱ=limp+gϱmaxtIH2tpqp+q,H1,ϱ=limp+qϱinftIH1tpqp+q,H2,ϱ=limp+qϱinftIH2tpqp+q,whereϱ0.δ1=maxtI01G1tsds,δ2=maxtI01G2tsds.

Assume that the conditions C1C3 together with given assumptions are satisfied:

(H5) H1,0γ12θ1θG11sds>1,H1,γ12θ1θG11sds>1 and

H2,0γ22θ1θG21sds>1,H2,γ22θ1θG21sds>1.

Moreover, H1,0=H2,0=H1,=H2,= also hold:

(H6) There exists constant α>0 such that

maxtI,pqCαH1tpq<α2δ1

and

maxtI,pqCαH2tpq<α2δ2.

Then, the system (1) of BVPs has at least two positive solutions pq,p¯q¯ which obeying

0<pqE×E<α<p¯q¯E×E.E22

Proof. Assume that H5 holds, and consider ϵ,α,λ such that 0<ϵ<α<λ. Further we define a set by

Ωr=uvE×E:uvE×E<r,whererϵαλ.

Now, if

H1,0γ12θ1θG11sds>1andH2,0γ22θ1θG21sds>1.

Then, obviously, we can obtain that

SpqE×EpqE×E,forpqCΩε.E23

Now, if H1,γ12θ1θG11sds>1andH2,γ22θ1θG21sds>1.

Then, like the proof of Eq. (23), we have

SpqE×EpqE×E,forpqCΩλ.E24

Also, from H5 and pqCΩα, we get

S1pqt=01G1tsH(susvs)ds01G11sH1spsqsds.

From which we have

S1pqE×E<α2ϱ101G11sds=α2.

Similarly, we have S1pqE×E<α2 as pqCΩα. Hence, we have

SpqE×E<pqE×E,forpqCΩα.E25

Now, applying Lemma 2.11 to Eqs. (23) and (25) yields that S has a fixed point pqCΩ¯α\Cε and a fixed point in p¯q¯CΩ¯λ\Ωα. Hence, we conclude that the system of BVPs (1) has at least two positive solutions pq,p¯q¯ such that pqE×Eα and p¯q¯E×Eα. Thus, relation (22) holds.

Consider that C1C3 together with the following hypothesis are satisfied:

(H7) δ1H1,0<1, δ1H1,<1;δ2H1,0<1, and δ2H2,<1;

(H8) There exist ρ>0 such that

maxtI,pqCαγ12H1tpq>α2θ1θG11sds1,maxtI,pqCαγ22H2tpq>α2θ1θG21sds1

such that

0<pqE×E<α<p¯q¯E×E.

Then, the proposed coupled system of BVPs (1) has at least two positive solutions.

Proof. Proof is like the proof of Theorem 3.4.

Analogously, we deduce from Theorem 3.5 and 3.6 the following results for multiplicity of solutions to the system (1) of BVPs.

Under the conditions C1C3, there exist 2k positive numbers aj,âj,j=1,2k with a1<γ1â1<â1<a2<γ1â2<â2ak<γ1âk<âk and a1<γ2â1<â1<a2<γ2â2<â2ak<γ2âk<âk such that.

(H9) H1tptqtγ101G11sdsaj, for tpqI×γ1ajaj×γ2ajaj, and

H1tptqtδ1âi,fortpqI×γ1âjâj×γ2ajai,j=1,2k;

(H10) H2tptqtγ201G21sdsaj, for tpqI×γ1ajaj×γ2ajaj, and

H1tptqtδ2âj,fortpqI×γ1ajaj×γ2âjâj,j=1,2k.

Then, system (1) of BVPs has at least k solutions pjqj, satisfying

ajpjqjE×Eâj,j=1,2k.

Further, if assumptions C1C3 hold such that there exist 2k positive numbers bj,b̂j,j=1,2k, with

b1<b̂1<b2<b̂2<bk<b̂k,

together with following hypothesis hold:

(H11) H1tpq and H2tpq are nondecreasing on 0b̂k for all tI;

(H11)H1tptqtγ1θ1θG11sdsbj,H1tptqtδ1b̂j,j=1,2k,E26
H2tutvtγ2θ1θG21sdsbj,H2tptqtδ2b̂j,j=1,2k.

Then, system (1) of BVPs has at least k solutions pjqj, satisfying

bjpjqjE×Eb̂j,j=1,2k.

## 4. Hyers-Ulam stability

([31, Definition 2]) Consider a Banach space E×E such that S1,S2:E×EE×E be the two operators. Then, the operator system provided by

pt=S1pqt,qt=S2pqtE27
is called Hyers-Ulam stability if we can find Cii=1,2,3,4>0, such that for each ρii=12>0 and for each solution pqE×E of the inequalities given by
pH1pqE×Eρ1,qH2pqE×Eρ2,E28
there exist a solution p¯q¯E×E of system (26) which satisfy
pp¯E×EC1ρ1+C2ρ2,qq¯E×EC3ρ1+C4ρ2.E29

If λi, for i=1,2,,n be the (real or complex) eigenvalues of a matrix MCn×n, then the spectral radius ρM is defined by

ρM=maxλifori=12n.

Further, the matrix will converge to zero if ρM<1..

([31, Theorem 4]) Consider a Banach space E×E with S1,S2:E×EE×E be the two operators such that

S1pqS1pqE×EC1ppE×E+C2qqE×E,S2pqS2pqE×EC3ppE×E+C4qqE×E,forallpq,pqE×E,E30
and if the matrix
M=C1C2C3C4
converges to zero ([31, Theorem 1]), then the fixed points corresponding to operatorial system (26) are Hyers-Ulam stable.

For the stability results, the following should be hold:

(H13) Under the continuity of Hi,i=1,2, there exist ai,bic01,i=1,2 and pq,p¯q¯ such that

HitpqHitp¯q¯aitpp¯+bitqq¯,i=1,2.

In this section, we study Hyers-Ulam stability for the solutions of our proposed system. Thanks to Definition 4.1 and Theorem 4.3, the respective results are received.

Suppose that the assumptions H13 along with condition that matrix

M=01G11sa1sds01G11sb1sds01G21sa2sds01G21sb2sds
is converging to zero. Then, the solutions of (1) are Hyers-Ulam stable.

Proof. In view of Theorem 4.3, we have

S1pqS1p¯q¯E×E01G11sa1spp¯E×Eds+01G11sb1sqq¯E×EdsS2pqS2p¯q¯E×E01G21sa2spp¯E×Eds+01G21sb2sqq¯E×Eds.

From which we get

S1pqS1p¯q¯E×E01G11sa1sds]pp¯E×E+01G11sb1sdsqq¯E×ES2pqS2p¯q¯E×E01G21sa2sdspp¯E×E+01G21sb2sdsqq¯E×E.E31

Hence, we get

SpqSp¯q¯E×EMpqp¯q¯E×E,E32

where M=01G11sa1sds01G11sb1sds01G21sa2sds01G21sb2sds. Hence, we received the required results.

## 5. Illustrative examples

Consider the given system of BVPs

D72pt+1t2+ptqt13=0,D113qt+1+t+ptqt14=0,t01,pt=p't=p''t=qt=q't=q''t=0,att=0,p1=p14,q1=q13.E33

Clearly, H1tpq0,H2tpq0, at pq=00, and H1tpq0,H2tpq0, at pq=11. Simple computation yields that H1,H2 are nondecreasing for every t01. Also, for τ,t01, and p,q0, one has max1413=13,

H1tτpτqτ13H1tpq,H2tτpτqτ13H2tpq.

Thus, all the conditions of Theorem 3.4 are fulfilled, so the system (32) of BVPs has unique positive solution in Bf where ft=t52t92.

Consider the following system of BVPs:

D92pt+1+t2+pt+qt3=0,D92qt+1+t+pt+qt2=0,t01,pjt=qjt=0,j=0,1,2,3,att=0,p1=p12,q1=q12.E34

It is obvious that H1tpq0,H2tpq0, at pq=00, and H1tpq0,H2tpq0, at pq=11. Also, an easy computation yields that H1,H2 are nondecreasing for each t01. Moreover, for τ,t01, and p,q0, we see that max32=3,

H1tτpτqτ3H1tpq,H2tτpτqτ3H2tpq.

Thus, all the assumption of Theorem 3.4 is fulfilled, so the coupled system (33) has a unique positive solution in Bf where ft=t34t43.

Consider the following system of BVPs:

D72pt=t40+t20cospt+t220sinqt,t01,D72qt=t250+t260sinpt+t60cosqt,t01,pjt=qjt=0,j=0,1,2,att=0,p1=p12,q1=q12.E35

From system (33), we see that

H1tpqt40+t20cospt+t220sinqt
and
H2tpqt250+t260sinpt+t60cosqt.
where φ1t=t40, φ2t=t250, ψ1t=t20, ψ2t=t260, σ1t=t220, σ2t=t60. Also, η=ξ=12,λ1=λ2=0.17677. Thus, by computation, we have
Gj1s=6.657101s52Γ72,forj=1,2.

Upon computation, we get

Δ1=01G11sφ1sds=0.003577<,Δ2=01G21sφ2sds=0.000924<.

Similarly, we can also compute.

Λ1=01G11sψ1s+σ1sds=0.03092853<1,Λ2=01G21sψ2s+σ2sds=0.00289<1.Further, we see that max0.0076260.00185=0.007626. So, all the conditions of Theorem 3.3 are satisfied. So, the BVP (34) has at least one solution and the solution lies in

E=pqC:pqE×E<0.007626.

Taking the following system of BVPs

D112pt+pt+qt2+115+t2δ1=0,D112qt+pt+qt2+t15+t2δ2=0,t01,pjt=qjt=0,j=0,1,2,3,4,att=0,p1=p14,q1=q14.E36

It is simple to check that H1,0=H2,0=H1,=H2,=. Also, for any tpqI×I×I, we see that

H1tpq13δ1H2tpq13δ2.

Thus, all the assumptions of Theorem 3.5 are satisfied with taking α=1, so the coupled system (35) has two solutions satisfying 0<pqE×E<1<pqE×E.

Consider the following coupled systems of boundary value problems:

D52pt+Γ52tpt16+t2qt32=0,t01,D52qt+Γ529t2cospt16π+9tcosqt32π=0,t01,pjt=qjt=0,j=0,1,2,att=0,p1=p12,q1=q12.E37

Here, a1t=Γ52t16,b1t=Γ52t232,a2t=Γ529t216π,b2t=Γ529t32π. Moreover

M=01G11sa1sds01G11sb1sds01G21sa2sds01G21sb2sds=0.04600.00070.00680.0058.

Here, ρM=4.61×102<1. Therefore, matrix M converges to zero, and hence the solutions of (36) are Hyers-Ulam stable by using Theorem 4.4.

## 6. Conclusion

We have developed a comprehensive theory on existence of solutions and its Hyers-Ulam stability for system of multipoint BVP of FDEs. The concerned theory has been enriched by providing suitable examples.

## Acknowledgments

We are very thankful to the reviewers for his/her careful reading and suggestion which improved this chapter very well.

## Conflict of interest

We declare the there is no conflict of interest regarding this chapter.

## Research funder

This work has been supported by the National Natural Science Foundation of China (11571378).

## References

1. 1. Hilfer R. Applications of Fractional Calculus in Physics. Singapore: World Scientific; 2000
2. 2. Kilbas AA, Srivastava H, Trujillo J. Theory and application of fractional differential equations, North Holland Mathematics Studies. Vol. 204. Amsterdam: Elseveir; 2006
3. 3. Miller KS, Ross B. An Introduction to the Fractional Calculus and Fractional Differential Equations. New York: Wiley; 1993
4. 4. Podlubny I. Fractional differential equations. New York: Academic press; 1993
5. 5. Podlubny I. Fractional differential equations. In: Mathematics in Science and Engineering. New York: Academic Press; 1999
6. 6. Jalili M, Samet B. Existence of positive solutions to a coupled system of fractional differential equation. Mathematical Methods in Applied Science. 2015;38:1014-1031
7. 7. Khan RA, Shah K. Existence and uniqueness of solutions to fractional order multi-point boundary value problems. Communications on Pure and Applied Analysis. 2015;19:515-526
8. 8. Shah K, Khan RA. Study of solution to a toppled system of fractional differential equations with integral boundary. International Journal of Applied and Computational Mathematics. 2016;2(3) 19 pages
9. 9. Shah K, Khalil H, Khan RA. Investigation of positive solution to a coupled system of impulsive boundary value problems for nonlinear fractional order differential equations. Journal of Chaos, Solitons and Fractals. 2015;77:240-246
10. 10. Ahmad B, Alsaedi A. Existence and uniqueness of solutions for coupled systems of higher-order nonlinear fractional differential equations. Fixed Point Theory Applications. 2010 Article ID 364560
11. 11. Su X. Boundary value problem for a coupled system of nonlinear fractional differential equations. Applied Mathematics Letters. 2009;22:64-69
12. 12. Ahmad B, Nieto JJ. Existence results for a coupled system of nonlinear fractional differential equations with three-point boundary conditions. Computers & Mathematcs with Applications. 2009;58:18381843
13. 13. Yang W. Positive solution to nonzero boundary values problem for a coupled system of nonlinear fractional differential equations. Computers & Mathematcs with Applications. 2012;63:288-297
14. 14. Bai Z, Lü H. Positive solutions for boundary value problem of nonlinear fractional differential equation. Journal of Mathematical Analysis and Applications. 2005;311:495-505
15. 15. Kaufmann ER, Mboumi E. Positive solutions for boundary value problem of a nonlinear fractional differential equations, Electronic Journal of Qualitative Theory of Differential Equations. 2008;8:1-11
16. 16. Shah K, Khan RA. Existence and uniqueness of positive solutions to a coupled system of nonlinear fractional order differential equations with anti periodic boundary conditions. Differential Equations and Applications. 2015;7(2):245-262
17. 17. Shah K, Khalil H, Khan RA. Upper and lower solutions to a coupled system of nonlinear fractional differential equations. Progress Fractional Differential Applications. 2016;2(1):1-10
18. 18. Zhang SQ. Positive solutions for boundary value problem problems of nonlinear fractional differential equations. Electronic Journal of Differential Equations. 2006;2006:1-12
19. 19. Wang J, Xiang H, Lu Z. Positive solutions to nonzero boundary value problem for a coupled system of nonlinear fractional differential equations, International Journal of Differential Equations. 2010. 12-pages
20. 20. Ulam SM. Problems in Modern Mathematics. New York, USA: John Wiley and Sons; 1940
21. 21. Ulam SM. A Collection of Mathematical Problems. New York: Interscience; 1960
22. 22. Obloza M. Hyers stability of the linear differential equation. Rocznik Nauk-Dydakt. Prace Mat. 1993;13:259-270
23. 23. Wang J, Li X. Ulam-Hyers stability of fractional Langevin equations. Applied Mathematics and Computation. 2015;258:72-83
24. 24. Stamova I. Mittag-Leffler stability of impulsive differential equations of fractional order. Quarterly of Applied Mathematics. 2015;73(3):525-535
25. 25. Kumama P, Ali A, Shah K, Khan RA. Existence results and Hyers-Ulam stability to a class of nonlinear arbitrary order differential equations. Journal of Nonlinear Science and Applications. 2017;10:2986-2997
26. 26. Wang J, Lv L, Zhou Y. Ulam stability and data dependence for fractional differential equations with Caputo derivative. Electronic Journal of Qualitative Theory of Differential Equations. 2011;(63):1-10
27. 27. Haq F, Shah K, Rahman G, Shahzad M. Hyers–Ulam stability to a class of fractional differential equations with boundary conditions. International Journal of Applied Computational Mathematics. 2017:1-13
28. 28. Deimling K. Nonlinear Functional Analysis. New York: Springer-Verlag; 1985
29. 29. Zeidler E. Non Linear Functional Analysis and Its Applications. New York, USA: Springer; 1986
30. 30. Agarwal R, Meehan M, Regan DO. Fixed Points Theory and Applications. Cambridge: Cambridge University Press; 2004
31. 31. Urs C. Coupled fixed point theorem and applications to periodic boundary value problem. Miskolic Mathematical Notes. 2013;14:323-333

Written By

Kamal Shah and Yongjin Li

Submitted: 01 November 2017 Reviewed: 16 February 2018 Published: 23 May 2018