Dynamical and computational analysis of a fractional predator-prey model with an infectious disease and harvesting policy

1.   Introduction and preliminaries

1.1.   Fractional calculus

Fractional calculus has come out as one of the most applicable subjects of mathematics ^[1]. Its importance is evident from the fact that many real-world phenomena can be best interpreted and modeled using this theory. It is also a fact that many disciplines of engineering and science have been influenced by the tools and techniques of fractional calculus. Its emergence can easily be traced and linked with the famous correspondence between the two mathematicians, L'Hospital and Leibnitz, which was made on 30th September 1695. After that, many researchers tried to explore the concept of fractional calculus, which is based on the generalization of nth order derivatives or n-fold integration ^[2,3,4].

Recently, Khan and Khan ^[5] have discovered novel definitions of fractional integral and derivative operators. These operators enjoy interesting properties such as continuity, boundedeness, linearity etc. The integral operators, they presented, are stated as under:

Definition 1 (^[5]). Let $h \in L_{\theta}[s, t]$(conformable integrable on $[s, t]\subseteq[0, \infty)$). The left-sided and right-sided generalized conformable fractional integrals ${}^{\tau}_{\theta}K_{s^{+}}^{\nu}$ and ${}^{\tau}_{\theta}K_{t^{-}}^{\nu}$ of order $\nu > 0$ with $\theta\in(0, 1]$, $\tau\in\mathbb{R}$, $\theta+\tau\neq0$ are defined by:

and

respectively, and ${}^{\tau}_{\theta}K_{s^{+}}^{0}h(r) = {}^{\tau}_{\theta}K_{t^{-}}^{0}h(r) = h(r)$. Here $\Gamma$ denotes the well-known Gamma function.

Here the integral $\int\limits_{s}^{t}d_{\theta}w$ represents the conformable integration, defined as:

The operators defined in Definition 1 are in generalized form and contain few important operators in themselves. Here, only the left-sided operators are presented, the corresponding right-sided operators may be deduced in the similar way. Moreover, to understand the theory of conformable fractional calculus, one can see ^[5,6,7]. Also, the basic theory of fractional calculus can be found in the books ^[1,8,9] and for the latest research in this field one can see ^{[3,4,10,11,12]} and the references there in.

Remark 1. 1) For $\theta = 1$ in the Definition 1, the following Katugampula fractional integral operator is obtained ^[13]:

2) For $\tau = 0$ in the Definition 1, the New Riemann Liouville type conformable fractional integral operator is obtained as given below:

3) Using the definition of conformable integral given in (1.3) and L'Hospital rule, it is straightforward that when $\theta\rightarrow 0$ in (1.5), we get the Hadamard fractional integral operator as follows:

4) For $\theta = 1$ in (1.5), the well-known Riemann-Liouville fractional integral operator is obtained as follows:

5) For the case $\nu = 1, \tau = 0$ in Definition 1, we get the conformable fractional integrals. And when $\theta = \nu = 1$, $\tau = 0$, we get the classical Riemann integrals.

1.2.   Hermite-Hadamard inequalities

This subsection is devoted to start with the definition of convex function, which plays a very important role in establishment of various kinds of inequalities ^[14]. This definition is given as follows ^[15]:

Definition 2. A function $h: I\subseteq \mathbb{R} \rightarrow \mathbb{R}$ is said to be convex on $I$ if the inequality

holds for all $s, t\in I$ and $0\leq\eta\leq 1$. The function $h$ is said to be concave on $I$ if the inequality given in (1.8) holds in the reverse direction.

Associated with the Definition 2 of convex functions the following double inequality is well-known and it has been playing a key role in various fields of science and engineering ^[15].

Theorem 1. Let $h: I\subseteq \mathbb{R} \rightarrow \mathbb{R}$ be a convex function and $s, t\in I$ with $s < t$. Then we have the following Hermite-Hadamard inequality:

This inequality (1.9) appears in a reversed order if the function $h$ is supposed to be concave. Also, the relation (1.9) provides upper and lower estimates for the integral mean of the convex function $h$. The inequality (1.9) has various versions (extensions or generalizations) corresponding to different integral operators ^{[16,17,18,19,20,21,22,23,24,25]} each version has further forms with respect to various kinds of convexities ^{[26,27,28,29,30,31,32]} or with respect to different bounds obtained for the absolute difference of the two leftmost or rightmost terms in the Hermite-Hadamard inequality.

By using the Riemann-Liouville fractional integral operators, Sirikaye et al. have proved the following Hermite-Hadamard inequality ^[33].

Theorem 2. (^[33]). Let $h:[s, t]\rightarrow\mathbb{R}$ be a function such that $0\leq s < t$ and $h\in L[s, t]$. If $h$ is convex on $[s, t]$, then the following double inequality holds:

For more recent research related to generalized Hermite-Hadamard inequality one can see ^{[34,35,36,37,38,39,40,41,42]} and the references therein.

Motivated from the Riemann-Liouville version of Hermite-Hadamard inequality (given above in (1.10)), we prove the same inequality for newly introduced generalized conformable fractional operators. As a result we get a more generalized inequality, containing different versions of Hermite-Hadamard inequality in single form. We also prove an identity for generalized conformable fractional operators and establish a bound for the absolute difference of two rightmost terms in the newly obtained Hermite-Hadamard inequality. We point out some relations of our results with those of other results from the past. At the end we present conclusion, where directions for future research are also mentioned.

2.   Main results

In the following theorem the well-known Hermite-Hadamard inequality for the newly defined integral operators is proved.

Theorem 3. Let $\nu > 0$ and $\tau \in \mathbb{R}, \theta\in (0, 1]$ such that $\tau+\theta > 0$. Let $h : [s, t] \subseteq [0, \infty)\rightarrow \mathbb{R}$ be a function such that $h \in L_{\theta}[s, t]$(conformal integrable on [s, t]). If $h$ is also a convex function on $[s, t]$, then the following Hermite-Hadamard inequality for generalized conformable fractional Integrals ${}^{\tau}_{\theta}K_{s^{+}}^{\nu}$ and ${}^{\tau}_{\theta}K_{t^{-}}^{\nu}$ holds:

where $H(x) = h(x)+\tilde{h}(x)$, $\tilde{h}(x) = h(s+t-x)$.

Proof. Let $\eta\in[0, 1]$. Consider $x, y\in[s, t]$, defined by $x = \eta s + (1 - \eta)t, y = (1 - \eta)s + \eta t$. Since $h$ is a convex function on $[s, t]$, we have

Multiplying both sides of (2.2) by

and integrating with respect to $\eta$, we get

Note that we have

Also, by using the identity $\tilde{h}((1 - \eta)s + \eta t) = h(\eta s+(1-\eta)t),$ and making substitution $(1-\eta)s+\eta t = w$, we get

Similarly

By substituting these values in (2.3), we get

Again, by multiplying both sides of (2.2) by

and then integrating with respect to $\eta$ and by using the same techniques used above, we can obtain:

Adding (2.7) and (2.6), we get:

Hence the left-hand side of the inequality (2.1) is established.

Also since $h$ is convex, we have:

Multiplying both sides

and integrating with respect to $\eta$ we get

that is,

Similarly multiplying both sides of (2.9) by

and integrating with respect to $\eta$, we can obtain

Adding the inequalities (2.11) and (2.12), we get:

Combining (2.8) and (2.13), we get the required result.

The inequality in (2.1) is in compact form containing few inequalities for different integrals in it. The following remark tells us about that fact.

Remark 2. 1) For $\theta = 1$ in (2.1), we get Hermite-Hadamard inequality for Katugampola fractional integral operators, as follows ^[38]:

where $H(x) = h(x)+\tilde{h}(x)$, $\tilde{h}(x) = h(s+t-x)$.

2) For $\tau = 0$ in (2.1), we get Hermite-Hadamard inequality for newly obtained Riemann Liouville type conformable fractional integral operators, as follows:

where $H(x) = h(x)+\tilde{h}(x)$, $\tilde{h}(x) = h(s+t-x)$.

3) For $\tau+\theta\rightarrow 0$, in (2.1), applying L'Hospital rule and the relation (1.3), we get Hermite-Hadamard inequality for Hadamard fractional integral operators, as follows:

4) For $\tau+\theta = 1$ in (2.1), the Hermite-Hadamard inequality is obtained for Riemann-Liouville fractional integrals ^[33]:

5) For the case $\nu = 1, \tau = 0$ in (2.1), the Hermite-Hadamard inequality is obtained for the conformable fractional integrals as follows:

6) When $\theta = \nu = 1$, $\tau = 0$ the Hermite-Hadamard inequality is obtained for classical Riemann integrals ^[15]:

To bound the difference of two rightmost terms in the main inequality (2.1), we need to establish the following Lemma.

Lemma 1. Let $\tau+\theta > 0$ and $\nu > 0.$ If $h \in L_{\theta}[s, t]$, then

where

Proof. With the help of integration by parts, we have

Similarly, we have

Using (2.21) and (2.22) we have

Also, we have

And

Also, we have

Using (2.23), (2.25) and (2.26) we get the required result.

Remark 3. When $\tau+\theta = 1$ in Lemma 1, we get the Lemma 2 in ^[33].

Definition 3. For $\nu > 0$, we define the operators

and

where $x, y\in[s, t]\subseteq [0, \infty)$ and $\tau+\theta > 0$.

Theorem 4. Let $h$ be a conformable integrable function over $[s, t]$ such that $|h'|$ is convex function. Then for $\nu > 0$ and $\tau+\theta > 0$ we have:

where $K_{\tau+\theta}^{\nu}(s, t) = \Omega_{1}^{\nu}(t, t, \tau+\theta)+\Omega_{2}^{\nu}(s, s, \tau+\theta)-\Omega_{2}^{\nu}(t, s, \tau+\theta)-\Omega_{1}^{\nu}(s, t, \tau+\theta)$.

Proof. Using Lemma 1 and convexity of $|h'|$, we have:

Here $\int\limits^{1}_{0}\eta\left|\Delta_{\tau+\theta}^{\nu}(\eta)\right|d\eta = \frac{1}{(t-s)^{2}}\int\limits^{t}_{s}|\psi(u)|(t-u)du$,

and $\psi(u) = (u^{\tau+\theta}-s^{\tau+\theta})^{\nu}-((t+s-u)^{\tau+\theta}-s^{\tau+\theta})^{\nu}+(t^{\tau+\theta}-(s+t-u)^{\tau+\theta})^{\nu} -(t^{\tau+\theta}-u^{\tau+\theta})^{\nu}$.

We observe that $\psi$ is a nondecreasing function on $[s, t]$. Moreover, we have:

and also $\psi\left(\frac{s+t}{2}\right) = 0$. As a consequence, we have

Thus we get

where

and

We can see here that $K_{1} = -\Omega_{2}^{\nu}(t, s, \tau+\theta)$, $K_{4} = \Omega_{1}^{\nu}(t, t, \tau+\theta)$.

Also, by using of change of the variables $v = s+t-u$, we get

By substituting these values in (2.31), we get

Similarly, we can find

Finally, by using (2.30), (2.37) and (2.38) we get the required result.

Remark 4. when $\tau+\theta = 1$ in (2.29), we obtain

which is Theorem 3 in ^[33].

3.   Conclusion and future works

A generalized version of Hermite-Hadamard inequality via newly introduced GC fractional operators has been acquired successfully. This result combines several versions (new and old) of the Hermite-Hadamard inequality into a single form, each one has been discussed by fixing parameters in the newly established version of the Hermite-Hadamard inequality. Moreover, an identity containing the GC fractional integral operators has been proved. By using this identity, a bound for the absolute of the difference between the two rightmost terms in the newly established Hermite-Hadamard inequality has been presented. Also, some relations of our results with those of already existing results have been pointed out. Since this is a fact that there exist more than one definitions for fractional derivatives ^[2] which makes it difficult to choose a convenient operator for solving a given problem. Thus, in the present paper, the GC fractional operators (containing various previously defined fractional operators into a single form) have been used in order to overcome the problem of choosing a suitable fractional operator and to provide a unique platform for researchers working with different operators in this field. Also, by making use of GC fractional operators one can follow the research work which has been performed for the two versions (1.9) and (1.10) of Hermite-Hadamard inequality.

Acknowledgments

This work was supported by the Natural Science Foundation of China (Grant Nos. 61673169, 11301127, 11701176, 11626101, 11601485).

Conflict of interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.