Asymptotic behavior of positive solutions of nonlinear fractional differential equations with Caputo-type Hadamard derivative

In this paper we are concerned with the problem of asymptotic integration of positive solutions of higher order fractional differential equations with Caputo-type Hadamard derivative of the form C,H Dr ax(t) = e(t) + f (t, x(t)), a > 1, where r = n + α − 1, α ∈ (0, 1), n ∈ Z+. We shall apply our technique to investigate the oscillatory and asymptotic behavior of all solutions of the integral equation x(t) = e(t) + ∫ t a (ln t s ) r−1k(t, s) f (s, x(s)) ds s , a > 1, r is as above.


Consider the initial value problem
C,H D r a x(t) = e(t) + f (t, x(t)), a > 1, and the integral equation where r = n + α − 1, α ∈ (0, 1), n ∈ Z + , C,H D r a x(t) is the Caputo-type Hadamard modification of fractional derivative of a C n − scalar valued function x(t) defined on the interval [a, ∞), which was recently proposed by Jarad et al. [1].
We only consider those solutions of equation ( 1) that are continuable and nontrivial in any neighborhood of ∞.Such a solution is said to be oscillatory if there exists a sequence t n ⊆ [c, ∞) t n → ∞ such that x(t n ) = 0, and it is nonoscillatory otherwise.The subject of fractional differential equations has recently evolved as an interesting and popular field of research.In fact, fractional derivatives provide an excellent tool for the description of memory and hereditary properties of various materials and processes.Many researchers have found that fractional differential equations play important roles in many research areas, such as physics, chemical technology, population dynamics, biotechnology, and economics [2][3][4][5][6][7][8][9][10][11][12][13][14].
However, it has been observed that most of the work on the topic involves either Riemann-Liouville or Caputo-type fractional derivative.Besides these derivatives, Hadamard fractional derivative is another kind of fractional derivatives that was introduced by Hadamard in 1892 [15].This fractional derivative differs from the other ones in the sense that the kernel of the integral (in the definition of Hadamard derivative) contains logarithmic function of arbitrary exponent.For background material of Hadamard fractional derivative and integral, we refer [4,[16][17][18][19][20][21][22][23][24][25][26][27].
The asymptotic behavior results for fractional differential equations with Caputo-type Hadamard derivative are scarce and it seems that there are no such results for such equations.The main objective of this paper is to establish some new criteria for the asymptotic behavior of all positive solutions of equation (1) as well as the oscillatory and asymptotic behavior of the integral equation (2).
We introduce some notations and definition of fractional calculus [3,4].
Definition 1.For at least n−times differentiable function f : (a, +∞) → R the Caputo-type Hadamard derivative of fractional order α is defined as where δ = t d dt and [α] denotes the integer part of the real number α.
Definition 2. The Hadamard fractional integral of order r ∈ R + of function f (t) for all t > a > 0 is defined by where Γ(α) is the Euler Gamma function.
Definition 3. The Hadamard derivative of order r ∈ [n − 1, n), n ∈ Z + of a function f (t) is given by The corresponding Volterra integral equation can be easily derived using Lemma 2.5 in [1] and it takes the form

Asymptotic Behavior of equation (1)
To obtain our main results of this paper, we need the following two lemmas.

Proof. Let ln
Lemma 5. [28].If X and Y are nonnogative, then where equality holds if and only if X = Y.
In what follows, we let and 0 < λ < 1, t ≥ t 1 for some t 1 ≥ a, where m : [a, ∞) → (0, ∞) is a continuous function.Now we give sufficient conditions under which any eventually positive solution x of equation ( 1) satisfies Theorem 6.Let 0 < λ < 1 and conditions ( 1), ( 2) hold and suppose that p > 1, p(r − 1) and If x is positive solution of equation (1), then Proof.Let x(t) be an eventually positive solution of equation ( 1), say x(t) > 0 for t ≥ t 1 for some t 1 ≥ a.Since equation ( 1) is equivalent to equation (3), we see that there exists a constant C 1 > 0 such that We let F(t) = e(t) + f (t, x(t)).In view of condition (2), we can write and so, or, where C is the upper bound of the function ln t a By applying the Holder inequality and Lemma 4, we obtain (m(s)) q x q (s)ds (m(s)) q x q (s)ds (m(s)) q x q (s)ds , where B = B[p(γ − 1) + 1, p(r − 1) + 1] and θ = p(r + γ − 2) + 1 and by using γ = (n − r) + 1 q , we have (m(s)) q x q (s)ds Using ( 13) in ( 12), we have Using this inequality and the elementary inequality (x + y) q ≤ 2 q−1 (x q + y q ), x, y ≥ 0 and q > 1, we obtain from ( 12) (n−1)q m q (s)z q (s)ds .
If we denote u(t) = z q (t), i.e. z(t) = u 1/q (t), The conclusion follows from Gronwall's inequality and we conclude that lim sup t→∞ ln t a Similar to the sublinear case, one can easily prove the following results.
Theorem 7. Let λ = 1 and the hypotheses of Theorem 6 hold with m(t) = h(t) and g(t) = 0. Then the conclusion of Theorem 6 holds.
Example 1.Let p > 1, r = α.i.e., n = 1, 0 < α = 1 − 1 2p < 1, α = γ and q = p p−1 .Clearly, Let the functions a(t) and b(t) be as in (i) and (ii) with b(t) be a bounded function and let and let e : R + → R be a continuous function with lim All conditions of Theorem 6 are satisfied and hence every positive solution x of equation ( 1) is bounded.

Oscillatory and asymptotic behavior of equation (2)
In this section, we present the following result on the asymptotic behavior of nonoscillatory solutions of equation ( 2).Theorem 8. Let 0 < λ < 1 and conditions (1), (3) hold and suppose that p > 1, q = p p−1 , γ = (n − r) + 1 q , p(r − 1) + 1 > 0, p(γ − 1) + 1 > 0, the functions b(t) is bounded, (7) Proof.Let x(t) be an eventually positive solution of equation ( 2).We may assume that for t ≥ t 1 for some We let F(t) = f (t, x(t)).In view of conditions ( 1) and (3) we may then write Proceeding exactly as in the proof of Theorem 6, we obtain Applying Holder's inequality and Lemma 4, we obtain where It follows from ( 18) and ( 20) that ln t a where C and c are the upper bound of the functions ln t a respectively.The rest of the proof is similar to that of Theorem 6 and conclude that lim sup t→∞ ln t a If x(t) is eventually negative, we can set y = −x to see that y satisfies equation ( 2) with e(t) be replaced by −e(t) and f (t, x) by − f (t, −y).It follows in a similar manner that lim sup t→∞ ln t a It follows from ( 22) and ( 23) that ( 16) holds.This completes the proof.
Next, by employing Theorem 6 we present the following oscillation result for equation (2).
Proof.Let x(t) be a nonoscillatory solution of equation ( 2), say x(t) > 0 for t ≥ t 1 for some t 1 ≥ 0. The proof when x(t) is eventually negative is similar.Proceeding as in the proof of Theorem 7 we arrive at (18) Clearly, the conclusion of Theorem 6 holds.This together with ( 7) and ( 8) imply that the second integral in this inequality is bounded and hence one can easily see that where M is a positive constant.Finally, taking lim inf in (25) as t → ∞ and using condition (22) result in a contradiction with the fact that x(t) is eventually positive.
Theorem 10.Let λ = 1 and the hypotheses of Theorems 8 or 9 hold with m(t) = h(t) and g(t) = 0. Then the conclusions of Theorems 8 and 9 hold.All conditions of Theorem 8 are satisfied and hence every nonocillatory solution x(t) of equation ( 2) satisfies (16).In addition, if the function e(t) satisfies condition (24), then equation ( 2) is oscillatory.We may take e(t) = (ln t) sin t.
Case (I) If γ > (n − r) + 1 q , then θ > p(r + γ − 2) + 1 and so σ > 0. In this case we need to impose the condition that the function b(t) ln t a σ is bounded and the results are valid.
and the results remain valid.The details are left to the reader.

General Remarks:
1.The results of this paper are presented in a form which is essentially new and of higher degree of generality.2. It would be of interest to study equations ( 1) and (2) when f satisfies condition (ii) with λ > 1.

Conclusion
We are concerned with the problem of asymptotic integration of positive solutions of higher order fractional differential equations with Caputo-type Hadamard derivative, and establish some new criteria for the asymptotic behavior of all positive solutions of equation ( 1) as well as the oscillatory and asymptotic behavior of the integral equation ( 2).