Computational Methods for Differential Equations http://cmde.tabrizu.ac.ir Vol. 5, No. 2, 217, pp. 158-169 Existence of triple positive solutions for boundary value problem of nonlinear fractional differential equations Kamal Shah Department of Mathematics, University of Malakand, Chakadara Dir(L), Khyber Pakhtunkhwa, Pakistan. E-mail: kamalshah48@gmail.com Salman Zeb Department of Mathematics, University of Malakand, Chakadara Dir(L), Khyber Pakhtunkhwa, Pakistan. E-mail: salmanzeb@uom.edu.pk Rahmat Ali Khan Department of Mathematics, University of Malakand, Chakadara Dir(L), Khyber Pakhtunkhwa, Pakistan. E-mail: rahmat alipk@yahoo.com Abstract This article is devoted to the study of existence and multiplicity of positive solutions to a class of nonlinear fractional order multi-point boundary value problems of the type D q + u(t) = f(t, u(t)), 1 < q 2, < t < 1, u() =, u(1) = δ i u(η i ), where D q + represents standard Riemann-Liouville fractional derivative, δ i, η i (, 1) with δ i η q 1 i < 1, and f : [, 1 [, ) [, ) is a continuous function. We use some classical results of fixed point theory to obtain sufficient conditions for the existence and multiplicity results of positive solutions to the problem under consideration. In order to show the applicability of our results, we provide some examples. Keywords. Fractional differential equations; Boundary value problems; Positive solutions; Green s function; Fixed point theorem. 21 Mathematics Subject Classification. 34A8 and 26A33. 1. Introduction This article is concerned with the existence triple positive solutions to multi-point boundary value problems with nonlinear fractional order differential equations of the Received: 14 January 217 ; Accepted: 29 April 217. Corresponding author. 158
CMDE Vol. 5, No. 2, 217, pp. 158-169 159 form D q + u(t) = f(t, u(t)), 1 < q 2, < t < 1 u() =, u(1) = δ i u(η i ), (1.1) where D q + is the standard Riemann-Liouville fractional derivative of order q and δ i, η i (, 1) with δ i η q 1 i < 1, and f : [, 1 [, ) [, ) is continuous. Differential equations of fractional order is one of the fast growing area of research in the field of mathematics and have recently been proved to be valuable tools in the modeling of many phenomena in various fields of science and engineering. Indeed, one can find numerous applications of fractional order differential equations in viscoelasticity, electro-chemistry, control theory, movement through porous media, electromagnetics, etc (see [2, 6, 1, 11, 18, 19, 21, 22). The theory of existence of solutions to boundary value problems associated with fractional differential equations have recently been attracted the attention of many researchers, see for example [1, 3, 4, 5, 8, 12, 13, 16, 23 and the references therein. In these cited references, existence of at least one solution is studied with the tools of classical fixed point theory. In [23, Rehman and Khan investigated multi-point boundary value problem for fractional order differential equation D α t y(t) = f(t, y(t), D β t y(t)); t (, 1), y() =, D β t y(1) ζ i D β t y(ξ i ) = y, where 1 < α 2, < β < 1, < ξ i < 1, ζ i [, + ) with ζ i ξ α β 1 < 1 and obtained sufficient conditions for existence and uniqueness of nontrivial solutions via Schauder fixed point theorem, and contraction mapping principle. As for the existence of multiple positive solutions are concerned, few papers can be found in the literature dealing with the existence and multiplicity of positive solutions to multi-point boundary value problems for fractional differential equations [7, 24, 29. Zhang et. al [29 studied existence of positive solutions for the eigenvalue problem corresponding to a class of fractional differential equation D α t x(t) = λf(t, x(t), D β t x(t)); t (, 1), p 2 D β t x() =, D γ t x(1) = a j D γ t x(ξ j ), where λ is a parameter, 1 < α 2, α β > 1, < β γ < 1, < ξ 1 < ξ 2 <... < ξ p 2 < 1, a j [, + ) with c = p 2 a j ξ α γ 1 j < 1, and D t is the standard Riemann-Liouville derivative.
16 K. SHAH, S. ZEB, AND R. A. KHAN Recently existence of positive solutions for single and system of boundary value problems of nonlinear fractional order differential equations corresponding to different boundary conditions are also studied by many authors for which we refer some recent works in [9, 17, 2, 26, 27, 28. Inspired from the above works, in this paper, we study a different problem and obtain sufficient conditions for existence, uniqueness as well as conditions for existence of triple positive solutions to the BVP (1.1). For the applicability of our results, we include some examples. 2. Preliminaries This section contains some necessary definitions and lemmas taken from fractional calculus [18, 21 and functional analysis books. These definitions and lemmas will be used frequently in the forthcoming section. Definition 2.1. The fractional integral of order q > of a function y : (, ) R is given by I q + y(t) = 1 t (t s) q 1 y(s) ds, provided that the integral converges. Definition 2.2. The fractional derivative of order q > of a continuous function y : (, ) R is given by ( ) n D q + y(t) = 1 d t (t s) n q 1 y(s) ds, Γ(n q) dt where n = [q + 1, provided that the right side is pointwise defined on (, ). Definition 2.3. A map θ is said to be a nonnegative continuous concave functional on a cone P of a real Banach space E provided that θ : P [, ) is continuous and θ(tx + (1 t)y) tθ(x) + (1 t)θ(y), for all x, y P and t 1. The next two lemmas provide an important rule for obtaining the equivalent integral equation of BVP (1.1). Lemma 2.4. [25 If we assume u C(, 1) L(, 1), then the fractional differential equation of order q > D q + u(t) =, has a unique solution of the form u(t) = C 1 + C 2 t q 2 +... + C N t q N, C i R, i = 1, 2,..., N. The following law of composition can be easily deduced from Lemma 2.4. Lemma 2.5. Assume that u C(, 1) L(, 1), with a fractional derivative of order q that belongs to C(, 1) L(, 1), then I q + Dq + u(t) = u(t) + C 1 + C 2 t q 2 +... + C N t q N, C i R, i = 1, 2,..., N.
CMDE Vol. 5, No. 2, 217, pp. 158-169 161 Lemma 2.6. [14 Let E be a Banach space, P E a cone, and Ω 1, Ω 2 be the two bounded open balls of E centered at the origin with Ω 1 Ω 2. Suppose that T : P (Ω 2 \ Ω 1 ) P is a completely continuous operator such that either (i) T u u, u P Ω 1 and T u u, u P Ω 2, or (ii) T u u, u P Ω 1 and T u u, u P Ω 2 holds. Then T has a fixed point in P (Ω 2 \ Ω 1 ). Lemma 2.7. [15 Let P be a cone in a real Banach space E, P c = {u P : u c}, θ a nonnegative continuous concave function on P such that θ(u) u for all u P c, and P (θ, b, d) = {u P : b θ(u), u d}. Suppose T : P c P c is completely continuous and there exist constants < a < b < d c such that (i) {u P (θ, b, d) θ(u) > b} =, and θ(t u) > b for u P (θ, b, d) (ii) T u < a for u a (iii) θ(t u) > b for u P (θ, b, c) with T u > d. Then T has at least three fixed points u 1, u 2, u 3 with u 1 < a, b < θ(u 2 ), a < u 3 with θ(u 3 ) < b. 3. Main Results In this section we develop sufficient conditions on the nonlinear function f, under which the BVP (1.1) has at least one solution and also we study sufficient conditions leading to multiplicity of positive solutions. Denote by E = C[, 1 the Banach space of all continuous real-valued functions on [,1 with norm u = cone by P such that P = {u E : u(t), t [, 1}. sup t 1 Define nonnegative continuous concave functional θ on the cone P as follow θ(u) = u(t) and a min u(t). (3.1) η i 2 t η i Lemma 3.1. For y(t) C[, 1, the linear BVP D q + u(t) + y(t) = ; < t < 1, 1 < q 2, u() =, u(1) = δ i u(η i ), (3.2) has a unique solution of the form u(t) = G(t, s)y(s) ds, where the Green function G(t, s) is given by [ (1 s) q 1 δ 1 λ j(η j s) q 1 (t s) q 1 ; s t, < s η i, G(t, s) = 1 [ i = 1, 2,..., m 1, (1 s) q 1 δ 1 λ j (η j s) q 1 ; t s, < s η i, i = 1, 2,..., m 1. (3.3)
162 K. SHAH, S. ZEB, AND R. A. KHAN Proof. In view of Lemma (2.5), we obtain u(t) = I q + y(t) + C 1 + C 2 t q 2, (3.4) for some C 1, C 2 R. The boundary condition u() [ = implies C 2 = and the condition u(1) = δ i u(η i ), yields C 1 = 1 1 λ I q + y(1) δ i I q + y(η i), where λ = δ i η q 1 i < 1. Hence, (3.4) takes the form u(t) = I q tq 1 + y(t) + 1 λ [ I q + y(1) δ i I q + y(η i). (3.5) We discuss several cases: for t η 1, we write (3.5) as u(t) = t [ (t s) q 1 ( + (1 s) q 1 δ j(η j s) q 1) y(s) ds + + i=2 ηi η1 t ( (1 s) q 1 ( (1 s) q 1 δ j (η j s) q 1) y(s) ds δ j (η j s) q 1) y(s) ds + For η l 1 t η l, 2 l m 2, we write (3.5) as u(t) = η 1 + i=2 [ (t s) ( q 1 + tq 1 (1 λ) (1 s) q 1 ηi [ (t s) q 1 + tq 1 (1 λ) η (1 s) q 1 y(s) ds. δ j (η j s) q 1) y(s) ds ( (1 s) q 1 δ j (η j s) q 1) + (1 s) q 1 y(s) ds + t [ η l 1 (t s) ( q 1 + tq 1 (1 λ) (1 s) q 1 δ j (η j s) q 1) y(s) ds [ + tq 1 ηl (1 λ) t + i=l+1 ηi ( (1 s) q 1 j=l j=l δ j (η j s) q 1) y(s) ds ( (1 s) q 1 δ j (η j s) q 1) y(s) ds + η (1 s) q 1 y(s) ds. For η t 1, we write (3.5) as u(t) = η1 + [ (t s) q 1 ( + (1 s) q 1 δ j (η j s) q 1) y(s) ds ηi i=2 [ (t s) q 1 + ( (1 s) q 1 δ j (η j s) q 1) y(s) ds [ t + (1 λ)(t s) q 1 + (1 s) q 1 y(s) ds + η t (1 s) q 1 y(s) ds.
CMDE Vol. 5, No. 2, 217, pp. 158-169 163 Therefore, the unique solution of the BVP (1.1) is given by u(t) = G(t, s)y(s) ds, where [ 1 λ (1 s) q 1 δ j (η j s) q 1 (t s) q 1 ; s t, < s η i, G(t, s) = 1 [ i = 1, 2,..., m 1, 1 λ (1 s) q 1 δ j (η j s) q 1 ; t s, < s η i, i = 1, 2,..., m 1. Lemma 3.2. The Green s function defined by (3.3) satisfies the following conditions: (i) G(t, s) >, for t, s (, 1); (ii) There exists a positive function γ C(, 1) such that min G(t, s) γ(s) max G(t, s) = γ(s)g(s, s), for < s < 1. (3.6) t η i t 1 Proof. (i): If η j 1 < s η j, For t < η j and for t η j G(t, s) = tq 1 [(1 s) q 1 δ j (η j s) q 1 (1 λ)(t s) q 1. G(t, s) > tq 1 [(1 s) q 1 δ j (η j s) q 1 >, G(t, s) tq 1 (1 s) q 1 (t s) q 1 For η j s t 1, >. G(t, s) = tq 1 (1 s) q 1 (1 λ)(t s) q 1 tq 1 (1 s) q 1 (t s) q 1 Therefore for each η j s t 1 and t s, s η j. Hence the expression for G(t, s) in (3.3) shows that G(t, s) >, for s, t (, 1). (ii) Let us denote g 1 (t, s) = g 2 (t, s) = ( (1 s) q 1 δ j (η j s) q 1) ( (1 s) q 1 δ j (η j s) q 1), (t s)q 1, then, g 1(t,s) t < and g 2(t,s) t > for s t, which implies that g 1 (t, s) is decreasing and g 2 (t, s) is an increasing function for s t. It follows that G(t, s) is decreasing >.
164 K. SHAH, S. ZEB, AND R. A. KHAN with respect to t for s t and increasing with respect to t for t s. Consequently, min G(t, s) = t η i which implies that min G(t, s) = t η i g 1 (η i, s), s (,, min{g 1 (η i, s), g 2 (, s)}, s [, η i, g 2 (, s), s [η i, 1), { g1 (η i, s), s (, r, g 2 (, s), s [r, 1), where < r < η i. Further, we note that s q 1 [ max G(t, s) = G(s, s) = (1 s) q 1 δ i (η i s) q 1, t 1 for s (, 1). As a result ( ηi ) q 1 s [ (η i s) q 1, s (, r, s γ(s) = q 1 (1 s) (1 λ) q 1 δ j (η j s) q 1 ) q 1, s [r, 1). ( ηi 1 s In view of Lemma (3.1), the BVP (1.1) is equivalent to the integral equation u(t) = G(t, s)f(s, u(s))ds, (3.7) and by a solution of the BVP (1.1), we mean a solution of the integral equation (3.7), that is a fixed point of the operator T : P P defined by T u(t) = G(t, s)f(s, u(s))ds. (3.8) Lemma 3.3. For nonnegative real-valued functions m, n L[, 1 such that f(t, u) m(t) + n(t)u, for almost every t [, 1, and all u [, ), (3.9) the operator T defined by (3.8) is completely continuous. Proof. Due to nonnegativity and continuity of G(t, s) and f(t, s), the operator T is continuous. For each u Ω = {u P : u M, M > }, we have T u(t) = G(t, s)f(s, u(s))ds G(t, s)(m(s) + n(s)u(s))ds G(s, s)m(s)ds + M which implies that T (Ω) is bounded. G(s, s)n(s)ds = l, For equicontinuity of T : P P, take t 1, t 2 [, 1 such that t 1 < t 2 with t 2 t 1 < δ. Then for u Ω, we have T u(t 2 ) T u(t 1 ) = [G(t 2, s)f(s, u(s)) G(t 1, s)f(s, u(s)) ds
CMDE Vol. 5, No. 2, 217, pp. 158-169 165 G(t 2, s) G(t 1, s) (m(s) + n(s)m)ds, which in view of the continuity of G(t, s) and K = max{m(s) + n(s)m : s [, 1} implies that 1 1 [ T u(t 2 ) T u(t 1 ) tq 1 2 1 1 (1 s) q 1 δ i (η i s) q 1 Kds 1 λ 1 tq 1 2 K 1 q(1 λ). Using mean value theorem on 2 1 and the choice t 2 t 1 < δ, we have K T u(t 2 ) T u(t 1 ) q(1 λ) (q K(q 1) 1)Cq 2 t 2 t 1 q(1 λ) δq 1 ϵ where C q 2 δ q 2, δ = ( ϵ Aq ) 1 q 1 and A = K q(1 λ). Hence T : P P is equicontinuous. By Arzela-Ascoli theorem, we conclude that the operator T : P P is completely continuous. Now, we show existence of at least one solution of the BVP (1.1). Fix M = ( ) 1 1 ( G(s, s)ds ηi 1 and N = γ(s)g(s, s)ds). Theorem 3.4. Assume that the condition (3.9) of Lemma (3.3) holds and there exist two positive constants r 2 > r 1 > such that f(t, u) Mr 2, for (t, u) [, 1 [, r 2 and f(t, u) Nr 1, for (t, u) [, 1 [, r 1. Then BVP (1.1) has at least one positive solution u such that r 1 u r 2. (3.1) Proof. In view of (3.8) and Lemma (3.9), T is completely continuous. By Schauder fixed point theorem, T has a fixed point, T u = u. Define Ω 2 = {u P : u < r 2 }. For u Ω 2, we have u(t) r 2 for all t [, 1 and using (3.1) and Lemma (3.2), it follows that T u = max t 1 G(t, s)f(s, u(s))ds Mr 2 G(s, s)ds = r 2 = u. Similarly, define Ω 1 = {u P : u < r 1 }. For u Ω 1, we have u(t) r 1 for all t [, 1 and for t [η i i, η i, we have T u = G(t, s)f(s, u(s))ds Nr 1 ηi γ(s)g(s, s)f(s, u(s))ds γ(s)g(s, s)ds = r 1 = u.
166 K. SHAH, S. ZEB, AND R. A. KHAN Example 1. For the BVP ( M = D 3 2 + u(t) + u 2 + cos t 4 u() =, u(1) = + 1 =, < t < 1, 5 ) 1 ( G(s, s)ds 3 1.88 and N = 4 1 4, we have r 1 = 1 25, r 2 = 3 4 δ i u(η i ) = 1 5, (3.11) γ(s)g(s, s) ds) 1 = 4.914. Choosing f(t, u) = u 2 + cos t + 1 4 5 1.151 Mr 2, for (t, u) [, 1 [, 4 3, f(t, u) = u 2 + cos t + 1 4 5 1 5 Nr 1, for (t, u) [, 1 [, 25 1. By Theorem 3.2, the BVP (3.11) has at least one solution u such that 1 25 u 3 4. Theorem 3.5. Let there exists h(t) L[, 1 such that f(t, u) f(t, v) h(t) u v, 1 G(s, s)h(s)ds < 1 and for almost t L[, 1, and u, v [, ), then the BVP (1.1) has a unique positive solution. Proof. For u, v P, we have T u(t) T v(t) = where α = G(t, s)(f(s, u(s) f(s, v(s))ds G(t, s) (f(s, u(s)) f(s, v(s)) ds, G(s, s)h(s) u(s) v(s) ds = α u v, 1 G(s, s)h(s)ds < 1. BVP has a unique solution. Example 2. Consider the BVP D 3 2 + u(t) + u() =, u(1) = G(s, s)h(s)ds u v Therefore, by Banach contraction principle, The e 2t u 4(1 + e 2t )(1 + u) + cos2 t + 1 =, < t < 1, δ i u(η i ) = 1 5. (3.12)
CMDE Vol. 5, No. 2, 217, pp. 158-169 167 e 2t u Here, f(t, u) = 4(1+e 2t )(1+u) + cos2 t + 1 and taking h(t) = e2t 4(1+e 2t ), then G(s, s)h(s)ds 1 s 5 + 2 s(1 s) e 2s 4(1+e 2s ) ds and π s 5 + 2 s(1 s) π ds =.576446 f(t, u) f(t, v) h(t) u v, for (t, u), (t, v) [, 1 [, ). By Theorem (3.5), the BVP has a unique positive solution. Theorem 3.6. Assume that the condition (3.9) of Lemma (3.3) holds and there exist positive constants < a < b < c such that (i) f(t, u) < Ma, for (t, u) [, 1 [, a (ii) f(t, u) Nb, for (t, u) [, η i [b, c (iii) f(t, u) Mc, for (t, u) [, 1 [, c, then the BVP (1.1) has at least three positive solutions u 1, u 2, and u 3 such that max u 1(t) < a, b < t 1 min u 2 (t) c, a < max u 3(t) < c, t η i t 1 Proof. For u P c, the relation T u = max G(t, s)f(s, u(s))ds G(s, s)f(s, u)s))ds t 1 min u 3 (t) < b. t η i (3.13) G(s, s)mcds c follows from item (iii) and the completely continuity of T : P c P c follows from Lemma (3.3). Choose u(t) = b+c b+c 2, t 1. Then using (3.1), we have u(t) = 2 P (θ, b, c), θ(u) = θ( (b+c) 2 ) > b implies that {u P (θ, b, c) θ(u) > b} =. Hence, if u P (θ, b, c), then b u(t) c for t η i. Also, from assumption (ii), we have f(t, u(t)) Nb, for t η i and θ(t u) = which implies that min (T (u)) t η i θ(t u) > b, for all u P (θ, b, c). γ(s)g(s, s)f(s, u(s))ds > ηi γ(s)g(s, s)nbds = b, Hence by Lemma (2.7), the BVP (1.1) has at least three positive solutions u 1, u 2,and u 3 satisfying (3.13). Example 3. For the problem where D 3 2 + u(t) + f(t, u) =, < t < 1, u() =, u(1) = f(t, u) = { cos t 1 + u2 ; u 1, 3 + cos t 1 + u; u > 1, δ i u(η i ) = 1 5, (3.14)
168 K. SHAH, S. ZEB, AND R. A. KHAN we find that M 1.88 and N = 4.914. Choosing a = 1 1, b = 3 4 and c = 44, we have f(t, u) = cos t 1 + u2.2 < Ma.188, for (t, u) [, 1 [, 1, 1 f(t, u) = 3 + cos t 1 + u 47.1 Nb 3.6855, for (t, u) [ 1 4, 3 [3 4 4, 44, f(t, u) = 3 + cos t + u 47.1 Mc 47.872, for (t, u) [, 1 [, 44. 1 Hence, by Theorem (3.6), the BVP (3.14) has at least three positive solutions u 1, u 2, and u 3 with < max t 1 u 1 (t) < 1 1, 1 < min 1 4 t 3 4 u 2(t) 44, 1 1 < max t 1 u 3 (t) < 44, min 1 4 t 3 4 u 3(t) < 3 4. Acknowledgements We are really thankful to the anonymous referees for their useful comments and suggestions which improved the quality of the paper. References [1 B. Ahmad and J. J. Nieto, Existence results for a coupled system of nonlinear fractional differential equations with three-point boundary conditions, Comput. Math. Appl., 58 (29), 1838 1843. [2 A. Bakakhani and V. D. Gejji, Existence of positive solutions of nonlinear fractional differential equations, J. Math. Anal. Appl., 278 (23), 434 442. [3 Z. Bai, On positive solutions of a nonlocal fractional boudary value problem, Nonlinear Anal., 72 (21), 916 924. [4 Z. Bai and H. Lu, Positive solutions for boundary value problem of nonlinear fractional differential equation, J. Math. Anal. Appl., 311 (25), 495 55. [5 M. Benchohra, S. Hamani, and S. K. Ntouyas, Boundary value problems for differential equations with fractional order and nonlocal conditions, Nonlinear Anal., 71 (29), 2391 2396. [6 D. Delbosco, Fractional calculus and function spaces, J. Fract. Calc., 6 (1996), 45 53. [7 M. El-Shahed and J. J. Nieto, Nontrivial solutions for a nonlinear multi-point boundary value problem of fractional order, Comput. Math. Appl., 59(11) (21), 3438 3443. [8 C. S. Goodrich, Existence of a positive solution to a class of fractional differential equations, Appl. Math. Lett., 23 (21), 15 155. [9 Z. Hu, W. Liu, and J. Liu, Boundary value problems for fractional differential equations, Boundary Value Problems., 214 (214), 12 pages. [1 A. A. Kilbas, H.M. Srivastava, and J. J. Trujillo, Theory and Applications of Fractional Differential Equations, North-Holland athematics Studies, vol. 24, Elsevier, Amsterdam, 26. [11 A. A. Kilbas, O.I. Marichev, and S. G. Samko, Fractional Integral and Derivatives (Theory and Applications), Gordon and Breach, Switzerland, 1993. [12 R. A. Khan, Existence and approximation of solutions to three-point boundary value problems for fractional differential equations, Elec. J. Qual. Theory Differ. Equ., 58 (211), 8 pages. [13 R. A. Khan and M. Rehman, Existence of multiple positive solutions for a general system of fractional differential equations, Commun. Appl. Nonlinear Anal., 18 (211), 25 35. [14 M. A. Krasnoselskii, Positive solutions of operator equations, Noordhoff, Groningen, 1964. [15 R. W. Leggett and L.R Williams, Multiple positive fixed points of nonlinear operators on ordered Banach spaces, Indiana Univ. Math. J., 28 (1979), 673 688.
CMDE Vol. 5, No. 2, 217, pp. 158-169 169 [16 C. F. Li, X. N. Luo, and Y. Zhou, Existence of positive solutions of the boundary value problem for nonlinear fractional differential equations, Comput. Math. Appl., 59 (21), 1363 1375. [17 N. Li and C. Y. Wang, New existence results of positive solution for a class of nonlinear fractional differential equations, Acta. Math. Scientia., 33 (213), 847 854. [18 K. S. Miller and B. Ross, An Introduction to the Fractional Calculus and Fractional Differential Equations, Wiley, New York, 1993. [19 K. S. Miller, Fractional differential equations, J. Fract. Calc., 3 (1993), 49 57. [2 A. Nanware and D. B. Dhaigude, Existence and uniqueness of solutions of differential equations of fractional order with integral boundary conditions, J. Nonlinear Sci. Appl., 7 (214), 246 254. [21 I. Podlubny, Fractional Differential Equations, Mathematics in Science and Engineering, Academic Press, New York, 1999. [22 I. Podlubny, The Laplace transform method for linear differential equations of the fractional order, Inst. Expe. Phys., Slov. cad. Sci., UEF-2-94, Kosice, 1994. [23 M. Rehman and R. A. Khan, Existence and uniqueness of solutions for multi-point boundary value problems for fractional differential equations, Appl. Math. Lett., 23 (21), 138 144. [24 A. H. Salem, On the fractional m-point boundary value problem in reflexive Banach space and the weak topologies, J. Comput. Appl. Math., 224 (29), 565 572. [25 J. Sabatier, O. P. Agarwal, and J. A. Ttenreiro Machado, Advances in Fractional Calculus, Theoretical Developments and Applications in Physics and Engineering, Springer, 27. [26 K. Shah, H. Khalil, and R. A. Khan, Investigation of positive solution to a coupled system of impulsive boundary value problems for nonlinear fractional order differential equations, Chaos. Solit. Fract., 77(215), 24 246. [27 K. Shah and R. A. Khan, Existence and uniquness of positive solutions to a coupled system of nonlinear fractional order differential equations with anti-periodic boundary conditions, Diff. Equ. Appl., 7(2) (215), 245-262. [28 X. Wang, L. Wang, and Q. Zeng, Fractional differential equations with integral boundary conditions, J. Nonlinear Sci. Appl., 8 (215), 39 314. [29 X. Zhang, L. Liu, B. Wiwatanapataphee, and Y.H. Wu, Positive solutions of eigenvalue problems for a class of fractional differential equations with derivatives, Abstr. Appl. Anal., 212 (212), 14 pages. [3 W. Zhong and W. Lin, Nonlocal and multiple-point boundary value problem for fractional differential equations, Comput. Math. Appl., 59(3) (21), 1345 1351.