Existence, Uniqueness and Ulam's Stability of Solutions for a ... - MDPI

mathematics Article

Existence, Uniqueness and Ulam’s Stability of Solutions for a Coupled System of Fractional Differential Equations with Integral Boundary Conditions Dimplekumar Chalishajar 1, * and Avadhesh Kumar 2 1 2

*

Department of Mathematics and Computer Science, Mallory Hall, Virginia Military Institute, Lexington, VA 24450, USA School of Basic Sciences, Indian Institute of Technology Mandi, Kamand 175 005, H.P., India; [email protected] Correspondence: [email protected]

Received: 2 May 2018; Accepted: 29 May 2018; Published: 7 June 2018

Abstract: In this paper, the existence and uniqueness of the solutions to a fractional order nonlinear coupled system with integral boundary conditions is investigated. Furthermore, Ulam’s type stability of the proposed coupled system is studied. Banach’s fixed point theorem is used to obtain the existence and uniqueness of the solutions. Finally, an example is provided to illustrate the analytical findings. Keywords: coupled system; green’s function; integral boundary conditions; Ulam’s stability AMS Subject Classification: 34A08; 34B15; 34B27; 35B35

1. Introduction Fractional calculus is a branch of mathematical analysis, in which arbitrary order differential and integral operators are studied. It started with a correspondence between L’Hospital and Leibnitz in 1695. Presently, plenty of literature is available on theoretical as well as numerical work on this topic. It has application in numerous fields, for example, control theory, signal and image processing, aerodynamics and biophysics [1–3]. For the fundamental concepts of fractional calculus, books like Kilbas et al. [4], Miller and Ross [5] and Halfer [6] are referred. Existence and uniqueness of solutions for fractional order differential systems in finite dimensional as well infinite dimensional spaces were studied by several authors [7–11]. Ahmad et al. [8] established existence results for nonlinear boundary value fractional integro-differential equations with integral boundary conditions. Integral boundary conditions have several applications in real-life problems such as population dynamics, blood flow problems, underground water flow, and chemical engineering. For more details on integral boundary conditions, we refer the reader to [12]. Here we would like to consider a practical example of the integral boundary condition:

− x 00 = g(t) f (t, x ), x (0) = 0, βx 0 (1) = x (η ), where t ∈ (0, 1), η ∈ (0, 1] and β is a positive constant. This is the model for a thermostat. Solutions of the problem are stationary solutions for a one-dimensional heat equation, corresponding to a heated bar, with a controller at 1, which adds or removes heat, depending on the temperature detected by a sensor at η. This problem can be generalized. One can consider the heat equation with nonlinear gradient source terms that vary in time. Moreover, now, the heated bar, with a controller at 1, which adds Mathematics 2018, 6, 96; doi:10.3390/math6060096

www.mdpi.com/journal/mathematics

Mathematics 2018, 6, 96

2 of 12

or removes heat depending on the temperature detected by sensors located at any points of the bar (it depends on how we define the function h). This problem can be written in the form x 00 = f (t, x, x 0 ), x (0) = 0, x 0 (1) =

Z 1 0

x (s)dh(s)(η ).

Recently, Ulam’s type stability has been of great interest to many researchers. In 1940, the above mentioned stability was first introduced by Ulam [13]. Then, it was explained by Hyers [14] in the subsequent years. Ulam and Hyers studied it for various kinds of differential equations with integer order. Nowadays, we describe the result of Hyers simply saying that Cauchy functional equation is Hyers-Ulam stable (or has the Hyers-Ulam stability). Next, Hyers and Ulam published some further stability results for polynomial functions, isometries, and convex functions. The Hyers’ results are extended and generalized by many researchers for integer order differential equations. Plenty of significant results on Ulam’s type stability can be found in the literature, we refer [15–17] and references cited therein. To the best of our knowledge, there are only few manuscripts devoted to the study of Ulam’s type stability for coupled system of fractional differential equations. Further, there is no manuscript considering the Ulam’s stability for coupled system of fractional order α ∈ (1, 2] differential equations with integral boundary conditions. Motivated by this fact, in this paper, the existence, uniqueness of solutions as well as Ulam’s type stability for the considered coupled system involving Caputo derivative is studied. The proposed system is given as follows:  c D α x ( t ) = f ( t, y ( t )), α ∈ (1, 2], t ∈ J         c β    D y(t) = g(t, x (t)), β ∈ (1, 2], t ∈ J

R1    px (0) + qx 0 (0) = 0 a1 ( x (s))ds,        R1  ˜ (0) + qy ˜ 0 (0) = 0 a˜ 1 (y(s))ds, py

(1)

=

R1

˜ (1) + qy ˜ 0 (1) = py

R1

px (1) + qx 0 (1)

0

0

a2 ( x (s))ds, a˜ 2 (y(s))ds,

where J = [0, 1] and f , g : J × R → R are continuous functions. Here, a1 , a2 , a˜ 1 , a˜ 2 : R → R and p, p˜ > 0; q, q˜ ≥ 0 are real numbers. The plan of the paper is as follows. In the second section, some useful definitions, notations, lemmas and results are given which will be required for the later sections. In the third section, existence and uniqueness of the solutions for the coupled system (1) is studied. In the fourth section, Ulam’s type stability results are obtained. In the last section, a few examples are given to show the application of the obtained abstract results. 2. Preliminaries and Assumptions In this section, some useful definitions, notations and lemmas are briefly reviewed. Definition 1. [4] For any function z ∈ ((0, 1), R), the Caputo derivative of fractional order α ∈ R+ is defined as Z t 1 c α D z(t) = (t − s)n−α−1 z(n) (s)ds, n = [α] + 1, Γ(n − α) 0 where [α] denotes the integer part of α and Γ(·) is the gamma function. Definition 2. [4] The Riemann-Liouville fractional integral of order α ∈ R+ is defined by J α z(t) =

1 Γ(α)

Z t 0

(t − s)α−1 z(s)ds,


3 of 12

where z(t) ∈ L1 ([0, 1], R+ ). Lemma 1. [4] For any α > 0, then the differential equations c

D α z(t) = 0

has solution given by z(t) = c0 + c1 t + c2 t2 + · · · + cn−1 tn−1 , ci ∈ R, i = 0, 1, · · · , n − 1, where n = [α] + 1. Lemma 2. [4] For any α > 0, then the solution of the differential equations c

D α z(t) = h(t)

will be given by J α [c D α z(t)] = J α z(t) + c0 + c1 t + c2 t2 + · · · + cn−1 tn−1 , ci ∈ R, i = 0, 1, · · · , n − 1, where n = [α] + 1. Lemma 3. [8] For any h, γ1 , γ2 ∈ C ([0, 1], R), the unique solution of the boundary value problem (

c D α z(t)

= h(t), α ∈ (1, 2], t ∈ [0, 1] R1 R1 pz(0) + qz0 (0) = 0 γ1 (s)ds, pz(1) + qz0 (1) = 0 γ2 (s)ds,

(2)

is given by z(t) =

Z 1 0

Gα (t, s)h(s)ds +

1h ( p (1 − t ) + q ) p2

Z 1 0

γ1 (s)ds + (q + pt)

Z 1 0

i γ2 (s)ds ,

(3)

where Gα (t, s) is the Green’s function given by

Gα (t, s) =

   

p(t−s)α−1 +(q− pt)(1−s)α−1 pΓ(α)

  

(q− pt)(1−s)α−1 pΓ(α)

+

q(q− pt)(1−s)α−2 , p2 Γ ( α −1)

s ≤ t, (4)

+

q(q− pt)(1−s)α−2 , p2 Γ ( α −1)

t ≤ s.

Lemma 4. The space B = {z(t) | z ∈ C (J )} is a Banach space under the defined norm kzkB = maxt∈J |z(t)|. Similarly, the norm on product space is defined by k(z, z˜ )kB×B = kzkB + kz˜ kB . Obviously (B × B , k(·, ·)kB×B ) is a Banach space. Further, the cone C ⊂ B × B is defined by

C = {(z, z˜) ∈ B × B | z(t) ≥ 0, z˜(t) ≥ 0}. Here, the problem (1) is transformed into a fixed point problem. Let F : B × B → B × B be the operator defined as


4 of 12

R 1 0

F ( x, y)(t) =

=

1 p2

h

( p (1 − t ) + q ) i R1 +(q + pt) 0 a2 ( x (s))ds

Gα (t, s) f (s, y(s))ds +

R1 0

a1 ( x (s))ds



           R h R1   1 1  0 Gβ (t, s) g(s, x (s))ds + ˜ 2 ( p˜ (1 − t) + q˜) 0 a˜ 1 (y(s))ds p   i R1 ˜ ) 0 a˜ 2 (y(s))ds +(q˜ + pt   Fα (y, x )(t)    . F β ( x, y)(t)

(5)

Then the fixed point of the operator F coincides with the solution of coupled system (1). In order to prove the existence and uniqueness of solutions of coupled system (1), following assumptions are taken: (A1) For t ∈ J , there exist λ, µ ∈ C (J , R+ ), such that

| f (t, y(t))| ≤ λ(t) + µ(t)|y(t)|, ∀ y(t) ∈ C (J , R) with λ∗ = supt∈J λ(t), µ∗ = supt∈J µ(t). Similarly, for t ∈ J , there exist ν, ξ ∈ C (J , R+ ), such that

| g(t, x (t))| ≤ ν(t) + ξ (t)| x (t)|, ∀ x (t) ∈ C (J , R) with ν∗ = supt∈J ν(t), ξ ∗ = supt∈J ξ (t). (A2) For t ∈ J , there exist positive constant L a1 , L a2 , such that

| a1 ( x (t))| ≤ L a1 | x (t)| and | a2 ( x (t))| ≤ L a2 | x (t)|, ∀ x (t) ∈ C (J , R) Similarly, For t ∈ J , there exist positive constant L a˜ 1 , L a˜ 2 , such that

| a˜ 1 (y(t))| ≤ L a˜1 |y(t)| and | a˜ 2 ( x (t))| ≤ L a˜2 |y(t)|, ∀ y(t) ∈ C (J , R) (A3)

P1 =

Z 1 0

| Gα (t, s)|λ(s)ds < ∞, Q1 =

Z 1 0

| Gα (t, s)|µ(s)ds +

( p + q)( L a1 + L a2 ) 1 < 2 2 p

| Gβ (t, s)|ξ (s)ds +

( p˜ + q˜)( L a˜1 + L a˜2 ) 1 < . 2 2 p˜

and

P2 =

Z 1 0

| Gβ (t, s)|ν(s)ds < ∞, Q2 =

Z 1 0

(A4) For all y, y˜ ∈ C (J , R) and for each t ∈ J there exists a positive constant K f , such that

| f (t, y) − f (t, y˜ )| ≤ K f |y − y˜ |. Similarly, for all x, x˜ ∈ C (J , R) and for each t ∈ J there exists a positive constant K g , such that

| g(t, x ) − g(t, x˜ )| ≤ K g | x − x˜ |. (A5) For all x, x˜ ∈ C (J , R) and for each t ∈ J there exist positive constants Ka1 , Ka2 , such that

| a1 ( x ) − a1 ( x˜ )| ≤ Ka1 | x − x˜ | and | a2 ( x ) − a2 ( x˜ )| ≤ Ka2 | x − x˜ |.


5 of 12

Similarly, for all y, y˜ ∈ C (J , R) and for each t ∈ J there exist positive constants Ka˜ 1 , Ka˜ 2 , such that | a˜ 1 (y) − a˜ 1 (y˜ )| ≤ Ka˜1 |y − y˜ | and | a˜ 2 (y) − a˜ 2 (y˜ )| ≤ Ka˜2 |y − y˜ |. (A6) Let (i) φ1 = K( β p˜ q˜) K g + K pq , where K pq =

K( β p˜ q˜)

( p+q)(Ka1 +Ka2 ) p2

and

Z 1 Gβ (t, s)ds max t∈[0,1] 0 Z " # t p˜ (t − s) β−1 + (q˜ − pt ˜ )(1 − s) β−1 ˜ )(1 − s) β−2 q˜(q˜ − pt + ds 0 ˜ ( β) pΓ p˜ 2 Γ( β − 1) Z " # 1 (q˜ − pt ˜ )(1 − s) β−1 ˜ )(1 − s) β−2 q˜(q˜ − pt + ds t ˜ ( β) pΓ p˜ 2 Γ( β − 1)

= = +

1 2( p˜ + q˜) 2(q˜2 + p˜ q˜) ; + + Γ ( β + 1) Γ ( β + 1) p˜ 2 Γ( β)

=

(ii) φ2 = K(αpq) K f + K p˜ q˜ , where K p˜ q˜ =

K(αpq)

( p˜ +q˜)(Ka˜ 1 +Ka˜ 2 ) p˜ 2

Z max

=

t∈[0,1]

1 0

and

Gα (t, s)ds

2( p + q ) 2(q2 + pq) 1 + + . Γ ( α + 1) Γ ( α + 1) p2 Γ ( α )

= 3. Existence and Uniqueness Analysis

Theorem 1. If all the assumptions (A1)–(A6) and φ = max{φ1 , φ2 } < 1 are fulfilled, then the fractional order coupled system (1) has a unique solution. Proof. For a positive number δ = max

2P1 2 P2 , , 1 − 2Q1 1 − 2Q2

we define a set

W = {( x, y) ∈ B × B : k( x, y)kB×B ≤ δ}. First, in order to prove that F maps W into itself, we have

|Fα (y, x )(t)|

≤

Z 1 0

+ ≤

1h |( p(1 − t) + q)| p2 Z 1 0

+ ≤

Z 1

Z 1 0

| Gα (t, s)|λ(s)ds +

1h ( p + q) p2 0

=

| Gα (t, s)|| f (s, y(s))|ds

Z 1 0

Z 1 0

Z 1 0

| a2 ( x (s))|ds

i

| Gα (t, s)|[µ(s)|y(s)|]ds

| a1 ( x (s))|ds + ( p + q) "Z

| Gα (t, s)|λ(s)ds + δ

δ P1 + δQ1 ≤ . 2

| a1 ( x (s))|ds + |(q + pt)|

1

0

Z 1 0

| a2 ( x (s))|ds

i

( p + q)( L a1 + L a2 ) | Gα (t, s)|µ(s)ds + p2

#

(6)


6 of 12

Now taking maximum on both side of the inequality (6) over J , we obtain δ kFα (y, x )kB ≤ . 2 Similarly, kF β ( x, y)kB ≤ 2δ . Hence, we can conclude that

kF ( x, y)kB×B ≤ δ.

(7)

Inequality (7) shows that F maps W into itself. Next, in order to show that F is the contraction operator when t ∈ J , we have ¯ x¯ )(t)| |Fα (y, x )(t) − Fα (y,

≤

Z 1 0

+

| Gα (t, s)|| f (s, y(s)) − f (s, y¯ (s))|ds

1h |( p(1 − t) + q)| p2

+ |(q + pt)| ≤ +

Z 1 0

Z 1 0

| a1 ( x (s)) − a1 ( x¯ (s))|ds

| a2 ( x (s)) − a2 ( x¯ (s))|ds

i

K(αpq) K f |y(t) − y¯ (t)| ( p + q)(Ka1 + Ka2 ) | x (t) − x¯ (t)|. p2

(8)

When we take maximum on both side of the inequality (8) over J , we obtain ¯ x¯ )kB kFα (y, x ) − Fα (y,

≤ K(αpq) K f ky − y¯ kB +

( p + q)(Ka1 + Ka2 ) k x − x¯ kB . p2

(9)

( p˜ + q˜)(Ka˜1 + Ka˜2 ) ky − y¯ kB . p˜ 2

(10)

Similarly, the following can be obtained ¯ y¯ )kB kF β ( x, y) − F β ( x,

≤ K( β p˜ q˜) K g k x − x¯ kB +

From (9) and (9), we get ¯ y¯ )kB×B kF ( x, y) − F ( x,

¯ y¯ )kB×B . ≤ φk( x, y) − ( x,

(11)

Thus, the operator F is strict contraction. By Banach’s fixed point method, it has a unique fixed point which is the unique solution of the considered coupled system (1). 4. Ulam’s Stability Analysis In this section, we study Ulam’s type stability for the coupled system (1). For some e = (eα , e β ) > 0, we consider the following inequality  c α   D x (t) − f (t, y(t)) ≤ eα , t ∈ J ,   c D β y(t) − g(t, x (t)) ≤ e , t ∈ J . β

(12)

The following definitions are inspired by Rus [18]. Definition 3. The coupled system (1) is said to be Ulam-Hyers stable, if there exist K(αβpq p˜ q˜) = (K(αpq) , K( β p˜ q˜) ) > 0 such that for every solution ( x, y) ∈ B × B of the inequality (12), there exists a unique solution (ϑ, κ ) ∈ B × B with ( x, y)(t) − (ϑ, κ )(t) ≤ K (13) (αβpq p˜ q˜) e, t ∈ J .


7 of 12

Definition 4. The coupled system (1) is said to be generalized Ulam-Hyers stable, if there exist Ψ ∈ C (R+ , R+ ) with Ψ(0) = 0, such that for every solution ( x, y) ∈ B × B of the inequality (12), there exist a unique solution (ϑ, κ ) ∈ B × B of the system (1) which satisfies ( x, y)(t) − (ϑ, κ )(t) ≤ Ψ(e), t ∈ J .

(14)

Remark 1. Let ( x, y) ∈ B × B is a solution of the system of inequality (12) if there exist functions φ, ψ ∈ C (J , R) which depend upon x and y respectively, such that (R1) |φ(t)| ≤ eα , |ψ(t)| ≤ e β , t ∈ J ; (R2) and  c α   D x (t) = f (t, y(t)) + φ(t), t ∈ J ,   c D β y(t) = g(t, x (t)) + ψ(t), t ∈ J . Lemma 5. Let ( x, y) ∈ B × B be the solution of the inequality (12), then the following inequality will be satisfied: h  R1 R1 1  ( t ) − G ( t, s ) f ( s, y ( s )) ds − x α 2 ( p (1 − t ) + q ) 0 a1 ( x ( s )) ds  0 p  i  R1    +( q + pt ) a ( x ( s )) ds ≤ K(αpq) eα , t ∈ J ,  2 0  h  R1 R1    (t) − 0 Gβ (t, s) g(s, x (s))ds − p˜12 ( p˜ (1 − t) + q˜) 0 a˜ 1 (y(s))ds y   i  R1   ˜ ) e , t ∈ J. +(q˜ + pt a˜ (y(s))ds ≤ K 0

2

( β p˜ q˜) β

Proof. By Remark 1 (R2), we have  c D α x ( t ) = f ( t, y ( t )) + φ ( t ), t ∈ J ,         c β    D y(t) = g(t, x (t)) + ψ(t), t ∈ J ,

R1    px (0) + qx 0 (0) = 0 a1 ( x (s))ds,        R1  ˜ (0) + qy ˜ 0 (0) = 0 a˜ 1 (y(s))ds, py

(15)

px (1) + qx 0 (1) =

R1

˜ (1) + qy ˜ 0 (1) = py

R1

0

0

a2 ( x (s))ds a˜ 2 (y(s))ds

By Applying Lemma 3, the solution of (15) will be as follows:  R1 R1  x (t) = 0 Gα (t, s) f (s, y(s))ds + 0 Gα (t, s)φ(s)ds +   i  R1    +( q + pt ) a ( x ( s )) ds , t ∈ J,  2 0   R1 R1    y(t) = 0 Gβ (t, s) g(s, x (s))ds + 0 Gα (t, s)ψ(s)ds +   i  R1   ˜ ) 0 a˜ 2 (y(s))ds , t ∈ J . +(q˜ + pt Considering first equation of the system (16), we have

1 p2

h

1 p˜ 2

h

( p (1 − t ) + q )

R1 0

a1 ( x (s))ds (16)

( p˜ (1 − t) + q˜)

R1 0

a˜ 1 (y(s))ds


Z x (t) −

1 0

8 of 12

Gα (t, s) f (s, y(s))ds −

1h ( p (1 − t ) + q ) p2

+(q + pt)

Z 1 0

Z 1 0

a1 ( x (s))ds Z 1 Gα (t, s)φ(s)ds ≤

i a2 ( x (s))ds

0

≤

Z 1 0

| Gα (t, s)||φ(s)|ds.

By Remark 1 (R1) and using condition of (A6), we get Z x (t) −

1 0

Gα (t, s) f (s, y(s))ds −

1h ( p (1 − t ) + q ) p2

+(q + pt)

Z 1 0

Z 1

a1 ( x (s))ds

0

i a2 ( x (s))ds ≤ K(αpq) eα .

(17)

Repeating the same procedure for second equation of the system (16), we get Z y(t) −

1 0

1 1h ˜ ˜ ( p ( 1 − t ) + q ) a˜ 1 (y(s))ds p˜ 2 0 Z 1 i ˜ ) +(q˜ + pt a˜ 2 (y(s))ds ≤ K( β p˜ q˜) e β .

Z

Gβ (t, s) g(s, x (s))ds −

0

(18)

Theorem 2. If all the assumptions (A4)–(A6) are fulfilled, then the fractional order coupled system (1) is Ulam-Hyers stable and consequently generalized Ulam-Hyers stable provided that

(1 − K pq )(1 − K p˜ q˜ ) − Kαpq K β p˜ q˜ 6= 0. Proof. Let ( x, y) ∈ B × B be the solution of the system (15) and (ϑ, κ ) ∈ B × B be the unique solution to the following considered system:  c D α ϑ ( t ) = f ( t, κ ( t )), t ∈ J ,         c β    D κ (t) = g(t, ϑ (t)), t ∈ J ,

R1    pϑ (0) + qϑ0 (0) = 0 a1 (ϑ (s))ds,        R1  ˜ (0) + qκ ˜ 0 (0) = 0 a˜ 1 (κ (s))ds, pκ

(19) pϑ (1) + qϑ0 (1) =

R1

˜ (1) + qκ ˜ 0 (1) = pκ

R1

0

0

a2 (ϑ (s))ds

a˜ 2 (κ (s))ds

Using Lemma 3, the solution of (19) h  R1 R1 1  ( p(1 − t) + q) 0 a1 (ϑ(s))ds ϑ ( t ) = G ( t, s ) f ( s, κ ( s )) ds + α 2  0 p  i  R1    +( q + pt ) a ( ϑ ( s )) ds , t ∈ J,  2 0  h  R1 R1    κ (t) = 0 Gβ (t, s) g(s, ϑ (s))ds + p˜12 ( p˜ (1 − t) + q˜) 0 a˜ 1 (κ (s))ds   i  R1   ˜ ) 0 a˜ 2 (κ (s))ds , t ∈ J . +(q˜ + pt We have

(20)


9 of 12

| x (t) − ϑ(t)|

=

Z x (t) −

1 0

Gα (t, s) f (s, κ (s))ds −

1h ( p (1 − t ) + q ) p2

Z 1 0

a1 (ϑ (s))ds

Z 1 i + (q + pt) a2 (ϑ (s))ds 0 Z 1 Z 1 1h Gα (t, s) f (s, y(s))ds − 2 ( p(1 − t) + q) ≤ x (t) − a1 ( x (s))ds p 0 0 Z 1 i + (q + pt) a2 ( x (s))ds 0 Z 1 Z 1 Gα (t, s) f (s, y(s))ds − Gα (t, s) f (s, κ (s))ds + 0 0 Z 1 Z 1 1h i + 2 ( p (1 − t ) + q ) a1 ( x (s))ds + (q + pt) a2 ( x (s))ds p 0 0 Z 1 Z 1 h i 1 − ( p ( 1 − t ) + q ) a ( ϑ ( s )) ds + ( q + pt ) a ( ϑ ( s )) ds 2 1 p2 0 0 ≤ Kαpq eα + Kαpq |y(t) − κ (t)| + K pq | x (t) − ϑ(t)|,

where Kαpq = K(αpq) K f . Hence, we get

(1 − K pq )k x − ϑkB ≤ K(αpq) eα + Kαpq ky − κ kB .

(21)

(1 − K p˜ q˜ )ky − κ kB ≤ K( β p˜ q˜) eβ + K β p˜ q˜ k x − ϑkB ,

(22)

Similarly, we have

where K β p˜ q˜ = K( β p˜ q˜) K g . From (21) and (22), it can be written as    (1 − K pq )k x − ϑ kB − Kαpq ky − κ kB ≤ K(αpq) eα

(23)

  (1 − K )ky − κ k − K k x − ϑ k ≤ K p˜ q˜ B B β p˜ q˜ ( β p˜ q˜) e β The matrix representation of (23) is as follows   

(1 − K pq ) −K β p˜ q˜



−Kαpq (1 − K p˜ q˜ )

k x − ϑ kB

 

ky − κ kB





  ≤

K(αpq) eα K( β p˜ q˜) eβ

  .

After simplification of the above inequality, we have   

k x − ϑ kB ky − κ kB

 (1− K ˜ ˜ )



  ≤

∆

pq

K β p˜ q˜ ∆

Kαpq ∆

(1−K pq ) ∆

  

K(αpq) eα K( β p˜ q˜) eβ

  ,

where ∆ = (1 − K pq )(1 − K p˜ q˜ ) − Kαpq K β p˜ q˜ 6= 0. Further simplification gives

k x − ϑ kB

≤

ky − κ kB

≤

(1 − K p˜ q˜ )K(αpq) eα Kαpq K( β p˜ q˜) eβ + ∆ ∆ K β p˜ q˜ K(αpq) eα ∆

+

(1 − K pq )K( β p˜ q˜) eβ ∆

(24)

(25)


10 of 12

From inequalities (24) and (25), we have

k x − ϑ kB + ky − κ kB

≤ +

(1 − K p˜ q˜ )K(αpq) eα Kαpq K( β p˜ q˜) eβ K β p˜ q˜ K(αpq) eα + + ∆ ∆ ∆ (1 − K pq )K( β p˜ q˜) eβ ∆

Therefore, we have

k( x, y) − (ϑ, κ )kB×B ≤ K(αβpq p˜ q˜) e,

(26)

where e = max{eα , e β } and

K(αβpq p˜ q˜) =

(1 − K p˜ q˜ )K(αpq) Kαpq K( β p˜ q˜) K β p˜ q˜ K(αpq) (1 − K pq )K( β p˜ q˜) + + + . ∆ ∆ ∆ ∆

Hence, by inequality (26), we can conclude that the coupled system (1) is Ulam-Hyers stable. Further, inequality (26) can be written as

k( x, y) − (ϑ, κ )kB×B ≤ Ψ(e), where Ψ(0) = 0.

(27)

By inequality (27), we further conclude that the coupled system (1) is generalized Ulam-Hyers stable.

5. Application Example 1. We consider the following fractional order coupled system:  c D α x (t) =           c β   D y(t) =

|y(t)| 1 , (t+7)2 1+|y(t)| 1 100

α ∈ (1, 2], t ∈ J = [0, 1]

t cos x (t) + x (t) sin t , β ∈ (1, 2], t ∈ J = [0, 1] (28)

 R 1 | x(s)| R 1 | x(s)|   x (0) + x 0 (0) = 0 13+| x(s)| ds, x (1) + x 0 (1) = 0 15+| x(s)| ds         R1 1 R1  y(0) + y0 (0) = 0 26 cos y(s) + sin y(s) ds, y(1) + y0 (1) = 0

1 30

cos y(s) + sin y(s) ds

By comparing the coupled systems (28) to (1), the following values are derived: p = q = p˜ = q˜ = 1, Ka1 = Ka˜ 1 =

1 1 and Ka2 = Ka˜ 2 = . 13 15

Here, f (t, y(t)) =

1 |y(t)| 1 and g(t, x (t)) = t cos x (t) + x (t) sin t . 2 100 (t + 7) 1 + |y(t)|

1 As, | f (t, y) − f (t, y˜ )| ≤ 49 |y − y˜ | and | g(t, x ) − g(t, x˜ )| ≤ 1 1 K f = 49 and K g = 50 . Further, we have

φ1

= = < 1.

1 50 | x − x˜ |,

therefore (A4) is satisfied with

! ( p + q)(Ka1 + Ka2 ) 1 2( p˜ + q˜) 2(q˜2 + p˜ q˜) + + Kg + 2 Γ ( β + 1) Γ ( β + 1) p˜ Γ( β) p2 ! 1 1 4 4 1 1 + + +2 + Γ( β + 1) Γ( β + 1) Γ( β) 50 13 15


11 of 12

Similarly, φ2 < 1. Hence, the coupled system (28) has a unique solution. Moreover, the condition (1 − K pq )(1 − K p˜ q˜ ) − Kαpq K β p˜ q˜ 6= 0 in Theorem 2 is also satisfied. Therefore, coupled system (28) is Ulam-Hyers stable as well as generalized Ulam-Hyers stable. 6. Conclusions Here we have studied the existence and uniqueness of the solutions as well as the stability for a coupled system of fractional order α ∈ (1, 2] differential equation with integral boundary conditions. We have discussed two types of stability, called Ulam-Hyers stability and generalized Ulam-Hyers stability. As a future work, one can generalize the same concept of stability to a neutral time delay system/inclusion as well as state delay system/inclusion (finite and infinite delay), which have some useful scientific applications. This will enhance a new direction of research: a special kind of phase space to be used for the study of controllability and stability of an infinite delay system/inclusion. Author Contributions: Both the authors have contributed equally in the article. Acknowledgments: We authors are thankful to referees for their valuable suggestions and comments which greatly improved the manuscript. Conflicts of Interest: The authors declare no conflicts of interest.

References 1. 2. 3. 4. 5. 6. 7. 8. 9. 10.

11.

12. 13. 14. 15. 16.

Muslim, M.; Kumar, A. Controllability of fractional differential equation of order α ∈ (1, 2] with non-instantaneous impulses. Asian J. Control 2017. [CrossRef] Yang, Q.; Chen, D. Fractional calculus in image processing: A review. Fract. Calcul. Appl. Anal. 2016, 19, 1222–1249. [CrossRef] Magin, R.L. Fractional Calculus in Bioengineering; Begell House: Redding, CA, USA, 2006. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematics Studies; Elsevier Science Inc.: New York, NY, USA, 2006; Volume 204. Miller, K.S.; Ross, B. An Introduction to Fractional Calculus and Fractional Differential Equations; A Wiley Interscience Publication, Wiley: New York, NY, USA, 1993. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. Khan, A.; Shah, K.; Li, Y.; Khan, T.S. Ulam type stability for a coupled system of boundary value problems of nonlinear fractional differential. Equ. J. Funct. Spaces 2017, 2, 1–8. [CrossRef] Ahmad, B.; Nieto, J.J. Existence results for nonlinear boundary value problems of fractional integro-differential equations with integral boundary conditions. Bound. Value Probl. 2009, 2009. [CrossRef] Muslim, M.; Kumar, A.; Agarwal, R.P. Exact controllability of fractional integro-differential systems of order α ∈ (1, 2] with deviated argument. Analele Universitˇatii Oradea 2017, XXIV, 185–194. Chalishajar, D.N.; Karthikeyan, K. Boundary value Problems for Impulsive Fractional Evolution Integrodifferential Equations with Gronwall’s inequality in Banach spaces. J. Discontinuity Nonlinearity Complex. 2014, 3, 33–48. [CrossRef] Chalishajar, D.N.; Karthikeyan, K. Existence and Uniqueness Results for Boundary Value Problems of Higher order Fractional Integro-differential Equations Involving Gronwall’s Inequality in Banach spaces. Acta Math. Sci. Ser. A 2013, 33B, 758–772. [CrossRef] Ahmad, B.; Alsaedi, A.; Alghamdi, B.S. Analytic approximation of solutions of the forced Duffing equation with integral boundary conditions. Nonlinear Anal. Real World Appl. 2008, 9, 1727–1740. [CrossRef] Ulam, S.M. A Collection of the Mathematical Problems; Interscience: New York, NY, USA, 1960. Hyers, D.H. On the stability of the linear functional equation. Proc. Natl. Acad. Sci. USA 1941, 27, 222–224. [CrossRef] [PubMed] Jung, S.M. Hyers-Ulam stability of linear differential equations of first order. Appl. Math. Lett. 2006, 19, 854–858. [CrossRef] Wang, J.; Feckan, M. A general class of impulsive evolution equations. Topol. Methods Nonlinear Anal. 2015, 46, 915–933. [CrossRef]


17. 18.

12 of 12

Muslim, M.; Kumar, A.; Feˇckan, M. Existence, uniqueness and stability of solutions to second order nonlinear differential equations with non-instantaneous impulses. J. King Saud Univ. Sci. 2016. [CrossRef] Rus, I.A. Ulam stabilities of ordinary differential equations in a Banach space. Carpath. J. Math. 2010, 26, 103–107. c 2018 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access

article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).