Stability analysis for boundary value problems with generalized nonlocal condition via Hilfer–Katugampola fractional derivative

*Correspondence: poom.kum@kmutt.ac.th 2Center of Excellence in Theoretical and Computational Science (TaCS-CoE), Science Laboratory Building, King Mongkut’s University of Technology Thonburi (KMUTT), Bangkok, Thailand 4Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, Taiwan Full list of author information is available at the end of the article Abstract


Introduction
Fractional calculus is a branch of mathematics that examines the properties of derivatives and integrals of arbitrary non-integer order. This theory has been studied over the years and was initiated in a letter written to Guillaume de l'Hopital by Gottfried Wilhelm Leibniz in 1965. Subsequently, the area has become a great area of interest to many famous mathematicians such as Fourier, Laplace, Abel, Liouville, Riemann, Letnikov etc. For more interpretation on the fundamentals ideas of fractional calculus, physical and geometric meanings we refer the reader to see the books, [29,37,39,54] etc.
There are various definitions of fractional order derivatives and integrals, for example, we have the Riemann, Liouville, Caputo, Riemann-Liouville, Weyl, Hadamard, and Grunwald-Letnikov versions. These derivatives do not coincide in general except for some special cases. This is due to the fact that the authors try to preserve some properties of these operators. However, the most famous definitions are the ones for Riemann, Liouville, Riemann-Liouville, and Caputo fractional order derivatives and integrals.
Fractional differential equations are an area that is related to the fractional calculus and it played an important role in numerous fields of engineering and scientific areas such as physics, chemistry, biology, economics, control theory, signal and image processing. Recently, it was proved that fractional derivatives and integrals of arbitrary order are more advantageous and pragmatic for characterizing the memory and hereditary properties of different materials and processes than the integer-order derivatives; see [14,18,32].
Existence, uniqueness and stability analyses of fractional differential equations with different types of initial and boundary conditions belong to the most essential parts in the theory of fractional differential equations to which this paper is devoted. However, many researchers studied the existence and uniqueness of solutions for boundary value problems (BVPs) of fractional differential equations with different types of fractional integrals and derivatives; see for example [5, 7, 13, 17, 25, 33, 36, 40-43, 45, 52, 55]. Hilfer initiated a new class of Riemann-Liouville fractional derivatives, called Hilfer fractional derivatives. These fractional derivatives coincide with the classical Riemann-Liouville and Caputo fractional derivatives. The main advantage of these operators is the degree of freedom on the initial conditions. Thus, based on this setting, it motivated a lot of researchers to initiate research and finding the existence, uniqueness and stability results; see [2-4, 8, 10, 15, 16, 21, 27, 48, 50].
More recently, motivated by Hilfer fractional derivatives and some well-known results in respect of the generalized classical definitions of Riemann-Liouville and Caputo fractional derivatives, Oliveira and Capelas [35] proposed new fractional derivatives called Hilfer-Katugampola fractional derivatives, which incorporate Hilfer and generalized fractional derivatives. These have been contemplated by some researchers in the field of fractional calculus, results can be found in [1,9,20,23,24,26,28,53] and the references therein. Stability analysis is one of the most important aspects of fractional differential equations. This notion was initiated by Ulam [22], which motivated many researchers to make contributions in this area [6,8,10,11,21,28,31,38,46,47,51]. For instance, Ahmed et al. [9] developed the existence and uniqueness criteria for solutions of boundary value problems (BVPs) with a generalized nonlocal integral and multi-point conditions described by (1.1) and the inclusion problem where ρ D ν 0 + (·) is the generalized fractional derivatives of order (1 < ν ≤ 2) and ρ I q 0 + (·) is a generalized fractional integral of order q > 0, δ = s 1-ρ d ds , ρ > 0, g : (0, T] × R → R is a continuous function, T > 0 and G : J × R → P(R) is a multivalued function. Recently, Suphawat et al. [12] initiated research and found existence and uniqueness results for the boundary value problems (BVPs) which involve Hilfer fractional derivatives and a nonlocal condition described by where H D p,q a + (·) is the Hilfer fractional derivatives of order (1 < p < 2) and the parameter 0 ≤ q ≤ 1, I r k a + (·) is the Riemann-Liouville fractional integral of order r k > 0, Θ : I × R → R is a continuous function. The authors employed a fixed point approach of the Banach contraction principle, and the Hölder, Boyd and Wang, Leray-Schauder and Krasnoselskii's approaches to obtain their results.
Motivated by the above series of papers and some basic ideas of the fractional calculus, the main objectives of this paper is to extend the results of [9,12,19] and also to discuss the stability analysis in the framework of Ulam-Hyers and generalized Ulam-Hyers theory. More specifically, we consider the fractional differential equation given by where ρ D p,q a + (·) is the Hilfer-Katugampola fractional derivative of order (1 < p ≤ 2) and parameter 0 ≤ q ≤ 1, ρ I ξ i a + (·) is the generalized fractional integral of order ξ i > 0, f : [a, b] × R → R is a continuous function, μ i ∈ R and δ i ∈ [a, b] satisfying 0 < a < δ 1 ≤ δ 2 ≤ · · · ≤ δ m < b for i = 1, 2, . . . , m. Problem (1.6) is simply called a boundary value problem (BVP).
Evidently, to the best of our knowledge, stability analysis of boundary value problems (BVPs) is still in the initial stage. Thus the main contribution of this paper is to investigate the stability analysis of Ulam-Hyers and generalized Ulam-Hyers theory, which is not established in [9,12,19]. Besides, we also show the equivalence between problem (1.6) and the Volterra integral equation and apply the Banach and Krasnoselskii's fixed point theorems to establish the existence and uniqueness results. Moreover, problem (1.6) coincides with: • the Riemann-Liouville and Hadamard fractional derivative when q = 0, ρ → 1; • the Hilfer and Hilfer-Hadamard fractional derivative (ρ → 1, and ρ → 0 + ); • the generalized (defined by Katugampola) and Caputo-type versions (q = 0) and (q = 1). We organize the remaining part of the paper as follows: Some preliminary facts, vital definitions, and theoretical results are recalled in Sect. 2. In Sect. 3, we derive the relation between the Volterra integral equation and the proposed problem (1.6). Also, with the help of these equations, the existence and uniqueness of solutions are proven using the concepts of the Banach and Krasnoselskii fixed point theorems. In Sect. 4, we present certain criteria under which the proposed problem (1.6) is both Ulam-Hyers and generalized Ulam-Hyers stable. Furthermore, as an application, two examples are given to show the applicability of the obtained results. Finally, we summarize the paper in the last section.

Preliminaries and theoretical results
This section summarizes the basic definitions and principles of fractional derivatives and integrals and describes the theoretical findings that are useful in this paper. For further information, see [29,35,37].
which is referred to as a left-sided generalized fractional integral of order p provided the right side is point-wise defined on (a, ∞); Γ (·) denotes the gamma function.
Likewise, the generalized fractional derivative is defined in terms of the generalized fractional integral (2.1), Definition 2.2 Suppose n -1 < p ≤ n, 0 ≤ q ≤ 1 with n ∈ N and ρ > 1. Then the fractional operator is called the left-sided Hilfer generalized (or simply Hilfer-Katugampola) fractional derivative of order p and parameter q.
Next, we consider some important properties of these operators as follows.
, the following relations hold:

Main results
In this section, firstly, we establish the relationship between the proposed problem (1.6) and the Volterra integral equation. Secondly, with the help of these equation, we derive the existence and uniqueness of solutions based on Banach's and Krasnoselskii's fixed point theorems.

Equivalence with the Volterra integral equation
We consider the generalized boundary value problem described by where p ∈ (1, 2], q ∈ [0, 1], ρ > 0, subject to the boundary conditions: 2) The following lemma shows the relationship between problem (3.1)-(3.2) and the Volterra integral equation of the form .

Then x is a solution of the problems (3.1), (3.2) if and only if it satisfies Eq. (3.3).
Proof Indeed, in view of (2.4), Eq. (3.1) can be written as Applying the operator ρ I p a + to both sides of Eq. (3.5) and using Lemma 2.6, we get Denoting e 1 = ( ρ I 2-r a + x)(a) and e 2 = ( ρ I 1-r a + x)(a) in Eq. (3.6) yields By substituting the boundary condition x(a) = 0 in Eq. (3.7), we get e 1 = 0. Thus Substituting t = δ i and multiplying through by μ i in Eq. (3.8), we have Applying ρ I ξ i a + to both sides of Eq. (3.9) and using Lemmas 2.5 and 2.4 yield which implies that Also, From the boundary condition a + x(δ i ) + and in view of Eqs. (3.11), (3.12), we obtain which implies that (3.14) Therefore, the results follow by substituting Eq. (3.14) in Eq. (3.8). The converse follows by simple calculation as above.
Denote by C = C([a, b], R) the Banach space of all continuous functions from [a, b] to R defined by and define the operator F : We show that the operator F has a fixed point which is the solutions of the problem (1.6). For simplicity, we denote

Existence result via Banach contraction principle
We prove first the existence and uniqueness result of the proposed problem (1.6) by applying the concept of the Banach contraction principle in this subsection. Thus, before we state the main theorem, we make the following assumptions: (T 1 ) Suppose that f : [a, b] × R → R is a continuous functions such that f (·, x(·)) ∈ C for any x ∈ C. (T 2 ) There exists a constant A > 0 such that where Φ is defined by (3.16).
Then we show that FB χ ⊂ B χ . Now, for any x ∈ B χ , we have This shows that (F)(x) ≤ χ and hence FB χ ⊂ B χ . Finally, we show that the operator F defined by (3.15) is contractive. Indeed, for any x, y ∈ B χ , we have Therefore, Fx -Fy ≤ AΦ xy . Thus in view of (3.17), the operator F is contractive and hence by the Banach contraction principle, problem (1.6) has a unique solution on [a, b]. Remark 3. 3 We now present some special cases of our Theorem 3.2. For instance, to obtain the result of [12] (see Theorem 3.2), we set the parameters ρ = 1 and = 0, then we obtain (3.20) where We now state the following corollary.
Our next existence results are established using the concepts of Krasnoselskii's fixed point theorem.

Existence result via Krasnoselskii's fixed point theorem
In this subsection, we establish the existence results via Kransnoselkii's fixed point theorem [30].

Theorem 3.5 ([30])
Let P be a closed, convex, and nonempty subset of a Banach space X, and let F 1 , F 2 be operators such that: (i) F 1 y + F 2 z ∈ M whenever y, z ∈ P.
(ii) F 1 is compact and continuous.
(iii) F 2 is a contraction mapping. Then there exists w ∈ P such that w = F 1 w + F 2 w. and Step 1. We show that Fx + Gy ∈ B d . Indeed, for any x, y ∈ B d and t ∈ [a, b], we get Step 2. We show that the operator G is a contraction. Let x, y ∈ B d and t ∈ [a, b]. We have Thus, in view of condition (3.23), we conclude that the operator G is a contraction.
Step 3. We show that the operator F is compact and continuous.
Since the function f ∈ C([a, b], R), the operator F is also continuous. Moreover, it is not difficult to see that the operator F is uniformly bounded. For any x ∈ B d , we have In order to show the compactness of the operator F , we set → ∞, as t 2 → t 1 . (3.26) This shows that the operator F is equicontinuous and, therefore, is relatively compact on B d . Thus, as a result of the Arzelá-Ascoli theorem, it follows that F is compact on Remark 3.7 We now present a special case of Theorem 3.2. For instance, to obtained the result of [12] (see Theorem 3.17), we set the parameters ρ = 1 and = 0, then we obtain where We now state the following corollary.
Example 5.2 Consider the generalized nonlocal boundary value problem which involves the Hilfer-Katugampola fractional derivative described by

Conclusions
In this article, we have shown that the proposed problem (1.6) is equivalent to the Volterra integral equation. Based on this equation, we have provided some sufficient conditions for the existence and uniqueness of solutions by utilizing the concepts of the Banach and Krasnoselskii's fixed point theorems. Also, we derive stability in the case of Ulam-Hyers and generalized Ulam-Hyers, respectively. Furthermore, two examples are given as an application to validate the results obtained. Also, it enables us to make the accompanying comments: • If ρ → 1, = 0 problem (1.6) reduces to the problem considered in [12].
• If ξ i = 1 for all i = 1, . . . , m, the generalized nonlocal Riemann-Liouville integral condition takes the form which plays a significant role in computational fluid dynamics, ill-posed problems and mathematical models and yields a better result than the initial condition x(b) = x b [14,34]. Therefore, based on this setting, we conclude that our results are new and extend the above cited results which can be considered as a further development of the qualitative analysis of fractional differential equations. The analysis of the Ulam-Hyers-Mittag-Leffler function is to be discussed shortly in the context of Hilfer-Katugampola fractional derivative with generalized nonlocal conditions proposed in this paper or another mixed boundary condition.