Ji et al. Boundary Value Problems (2015) 2015:214 DOI 10.1186/s13661-015-0485-3
RESEARCH
Open Access
Positive solutions for higher order differential equations with integral boundary conditions Yude Ji1* , Yanping Guo2 and Yukun Yao3 *
Correspondence:
[email protected] 1 College of Sciences, Hebei University of Science and Technology, Shijiazhuang, Hebei 050018, P.R. China Full list of author information is available at the end of the article
Abstract In this paper, we consider the existence of at least three positive solutions for the 2nth order differential equations with integral boundary conditions
x(2n) (t) = f (t, x(t), x (t), . . . , x(2(n–1)) (t)), 0 ≤ t ≤ 1, 1 x(2i) (0) = 0 ki (s)x(2i) (s) ds, x(2i) (1) = 0, 0 ≤ i ≤ n – 1,
where (–1)n f > 0 is continuous, and ki (t) ∈ L1 [0, 1] (i = 0, 1, . . . , n – 1) are nonnegative. The associated Green’s function for the higher order differential equations with integral boundary conditions is first given, and growth conditions are imposed on f which yield the existence of multiple positive solutions by using the Leggett-Williams fixed point theorem. Keywords: boundary value problem; integral boundary conditions; positive solution; fixed point theorem
1 Introduction The multi-point boundary value problems (BVPs) for ordinary differential equations arise in a variety of different areas of applied mathematics and physics. The study of nonlocal BVPs for second order ordinary differential equations has been widely investigated in [–]. Since then, nonlinear high order nonlocal BVPs have been studied by many authors. We refer the reader to [–] and references therein. Recently, Guo et al. [] used LeggettWilliams fixed point theorem to obtain the existence of at least three positive solutions for the nth order m-point BVP
y(n) (t) = f (t, y(t), y (t), . . . , y((n–)) (t)), ≤ t ≤ , (i) y(i) () = m– (ξj ), ≤ i ≤ n – , y(i) () = , j= kij y
where kij > (i = , , . . . , n – ; j = , , . . . , m – ), < ξ < ξ < · · · < ξm– < , and (–)n f : [, ] × Rn → [, +∞) is continuous. BVPs with integral boundary conditions for ordinary differential equations represent a very interesting and important class of problems and arise in the study of various physical, biological and chemical processes, such as heat conduction, chemical engineering, © 2015 Ji et al. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.
Ji et al. Boundary Value Problems (2015) 2015:214
Page 2 of 11
thermo-elasticity, underground water flow, population dynamics, and plasma physics. Such problems include two-, three-, multi-point and nonlocal BVPs as special cases. The existence and multiplicity of positive solutions for such problems have received a great deal of attention; see [–] and references therein. In particular, we would mention the result of [], Zhang and Ge investigated the existence and nonexistence of positive solutions of the following fourth-order BVP with integral boundary conditions ⎧ () ⎪ ⎨ x (t) = ω(t)f (t, x(t), x (t)), < t < , x() = , x() = g(s)x(s) ds, ⎪ ⎩ x () = , x () = h(s)x (s) ds, where ω may be singular at t = and (or) t = , f ∈ C([, ] × [, +∞) × (–∞, ], [, +∞)), and g, h ∈ L [, ] are nonnegative. Motivated by [, ], in this paper, we consider the existence of at least three positive solutions for the nth order differential equations with integral boundary conditions
x(n) (t) = f (t, x(t), x (t), . . . , x((n–)) (t)), ≤ t ≤ , x(i) () = , ≤ i ≤ n – , x(i) () = ki (s)x(i) (s) ds,
(.)
where (–)n f > is continuous, and ki (t) ∈ L [, ] (i = , , . . . , n – ) are nonnegative. For more precise conditions on f , let (–)j [a, b] = [a, b] if j is even and (–)j [a, b] = [–b, –a] if j is odd. Let n–
[aj , bj ] = [a , b ] × · · · × [an– , bn– ]. j=
We shall require that (–)n f : [, ] ×
n–
(–)j [, +∞) → [, +∞).
j=
We shall suppose the following conditions are satisfied: (H ) ki (t) ∈ L [, ] are nonnegative, and Ki ∈ [, ), where Ki =
( – s)ki (s) ds,
≤ i ≤ n – ;
(H ) (–)n f : [, ] ×
n–
j= (–)
j
[, +∞) → [, +∞) is continuous.
2 Preliminary results Definition . Let E be a Banach space over R. A nonempty convex closed set K ⊂ E is said to be a cone provided that (i) au ∈ K for all u ∈ K and all a ≥ ; (ii) u, –u ∈ K implies u = . Definition . The map α is said to be a nonnegative continuous concave functional on K provided that α : K → [, ∞) is continuous and
α tx + ( – t)y ≥ tα(x) + ( – t)α(y)
Ji et al. Boundary Value Problems (2015) 2015:214
Page 3 of 11
for all x, y ∈ K and ≤ t ≤ . Similarly, we say the map γ is a nonnegative continuous convex functional on K provided that γ : K → [, ∞) is continuous and
γ tx + ( – t)y ≤ tγ (x) + ( – t)γ (y) for all x, y ∈ K and ≤ t ≤ . Definition . Let < a < b be given and let α be a nonnegative continuous concave functional on K . Define the convex sets Pr and P(α, a, b) by and Pr = x ∈ K| x < r
P(α, a, b) = x ∈ K|a ≤ α(x), x ≤ b .
Theorem . (Leggett-Williams fixed point theorem []) Let A : Pc → Pc be a completely continuous operator and let α be a nonnegative continuous concave functional on K such that α(x) ≤ x for all x ∈ Pc . Suppose there exist < a < b < d ≤ c such that (C ) {x ∈ P(α, b, d)|α(x) > b} = ∅, and α(Ax) > b for x ∈ P(α, b, d), (C ) Ax < a for x ≤ a, and (C ) α(Ax) > b for x ∈ P(α, b, c), with Ax > d. Then A has at least three fixed points x , x , and x such that x < a,
b < α(x ) and x > a
with α(x ) < b.
Remark . If we have d = c, then condition (C ) of Theorem . implies condition (C ) of Theorem ..
3 Preliminary lemmas Lemma . Suppose (H ) holds. Then gi (t, s) ≤ ( ≤ i ≤ n – ), where gi (t, s) is the Green’s function for the problem
x (t) = , ≤ t ≤ , x() = ki (s)x(s) ds,
x() = .
Proof It is easy to see that gi (t, s) ≤ ( ≤ i ≤ n – ) by using Lemma . of [].
Let G (t, s) = gn– (t, s), then for ≤ j ≤ n – , we recursively define Gj (t, s) =
gn–j– (t, τ )Gj– (τ , s) dτ .
Lemma . Suppose (H ) holds. If y ∈ C[, ], then the BVP
x(l) (t) = y(t), ≤ t ≤ , x(i) () = kn–l+i– (s)x(i) (s) ds,
x(i) () = ,
≤ i ≤ l – ,
(.)
has a unique solution for each ≤ l ≤ n – , where Gl (t, s) is the associated Green’s function for the BVP (.).
Ji et al. Boundary Value Problems (2015) 2015:214
Page 4 of 11
Proof We will prove the result by using mathematical induction. When l = , which implies that i = , then the BVP (.) reduces to
x (t) = y(t), ≤ t ≤ , x() = kn– (s)x(s) ds,
(.)
x() = .
By using Lemma ., it is easy to see that the BVP (.) has a unique solution
G (t, s)y(s) ds.
x(t) =
Therefore, the result holds for l = . We assume that the result holds for l – . Now, we deal with the case for l. Let x (t) = u(t), then the BVP (.) is equivalent to the following BVPs:
x (t) = u(t), ≤ t ≤ , x() = kn–l– (s)x(s) ds,
(.)
x() =
and
u((l–)) (t) = y(t), ≤ t ≤ , u(i) () = kn–l+i (s)u(i) (s) ds,
u(i) () = ,
≤ i ≤ l – .
(.)
By applying Lemma ., the BVP (.) has a unique solution
gn–l– (t, r)u(r) dr.
x(t) =
(.)
Replacing l by l – and x by u in (.), by applying the inductive hypothesis, the BVP (.) has also a unique solution
Gl– (t, s)y(s) ds.
u(t) =
(.)
Substituting (.) into (.), we see that the BVP (.) has a unique solution
x(t) =
gn–l– (t, r)
=
Gl– (r, s)y(s) ds dr
gn–l– (t, r)Gl– (r, s) dr y(s) ds
=
Gl (t, s)y(s) ds.
Therefore, the result holds for l. Lemma . is now completed. For each ≤ l ≤ n – , we define Al : C[, ] → C[, ] by Al u(t) =
Gl (t, τ )u(τ ) dτ .
Ji et al. Boundary Value Problems (2015) 2015:214
Page 5 of 11
With the use of Lemma ., for each ≤ l ≤ n – , we have
(Al u)(l) (t) = u(t), ≤ t ≤ , (Al u)(i) () = kn–l+i– (s)(Al u)(i) (s) ds,
(Al u)(i) () = ,
≤ i ≤ l – .
Therefore (.) has a solution if and only if the BVP
u (t) = f (t, An– u(t), An– u(t), . . . , A u(t), u(t)), u() = u() = kn– (s)u(s) ds,
≤ t ≤ ,
(.)
has a solution. If x is a solution of (.), then u = x((n–)) is a solution of (.). Conversely, if u is a solution of (.), then x = An– u is a solution of (.). In addition if (–)n– u(t) ≥ ( = ) on [, ], then x = An– u is a positive solution of (.). For ≤ i ≤ n – , let mi = min
–τ
t∈[τ ,–τ ] τ
gi (t, s) ds and
Mi = max
t∈[,]
gi (t, s) ds.
(.)
Obviously, < mi < Mi . Let E denote the Banach space C[, ] with the maximum norm u = max u(t), ≤t≤
and define the cone K ⊂ E by K = u ∈ E(–)n– u(t) ≥ , (–)n– u(t) is concave on [, ], and min (–)n– u(t) ≥ τ u . t∈[τ ,–τ ]
Finally, we define the nonnegative continuous concave functional α on K by α(u) = min u(t) t∈[τ ,–τ ]
for each u ∈ K and it is easy to see that α(u) ≤ u .
4 Main results Theorem . Suppose conditions (H ), (H ) hold. In addition assume there exist nonnegmn– }c and f (t, un– , un– , . . . , u , u ) ative numbers a, b, and c such that < a < b ≤ min{τ , M n– satisfies the following growth conditions: c , for (t, un– , un– , . . . , u , u ) ∈ [, ] × (H ) (–)n f (t, un– , un– , . . . , u , u ) ≤ Mn– j+ n––j n– [, i= Mn–i c] × (–) [, c]; j=n– (–) (H ) (–)n f (t, un– , un– , . . . , u , u ) < Man– , for (t, un– , un– , . . . , u , u ) ∈ [, ] × j+ n––j [, i= Mn–i a] × (–)n– [, a]; j=n– (–) b , for (t, un– , un– , . . . , u , u ) ∈ [τ , – τ ] × (H ) (–)n f (t, un– , un– , . . . , u , u ) ≥ mn– j+ j+ n––j [ i= mn–i b, i= Mn–i τb ] × (–)n– [b, τb ]. j=n– (–)
Ji et al. Boundary Value Problems (2015) 2015:214
Page 6 of 11
Then the BVP (.) has at least three positive solutions x , x , and x , which satisfy ((n–))
x (t) < a, b < α x((n–)) (t) and ((n–))
x (t) > a with α x((n–)) (t) < b. Proof Define the completely continuous operator A by
Au(t) =
gn– (t, s)f s, An– u(s), . . . , A u(s), u(s) ds.
We will first verify that A : K → K . Let u ∈ K , then (–)n– Au(t) ≥ , ((–)n– Au) (t) = (–)n– f (t, An– u(t), . . . , A u(t), u(t)) ≤ , ≤ t ≤ , and by Proposition . of [], we have Au = max Au(t) t∈[,]
= max gn– (t, s)f s, An– u(s), . . . , A u(s), u(s) ds t∈[,] –t h(s, τ )kn– (τ ) dτ h(t, s) + = max t∈[,] – Kn–
× f s, An– u(s), . . . , A u(s), u(s) ds h(τ , τ )kn– (τ ) dτ h(s, s) + ≤ – Kn– × f s, An– u(s), . . . , A u(s), u(s) ds. On the other hand, by Proposition . of [], we obtain min (–)n– Au(t)
t∈[τ ,–τ ]
= min (–)
n–
t∈[τ ,–τ ]
gn– (t, s)f s, An– u(s), . . . , A u(s), u(s) ds
–t h(t, s) + h(s, τ )kn– (τ ) dτ = min (–) t∈[τ ,–τ ] – Kn–
× f s, An– u(s), . . . , A u(s), u(s) ds h(s, s) + h(τ , τ )kn– (τ ) dτ ≥τ – Kn– × f s, An– u(s), . . . , A u(s), u(s) ds n
≥ τ Au . Consequently, A : K → K . It is a standard argument to show that the operator A is completely continuous. Equicontinuity and uniform boundedness follow readily from the properties of Gl , ≤ l ≤ n – . If u ∈ Pc , then u ≤ c. For each ≤ j ≤ n – , note that inductively (using (.)) we have
j+ Mn–i c. Aj u = max Gj (t, s)u(s) ds ≤ t∈[,]
i=
Ji et al. Boundary Value Problems (2015) 2015:214
Page 7 of 11
From the condition (H ) and (.), we obtain Au = max Au(t) t∈[,]
= max gn– (t, s)f s, An– u(s), An– u(s), . . . , A u(s), u(s) ds t∈[,] c gn– (t, s) ds = c. max ≤ Mn– t∈[,] So, A : Pc → Pc . In a completely analogous argument, the condition (H ) implies the condition (C ) of Theorem . is satisfied. We now show that condition (C ) is satisfied. Note that, for ≤ t ≤ , u(t) = (–)n–
b b ∈ P α, b, τ τ
and
α(u) =
b > b. τ
Thus, b u ∈ P α, b, α(u) > b = ∅. τ Also, if u ∈ P(α, b, τb ), then b ≤ (–)n– u(t) ≤ n – , t ∈ [τ , – τ ], inductively, (–)n––j Aj u(t) ≤
j+
Mn–i
i=
b τ
for t ∈ [τ , – τ ], implies for each ≤ j ≤
b , τ
n––j (–) Aj u(t) = Gj (t, s)u(s) ds ≥ b
j+
Gj (t, s) ds ≥ mn–i b.
–τ
τ
i=
With the use of condition (H ) and (.), we get
α(Au) = min gn– (t, s)f s, An– u(s), . . . , A u(s), u(s) ds t∈[τ ,–τ ] –τ
> min gn– (t, s)f s, An– u(s), . . . , A u(s), u(s) ds t∈[τ ,–τ ]
≥
τ
b min mn– t∈[τ ,–τ ]
gn– (t, s) ds = b.
Therefore, condition (C ) is satisfied. Finally, we show that condition (C ) is also satisfied. That is, we show that if u ∈ P(α, b, c) and Au > d = τb , then α(Au) > b. This follows since A : K → K . In particular, since (–)n– (Au) is concave and mint∈[τ ,–τ ] (–)n– (Au)(t) ≥ τ Au . That is, α(Au) = min Au(t) ≥ τ Au > b. t∈[τ ,–τ ]
Ji et al. Boundary Value Problems (2015) 2015:214
Page 8 of 11
Therefore, condition (C ) is also satisfied. By Theorem ., there exist three solutions u , u , u ∈ K for the BVP (.). Moreover, let
xi (t) = An– ui (t) =
Gn– (t, s)ui (s) ds,
i = , , ,
then x , x , x are three positive solutions for the BVP (.) and satisfy ((n–))
x (t) < a, b < α x((n–)) (t) and ((n–))
x (t) > a with α x((n–)) (t) < b.
Consider the following nth order differential equations with integral boundary conditions:
x(n) (t) = f (t, x(t), x (t), . . . , x((n–)) (t)), ≤ t ≤ , x(i) () = , x(i) () = ki∗ (s)x(i) (s) ds, ≤ i ≤ n – ,
(.)
∗ j where (–)n f ∈ C([, ] × n– j= (–) [, +∞) → [, +∞)) and ki (t) ∈ L [, ] (i = , , . . . , n – ) are nonnegative. Now we deal with problem (.). The method is just similar to what we have done for the problem (.), so we omit the proof of main results in this section. For convenience, we list the following assumptions: (H∗ ) ki∗ (t) ∈ L [, ] are nonnegative, and Ki∗ ∈ [, ), where Ki∗ =
ski∗ (s) ds,
≤ i ≤ n – .
By analogous methods, we have the following results. Lemma . Suppose (H∗ ) holds. Then gi∗ (t, s) ≤ ( ≤ i ≤ n – ), where gi∗ (t, s) is the Green’s function for the problem
x (t) = , x() = ,
≤ t ≤ , x() = ki∗ (s)x(s) ds.
For ≤ i ≤ n – , let m∗i = min
t∈[τ ,–τ ] τ
–τ
g ∗ (t, s) ds i
and Mi∗ = max
t∈[,]
g ∗ (t, s) ds. i
Theorem . Suppose condition (H∗ ) holds. In addition assume there exist nonnegative m∗ numbers a, b, and c such that < a < b ≤ min{τ , Mn– ∗ }c and f (t, un– , un– , . . . , u , u ) satn– isfies the following growth conditions: (H∗ ) (–)n f (t, un– , un– , . . . , u , u ) ≤ M∗c , for (t, un– , un– , . . . , u , u ) ∈ [, ] × n– j+ ∗ n––j n– (–) [, M c] × (–) [, c]; n–i i= j=n– a ∗ n (H ) (–) f (t, un– , un– , . . . , u , u ) < M∗ , for (t, un– , un– , . . . , u , u ) ∈ [, ] × n– j+ ∗ n––j n– (–) [, M a] × (–) [, a]; n–i i= j=n–
Ji et al. Boundary Value Problems (2015) 2015:214
Page 9 of 11
(H∗ ) (–)n f (t, un– , un– , . . . , u , u ) ≥ m∗b , for (t, un– , un– , . . . , u , u ) ∈ [τ , – τ ] × j+ ∗ j+ n– ∗ b n––j n– (–) [ m b, [b, τb ]. i= n–i i= Mn–i τ ] × (–) j=n– Then the BVP (.) has at least three positive solutions x , x , and x , which satisfy ((n–))
x (t) < a, b < α x((n–)) (t) and ((n–))
x (t) > a with α x((n–)) (t) < b.
5 Example Example . As an example of problem (.), consider the following sixth order BVP: ⎧ () x (t) = f (t, x(t), x (t), x() (t)), ≤ t ≤ , ⎪ ⎪ ⎪ ⎨ x() = , x() = x(s) ds, ⎪ () = sx (s) ds, x () = , x ⎪ ⎪ () ⎩ () x () = s x (s) ds, x() () = ,
(.)
where f (t, u , u , u ) ⎧ , u · +sin(t+u +u ) , – ⎪ ⎪ ⎪ , , ⎪ +sin(t+u +u ) , ⎪ , ⎪ ⎪ –[( × , – )(u – ) + ] · ⎪ ⎨ – × , u · +sin(t+u +u ) , , = , ⎪ –[ × , + ( × – × , )(u – )] ⎪ , , ⎪ ⎪ +sin(t+u +u ) ⎪ ⎪ · , ⎪ ⎪ ⎩ +u ) · +sin(t+u , – × ,
≤ u ≤ , ≤ u ≤ , ≤ u ≤ , ≤ u ≤ , u ≥ .
We notice that n = , ki (s) = si (i = , , ) and K = , K = , K = If we take τ = , by calculation we obtain m = min
t∈[ , ]
m = min
t∈[ , ]
m = min
t∈[ , ]
g (t, s) ds = , g (t, s) ds = , ,
M = max
g (t, s) ds = ,
t∈[,]
M = max
g (t, s) ds = , , ,
.
t∈[,]
g (t, s) ds = ,
M = max
t∈[,]
In addition, if we take a = , b = , c = , then m < a = < b = ≤ min τ , c M , × , = min × = , , , × and f (t, u , u , u ) satisfies the growth conditions (H )-(H ).
g (t, s) ds = . ,
Ji et al. Boundary Value Problems (2015) 2015:214
Page 10 of 11
Therefore all the conditions of Theorem . are satisfied. Hence, the problem (.) has at least three positive solutions x , x , and x , which satisfy max x() (t) < ,
≤t≤
< min x() (t) t∈[ , ]
and
() max x() (t) > with min x (t) < . t∈[ , ]
≤t≤
Example . As another example of problem (.), consider the following sixth order BVP: ⎧ () x (t) = f (t, x(t), x (t), x() (t)), ≤ t ≤ , ⎪ ⎪ ⎪ ⎨ x() = , x() = x(s) ds, ⎪ () = , x () = x ⎪ sx (s) ds, ⎪ ⎩ () () x () = s x() (s) ds, x () = ,
(.)
where f (t, u , u , u ) ⎧ +sin(t+u +u ) ⎪ , ⎪ – u · ⎪ ⎪ +sin(t+u +u ) , ⎪ ⎪ , ⎨ –[( × , – )(u – ) + ] · +sin(t+u +u ) = – × , u · , , ⎪ ⎪ , ⎪ + ( × – × , )(u – )] · –[ × ⎪ , , ⎪ ⎪ ⎩ – × · +sin(t+u +u ) ,
≤ u ≤ , ≤ u ≤ , ≤ u ≤ , +sin(t+u +u ) , ≤ u ≤ , u ≥ .
We notice that n = , ki∗ (s) = si (i = , , ) and K∗ = , K∗ = , K∗ = . If we take τ = , by calculation we obtain m∗ m∗ m∗
= min
t∈[ , ]
= min
t∈[ , ]
= min
t∈[ , ]
∗ g (t, s) ds = ,
M∗
∗ g (t, s) ds = , ,
M∗
∗ g (t, s) ds = , , ,
g ∗ (t, s) ds = ,
= max
t∈[,]
= max
M∗
t∈[,]
g ∗ (t, s) ds = ,
= max
t∈[,]
g ∗ (t, s) ds = .
In addition, if we take a = , b = , c = , then ∗ , × m < a = < b = ≤ min τ , ∗ c = min × = , , M , × and f (t, u , u , u ) satisfies the growth conditions (H∗ )-(H∗ ). Therefore all the conditions of Theorem . are satisfied. Hence, the problem (.) has at least three positive solutions x , x , and x , which satisfy max x() (t) < ,
≤t≤
< min x() (t) t∈[ , ]
and
() max x() (t) > with min x (t) < .
≤t≤
t∈[ , ]
Ji et al. Boundary Value Problems (2015) 2015:214
Page 11 of 11
Competing interests The authors declare that there is no conflict of interests regarding the publication of this article. Authors’ contributions The authors declare that the study was realized in collaboration with the same responsibility. All authors read and approved the final manuscript. Author details 1 College of Sciences, Hebei University of Science and Technology, Shijiazhuang, Hebei 050018, P.R. China. 2 School of Electrical Engineering, Hebei University of Science and Technology, Shijiazhuang, Hebei 050018, P.R. China. 3 Hebei Medical University, Shijiazhuang, Hebei 050200, P.R. China. Acknowledgements The authors express their sincere thanks to the referees for the careful and details reading of the manuscript and very helpful suggestions. The project was supported by the Natural Science Foundation of China (11371120), the Natural Science Foundation of Hebei Province (A2013208147, A2015208051). Received: 13 July 2015 Accepted: 12 November 2015 References 1. Chasreechai, S, Tariboon, J: Positive solutions to generalized second-order three-point integral boundary-value problems. Electron. J. Differ. Equ. 2011, 14 (2011) 2. Feng, M: Existence of symmetric positive solutions for a boundary value problem with integral boundary conditions. Appl. Math. Lett. 24, 1419-1427 (2011) 3. Galvis, J, Rojas, EM, Sinitsyn, AV: Existence of positive solutions of a nonlinear second-order boundary-value problem with integral boundary conditions. Electron. J. Differ. Equ. 2015, 236 (2015) 4. Wu, H, Zhang, J: Positive solutions of higher-order four-point boundary value problem with p-Laplacian operator. J. Comput. Appl. Math. 233, 2757-2766 (2010) 5. Jankowski, T: Positive solutions for fourth-order differential equations with deviating arguments and integral boundary conditions. Nonlinear Anal., Theory Methods Appl. 73, 1289-1299 (2010) 6. Zhao, J, Wang, P, Ge, W: Existence and nonexistence of positive solutions for a class of third order BVP with integral boundary conditions in Banach spaces. Commun. Nonlinear Sci. Numer. Simul. 16, 402-413 (2011) 7. Zhang, X, Ge, W: Symmetric positive solutions of boundary value problems with integral boundary conditions. Appl. Math. Comput. 219, 3553-3564 (2012) 8. Yang, W: Positive solutions for a coupled system of nonlinear fractional differential equations with integral boundary conditions. Comput. Math. Appl. 63, 288-297 (2012) 9. Cabada, A, Wang, G: Positive solutions of nonlinear fractional differential equations with integral boundary value conditions. J. Math. Anal. Appl. 389, 403-411 (2012) 10. Agarwal, RP, Wang, G, Ahmad, B, Zhang, L, Hobiny, A, Monaquel, S: On existence of solutions for nonlinear q-difference equations with nonlocal q-integral boundary conditions. Math. Model. Anal. 20, 604-618 (2015) 11. Pang, H, Xie, W, Cao, L: Successive iteration and positive solutions for a third-order boundary value problem involving integral conditions. Bound. Value Probl. 2015, 139 (2015) 12. Wang, Q, Guo, Y, Ji, Y: Positive solutions for fourth-order nonlinear differential equation with integral boundary conditions. Discrete Dyn. Nat. Soc. 2013, Article ID 684962 (2013) 13. Jiang, W, Zhang, Y, Qiu, J: The existence of solutions for p-Laplacian boundary value problems at resonance on the half-line. Bound. Value Probl. 2015, 179 (2015) 14. Guo, Y, Liu, X, Qiu, J: Three positive solutions for higher order m-point boundary value problems. J. Math. Anal. Appl. 289, 545-553 (2004) 15. Ahmad, B, Alsaedi, A, Alghamdi, BS: Analytic approximation of solutions of the forced Duffing equation with integral boundary conditions. Nonlinear Anal., Real World Appl. 9, 1727-1740 (2008) 16. Ahmad, B, Ntouyas, SK, Alsulami, HH: Existence theory for nth order nonlocal integral boundary value problems and extension to fractional case. Abstr. Appl. Anal. 2013, Article ID 183813 (2013) 17. Hao, X, Liu, L: Multiple monotone positive solutions for higher order differential equations with integral boundary conditions. Bound. Value Probl. 2014, 74 (2014) 18. Li, Y, Zhang, X: Positive solutions for systems of nonlinear higher order differential equations with integral boundary conditions. Abstr. Appl. Anal. 2014, Article ID 591381 (2014) 19. Zhang, X, Ge, W: Positive solutions for a class of boundary-value problems with integral boundary conditions. Comput. Math. Appl. 58, 203-215 (2009) 20. Leggett, RW, Williams, LR: Multiple positive fixed points of nonlinear operators on ordered Banach spaces. Indiana Univ. Math. J. 28, 673-688 (1979)