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

1. Introduction

Consider the initial value problem and the integral equation where \(r = n +\alpha -1, \;\alpha\in (0, 1), \;n \in \mathbb{Z}^+\), \(^{C,H}D_{a}^r 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,\infty)\), which was recently proposed by Jarad et al. [1]. $$ ^{C,H}D_{a}^r x(t)=^{C,H}J_a ^{n-r}\delta ^{n}x(t), \; \delta: =t \frac{d}{dt}, $$ where \(r = n +\alpha -1, \;\alpha\in (0, 1), \;n \in\mathbb{Z}^+\). In the sequel we assume that:

1. \(e:[a, \infty)\rightarrow R^+=(0, \infty)\) is a continuous function;
2. \(k: [a, \infty)\times [c, \infty)\rightarrow \mathbb{R}\) is a continuous function, and assume that there exists a continuous function, \(b:[a, \infty)\rightarrow \mathbb{R}^+\) such that \(|k(t, s)| \leq b(t)\) for all \(t>s>a\);
3. \(f: [a, \infty)\times \mathbb{R} \rightarrow \mathbb{R}\) is continuous and assume that there exists a continuous function, \(h:[a, \infty)\rightarrow (0, \infty)\) and a real number \(\lambda\), \(0< \lambda \leq 1\) and \(\gamma\) is a real number such that \(0\leq xf(t,x)\leq t^{\gamma-1} h(t)|x|^{\lambda+1}, \;\text{for}\; x\neq 0, \;\text{and}\; t\geq a \).

We only consider those solutions of equation (1) that are continuable and nontrivial in any neighborhood of \(\infty\). Such a solution is said to be oscillatory if there exists a sequence \({t_n} \subseteq [c, \infty)\; t_n\rightarrow \infty\) 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, , 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,+\infty)\to \mathbb{R}\) the Caputo-type Hadamard derivative of fractional order \(\alpha\) is defined as $$^{C,H}D_a ^r f(t)=\frac{1}{\Gamma{(n-\alpha)}}\int_{a}^{t}(\ln\frac{t}{s} )^{n-\alpha-1} \delta^{n}f(s)\frac{d s}{s},\; n-1< \alpha< n, \;n=[\alpha],$$ where \(\delta=t\frac{d}{dt}\) and \([\alpha]\) denotes the integer part of the real number \(\alpha\).

Definition 2. The Hadamard fractional integral of order \(r\in \mathbb{R}^+\) of function \(f(t)\) for all \(t>a>0\) is defined by $$^{H}D_a ^{-\alpha} f(t)=\frac{1}{\Gamma{(\alpha)}}\int_{a}^{t}(\ln\frac{t}{s} )^{\alpha-1} f(s)\frac{d s}{s}, $$ where \(\Gamma (\alpha)\) is the Euler Gamma function.

Definition 3. The Hadamard derivative of order \(r\in[n-1,n), \;n\in \mathbb{Z}^+\) of a function \(f(t)\) is given by $$^{H}D_a ^\alpha f(t)=\frac{1}{\Gamma{(n-\alpha)}}(t\frac{d}{dt})^n\int_{a}^{t}(\ln\frac{t}{s} )^{n-\alpha-1} f(s)\frac{d s}{s}. $$

The corresponding Volterra integral equation can be easily derived using Lemma 2.5 in [1] and it takes the form

2. Asymptotic Behavior of equation (1)

To obtain our main results of this paper, we need the following two lemmas.

Lemma 4. Let \(\beta\), \(\gamma\) and \(p\) be positive constants such that \([p(\beta-1)+1]>0\), \(p(\gamma-1)+1>0\). Then

Where \(B:=B[p(\gamma-1)+1,p(\beta-1)+1]\), \(B[\xi, \eta]=\int_0^1s^{\xi-1}(1-s)^{\eta-1}ds\), \((\xi>0,\;\eta>0)\) and \(\theta=p(\beta+\gamma-2)+1\).

Proof. Let \(\ln\frac{s}{a}=\tau\ln\frac{t}{a}\). Then \(s=a\left(\frac{t}{a}\right)^\tau\), \(ds=a\left(\frac{t}{a}\right)^\tau\ln\left(\frac{t}{a}\right)d\tau\) and \(\frac{ds}{s}=\ln\left(\frac{t}{a}\right)d\tau\). Therefore we obtain \begin{equation*}\begin{array}{ll} \int_a^t\left(\ln\frac{t}{s}\right)^{p(\beta-1)}\left(\ln\frac{s}{a}\right)^{p(\gamma-1)}\frac{ds}{s}&=\int_a^1\left(\ln\left(\frac{t}{a}\left (\frac{t}{a}\right)^{-\tau}\right)\right)^{p(\beta-1)}\left(\tau\ln\left(\frac{t}{a}\right)\right)^{p(\gamma-1)}\ln\left(\frac{t}{a}\right)d\tau\\\\ &=\left(\ln\frac{t}{a}\right)^{p(\beta+\gamma-2)+1}\int_0^1(1-\tau)^{p(\beta-1)}\tau^{p(\gamma-1)}d\tau. \end{array} \end{equation*}

Lemma 5. [1]. If \(X\) and \(Y\) are nonnogative, then

where equality holds if and only if \(X=Y\).

In what follows, we let and \(0< \lambda< 1 \), \(t\geq t_1\) for some \(t_1\geq a\), where \(m:[a,\infty)\rightarrow(0,\infty)\) is a continuous function. Now we give sufficient conditions under which any eventually positive solution \(x\) of equation (1) satisfies $$x(t)=O\left(\ln\frac{t}{a}\right)^{n-1}\;\text{as}\;t\rightarrow\infty.$$

Theorem 6. Let \(0< \lambda< 1\) and conditions (1), (2) hold and suppose that \(p>1\), \(p(r-1)+1>0\), \(p(\gamma-1)+1>0\), \(q=\frac{p}{p-1}\), \(\gamma=(n-r)+\frac{1}{q}\),

and If \(x\) is positive solution of equation (1), then $$\limsup_{t\rightarrow\infty}\left(\ln\frac{t}{a}\right)^{1-n}x(t)< \infty.$$

Proof. Let \(x(t)\) be an eventually positive solution of equation (1), say \(x(t)>0\) for \(t\geq t_1\) for some \(t_1\geq 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 \begin{equation*}\begin{array}{ll} x(t)&\leq \ \ \ \ C_1\left(\ln\frac{t}{a}\right)^{n-1}+\frac{1}{\Gamma(r)}\int_a^{t_1}\left(\ln\frac{t}{s}\right)^{r-1}|F(s)|\frac{ds}{s}+\frac{1}{\Gamma(r)} \int_{t_1}^t\left(\ln\frac{t}{s}\right)^{r-1}e(s)\frac{ds}{s}\\\\ &~~+\frac{1}{\Gamma(r)}\int_{t_1}^t\ln\left(\frac{t}{s}\right)^{r-1}\left(\ln\frac{s}{a}\right)^{\gamma-1}[h(s)x^\lambda(s)-m(s)x(s)]\frac{ds}{s}\\\\ &~~+\frac{1}{\Gamma(r)}\int_{t_1}^t\left(\ln\frac{t}{s}\right)^{r-1}\left(\ln\frac{s}{a}\right)^{\gamma-1}m(s)x(s)\frac{ds}{s}. \end{array} \end{equation*} Applying (5) of Lemma 5 to \([h(s)x^\lambda(s)-m(s)x(s)]\) with $$X=(h)^{1/\lambda}x\;\;\text{and}\;\;Y=\left(\frac{1}{\lambda}m(h^{-1/\lambda})\right)^{1/(\lambda-1)},$$ we have $$h(s)x^\lambda(s)-m(s)x(s)\leq(1-\lambda)\lambda^{\lambda/(1-\lambda)}m^{\lambda/(\lambda-1)}(s)h^{1/(1-\lambda)}(s),$$ and so, \begin{equation*}\begin{array}{ll} x(t)&\leq \ \ \ \ C_1\left(\ln\frac{t}{a}\right)^{n-1}+\frac{1}{\Gamma(r)}\int_a^{t_1}\left(\ln\frac{t}{s}\right)^{r-1}|F(s)|\frac{ds}{s}+\frac{1}{\Gamma(r)} \int_{t_1}^t\left(\ln\frac{t}{s}\right)^{r-1}|e(s)|\frac{ds}{s}\\\\ &~~+\frac{(1-\lambda)\lambda^{\lambda/(1-\lambda)}}{\Gamma(r)}\int_{t_1}^t\ln\left(\frac{t}{s}\right)^{r-1}\left(\ln\frac{s}{a}\right)^{r-1} [m^{\lambda/(\lambda-1)}(s)h^{1/(1-\lambda)}(s)]\frac{ds}{s}\\\\ &~~+\frac{1}{\Gamma(r)}\int_{t_1}^t\left(\ln\frac{t}{s}\right)^{r-1}\left(\ln\frac{s}{a}\right)^{\gamma-1}m(s)x(s)\frac{ds}{s}. \end{array} \end{equation*} or, or, where \(C\) is the upper bound of the function $$\left(\ln\frac{t}{a}\right)^{1-n}\left[\frac{1}{\Gamma(r)}\left(\int_a^{t_1}\left(\ln\frac{t}{s}\right)^{r-1}|F(s)|\frac{ds}{s}\right)+g(t)+ \frac{1}{\Gamma(r)}\int_{t_1}^t\left(\ln\frac{t}{s}\right)^{r-1}|e(s)|\frac{ds}{s}\right].$$ By applying the Holder inequality and Lemma 4, we obtain \begin{equation*}\begin{array}{ll} \int_{t_1}^t\left(\ln\frac{t}{s}\right)^{r-1}\left(\ln\frac{t}{a}\right)^{\gamma-1}m(s)x(s)\frac{ds}{s}&\leq\left(\int_{t_1}^t\left(\ln\frac{t}{s}\right) ^{p(r-1)}\left(\ln\frac{t}{a}\right)^{p(\gamma-1)}m(s)x(s)\frac{ds}{s}\right)^{1/p}\left(\int_{t_1}^t(m(s))^qx^q(s)ds\right)^{1/q}\\\\ &\leq\left(\int_0^t\left(\ln\frac{t}{s}\right) ^{p(r-1)}\left(\ln\frac{t}{a}\right)^{p(\gamma-1)}\frac{ds}{s}\right)^{1/p}\left(\int_{t_1}^t(m(s))^qx^q(s)ds\right)^{1/q}\\\\ &\leq\left(B\left(\ln\frac{t}{a}\right)^\theta\right)^{1/p}\left(\int_{t_1}^t(m(s))^qx^q(s)ds\right)^{1/q}, \end{array} \end{equation*} where \(B=B[p(\gamma-1)+1,p(r-1)+1]\) and \(\theta=p(r+\gamma-2)+1\) and by using \(\gamma=(n-r)+\frac{1}{q}\), we have Using (13) in (12), we have $$z(t)\leq 1+C+\frac{B^{1/p}}{\Gamma(r)}\int_{t_1}^t\left(\ln\frac{s}{a}\right)^{(n-1)q}(m(s)z(s))^qds.$$ Using this inequality and the elementary inequality $$(x+y)^q\leq 2^{q-1}(x^q+y^q),\;x,y\geq 0\;\text{and}\;q>1,$$ we obtain from (12) $$z^q(t)\leq 2^{q-1}\left((1+C)^q+\left(\frac{B^{1/p}}{\Gamma(r)}\right)^q\int_{t_1}^t\left(\ln\frac{s}{a}\right)^{(n-1)q}m^q(s)z^q(s)ds\right).$$ If we denote \(u(t)=z^q(t)\), i.e. \(z(t)=u^{1/q}(t)\), \(P=2^{q-1}[(1+C)^q]\) and \(Q=2^{q-1}\left(\frac{B^{1/p}}{\Gamma(r)}\right)^q\) then The conclusion follows from Gronwall's inequality and we conclude that $$\limsup_{t\rightarrow\infty}\left(\ln\frac{t}{a}\right)^{1-n}x(t)< \infty.$$ Similar to the sublinear case, one can easily prove the following results.

Theorem 7. Let \(\lambda=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=\alpha\). i.e., \(n=1\), \(0< \alpha=1-\frac{1}{2p}< 1\), \(\alpha=\gamma\) and \(q=\frac{p}{p-1}\). Clearly, $$p(\alpha-1)+1=p(\gamma-1)+1=p(1-\frac{1}{2p}-1)+1=\frac{1}{2}>0\;\text{and}\;\theta=p(\alpha+\gamma-2)+1=0.$$ Let the functions \(a(t)\) and \(b(t)\) be as in (i) and (ii) with \(b(t)\) be a bounded function and let \(f(t,x)=t^{\gamma-1}h(t)x^\lambda\), where \(0< \lambda< 1\), \(h:\mathbb{R}^+\rightarrow\mathbb{R}^+\) is a continuous function with \(h(t)=m(t)\), $$\int^\infty h^q(s)ds< \infty\;\text{and}\;\limsup_{t\rightarrow\infty}\left(\int_{t_1}^t\left(\ln\frac{t}{s}\right)^{\alpha-1}s^{\gamma-1}h(s)ds\right)< \infty,$$ and let \(e:\mathbb{R}^+\rightarrow\mathbb{R}\) be a continuous function with \(\lim\limits_{t\rightarrow\infty}\int_c^t\left(\ln\frac{t}{s}\right)^{\alpha-1}|e(s)|ds< \infty\). All conditions of Theorem 6 are satisfied and hence every positive solution \(x\) of equation (1) is bounded.

3. Oscillatory and asymptotic behavior of equation

In this section, we present the following result on the asymptotic behavior of nonoscillatory solutions of equation (2).

Theorem 8. Let \(0< \lambda< 1\) and conditions (1), (3) hold and suppose that \(p > 1\), \(q=\frac{p}{p-1}\), \(\gamma=(n-r)+\frac{1}{q}\), \(p(r-1)+1>0\), \(p(\gamma-1)+1>0\), the functions \(b(t)\) is bounded, (7) and (9) hold and

If \(x(t)\) is any nonoscillatory solution of equation (2), then

Proof. Let \(x(t)\) be an eventually positive solution of equation (2). We may assume that for \(t\geq t_1\) for some \(t_1\geq a\), $$x(t)\leq e(t)+\int_a^{t_1}\left(\frac{t}{s}\right)^{r-1}|k(t,s)|f(s,x(s))\frac{ds}{s}+\int_{t_1}^t\left(\frac{t}{s}\right)^{r-1}|k(t,s)|f(s,x(s))\frac{ds}{s},$$

We let \(F(t)=f(t,x(t))\). In view of conditions (1) and (3) we may then write \begin{equation*}\begin{array}{ll} x(t)&\leq\ \ \ \ e(s)+\int_a^{t_1}\left(\ln\frac{t}{s}\right)^{r-1}|F(s)|\frac{ds}{s}\\\\ &~~+b(t)\int_{t_1}^t\left(\frac{t}{s}\right)^{r-1}\left(\ln\frac{s}{a}\right)^{\gamma-1}h(s)[h(s)x^\lambda(s)-m(s)x(s)]\frac{ds}{s}\\\\ &~~+b(t)\int_{t_1}^t\left(\frac{t}{s}\right)^{r-1}\left(\ln\frac{s}{a}\right)^{\gamma-1}m(s)x(s)\frac{ds}{s}. \end{array} \end{equation*} Proceeding exactly as in the proof of Theorem 6, we obtain \begin{equation*}\begin{array}{ll} x(t)&\leq\ \ \ \ e(s)+b(t)\int_{t_1}^t\left(\frac{t}{s}\right)^{r-1}|F(s)|\frac{ds}{s}\\\\ &~~+\left((1-\lambda)\lambda^{\lambda/(1-\lambda)}\right)b(t)\int_{t_1}^t\left(\frac{t}{s}\right)^{r-1}\left(\ln\frac{s}{a}\right)^{\gamma-1} [m^{\lambda/(\lambda-1)}(s)h^{1/(1-\lambda)}(s)]\frac{ds}{s}\\\\ &~~+b(t)\int_{t_1}^t\left(\frac{t}{s}\right)^{r-1}\left(\ln\frac{s}{a}\right)^{\gamma-1}m(s)x(s)\frac{ds}{s}, \end{array} \end{equation*} or, Applying Holder's inequality and Lemma 4, we obtain where \(B=B[p(\gamma-1)+1,p(r-1)+1]\) and \(\theta=p(r+\gamma-2)+1\). Using \(\gamma=(n-r)+\frac{1}{q}\), we have It follows from (18) and (20) that where \(C\) and \(c\) are the upper bound of the functions $$\left(\ln\frac{t}{a}\right)^{1-n}\left[|e(s)|+b(t)\int_a^{t_1}\left(\ln\frac{t}{s}\right)^{r-1}|F(s)|ds+b(t)g(t)\right]\;\text{and}\;B^{1/p}b(t),$$ respectively. The rest of the proof is similar to that of Theorem 6 and conclude that 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 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).

Theorem 9. Let \(0< \lambda< 1 \), conditions (1)-(3) hold and suppose that \(p > 1\), \(q=\frac{p}{p-1}\), \(\gamma=n-r-\frac{1}{q}\), \(p(r-1)+1>0\), and \(p(\gamma-1)+1>0\). In addition we assume that the function \(b(t)\left(\ln\frac{t}{a}\right)^{n-1}\) is bounded and conditions (7), (9) and (15) hold. If

for all \(t_1\geq a\), then equation (2) is oscillatory.

Proof. Let \(x(t)\) be a nonoscillatory solution of equation (2), say \(x(t)>0\) for \(t\geq t_1\) for some \(t_1\geq 0\). The proof when \(x(t)\) is eventually negative is similar. Proceeding as in the proof of Theorem 7 we arrive at (18). Therefore, \begin{equation*}\begin{array}{ll} x(t)&\leq\ \ \ \ e(t)+b(t)\int_a^{t_1}\left(\ln\frac{t}{s}\right)^{r-1}|F(s)|\frac{ds}{s}+b(t)g(t)\\\\ &~~+b(t)\int_{t_1}^t\left(\ln\frac{t}{s}\right)^{r-1}\left(\ln\frac{s}{a}\right)^{\gamma-1}m(s)x(s)\frac{ds}{s}. \end{array} \end{equation*} Using (19) in the above inequality we have \begin{equation*}\begin{array}{ll} x(t)&\leq\ \ \ \ e(t)+b(t)\int_a^{t_1}\left(\ln\frac{t}{s}\right)^{r-1}|F(s)|\frac{ds}{s}+b(t)g(t)+B^{1/p}b(t)\left(\ln\frac{t}{a}\right)^{n-1} \left(\int_{t_1}^tm^q(s)x^q(s)\frac{ds}{s}\right)^{1/q}\\\\ &\leq e(t)+b(t)\int_a^{t_1}\left(\ln\frac{t}{s}\right)^{r-1}|F(s)|\frac{ds}{s}+b(t)g(t)\\\\ &~~+B^{1/p}b(t)\left(\ln\frac{t}{a}\right)^{n-1}\left(\int_{t_1}^t\left(\ln\frac{s}{a}\right)^{(n-1)q}m^q(s)\left(\left(\ln\frac{s}{a}\right) ^{1-n}x(s)\right)^qds\right)^{1/q}. \end{array} \end{equation*} 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 \(\liminf\) in (25) as \(t\rightarrow\infty\) and using condition (22) result in a contradiction with the fact that \(x(t)\) is eventually positive.

Theorem 10. Let \(\lambda=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.

Example 2. Let \(p>1\), \(0< r=2-\frac{1}{p}\in(1,2)\), \(q=\frac{p}{p-1}\) and \(\gamma=2-r+\frac{1}{q}\). Clearly $$p(r-1)+1=p(2-\frac{1}{p}-1)+1=p>1,\;p(\gamma-1)+1=p(2-\gamma+\frac{1}{q}-1)+1>0\;\text{and}\;\theta=p(r+\gamma-2)+1=p.$$ Let \(f(t,x)=\left(\ln\frac{t}{a}\right)^{\gamma-1}h(t)x^\lambda\), where \(0< \lambda< 1\), \(h:\mathbb{R}^+\rightarrow\mathbb{R}^+\) is a continuous function with \(h(t)=m(t)\), $$\int^\infty h^q(s)ds< \infty\;\text{and}\;\limsup_{t\rightarrow\infty}\left(\ln\frac{t}{a}\right)^{-1}\left(\int_{t_1}^t\left(\ln\frac{t}{a}\right)^{r-1} s^{\gamma-1}h(s)ds\right)< \infty,$$ and let \(e:\mathbb{R}^+\rightarrow\mathbb{R}\) be a continuous function with $$\limsup_{t\rightarrow\infty}\left(\ln\frac{t}{a}\right)^{-1}e(t)< \infty\;\text{and}\;\liminf_{t\rightarrow\infty}\left(\ln\frac{t}{a}\right)^{-1}e(t)>-\infty.$$ 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\).

Remark 1. We note that Theorems 8 and 9 are valid when \(\frac{\theta}{p}=n-1\). Here, we study the cases when $$(I)\;\sigma:=\frac{\theta}{p}-n+1>0\;\;\text{or}\;\;(II)\;\sigma< 0.$$ Case (I) If \(\gamma>(n-r)+\frac{1}{q}\), then \(\theta >p(r+\gamma-2)+1\) and so \(\sigma >0\). In this case we need to impose the condition that the function \(b(t)\left(\ln\frac{t}{a}\right)^\sigma\) is bounded and the results are valid.
Case (II) If \(\gamma< (n-r)+\frac{1}{q}\), then \(\theta < p(r+\gamma-2)+1\) and so \(\sigma < 0\). In this case condition (9) may be replaced by $$\int^\infty\left(\ln\frac{s}{a}\right)^{\sigma q}\left(\ln\frac{s}{a}\right)^{(n-1)q}m^q(s)ds< \infty.$$ and the results remain valid. The details are left to the reader.

4. 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 \(\lambda>1\).

5. 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).

Acknowledgments

The authors sincerely thank the reviewers for their valuable suggestions and useful comments that have led to the present improved version of the original manuscript.

Competing Interests

The author(s) do not have any competing interests in the manuscript.