Next Article in Journal
Solving Non-Linear Fractional Variational Problems Using Jacobi Polynomials
Next Article in Special Issue
Existence and Stability Results for a Fractional Order Differential Equation with Non-Conjugate Riemann-Stieltjes Integro-Multipoint Boundary Conditions
Previous Article in Journal
A Novel Hybrid Algorithm for Minimum Total Dominating Set Problem
Previous Article in Special Issue
Nonlocal Fractional Evolution Inclusions of Order α ∈ (1,2)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Ulam Stability and Multiplicity Results to a Nonlinear Coupled System with Integral Boundary Conditions

1
Department of Mathematics, University of Malakand, Chakdara Dir(L), Khyber Pakhtunkhwa 18800, Pakistan
2
KMUTT Fixed Point Research Laboratory, KMUTT-Fixed Point Theory and Applications Research Group, Theoretical and Computational Science Center (TaCS), Science Laboratory Building, Faculty of Science, King Mongkut’s University of Technology Thonburi (KMUTT), 126 Pracha-Uthit Road, Bang Mod, Thrung Khru, Bangkok 10140, Thailand
*
Author to whom correspondence should be addressed.
Mathematics 2019, 7(3), 223; https://doi.org/10.3390/math7030223
Submission received: 1 February 2019 / Revised: 16 February 2019 / Accepted: 17 February 2019 / Published: 27 February 2019

Abstract

:
This manuscript is devoted to establishing existence theory of solutions to a nonlinear coupled system of fractional order differential equations (FODEs) under integral boundary conditions (IBCs). For uniqueness and existence we use the Perov-type fixed point theorem. Further, to investigate multiplicity results of the concerned problem, we utilize Krasnoselskii’s fixed-point theorems of cone type and its various forms. Stability analysis is an important aspect of existence theory as well as required during numerical simulations and optimization of FODEs. Therefore by using techniques of functional analysis, we establish conditions for Hyers-Ulam (HU) stability results for the solution of the proposed problem. The whole analysis is justified by providing suitable examples to illustrate our established results.

1. Introduction

Fractional order differential equations (FODEs) emerge in the scientific demonstration of numerous frameworks and different fields of science such as physics, chemistry, economics, polymer rheology, aerodynamics, electrodynamics of complicated medium, blood flow phenomena, biophysics, etc. (see [1,2,3,4,5]). Recently, many authors have studied FODEs from different aspects, one is the numerical and scientific techniques for finding solutions and the other is the theoretical perspective of uniqueness and existence of solutions. The interest of the researchers in the investigation of FODEs lies in the incontrovertible fact that fractional-order models (FOM) are found to be highly realistic and practical, compared to the integer order models. Because there are additional degrees of opportunity in the FOM, in consequence, the subject of FODEs is gaining more attention from researchers. Another facet of research, which has been completely studied for integer order differential equations is devoted to uniqueness and existence of solutions to boundary value problems (BVPs). The mentioned aspect has been very well studied for FODEs, we refer the readers [6,7,8,9,10]. Uniqueness and existence results of solutions to multi-point BVPs have been studied via classical fixed point theorems such as the Schauder fixed point theorem and the Banach contraction principle, see [11,12,13,14,15,16,17].
FODEs under integral boundary conditions (IBCs) have been investigated very well because these type of equations are increasingly used in fluid-mechanics and dynamical problems. Jankowski [18] studied the ordinary differential equation under IBCs given by
y ( ϑ ) = F ( ϑ , y ( ϑ ) ) , ϑ [ 0 , T ] , T > 0 , y ( ϑ ) | ϑ = 0 = δ 0 T y ( s ) d s + d 0 , d 0 R ,
where F C ( [ 0 , T ] × R , R ) and δ = 1 or 1 . He developed a sufficient condition for iterative approximate solutions to the above problem.
Nanware and Dhaigude [19] have investigated the aforementioned BVP under the IBCs for FODE as given by
D + 0 σ y ( ϑ ) = F ( ϑ , y ( ϑ ) ) , ϑ [ 0 , T ] , T > 0 , y ( ϑ ) | ϑ = 0 = δ 0 T y ( s ) d s + d 0 , d 0 R ,
where 0 < σ 1 , δ is 1 or 1 and F C ( [ 0 , T ] × R , R ) , D + 0 σ is Riemann-Liouville fractional derivative of order σ is defined in (2). The aforementioned author also studied the iterative approximate solution to the above FODEs.
In the same line Cabada and Wang [20] studied the following problem under IBCs as
C D + 0 σ y ( ϑ ) + φ ( ϑ , y ( ϑ ) ) = 0 ; ϑ ( 0 , 1 ) , y ( 0 ) = y ( 0 ) = 0 , y ( 1 ) = δ 0 1 y ( s ) d s ,
where σ ( 2 , 3 ] , δ ( 0 , 2 ) and y : [ 0 , 1 ] × [ 0 , ] [ 0 , ] are the continuous functions. Also we remark that C D + 0 σ stands for Caputo’s fractional derivative.
Inspired from the aforementioned work, in this article we investigate a system of nonlinear FODEs with IBCs as
D + 0 σ y ( ϑ ) + φ ( ϑ , y ( ϑ ) , z ( ϑ ) ) = 0 ; ϑ ( 0 , 1 ) ; m 1 < σ m , D + 0 æ z ( ϑ ) + χ ( ϑ , y ( ϑ ) , z ( ϑ ) ) = 0 ; ϑ ( 0 , 1 ) ; m 1 < æ m , y ( 0 ) = y ( 0 ) = y ( 0 ) = = y ( m 2 ) ( 0 ) = 0 , y ( 1 ) = δ 0 1 y ( s ) d s , z ( 0 ) = z ( 0 ) = z ( 0 ) = 0 = z ( m 2 ) ( 0 ) = 0 , z ( 1 ) = ϱ 0 1 z ( s ) d s ,
such that m 3 , δ , ϱ ( 0 , 2 ) , the functions φ , χ : [ 0 , 1 ] × [ 0 , ] × [ 0 , ] [ 0 , ] are continuous functions and D + 0 σ , D + 0 æ stand for Riemann-Liouville fractional derivatives is defined in (2). We claim that such a system of FODEs are very rarely considered for stability as well as multiplicity results. Our analysis is devoted to the existence theory of a solution, multiplicity results and stability analysis of the suggested problem.
During the last few decades another part of research, which has been considered for FODEs and got much attention from the researchers is stability analysis. Numerous forms of stabilities have been studied in literature which are Mittag-Leffer stability, exponential stability, Lyapunov stability etc., we refer [21,22,23].
The Ulam stability was first presented by Ulam in 1940 and then brilliantly explained by Hyers in 1941. For more information about HU stability, we refer [24,25]. The HU stability results were generalized and extended by many researchers for FODEs under IBCs. In 1978, Jung studied the said stability for ODEs. Oblaz, Benchohra, etc., have studied the said stability for FODEs but their investigation was limited to initial value problems, we refer to [26,27,28]. To the best of our information and knowledge, the HU stability has been very rarely studied for coupled system of FODEs under IBCs. Therefore in this article we investigate HU stability to the considered problem. Here we remark that we also provide some necessary results for nonexistence of solution. Finally a series of examples are provided to support our analysis.

2. Axillary Results

In the current section, we review some fundamental definitions and useful results of functional analysis, fractional calculus and fixed point theory (see reference [1,2,8,29,30,31,32]). Here, first of all, we define the Banach space which is utilized throughout in this article.
Let us define E = { y ( ϑ ) | y C [ 0 , 1 ] } with the norm | | y | | = max ϑ [ 0 , 1 ] | y ( ϑ ) | . We define the norm for the product space as | | ( y , z ) | | = | | y | | + | | z | | . Obviously ( E × E , | | · | | ) is a Banach space.
Let K = [ θ , 1 θ ] for each θ ( 0 , 1 ) , then, we define the cone C E × E by
C = { ( y , z ) E × E : min ϑ K [ y ( ϑ ) + z ( ϑ ) ] λ | | ( y , z ) | | } . C r = { ( y , z ) C : ( y , z ) r } , C r = { ( y , z ) K : ( y , z ) = r } .
As in [31], we define positive solution as follows.
Definition 1.
A pair of functions ( y , z ) E × E is called a positive solution of problem (1) under the given IBCs if D + 0 σ y , D + 0 ρ z L 1 [ 0 , 1 ] with ( y , z ) > ( 0 , 0 ) on ( 0 , 1 ] × ( 0 , 1 ] , where the functions y , z satisfy the IBCs given in (1) respectively, for all ϑ [ 0 , 1 ] .
Definition 2.
The Riemann-Liouville fractional derivative of order σ > 0 of a continuous function y : ( 0 , ) R is defined as
D + 0 σ y ( ϑ ) = 1 Γ ( m σ ) d d ϑ m 0 ϑ ( ϑ s ) m σ 1 y ( s ) d s ,
where m = [ σ ] + 1 and [ σ ] denotes the integer part of σ.
Definition 3.
The Riemann-Liouville fractions of integration of order σ > 0 of a continuous function y : ( 0 , ) R is defined by
I + 0 σ y ( ϑ ) = 1 Γ ( σ ) 0 ϑ ( ϑ s ) σ 1 y ( s ) d s ,
where the integral is point-wise defined on ( 0 , ) .
Lemma 1.
Let σ > o , then the FODE
D + 0 σ y ( ϑ ) = 0
has a solution given by
y ( ϑ ) = i = 1 m y i ( 0 ) i ! ϑ i .
Lemma 2.
Let σ > 0 . Then we have
I + 0 σ [ D + 0 σ y ( ϑ ) ] = y ( ϑ ) i = 0 m y i ( 0 ) i ! ϑ i .
Lemma 3.
[2] Let σ > o and ϑ C ( 0 , 1 ) L ( 0 , 1 ) , then the FODE
D + 0 σ y ( ϑ ) = h ( ϑ )
has a solution given by
y ( ϑ ) = c 1 ϑ σ 1 + c 2 ϑ σ 2 + + c m ϑ σ m + I + 0 σ h ( ϑ ) ,
where c i R for i = 0 , 1 , 2 , , m and m = [ σ ] + 1 .
Definition 4.
[32,33] On the Banach space E defined afore, the mapping d : E × E R n is called a generalized metric on E if x , y , and y , z E with y x , z y , z y , then the following hold
( A 1 )
d ( y , z ) = 0 y = z , y , z E
( A 2 )
d ( y , z ) = d ( z , y ) , y , z E
( A 3 )
d ( x , y ) = d ( x , z ) + d ( z , y ) + d ( y , y ) , x , y , y , z , E .
Further the pair ( E , d ) is called a generalized metric space.
Definition 5.
[32,33] Let M = { M m , m R + m × m } , for any matrix B M the spectral radius is defined by æ ( B ) = sup { | λ ^ i | , i = 1 , 2 , , m } , where λ ^ i , for i = 1 , 2 , , m are the eigenvalues of the matrix B and the matrix will converge to zero if æ ( B ) < 1 .
Lemma 4.
[32,33] A complete generalized metric space ( M , d ) , with operator B : M M such that there a matrix B M with
d ( B y , B z ) B d ( y , z ) , for all y , z M ,
if æ ( B ) < 1 , then B has a fixed point in M .
Lemma 5.
[32,33] Consider a Banach space E with cone C E and y C is relatively open set with 0 y and B : y ¯ y be a completely continuous mapping. Then one of the following hold
( A 1 )
The mapping B has a fixed point in y ¯
( A 2 )
There exist y y and η ( 0 , 1 ) with y = η B y .
Lemma 6.
[33,34] Consider a cone C in the Banach space E and if A 1 and A 2 be two bounded open sets in E , such that 0 A 1 A ¯ 1 A 2 . Let B : C ( A ¯ 2 A 1 ) C be completely continuous operator and one of the following satisfied:
( 1 )
B y y y C A 1 ; B y , y C A 2
( 2 )
B y y y C A 1 ; B y , y C A 2
Then B has at least one fixed point in C ( A ¯ 2 A 1 ) .

3. Existence of at Least One Solution

Lemma 7.
Let h C [ 0 , 1 ] , then the BVP
D + 0 σ y ( ϑ ) + h ( ϑ ) = 0 ; ϑ ( 0 , 1 ) ; m 1 < σ m , y ( 0 ) = y ( 0 ) = y ( 0 ) = = y ( m 2 ) ( 0 ) = 0 , y ( 1 ) = δ 0 1 y ( s ) d s , ,
where δ ( 0 , 2 ) , has the following unique solution
y ( ϑ ) = 0 1 H σ ( ϑ , s ) h ( s ) d s ,
where H σ is the Green’s function given by
H σ ( ϑ , s ) = ϑ σ 1 ( 1 s ) σ 1 ( σ δ + δ s ) ( σ δ ) ( ϑ s ) σ 1 ( σ δ ) Γ ( δ ) , 0 s ϑ 1 , ϑ σ 1 ( 1 s ) σ 1 ( σ δ + δ s ) ( σ δ ) Γ ( δ ) , 0 ϑ s 1 .
Proof. 
Thanks to Lemma 3 for (7), one has
y ( ϑ ) = I + 0 σ h ( ϑ ) + c 1 ϑ σ 1 + c 2 ϑ σ 2 + + c m ϑ σ m .
By using initial condition y ( 0 ) = y ( 0 ) = y ( 0 ) = = y ( m 2 ) ( 0 ) = 0 , we get c 2 = c 3 = = c m = 0 . Therefore (9) implies that
y ( ϑ ) = c 1 ϑ σ 1 I + 0 σ h ( ϑ ) .
By using boundary condition y ( 1 ) = δ 0 1 y ( s ) d s in (10), we get
c 1 = 0 1 ( ϑ s ) σ 1 Γ ( σ ) h ( s ) d s + δ 0 1 y ( s ) d s .
Hence we have the following solution to (1)
y ( ϑ ) = 0 ϑ ( ϑ s ) σ 1 Γ ( σ ) h ( s ) d s + ϑ σ 1 0 1 ( 1 s ) σ 1 Γ ( σ ) h ( s ) d s + δ ϑ σ 1 0 1 y ( s ) d s .
Let B = 0 1 y ( s ) d s , then from Equation (11), we have
B = 0 1 0 ϑ ( ϑ s ) σ 1 Γ ( σ ) h ( s ) d s + 0 1 0 1 ϑ σ 1 ( ϑ s ) σ 1 Γ ( σ ) h ( s ) d s + 0 1 δ B ϑ σ 1 d s
B = 0 1 ( 1 s ) σ σ Γ ( σ ) h ( s ) d s + 0 1 ( 1 s ) σ 1 σ Γ ( σ ) h ( s ) d s + 1 σ δ B
implies Equation (12), so we get
B = 1 σ δ 0 1 ( 1 s ) σ Γ ( σ ) h ( s ) d s + 1 σ δ 0 1 ( 1 s ) σ 1 Γ ( σ ) h ( s ) d s .
Replacing this valve in (11), we get
y ( ϑ ) = 0 t ( ϑ s ) σ 1 Γ ( σ ) h ( s ) d s + ϑ σ 1 0 1 ( 1 s ) σ 1 Γ ( σ ) h ( s ) d s δ σ δ 0 1 ϑ σ 1 ( 1 s ) σ σ Γ ( σ ) h ( s ) d s + δ σ δ 0 1 ϑ σ 1 ( 1 s ) σ 1 Γ ( σ ) h ( s ) d s . = 0 ϑ ( ϑ s ) σ 1 Γ ( σ ) h ( s ) d s + 0 1 ϑ σ 1 ( 1 s ) σ 1 ( σ δ + δ s ) ( σ δ ) Γ ( σ ) h ( s ) d s = 0 ϑ ϑ σ 1 ( 1 s ) σ 1 ( σ δ + δ s ) ( σ δ ) ( ϑ s ) σ 1 ( σ δ ) Γ ( σ ) h ( s ) d s + ϑ 1 ϑ σ 1 ( 1 s ) σ 1 ( σ δ + δ s ) ( σ δ ) Γ ( σ ) h ( s ) d s = 0 1 H σ ( ϑ , s ) h ( s ) d s ,
where H σ ( ϑ , s ) is the Green’s function of BVP (7). Similarly we can obtain z ( ϑ ) = 0 1 H æ ( ϑ , s ) h ( s ) d s , where H æ ( ϑ , s ) is the Green’s function for the second equation of the system (1) and is given by
H æ ( ϑ , s ) = ϑ æ 1 ( 1 s ) æ 1 ( æ ϱ + ϱ s ) ( æ ϱ ) ( ϑ s ) æ 1 ( æ ϱ ) Γ ( æ ) , 0 s ϑ 1 , ϑ æ 1 ( 1 s ) æ 1 ( æ ϱ + ϱ s ) ( æ ϱ ) Γ ( æ ) , 0 ϑ s 1 .
 □
Lemma 8.
Let H ( ϑ , s ) = ( H σ ( ϑ , s ) , H æ ( ϑ , s ) ) be the Green’s function of (1) defined in Equations (8) and (13). This H ( ϑ , s ) has the given properties
(F1) 
H ( ϑ , s ) is continuous function on the unit square ( ϑ , s ) [ 0 , 1 ] × [ 0 , 1 ]
(F2) 
H ( ϑ , s ) 0 ϑ , s [ 0 , 1 ] and H ( ϑ , s ) > 0 ϑ , s ( 0 , 1 )
(F3) 
max 0 ϑ 1 H ( ϑ , s ) = H ( 1 , s ) , s [ 0 , 1 ]
(F4) 
min ϑ [ θ , 1 θ ] H ( ϑ , s ) λ ( s ) H ( 1 , s ) for each θ , s ( 0 , 1 ) ,
where λ = min { λ σ = θ σ 1 , λ æ = θ æ 1 } .
Now according to Lemma 7, we can write system (1) as follows
y ( ϑ ) = 0 1 H σ ( ϑ , s ) φ ( s , y ( s ) , z ( s ) ) d s , z ( ϑ ) = 0 1 H æ ( ϑ , s ) χ ( s , y ( s ) , z ( s ) ) d s .
Let B : E × E E × E be the operator defined as
B ( y , z ) ( ϑ ) = 0 1 H σ ( ϑ , s ) φ ( s , y ( s ) , z ( s ) ) d s , 0 1 H æ ( ϑ , s ) χ ( s , y ( s ) , z ( s ) ) d s = B 1 ( y , z ) ( ϑ ) , B 2 ( y , z ) ( ϑ ) . .
Then the fixed point of operator B coincides with the solution of the coupled system (1).
Theorem 1.
Consider that u , v : [ 0 , 1 ] × [ 0 , ) × [ 0 , ) [ 0 , ) are continuous. Then B ( C ) C and B : C C is completely continuous, where B is defined in (15).
Proof. 
To prove that B ( C ) C , let ( y , z ) C , then by Lemma 8, we have B ( y , z ) C and from ( F 4 ) and ϑ K , we obtain
B 1 ( y ( ϑ ) , z ( ϑ ) ) = 0 1 H σ ( ϑ , s ) φ ( s , y ( s ) , z ( s ) ) d s λ σ 0 1 H σ ( 1 , s ) φ ( s , y ( s ) ) , z ( s ) d s .
Also from ( F 3 ) , we obtain
B 1 ( y ( ϑ ) , z ( ϑ ) ) = 0 1 H σ ( ϑ , s ) φ ( s , y ( s ) , z ( s ) ) d s 0 1 H σ ( 1 , s ) φ ( s , y ( s ) ) , z ( s ) d s .
Thus from (16) and (17), we have
B 1 ( y ( ϑ ) , z ( ϑ ) ) λ B 1 ( y , z ) , for all ϑ K .
Similarly, one can write that
B 2 ( y ( ϑ ) , z ( ϑ ) ) λ B 2 ( y , z ) , for all ϑ K .
Thus
B 1 ( y ( ϑ ) , z ( ϑ ) ) + B 2 ( y ( ϑ ) , z ( ϑ ) ) λ B ( y , z ) , for all ϑ K ,
min ϑ K [ B 1 ( y ( ϑ ) , z ( ϑ ) ) + B 2 ( y ( ϑ ) , z ( ϑ ) ) ] λ B ( y , z ) .
Hence we have B ( y , z ) C B ( C ) C . Next, like the proof of Theorem 1 of [35], and applying the Arzelà-Ascoli’s theorem, it can be easily proven that B : C C is completely continuous □
Theorem 2.
Consider that φ and χ are continuous on [ 0 , 1 ] × [ 0 , ) × ( 0 , ) [ 0 , ) , and there exist f i ( ϑ ) , H i ( ϑ ) , ( i = 1 , 2 ) : ( 0 , 1 ) [ 0 , ) that satisfy
( A 1 )
| φ ( ϑ , y , z ) φ ( ϑ , y ¯ , z ¯ ) | u 1 ( ϑ ) | y y ¯ | + v 1 ( ϑ ) | z z ¯ | , for ϑ ( 0 , 1 ) and y , z , y ¯ , z ¯ 0
( A 2 )
| χ ( ϑ , y , z ) χ ( ϑ , y ¯ , z ¯ ) | u 2 ( ϑ ) | y y ¯ | + v 2 ( ϑ ) | z z ¯ | , for ϑ ( 0 , 1 ) and y , z , y ¯ , z ¯ 0
( A 3 )
æ ( B ) < 1 , where B { M 2 , 2 R + 2 × 2 } is a matrix given by
0 1 H σ ( 1 , s ) u 1 ( s ) d s 0 1 H σ ( 1 , s ) v 1 ( s ) d s 0 1 H æ ( 1 , s ) u 2 ( s ) d s 0 1 H æ ( 1 , s ) v 2 ( s ) d s .
Then the system (1) has a unique positive solution ( y , z ) C .
Proof. 
Let us define a generalized metric d : E 2 × E 2 R 2 by
d ( ( y , z ) , ( y ¯ , z ¯ ) ) = | | y y ¯ | | | | z z ¯ | | , for all ( y , z ) , ( y ¯ , z ¯ ) E × E .
Obviously ( E × E , d ) is a generalized complete metric space. Then for any ( y , z ) , ( y ¯ , z ¯ ) E × E and using property ( F 3 ) we get
| B 1 ( y , z ) ( ϑ ) B 1 ( y ¯ , z ¯ ) ( ϑ ) | max ϑ [ 0 , 1 ] 0 1 | H σ ( ϑ , s ) | [ | φ ( s , y ( s ) , z ( s ) ) φ ( s , y ¯ ( s ) , z ¯ ( s ) ) | ] d s 0 1 H σ ( 1 , s ) [ u 1 ( s ) y y ¯ + v 1 ( s ) z z ¯ ] d s 0 1 u 1 ( s ) H σ ( 1 , s ) d s y y ¯ + 0 1 v 1 ( s ) H σ ( 1 , s ) d s z z ¯ .
Similarly we can show that
| B 2 ( y , z ) B 2 ( y ¯ , z ¯ ) | 0 1 u 2 ( s ) H æ ( 1 , s ) d s y y ¯ + 0 1 v 2 ( s ) H æ ( 1 , s ) d s z z ¯ .
Thus we have
| B ( y , z ) B ( y ¯ , z ¯ ) | B d ( y , z ) , ( y ¯ , z ¯ ) , ( y , z ) , ( y ¯ , z ¯ ) E × E ,
where
B = 0 1 H σ ( 1 , s ) u 1 ( s ) d s 0 1 H σ ( 1 , s ) v 1 ( s ) d s 0 1 H æ ( 1 , s ) u 2 ( s ) d s 0 1 H æ ( 1 , s ) v 2 ( s ) d s .
As æ ( B ) < 1 , in the light of Lemma 4, system (1) has a unique positive solution. □
Theorem 3.
Consider that φ and χ are continuous on [ 0 , 1 ] × [ 0 , ) × ( 0 , ) [ 0 , ) and there exist a i , b i , c i ( i = 1 , 2 ) : ( 0 , 1 ) [ 0 , ) satisfying:
( A 4 )
φ ( ϑ , y ( ϑ ) , z ( ϑ ) ) a 1 ( ϑ ) + b 1 ( ϑ ) y ( ϑ ) + c 1 ( ϑ ) z ( ϑ ) , ϑ ( 0 , 1 ) , y , z 0
( A 5 )
χ ( ϑ , y ( ϑ ) , z ( ϑ ) ) a 2 ( ϑ ) + b 2 ( ϑ ) y ( ϑ ) + c 2 ( ϑ ) z ( ϑ ) , ϑ ( 0 , 1 ) , y , z 0
( A 6 )
Λ 1 = 0 1 H σ ( 1 , s ) a 1 ( s ) d s < , Δ 1 = 0 1 H σ ( 1 , s ) [ b 1 ( s ) + c 1 ( s ) ] d s < 1 2
( A 7 )
Λ 2 = 0 1 H æ ( 1 , s ) a 2 ( s ) d s < , Δ 2 = 0 1 H æ ( s , s ) [ b 2 ( s ) + c 2 ( s ) ] d s < 1 2 .
Then the system (1) has at least one positive solution in
( y , z ) C : ( y , z ) r , where max Λ 1 1 2 Δ 1 , Λ 2 1 2 Δ 2 < r .
Proof. 
Define Ω = ( y , z ) C : ( y , z ) < r with max Λ 1 1 2 Δ 1 , Λ 2 1 2 Δ 2 < r .
According to the Theorem 1, the operator B : Ω ¯ C is completely continuous. Let ( y , z ) Ω , such that ( y , z ) < r . Then, we have
B 1 ( y , z ) = max ϑ [ 0 , 1 ] | 0 1 H σ ( ϑ , s ) φ ( s , y ( s ) , z ( s ) ) | d s 0 1 H σ ( 1 , s ) a 1 ( s ) d s + 0 1 H σ ( 1 , s ) b 1 ( s ) | y ( s ) | d s + 0 1 H σ ( 1 , s ) c 1 ( s ) | z ( s ) | d s 0 1 H σ ( 1 , s ) a 1 ( s ) d s + r 0 1 H σ ( 1 , s ) [ b 1 ( s ) + c 1 ( s ) ] d s Λ 1 + r Δ 1 < r 2 ,
Similarly, B 2 ( y , z ) < r 2 , thus B ( y , z ) < r . Therefore, thanks to Lemma 5, we have B ( y , z ) Ω ¯ , thus B : Ω ¯ Ω ¯ . Let there exist ς ( 0 , 1 ) and ( y , z ) Ω such that ( y , z ) = ς B ( y , z ) . Then in the light of assumptions ( A 4 ) , ( A 5 ) and by ( F 4 ) of Lemma 8, we get ϑ [ 0 , 1 ]
| z ( ϑ ) | ς 0 1 H σ ( ϑ , s ) | φ ( s , y ( s ) , z ( s ) ) | d s ς 0 1 H σ ( 1 , s ) a 1 ( s ) d s + 0 1 H σ ( 1 , s ) b 1 ( s ) y ( s ) d s + 0 1 H σ ( 1 , s ) c 1 ( s ) z ( s ) d s ζ Δ 1 + r Λ 1 < ς r 2
which implies that y < ς r 2 . Similarly, it can be proved that z < ς r 2 . From which, we have ( y , z ) < ς r , with ς ( 0 , 1 ) which is a contradiction that ( y , z ) Ω as r = ( y , z ) . Thus, according to Lemma 5, B has at least one fixed point ( y , z ) Ω ¯ . □
Next the following assumptions and notations will be used:
( C 1 )
φ , χ : [ 0 , 1 ] × [ 0 , ) × [ 0 , ) [ 0 , ) are continuous and φ ( ϑ , 0 , 0 ) = χ ( ϑ , 0 , 0 ) = 0 uniformly with respect to ϑ on [ 0 , 1 ]
( C 2 )
H σ ( 1 , s ) , H æ ( 1 , s ) defined in Lemma 8 satisfy
0 < 0 1 H σ ( 1 , s ) d s < , 0 < 0 1 H æ ( 1 , s ) d s <
( C 3 )
Let these limits hold
φ α = lim ( y , z ) ( α , α ) sup ϑ [ 0 , 1 ] φ ( ϑ , y , z ) y + z , χ α = lim ( y , z ) ( α , α ) sup ϑ [ 0 , 1 ] χ ( ϑ , y , z ) y + z , φ α = lim ( y , z ) ( α , α ) f ϑ [ 0 , 1 ] φ ( ϑ , y , z ) y + z , χ æ = lim ( y , z ) ( α , α ) f ϑ [ 0 , 1 ] χ ( ϑ , y , z ) y + z , where α { 0 , } .
α σ = max ϑ [ 0 , 1 ] 0 1 H σ ( ϑ , s ) d s , α æ = max ϑ [ 0 , 1 ] 0 1 H æ ( ϑ , s ) d s .
Theorem 4.
If the assumptions ( C 1 ) ( C 2 ) hold and one of the following conditions is also satisfied:
( D 1 )
φ 0 λ 2 θ 1 θ H σ ( 1 , s ) d s > 1 , φ α σ < 1 and χ 0 λ 2 θ 1 θ H æ ( 1 , s ) d s > 1 , χ α æ < 1 .
Moreover, φ 0 = χ 0 = and φ = χ = 0
( D 2 )
There exist two constants η 1 , η 2 with 0 < η 1 η 2 such that φ ( ϑ , · , · ) and χ ( ϑ , · , · ) are nondecreasing on [ 0 , η 2 ] ϑ [ 0 , 1 ] ,
φ ( ϑ , λ σ η 1 , λ æ η 1 ) η 1 2 λ σ θ 1 θ H σ ( 1 , s ) d s 1 , χ ( ϑ , λ σ η 1 , λ æ η 1 ) η 1 2 λ æ θ 1 θ H æ ( 1 , s ) d s 1 and φ ( ϑ , η 2 , η 2 ) η 2 2 α σ , χ ( ϑ , η 2 , η 2 ) η 2 2 α æ , for all ϑ [ 0 , 1 ] ,
where λ , H σ ( 1 , s ) , H æ ( 1 , s ) defined in Lemma 8 and φ 0 , χ 0 , φ , χ , α σ , σ α defined in Equations (18) and (19). Then the coupled system (1) has at least one positive solution.
Proof. 
B as defined in (15) is completely continuous.
Case I. Let the condition ( D 1 ) hold. Taking φ 0 λ σ 2 θ 1 θ H σ ( 1 , s ) d s > 1 , then there exists a constant κ 1 > 0 such that
φ ( ϑ , y , z ) ( φ 0 r 1 ) ( y ( ϑ ) + z ( ϑ ) ) , χ ( ϑ , y , z ) ( χ 0 r 2 ) ( y ( ϑ ) + z ( ϑ ) ) , for all ϑ [ 0 , 1 ] , y , z [ 0 , κ 1 ] ,
where r 1 > 0 , and satisfies the conditions
( φ 0 r 1 ) λ σ 2 2 θ 1 θ H σ ( 1 , s ) d s 1 , ( χ 0 r 1 ) λ æ 2 2 θ 1 θ H æ ( 1 , s ) d s 1 .
So for ϑ [ 0 , 1 ] , ( y , z ) C κ 1 , we have
B 1 ( y , z ) ( ϑ ) = 0 1 H σ ( ϑ , s ) φ ( s , y ( s ) , z ( s ) ) d s λ σ 0 1 H σ ( 1 , s ) φ ( s , y ( s ) , z ( s ) ) d s ( φ 0 r 1 ) λ σ 2 2 0 1 H σ ( 1 , s ) d s ( y , z ) ( y , z ) 2 .
Analogously
B 2 ( y , z ) ( ϑ ) = 0 1 H æ ( ϑ , s ) χ ( s , y ( s ) , z ( s ) ) d s λ æ 0 1 H æ ( 1 , s ) φ ( s , y ( s ) , z ( s ) ) d s ( χ 0 r 2 ) λ æ 2 2 0 1 H æ ( 1 , s ) d s ( y , z ) ( y , z ) 2 .
Therefore, we have
B ( y , z ) B 1 ( y , z ) + B 2 ( y , z ) ( y , z ) .
Also for φ α σ < 1 and χ α æ < 1 , there exists a constant say κ 2 ¯ > 0 such that φ ( ϑ , y , z ) ( φ + r 2 ) ( y + z ) , χ ( ϑ , y , z ) ( χ + r 2 ) ( y + z ) , for ϑ [ 0 , 1 ] , y , z ( κ 2 ¯ , ) , where r 2 > 0 satisfies the conditions α σ ( φ + r 2 ) 1 , α æ ( χ + r 2 ) 1 . Let J = max ϑ [ 0 , 1 ] , y , z [ 0 , κ 2 ¯ ] φ ( ϑ , y , z ) , L = max ϑ [ 0 , 1 ] , y , z [ 0 , κ 2 ¯ ] χ ( ϑ , y , z ) , then φ ( ϑ , y , z ) J + ( φ + r 2 ) ( y , z ) , χ ( ϑ , y , z ) L + ( χ + r 2 ) ( y , z ) . Now setting max { κ 1 , κ 2 ¯ , J α σ ( 1 α σ ( φ + r 2 ) ) 1 } κ 2 2 , max { κ 1 , κ 2 ¯ , L α æ ( 1 α æ ( χ + r 2 ) ) 1 } κ 2 2 .
So for any ϑ [ 0 , 1 ] , ( y , z ) C κ 2 , we obtain
B 1 ( y , z ) ( ϑ ) = 0 1 H σ ( ϑ , s ) φ ( s , y ( s ) , z ( s ) ) d s λ σ 0 1 H σ ( 1 , s ) φ ( s , y ( s ) , z ( s ) ) d s 0 1 H σ ( 1 , s ) ( J + ( φ + r 2 ) [ u ( s ) + z ( s ) ] d s J 0 1 H σ ( 1 , s ) d s + ( φ + r 2 ) 0 1 H σ ( 1 , s ) d s ( y , z ) < κ 2 2 α σ ( φ + r 2 ) κ 2 2 + ( φ + r 2 ) α σ ( y , z ) < κ 2 2 .
Similarly B 2 ( y , z ) ( ϑ ) < κ 2 2 , as ( y , z ) C κ 2 , thus we have
B ( y , z ) < ( y , z ) .
Case II. If assumptions in ( D 2 ) hold, then in light of the definition of C for ( y , z ) C η 1 , we have ( y , z ) = η 1 , for ϑ K . Then from ( D 2 ) , we have
B 1 ( y , z ) ( ϑ ) = 0 1 H σ ( ϑ , s ) φ ( s , y ( s ) , z ( s ) ) d s λ σ θ 1 θ H σ ( 1 , s ) φ ( s , y ( s ) , z ( s ) ) d s λ σ θ 1 θ H σ ( 1 , s ) d s η 1 2 λ σ θ 1 θ H σ ( 1 , s ) d s 1 = η 1 2 .
Similarly it can also be obtained that B 2 ( y , z ) ( ϑ ) η 1 2 , for ( y , z ) C η 1 , and we get
B ( y , z ) = B 1 ( y , z ) + B 2 ( y , z ) ( y , z ) .
Also for ( y , z ) C η 2 , we get that ( y , z ) = η 2 for ϑ [ 0 , 1 ] . Then from ( D 2 ) , one can get
B 1 ( y , z ) ( ϑ ) = 0 1 H σ ( ϑ , s ) φ ( s , y ( s ) , z ( s ) ) d s 0 1 H σ ( 1 , s ) φ ( s , y ( s ) , z ( s ) ) d s η 2 2 α σ 0 1 H σ ( 1 , s ) d s = η 2 2 .
Similarly, it can also obtained that B 2 ( y , z ) ( ϑ ) η 2 2 , ( y , z ) C η 2 . Hence, we have
B ( y , z ) = B 1 ( y , z ) + B 2 ( y , z ) ( y , z ) .
Now according to the application of Lemma 6 to (20) and (21) or (22) and (23) implies that B has a fixed point ( y 1 , z 1 ) C ¯ κ , η or ( y 1 , z 1 ) C ¯ κ i , η i ( i = 1 , 2 ) such that y 1 ( ϑ ) λ σ y 1 > 0 and z 1 ( ϑ ) λ æ z 1 > 0 , ϑ [ 0 , 1 ] . From which it follows that the coupled system (1) has at least one positive solution. □
Theorem 5.
Under the conditions ( C 1 ) ( C 3 ) and if the following assumptions hold
( D 3 )
If φ 0 α σ < 1 ; φ λ σ 2 θ 1 θ H σ ( 1 , s ) d s > 1 and χ 0 α æ < 1 ; χ λ æ 2 θ 1 θ H σ ( 1 , s ) d s > 1 ,
then the coupled system (1) has at least one positive solution. Further, if φ 0 = χ 0 = 0 and φ = χ = , where λ , H σ ( 1 , s ) , H æ ( 1 , s ) defined in Lemma 8 and φ 0 , χ 0 , φ , χ , α σ , σ α defined in Equations (18) and (19), then the the considered system (1) has at least one positive solution.
Proof. 
Proof can be obtained as proof of Theorem 4. □

4. Existence of More Than One Solutions

Theorem 6.
Consider that ( C 1 ) ( C 3 ) hold and the following conditions are satisfied:
( D 4 )
If φ 0 λ σ 2 θ 1 θ H σ ( 1 , s ) d s > 1 , φ λ σ 2 θ 1 θ H σ ( 1 , s ) d s > 1 and
χ 0 λ æ 2 θ 1 θ H æ ( 1 , s ) d s > 1 , χ λ æ 2 θ 1 θ H σ ( 1 , s ) d s > 1 .
Moreover, φ 0 = χ 0 = φ = χ = also hold:
( D 5 )
there exists a > 0 such that
max ϑ [ 0 , 1 ] , ( y , z ) C a φ ( ϑ , y , z ) < a 2 α σ and max ϑ [ 0 , 1 ] , ( y , z ) C a χ ( ϑ , y , z ) < a 2 α æ .
Then the coupled system (1) has at least two positive solutions ( y , z ) , ( y ¯ , z ¯ ) such that
0 < ( y , z ) < a < ( y ¯ , z ¯ ) .
Where λ , H σ ( 1 , s ) , H æ ( 1 , s ) are defined in Lemma 8 and φ 0 , χ 0 , φ , χ , α σ , σ α defined in Equations (18) and (19)
Proof. 
Let ( D 4 ) hold. Select κ , η such that 0 < κ < μ < η . Now if φ 0 λ σ 2 θ 1 θ H σ ( 1 , s ) d s > 1 and χ 0 λ æ 2 θ 1 θ H æ ( 1 , s ) d s > 1 , then like the proof of Theorem 4, we have
B ( y , z ) | ( y , z ) , for ( y , z ) C κ .
Now, if φ λ σ 2 θ 1 θ H σ ( 1 , s ) d s > 1 and χ λ æ 2 θ 1 θ H æ ( 1 , s ) d s > 1 , then like the proof of Theorem 4, we have
B ( y , z ) ( y , z ) , for ( y , z ) C η .
Also from ( D 5 ) , ( y , z ) C μ , we get
B 1 ( y , z ) ( ϑ ) = 0 1 H σ ( ϑ , s ) φ ( s , y ( s ) , z ( s ) ) d s 0 1 H σ ( 1 , s ) φ ( s , y ( s ) , z ( s ) ) d s < μ 2 α σ 0 1 H σ ( 1 , s ) d s = μ 2 .
Similarly, we have B 1 ( y , z ) ( ϑ ) < μ 2 as ( y , z ) C μ . Hence, we have
B ( y , z ) < | ( y , z ) , for ( y , z ) C μ .
Now according to Lemma 6 for (25) and (27), we have gives that B has a fixed point ( y , z ) C ¯ κ , μ and a fixed point in ( y ¯ , z ¯ ) C ¯ μ , η . Therefore system (1) has at least two positive solutions ( y , z ) , ( y ¯ , z ¯ ) such that ( y , z ) μ and ( y ¯ , z ¯ ) μ . Thus the relation (24) holds. □
Theorem 7.
Consider that ( C 1 ) ( C 3 ) hold together with the given conditions
( D 6 )
α σ φ 0 < 1 and φ α σ < 1 ; α æ χ 0 < 1 , and χ α æ < 1
( D 7 )
there exist μ > 0 such that
max ϑ K , ( y , z ) C μ φ ( ϑ , y , z ) > μ 2 λ σ 2 θ 1 θ H σ ( 1 , s ) d s 1 , max ϑ K , ( y , z ) C μ χ ( ϑ , y , z ) > μ 2 λ æ 2 θ 1 θ H æ ( 1 , s ) d s 1 ,
such that
0 < ( y , z ) < μ < ( y ¯ , z ¯ ) ,
where λ , H σ ( 1 , s ) , H æ ( 1 , s ) defined in Lemma 8 and φ 0 , χ 0 , φ , χ , α σ , σ α defined in Equations (18) and (19). Thus the system (1) has at least two positive solutions.
Proof. 
We left the proof out, as it similar to the proof of Theorem 6. □
In same line for multiple solutions we give the following results.
Theorem 8.
Let ( C 1 ) ( C 3 ) hold. If there exist 2 m positive numbers u L , u ^ L , L = 1 , 2 m with u 1 < λ σ u ^ 1 < u ^ 1 < u 2 < λ σ u ^ 2 < u ^ 2 u m < λ σ u ^ m < u ^ m and u 1 < λ æ u ^ 1 < u ^ 1 < u 2 < λ æ u ^ 2 < u ^ 2 u m < λ æ u ^ m < u ^ m , such that
( D 8 )
φ ( ϑ , y ( ϑ ) , z ( ϑ ) ) u L λ σ 0 1 H σ ( 1 , s ) d s 1 , for ( ϑ , y , z ) [ 0 , 1 ] × [ λ σ u L , u L ] × [ λ æ u L , u L ] , and
φ ( ϑ , y ( ϑ ) , z ( ϑ ) ) α σ 1 u ^ L , for ( ϑ , y , z ) [ 0 , 1 ] × [ λ σ u ^ L , u ^ L ] × [ λ æ u L , u L ] , L = 1 , 2 m ,
( D 9 )
χ ( ϑ , y ( ϑ ) , z ( ϑ ) ) u L λ æ 0 1 H æ ( 1 , s ) d s ) 1 , for ( ϑ , y , z ) [ 0 , 1 ] × [ λ æ u L , u L ] × [ λ σ u L , u L ] , and
χ ( ϑ , y ( ϑ ) , z ( ϑ ) ) α æ 1 u ^ L , for ( ϑ , y , z ) [ 0 , 1 ] × [ λ σ u L , u L ] × [ λ æ u ^ L , u ^ L ] , L = 1 , 2 m .
where λ , H σ ( 1 , s ) , H æ ( 1 , s ) defined in Lemma 8.
Then the coupled system (1) has at least m -positive solutions ( y L , z L ) , satisfying
u L ( y L , z L ) u ^ L , L = 1 , 2 m .
Theorem 9.
Suppose that ( C 1 ) ( C 3 ) holds. If there exist 2 m positive numbers u L , u ^ L , L = 1 , 2 m , with u 1 < u ^ 1 < u 2 < u ^ 2 < u m < u ^ m such that
( D 10 )
φ and χ are non-decreasing on [ 0 , u ^ m ] ϑ [ 0 , 1 ] ;
( D 11 )
φ ( ϑ , y ( ϑ ) , z ( ϑ ) ) u L λ σ θ 1 θ H σ ( 1 , s ) d s ) 1 , φ ( ϑ , y ( ϑ ) , z ( ϑ ) ) u ^ L α σ , L = 1 , 2 m ,
χ ( ϑ , y ( ϑ ) , z ( ϑ ) ) u L λ æ θ 1 θ H æ ( 1 , s ) d s ) 1 , χ ( ϑ , y ( ϑ ) , z ( ϑ ) ) u ^ L α æ , L = 1 , 2 m .
Hence we conclude that there exist at least m positive solutions ( y L , z L ) , corresponding to coupled system (1) which satisfy
u L ( y L , z L ) u ^ L , L = 1 , 2 m .

5. Hyers-Ulam Stability

Definition 6.
[30] Let B 1 , B 2 : E × E E × E be the two operators. Then the system of operator equations
y ( ϑ ) = B 1 ( y , z ) ( ϑ ) z ( ϑ ) = B 2 ( y , z ) ( ϑ )
is called the HU stability if we can find J i ( i = 1 , 2 , 3 , 4 ) > 0 , with æ i ( i = 1 , 2 ) > 0 and for each solution ( y * , z * ) E × E of the inequalities given by
| | y * φ ( y * , z * ) | | E × E æ 1 , | | z * χ ( y * , z * ) | | E × E æ 2 ,
there exists a solution ( y ¯ , z ¯ ) E × E of system (28) such that
| | y * y ¯ | | E × E k 1 æ 1 + k 2 æ 2 , | | z * z ¯ | | E × E k 3 æ 1 + k 4 æ 2 ,
Theorem 10.
[30] Let B 1 , B 2 : E × E E × E be the two operators such that
| | B 1 ( y , z ) B 1 ( y * , z * ) | | E × E k 1 | | y y * | | E × E d s + k 2 | | z z * | | E × E d s , | | B 2 ( y , z ) B 2 ( y * , z * ) | | E × E k 3 | | y y * | | E × E d s + k 4 | | z z * | | E × E d s , f o r a l l ( y , z ) , ( y * , z * ) E × E ,
and if the matrix
B = k 1 k 2 k 3 k 3
converges to zero, then the fixed points corresponding to operator system (28) are HU-stable. Further, the given condition holds ( M 11 ) under the continuity of φ i , i = 1 , 2 , there exist f i , H i C ( 0 , 1 ) , i = 1 , 2 and ( y , z ) , ( y ¯ , z ¯ ) such that
| φ i ( ϑ , y , z ) φ i ( ϑ , y ¯ , z ¯ ) | f i ( ϑ ) | y y ¯ | + H i ( ϑ ) | z z ¯ | , i = 1 , 2 .
In this section, we study HU stability for the solutions of our proposed system.
Theorem 11.
Suppose that the assumption ( M 11 ) along with condition that matrix
B = 0 1 H σ ( 1 , s ) u 1 ( s ) d s 0 1 H σ ( 1 , s ) v 1 ( s ) d s 0 1 H æ ( 1 , s ) u 2 ( s ) d s 0 1 H æ ( 1 , s ) v 2 ( s ) d s .
is converging to zero. Then, the solutions of (1) are HU-stable.
Proof. 
Thanks to Theorem 2, we have
| | B 1 ( y , z ) B 1 ( y * , z * ) | | E × E 0 1 H σ ( 1 , s ) u 1 ( s ) | | y y * | | E × E d s + 0 1 H σ ( 1 , s ) v 1 ( s ) | | z z * | | E × E d s , | | B 2 ( y , z ) B 2 ( y * , z * ) | | E × E 0 1 H æ ( 1 , s ) u 2 ( s ) | | y y * | | E × E d s + 0 1 H æ ( 1 , s ) v 2 ( s ) | | z z * | | E × E d s .
From which we get
| | B 1 ( y , z ) B 1 ( y * , z * ) | | E × E 0 1 H σ ( 1 , s ) u 1 ( s ) d s | | y y * | | E × E + 0 1 H σ ( 1 , s ) v 1 ( s ) d s | | z z * | | E × E , | | B 2 ( y , z ) B 2 ( y * , z * ) | | E × E 0 1 H æ ( 1 , s ) u 2 ( s ) d s | | y y * | | E × E + 0 1 H æ ( 1 , s ) v 2 ( s ) d s | | z z * | | E × E .
Analogously one has
| | P ( y , z ) P ( y * , z * ) | | E × E B | | ( y , z ) ( y * , z * ) | | E × E ,
such that
B = 0 1 H σ ( 1 , s ) u 1 ( s ) d s 0 1 H σ ( 1 , s ) v 1 ( s ) d s 0 1 H æ ( 1 , s ) u 2 ( s ) d s 0 1 H æ ( 1 , s ) v 2 ( s ) d s .
Hence, we get the required results. □

6. Example

To verify the aforesaid established analysis we provide some test problems here in the given sequel.
Example 1.
Take the system of given BVPs with IBCs as
D + 0 7 2 y ( ϑ ) + ϑ + 1 4 [ Γ 5 2 | y ( ϑ ) | + cos | z ( ϑ ) | ] = 0 , ϑ [ 0 , 1 ] , y , z 0 D + 0 7 2 z ( ϑ ) + ϑ 2 + 1 4 [ sin | y ( ϑ ) | + | z ( ϑ ) | ] = 0 , ϑ [ 0 , 1 ] , y , z 0 y ( 0 ) = y ( 0 ) = y ( 0 ) = 0 y ( 1 ) = 1 2 0 1 y ( s ) d s , z ( 0 ) = z ( 0 ) = z ( 0 ) = 0 z ( 1 ) = 1 3 0 1 z ( s ) d s .
Since φ ( ϑ , y ( ϑ ) , z ( ϑ ) ) = ϑ + 1 4 [ Γ ( 5 2 ) | y ( ϑ ) | + cos | z ( ϑ ) | ] , χ ( ϑ , y ( ϑ ) , z ( ϑ ) ) = ϑ 2 + 1 4 [ sin | y ( ϑ ) | + | z ( ϑ ) | ] .
Also as m = [ 3.5 ] + 1 = 4 , δ = 1 2 and ϱ = 1 3 .
Then
| φ ( ϑ , y 2 , z 2 ) φ ( ϑ , y 1 , z 1 ) | Γ ( 5 2 ) ϑ + 1 4 | y 2 y 1 | + ϑ + 1 4 | z 2 z 1 | , | χ ( ϑ , y 2 , z 2 ) χ ( ϑ , y 1 , z 1 ) | ϑ 2 + 1 4 | y 2 y 1 | + ϑ 2 + 1 4 | z 2 z 1 | .
where u 1 ( ϑ ) = ϑ + 1 4 Γ ( 5 2 ) , v 1 ( ϑ ) = ϑ + 1 4 , u 2 ( ϑ ) = v 2 ( ϑ ) = ϑ 2 + 1 4 , so one can get
B = 0 1 H σ ( 1 , s ) u 1 ( s ) d s 0 1 H σ ( 1 , s ) v 1 ( s ) d s 0 1 H æ ( 1 , s ) u 2 ( s ) d s 0 1 H æ ( 1 , s ) v 2 ( s ) d s = 8 11 8 165 π 496 11583 π 496 11583 π .
d e t ( B λ ^ I ) = 8 11 λ ^ 8 165 π 496 11583 π 496 11583 π λ ^ .
We get λ ^ 1 = 0.728 and λ ^ 2 = 0.024 since æ ( B ) = sup { | λ ^ i | , i = 1 , 2 } = 0.728 < 1 . Therefore due to Theorem 2, BVPs (34) has a unique positive solution given by
y ( ϑ ) = 0 1 H 7 2 ( ϑ , s ) s + 1 4 [ Γ ( 5 2 ) | y ( s ) | + cos | z ( s ) | ] d s , z ( ϑ ) = 0 1 H 7 2 ( ϑ , s ) s 2 + 1 4 [ sin | y ( s ) | + | z ( s ) | ] d s ,
where H 7 2 ( ϑ , s ) and H 7 2 ( ϑ , s ) are the Green’s functions given by
H 7 2 ( ϑ , s ) = ϑ 5 2 ( 1 s ) 5 2 ( 3 + 1 2 s ) 3 ( ϑ s ) 5 2 3 Γ ( 7 2 ) , 0 s ϑ 1 , ϑ 5 2 ( 1 s ) 5 2 ( 3 + 1 2 s ) 3 Γ ( 7 2 ) , 0 ϑ s 1 .
H 7 2 ( ϑ , s ) = ϑ 5 2 ( 1 s ) 5 2 ( 19 6 + 1 3 s ) 19 6 ( ϑ s ) 5 2 19 6 Γ ( 7 2 ) , 0 s ϑ 1 , ϑ 5 2 ( 1 s ) 5 2 ( 19 6 + 1 3 s ) 19 6 Γ ( 7 2 ) , 0 ϑ s 1 .
Further, by the use of Theorem 11, the solution is HU-stable.
Example 2.
Taking a system of FODEs with IBCs as
D + 0 10 3 y ( ϑ ) + a ( ϑ ) y ( ϑ ) + z ( ϑ ) = 0 , D + 0 7 2 z ( ϑ ) + b ( ϑ ) y ( ϑ ) + z ( ϑ ) 3 = 0 , ϑ ( 0 , 1 ) , y ( 0 ) = y ( 0 ) = y ( 0 ) = y ( 0 ) = 0 y ( 1 ) = 0 1 y ( s ) d s , z ( 0 ) = z ( 0 ) = z ( 0 ) = z ( 0 ) = 0 z ( 1 ) = 0 1 z ( s ) d s .
From the given system one has
φ ( ϑ , y , z ) = a ( ϑ ) y ( ϑ ) + z ( ϑ )
and
χ ( ϑ , y , z ) = b ( ϑ ) y ( ϑ ) + z ( ϑ ) 3 , m = 4 , δ = ϱ = 1 .
Also a , b : [ 0 , 1 ] [ 0 , ) are continuous. Now φ 0 = lim ( y , z ) 0 φ ( ϑ , y , z ) y + z = , similarly χ 0 = .
Obviously we compute φ = 0 = χ . Hence due to Theorem 4, system (36) has at least one positive solution.
Example 3.
Taking another test problem with IBCs as
D + 0 9 2 y ( ϑ ) + ( 1 ϑ 2 ) [ y ( ϑ ) + z ( ϑ ) ] 2 = 0 , D + 0 14 3 z ( ϑ ) + [ y ( ϑ ) + z ( ϑ ) ] 3 = 0 , ϑ ( 0 , 1 ) , y ( 0 ) = y ( 0 ) = y ( 0 ) = y ( 0 ) = y ( 0 ) = 0 , y ( 1 ) = 3 2 0 1 y ( s ) d s , z ( 0 ) = z ( 0 ) = z ( 0 ) = z ( 0 ) = z ( 0 ) = 0 , z ( 1 ) = 3 2 0 1 z ( s ) d s .
From the considered problem (37), one has δ = ϱ = 3 2 , as m = 5 . It is easy to see that φ 0 = χ 0 = 0 and φ = χ = . Therefore thanks to Theorem 5, the given system (37) has a positive solution.
Example 4.
Further we take another system of FODEs with IBCs as
D + 0 11 2 y ( ϑ ) + ( 1 + ϑ 2 ) [ u 2 ( ϑ ) + z ( ϑ ) ] ( 4 ϑ 2 + 4 ) α σ = 0 , ϑ ( 0 , 1 ) , D + 0 16 3 z ( ϑ ) + ( ϑ 3 + 1 ) [ y ( ϑ ) + v 2 ( ϑ ) ] ( 4 ϑ 3 + 4 ) α æ = 0 , ϑ ( 0 , 1 ) , y ( 0 ) = y ( 0 ) = y ( 0 ) = y ( 0 ) = y ( 0 ) = 0 , y ( 1 ) = 3 2 0 1 y ( s ) d s , z ( 0 ) = z ( 0 ) = z ( 0 ) = z ( 0 ) = z ( 0 ) = 0 , z ( 1 ) = 3 2 0 1 z ( s ) d s .
where δ = ϱ = 3 2 and m = 6 . It is easy to obtain φ 0 = χ 0 = and φ = χ = .
Further ( ϑ , y , z ) [ 0 , 1 ] × [ 0 , 1 ] × [ 0 , 1 ] , we have
φ ( ϑ , y , z ) ( ϑ 2 + 1 ) 2 4 ( ϑ 2 + 1 ) α σ = α σ 1 2 , χ ( ϑ , y , z ) ( ϑ 3 + 1 ) 2 4 ( ϑ 3 + 1 ) α æ = α æ 1 2 .
Hence all the conditions of Theorem 6 hold. Thanks to Theorem 6, the given system (38) has at least two positive solutions ( y 1 , z 1 ) and ( y 2 , z 2 ) which satisfy
0 < ( y 1 , z 1 ) < 1 < ( y 2 , z 2 ) .

7. Non-Existence of Positive Solution

Here some conditions are developed under which the coupled system (1) with given IBCs has no solution.
Theorem 12.
Consider that ( C 1 ) ( C 3 ) hold and φ ( ϑ , y , z ) < ( y , z ) 2 α σ and χ ( ϑ , y , z ) < ( y , z ) 2 α æ for all ϑ [ 0 , 1 ] , y > 0 , z > 0 , then there is no positive solution for BVPs (1).
Proof. 
Consider ( y , z ) to be the positive solution of BVPs (1). Then, ( y , z ) C for 0 < ϑ < 1 and
( y , z ) = y + z = max ϑ [ 0 , 1 ] | y ( ϑ ) | + max ϑ [ 0 , 1 ] | z ( ϑ ) | max ϑ [ 0 , 1 ] 0 1 H σ ( ϑ , s ) | φ ( s , y ( s ) , z ( s ) ) | d s + max ϑ [ 0 , 1 ] 0 1 H æ ( ϑ , s ) | χ ( s , y ( s ) , z ( s ) | d s < 0 1 H σ ( 1 , s ) ( y , z ) 2 α æ d s + 0 1 H æ ( 1 , s ) ( y , z ) 2 α æ d s ( y , z ) < ( y , z ) ,
which is contradiction. Hence the considered system (1) has no solution. □
Theorem 13.
Let the hypothesis ( C 1 ) ( C 3 ) hold along with the conditions
φ ( ϑ , y ( ϑ ) , z ( ϑ ) ) > ( y , z ) 2 λ α 2 θ 1 θ H σ ( 1 , s ) d s 1 , χ ( ϑ , y ( ϑ ) , z ( ϑ ) ) > ( y , z ) 2 λ β 2 θ 1 θ H æ ( 1 , s ) d s 1 , for all ϑ [ 0 , 1 ] , y > 0 and z > 0 .
Then there does not exist positive solution to BVPs (1).
To demonstrate the results of Theorems 12 and 13 respectively, we give the following example.
Example 5.
Taking the given system of FODEs with given IBCs as
D + 0 5 2 y ( ϑ ) = 5 4 y + z + π 3 5 2 , ϑ [ 0 , 1 ] , D + 0 5 2 z ( ϑ ) = 30 + 30 y + z 3 2 + 1 50 , ϑ [ 0 , 1 ] , y ( 0 ) = y ( 0 ) = y ( 0 ) = 0 , y ( 1 ) = 1 2 0 1 y ( ϑ ) d ϑ , z ( 0 ) = z ( 0 ) = z ( 0 ) = 0 , z ( 1 ) = 1 2 0 1 z ( ϑ ) d ϑ .
Also as ( C 1 ) ( C 3 ) hold, where m = [ 2.5 ] + 1 = 3 and δ = ϱ = 1 2 . We calculate
φ 0 = 5 3 π 5 2 , χ 0 = 1 50 , φ = 5 , χ = 51 , 5 3 π 5 2 ( y , z ) < φ ( ϑ , y ( ϑ ) , z ( ϑ ) ) < 5 ( y , z ) , 1 50 ( y , z ) < χ ( ϑ , y ( ϑ ) , z ( ϑ ) ) < 51 ( y , z ) .
Therefore we have
5 3 π 5 2 ( y , z ) < φ ( ϑ , y ( ϑ ) , z ( ϑ ) ) < 5 ( y , z ) and χ ( ϑ , y ( ϑ ) , z ( ϑ ) ) < 5 ( y , z ) < ( y , z ) α σ ,
where α σ 0.32239 and α æ 0.32239 .
Case I:Now
φ ( ϑ , y ( ϑ ) , z ( ϑ ) ) < ( y , z ) α σ 1.1413 ( y , z )
yields that
φ ( ϑ , y ( ϑ ) , z ( ϑ ) ) < 5 ( y , z ) 3.1018 ( y , z )
and
χ ( ϑ , y ( ϑ ) , z ( ϑ ) ) < 51 ( y , z ) 3.1018 ( y , z ) .
Hence under the condition of Theorem 12, there is no solution corresponding to problem (39).
Case II:Also
φ ( ϑ , y ( ϑ ) , z ( ϑ ) ) > 5 3 π 5 2 ( y , z ) > ( y , z ) λ α 2 1 100 99 100 H α ( 1 , s ) d s 1 0.615 ( y , z )
and
χ ( ϑ , y ( ϑ ) , z ( ϑ ) ) > 1 50 ( y , z ) > ( y , z ) λ æ 2 1 100 99 100 H æ ( 1 , s ) d s 1 0.615 ( y , z ) .
Hence under the condition of Theorem 13, there is no solution corresponding to coupled system (39).

8. Conclusions

In the above research work we have successfully investigated a coupled system of nonlinear FODEs with IBCs for multiplicity results. Further, the aforesaid investigation has been strengthened by developing some conditions under which the solutions of the proposed system are HU-stable. Further some results which demonstrate the conditions of nonexistence of solutions have been established. The whole results have been verified by considering some examples where needed.

Author Contributions

All authors have equal contribution in this paper.

Funding

This research work has been supported financially by KMUTT Fixed Point Research Laboratory, KMUTT-Fixed Point Theory and Applications Research Group, Faculty of Science, King Mongkut’s University of Technology Thonburi (KMUTT), 126 Pracha-Uthit Road, Bang Mod, Thrung Khru, Bangkok 10140, Thailand.

Acknowledgments

All authors are very thankful to the referees for their useful comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley: New York, NY, USA, 1993. [Google Scholar]
  2. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations. In North-Holland Mathematics Studies; Elsevier: Amsterdam, The Nederlands, 2006; Volume 204. [Google Scholar]
  3. Podlubny, I. Fractional Differential Equations, Mathematics in Science and Engineering; Academic Press: New York, NY, USA, 1999. [Google Scholar]
  4. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  5. Lakshmikantham, V.; Vatsala, A.S. Basic theory of fractional differential equations. Nonlinear Anal. 2008, 69, 2677–2682. [Google Scholar] [CrossRef]
  6. Benchohra, M.; Graef, J.R.; Hamani, S. Existence results for boundary value problems with nonlinear fractional differential equations. Appl. Anal. 2008, 87, 851–863. [Google Scholar] [CrossRef]
  7. Shah, K.; Khan, R.A. Study of solution to a toppled system of fractional differential equations with integral boundary conditions. Int. J. Appl. Comput. Math. 2017, 3, 2369–2388. [Google Scholar] [CrossRef]
  8. Ahmad, B.; Nieto, J.J. Existence results for a coupled system of nonlinear fractional differential equations with three-point boundary conditions. Comput. Math. Appl. 2009, 58, 1838–1843. [Google Scholar] [CrossRef]
  9. Agarwal, R.P.; Benchohra, M.; Hamani, S. A survey on existence results for boundary value problems of nonlinear fractional differential equations and inclusions. Acta Appl. Math. 2010, 109, 973–1033. [Google Scholar] [CrossRef]
  10. Li, C.F.; Luo, X.N.; Zhou, Y. Existence of positive solutions of the boundary value problem for nonlinear fractional differential equations. Comput. Math. Appl. 2010, 59, 1363–1375. [Google Scholar] [CrossRef]
  11. Rehman, M.; Khan, R.A. Existence and uniqueness of solutions for multi-point boundary value problems for fractional differential equations. Appl. Math. Lett. 2010, 23, 1038–1044. [Google Scholar] [CrossRef]
  12. Zhong, W.; Lin, W. Nonlocal and multiple-point boundary value problem for fractional differential equations. Comput. Math. Appl. 2010, 59, 1345–1351. [Google Scholar] [CrossRef]
  13. El-Shahed, M.; Shammakh, W.M. Existence of positive solutions for m-point boundary value problem for nonlinear fractional differential equation. In Abstract and Applied Analysis; Hindawi: Cairo, Egypt, 2011; Volume 2011, p. 986575. [Google Scholar]
  14. Khan, R.A.; Shah, K. Existence and uniqueness of solutions to fractional order multi point boundary value problems. Commun. Appl. Anal. 2015, 19, 515–526. [Google Scholar]
  15. Khan, R.A. Three-point boundary value problems for higher order nonlinear fractional differential equations. J. Appl. Math. Inform. 2013, 31, 221–228. [Google Scholar] [CrossRef]
  16. EL-Sayed, A.M.A.; Taher, E.O.B. Positive solutions for a nonlocal multi-point boundary-value problem of fractional and second order. Electr. J. Differ. Equ. 2013, 64, 1–8. [Google Scholar]
  17. Yang, A.; Ge, W. Positive solutions of multi-point boundary value problems of nonlinear fractional differential equation at resonance. J. Korea Soc. Math. Educ. Ser. B Pure Appl. Math. 2009, 16, 181–193. [Google Scholar]
  18. Jankowski, T. Differential equations with integral boundary condition. J. Comput. Appl. Math. 2002, 147, 2274–2280. [Google Scholar] [CrossRef]
  19. Nanware, J.A.; Dhaigude, D.B. Existence and uniqueness of solution of differential equations of fractional order with integral boundary conditions. J. Nonlinear Sci. Appl. 2014, 7, 246–254. [Google Scholar] [CrossRef]
  20. Cabada, A.; Wang, G. Positive solution of nonlinear fractional differential equations with integral boundary value conditions. J. Math. Anal. Appl. 2012, 389, 403–411. [Google Scholar] [CrossRef]
  21. Deng, W. Smoothness and stability of the solutions for nonlinear fractional differential equations. Nonlinear Anal. TMA 2010, 72, 1768–1777. [Google Scholar] [CrossRef]
  22. Li, Y.; Chen, Y.; Podlubny, I. Mittag-Leffler stability of fractional order nonlinear dynamic systems. Automatica 2009, 45, 1965–1969. [Google Scholar] [CrossRef]
  23. Li, Y.; Chen, Y.; Podlubny, I. Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized Mittag-Leffler stability. Comput. Math. Appl. 2010, 59, 1810–1821. [Google Scholar] [CrossRef]
  24. Ulam, S.M. Problems in Modern Mathematics; John Wiley and Sons: New York, NY, USA, 1940; Chapter 6. [Google Scholar]
  25. Hyers, D.H. On the stability of linear functional equation. Proc. Natl. Acad. Sci. USA 1941, 27, 222–224. [Google Scholar] [CrossRef] [PubMed]
  26. Jung, S.M. On the Hyers-Ulam stability of the functional equations that have the quadratic property. J. Math. Anal. Appl. 1998, 222, 126–137. [Google Scholar] [CrossRef]
  27. Jung, S.M. Hyers-Ulam stability of linear differential equations of first order. Appl. Math. Lett. 2006, 19, 854–858. [Google Scholar] [CrossRef]
  28. Obloza, M. Hyers stability of the linear differential equation. Rocznik Nauk. Dydakt. Prace Mat. 1993, 13, 259–270. [Google Scholar]
  29. Krasnoselskii, M.A. Positive Solutions of Operator Equations; Noordhoff: Groningen, The Netherlands, 1964. [Google Scholar]
  30. Urs, C. Coupled fixed point theorem and applications to periodic boundary value problem. Miskolic Math. Notes 2013, 14, 323–333. [Google Scholar] [CrossRef]
  31. Stanêk, S. The existence of positive solutions of singular fractional boundary value problems. Comput. Math. Appl. 2011, 62, 1379–1388. [Google Scholar] [CrossRef]
  32. Shah, K.; Khan, R.A. Multiple positive solutions to a coupled systems of nonlinear fractional differential equations. SpringerPlus 2016, 5, 1–20. [Google Scholar] [CrossRef] [PubMed]
  33. Agarwal, R.; Meehan, M.; Regan, D.O. Fixed Points Theory and Applications; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  34. Guo, D.; Lakshmikantham, V. Nonlinear Problems in Abstract Cones; Academic Press: New York, NY, USA, 1988. [Google Scholar]
  35. Shah, K.; Khan, R.A. Existence and uniqueness of positive solutions to a coupled system of nonlinear fractional order differential equations with anti periodic boundary conditions. Differ. Equ. Appl. 2015, 7, 245–262. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Shah, K.; Kumam, P.; Ullah, I. On Ulam Stability and Multiplicity Results to a Nonlinear Coupled System with Integral Boundary Conditions. Mathematics 2019, 7, 223. https://doi.org/10.3390/math7030223

AMA Style

Shah K, Kumam P, Ullah I. On Ulam Stability and Multiplicity Results to a Nonlinear Coupled System with Integral Boundary Conditions. Mathematics. 2019; 7(3):223. https://doi.org/10.3390/math7030223

Chicago/Turabian Style

Shah, Kamal, Poom Kumam, and Inam Ullah. 2019. "On Ulam Stability and Multiplicity Results to a Nonlinear Coupled System with Integral Boundary Conditions" Mathematics 7, no. 3: 223. https://doi.org/10.3390/math7030223

APA Style

Shah, K., Kumam, P., & Ullah, I. (2019). On Ulam Stability and Multiplicity Results to a Nonlinear Coupled System with Integral Boundary Conditions. Mathematics, 7(3), 223. https://doi.org/10.3390/math7030223

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop