Open access peer-reviewed chapter

# Existence Theory of Differential Equations of Arbitrary Order

By Kamal Shah and Yongjin Li

Submitted: November 1st 2017Reviewed: February 16th 2018Published: May 23rd 2018

DOI: 10.5772/intechopen.75523

## 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 Dis 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 Dis 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,θ212and μ,νI,ξ01,Hi:01×RRfor i=1,2are 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θ2stand for Riemann-Liouville fractional derivative of order θ1,θ2in 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ψ:0Ris recalled as

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

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

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

[16] Letθ1>0,then for arbitraryCjR,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 spaceEwith a closed setCE. Then,Cis said to be partially ordered ifpqsuch thatqpC.Further,Cis said to be a cone if it holds the given conditions:

1. pCand for a real constantκ0the relationκpCholds.

2. pandpCyield that0C,where0is zero element of Banach spaceE

[17, 28] A closed and convex setCofEis said to be a normal cone if it obeys the given properties:

1. For0pqE,there existsβ>0,such thatpEβqE;

2. pq,for allp,qEyields that there exist constantsa,b>0such thatapqbq.

Asis an equivalence relation, therefore defines a setCf=pE:pfforfC. Obviously, one can derive thatCfCforf0.

The operatorS:CCis said to beλconcave for everyθ,λ01,pC,if and only ifSλpθλSp.

The operatorS:CCis said to be to be increasing ifp,qC,pqgives thatSpSq.

[17, 28] Assume thatS:CCis increasingλconcave operator for a normal coneCproduced by Banach spaceE,such that there existsp0withSfCf. Then,Shas a unique fixed pointpCf

[30] LetEbe a Banach space withCB, which is closed and convex. LetEbe a relatively open subset ofCwith0EandS:E¯Cbe a continuous and compact operator. Then.

1. The operatorShas a fixed point inE¯,

2. There existw∂ℰandλ01withw=λSw.

[30] For a Banach spaceEtogether with a coneC, there exist two relatively open subsetsA1andA2ofEsuch that0A1A¯1A2. Moreover, for a completely continuous operatorS:CA¯2\A1C,one of the given conditions holds:

1. Sppfor allpCA1;Spp, for allpCA2;

2. Sppfor allpCA1;Spp, for allpCA2

Then,Shas at least one fixed point inCA¯2\A1.

## 3. Main results

LetφC01R,η01andλ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
whereGtsis 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 Gtsis 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,G2tsare Green’s functions, which can be similarly computed like in Theorem 3.1. Further, they are continuous on I×Iand satisfy the following properties:

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

2. mintθ1θG1tsγ1s2G1sfor every θs01;

mintθ1θG2tsγ2s2G1sfor every θs01;

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

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

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

Consider an operator S:E×EE×Edefined by

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

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

Under the continuity ofH1,H2:I×R+0×R+0R+0, the operatorSsatisfies thatSCCandS:CCis 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 S2as

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 Sis equi-continuous. Therefore, thanks to the Arzelá-Ascoli theorem, we receive that S=S1S2:CCis completely continuous.

Due to continuity ofH1andH2onI×R+0×R+0R+, there existφj,ψj,σjj=12:01R+0fort01,p,q0such 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 solutionpqwhich lies in
E=pqC:(pq)E×E<min2Δ112Λ12Δ212Λ2.

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

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

Let pqEthat 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¯Eis completely continuous.

Consider the eigenvalue problem:

pq=ρSpq,withρ01.E18

Under the assumption that pqis 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 pqdoes not belong to Efor all ρ01.Therefore, due to Theorem 2.10, Shas a fixed point in E¯

Assume that the given hypothesis holds:

(H1) The nonlinear functions H1and H2are 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 tIsuch 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 assumptionsH1H4,the BVP(1) has a unique solution inCfwhereft=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×Einduced by the cone C, we get SτpτqτκSptqt,τ01,pqC.Which yields that Sis τconcave and nondecreasing operator with respect to the partial order by using hypothesis H4. Hence, taking fCfor each tIdefined 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,H2in 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 Sis concave; hence, it has at most one fixed point pqCfwhich is the corresponding solution of BVPs (1).

Now, we define the following:

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

(C2) Green’s functions G11s,G21ssatisfy

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 conditionsC1C3together with given assumptions are satisfied:

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

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

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

(H6) There exists constantα>0such that

maxtI,pqCαH1tpq<α2δ1

and

maxtI,pqCαH2tpq<α2δ2.

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

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

Proof.Assume that H5holds, 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 H5and pqCΩα,we get

S1pqt=01G1tsH(susvs)ds01G11sH1spsqsds.

From which we have

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

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

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

Now, applying Lemma 2.11 to Eqs. (23) and (25) yields that Shas 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 thatC1C3together with the following hypothesis are satisfied:

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

(H8) There existρ>0such 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 conditionsC1C3, there exist2kpositive numbersaj,âj,j=1,2kwitha1<γ1â1<â1<a2<γ1â2<â2ak<γ1âk<âkanda1<γ2â1<â1<a2<γ2â2<â2ak<γ2âk<âksuch that.

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

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

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

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

Then, system(1) of BVPs has at leastksolutionspjqj,satisfying

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

Further, if assumptionsC1C3hold such that there exist2kpositive numbersbj,b̂j,j=1,2k, with

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

together with following hypothesis hold:

(H11) H1tpqandH2tpqare nondecreasing on0b̂kfor alltI;

(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 leastksolutionspjqj,satisfying

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

## 4. Hyers-Ulam stability

([31, Definition 2]) Consider a Banach spaceE×Esuch thatS1,S2:E×EE×Ebe the two operators. Then, the operator system provided by

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

Ifλi,fori=1,2,,nbe the (real or complex) eigenvalues of a matrixMCn×n, then the spectral radiusρMis defined by

ρM=maxλifori=12n.

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

([31, Theorem 4]) Consider a Banach spaceE×EwithS1,S2:E×EE×Ebe 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,2and 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 assumptionsH13along 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,atpq=00,andH1tpq0,H2tpq0,atpq=11. Simple computation yields thatH1,H2are nondecreasing for everyt01.Also, forτ,t01,andp,q0, one hasmax1413=13,

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

Thus, all the conditions ofTheorem 3.4 are fulfilled, so the system(32) of BVPs has unique positive solution inBfwhereft=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 thatH1tpq0,H2tpq0,atpq=00,andH1tpq0,H2tpq0,atpq=11. Also, an easy computation yields thatH1,H2are nondecreasing for eacht01.Moreover, forτ,t01,andp,q0, we see thatmax32=3,

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

Thus, all the assumption ofTheorem 3.4 is fulfilled, so the coupled system(33) has a unique positive solution inBfwhereft=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 thatmax0.0076260.00185=0.007626. So, all the conditions ofTheorem 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 thatH1,0=H2,0=H1,=H2,=.Also, for anytpqI×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 satisfying0<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, matrixMconverges to zero, and hence the solutions of(36) are Hyers-Ulam stable by usingTheorem 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).

chapter PDF
Citations in RIS format
Citations in bibtex format

## More

© 2018 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## How to cite and reference

### Cite this chapter Copy to clipboard

Kamal Shah and Yongjin Li (May 23rd 2018). Existence Theory of Differential Equations of Arbitrary Order, Differential Equations - Theory and Current Research, Terry E. Moschandreou, IntechOpen, DOI: 10.5772/intechopen.75523. Available from:

### chapter statistics

1Crossref citations

### Related Content

Next chapter

#### An Extension of Massera’s Theorem for N-Dimensional Stochastic Differential Equations

By Boudref Mohamed Ahmed, Berboucha Ahmed and Osmanov Hamid Ibrahim Ouglu

First chapter

#### Platelets

By Gökhan Cüce and Tahsin Murad Aktan

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.