A study on the fractal-fractional tobacco smoking model

Hasib Khan; Jehad Alzabut; Anwar Shah; Sina Etemad; Shahram Rezapour; Choonkil Park; Hasib Khan; Jehad Alzabut; Anwar Shah; Sina Etemad; Shahram Rezapour; Choonkil Park

doi:10.3934/math.2022767

AIMS Mathematics

2022, Volume 7, Issue 8: 13887-13909. doi: 10.3934/math.2022767

Previous Article Next Article

Research article

A study on the fractal-fractional tobacco smoking model

1.
Department of Mathematics, Shaheed Benazir Bhutto Uniersity, Sheringal, Dir Upper, Khyber Pakhtunkhwa, Pakistan
2.
Department of Mathematics and Sciences, Prince Sultan University, 11586 Riyadh, Saudi Arabia
3.
Department of Industrial Engineering, OSTİM Technical University, 06374 Ankara, Türkiye
4.
Department of Mathematics, University of Malakand, Chakdara, Dir Lower, Khyber Pakhtunkhwa, Pakistan
5.
Department of Mathematics, Azarbaijan Shahid Madani University, Tabriz, Iran
6.
Department of Medical Research, China Medical University Hospital, China Medical University, Taichung, Taiwan
7.
Research Institute for Natural Sciences, Hanyang University, Seoul 04763, Korea

Received: 30 March 2022 Revised: 01 May 2022 Accepted: 08 May 2022 Published: 25 May 2022
MSC : 26A33, 34A08, 35R11

In this article, we consider a fractal-fractional tobacco mathematical model with generalized kernels of Mittag-Leffler functions for qualitative and numerical studies. From qualitative point of view, our study includes; existence criteria, uniqueness of solution and Hyers-Ulam stability. For the numerical aspect, we utilize Lagrange's interpolation polynomial and obtain a numerical scheme which is further illustrated simulations. Lastly, a comparative analysis is presented for different fractal and fractional orders. The numerical results are divided into four figures based on different fractal and fractional orders. We have found that the fractional and fractal orders have a significant impact on the dynamical behaviour of the model.

Keywords:

Citation: Hasib Khan, Jehad Alzabut, Anwar Shah, Sina Etemad, Shahram Rezapour, Choonkil Park. A study on the fractal-fractional tobacco smoking model[J]. AIMS Mathematics, 2022, 7(8): 13887-13909. doi: 10.3934/math.2022767

Related Papers:

[1]	Bandar Bin-Mohsin, Muhammad Uzair Awan, Muhammad Zakria Javed, Artion Kashuri, Muhammad Aslam Noor . Fractional integral estimations pertaining to generalized ${\gamma}$ -convex functions involving Raina's function and applications. AIMS Mathematics, 2022, 7(8): 13633-13663. doi: 10.3934/math.2022752
[2]	Thongchai Botmart, Soubhagya Kumar Sahoo, Bibhakar Kodamasingh, Muhammad Amer Latif, Fahd Jarad, Artion Kashuri . Certain midpoint-type Fejér and Hermite-Hadamard inclusions involving fractional integrals with an exponential function in kernel. AIMS Mathematics, 2023, 8(3): 5616-5638. doi: 10.3934/math.2023283
[3]	Shahid Mubeen, Rana Safdar Ali, Iqra Nayab, Gauhar Rahman, Kottakkaran Sooppy Nisar, Dumitru Baleanu . Some generalized fractional integral inequalities with nonsingular function as a kernel. AIMS Mathematics, 2021, 6(4): 3352-3377. doi: 10.3934/math.2021201
[4]	Yu-Pei Lv, Ghulam Farid, Hafsa Yasmeen, Waqas Nazeer, Chahn Yong Jung . Generalization of some fractional versions of Hadamard inequalities via exponentially $(\alpha, h-m)$ -convex functions. AIMS Mathematics, 2021, 6(8): 8978-8999. doi: 10.3934/math.2021521
[5]	Muhammad Imran Asjad, Waqas Ali Faridi, Mohammed M. Al-Shomrani, Abdullahi Yusuf . The generalization of Hermite-Hadamard type Inequality with exp-convexity involving non-singular fractional operator. AIMS Mathematics, 2022, 7(4): 7040-7055. doi: 10.3934/math.2022392
[6]	Shuhong Yu, Tingsong Du . Certain inequalities in frame of the left-sided fractional integral operators having exponential kernels. AIMS Mathematics, 2022, 7(3): 4094-4114. doi: 10.3934/math.2022226
[7]	Hu Ge-JiLe, Saima Rashid, Muhammad Aslam Noor, Arshiya Suhail, Yu-Ming Chu . Some unified bounds for exponentially $tgs$ -convex functions governed by conformable fractional operators. AIMS Mathematics, 2020, 5(6): 6108-6123. doi: 10.3934/math.2020392
[8]	Miguel Vivas-Cortez, Muhammad Aamir Ali, Artion Kashuri, Hüseyin Budak . Generalizations of fractional Hermite-Hadamard-Mercer like inequalities for convex functions. AIMS Mathematics, 2021, 6(9): 9397-9421. doi: 10.3934/math.2021546
[9]	Muhammad Tariq, Asif Ali Shaikh, Sotiris K. Ntouyas, Jessada Tariboon . Some novel refinements of Hermite-Hadamard and Pachpatte type integral inequalities involving a generalized preinvex function pertaining to Caputo-Fabrizio fractional integral operator. AIMS Mathematics, 2023, 8(11): 25572-25610. doi: 10.3934/math.20231306
[10]	Artion Kashuri, Rozana Liko, Silvestru Sever Dragomir . New generalized integral inequalities with applications. AIMS Mathematics, 2019, 4(3): 984-996. doi: 10.3934/math.2019.3.984

Abstract

1. Introduction and preliminaries

Fractional calculus has emerged as one of the most important interdisciplinary subjects. In recent past it experienced rapid development and consequently several new generalizations of classical concepts of fractional calculus have been obtained in the literature, for example, see ^[9].

The classical Riemann-Liouville fractional integrals are defined as:

Definition 1.1 (^[9]). Let ${\mathcal{F}}\in L_{1}[a, b]$ . Then the Riemann-Liouville integrals $J_{a^+}^{\alpha}{\mathcal{F}}$ and $J_{b^-}^{\alpha}{\mathcal{F}}$ of order $\alpha > 0$ with $a\geq 0$ are defined by

$\begin{equation*} J_{a^+}^{\alpha}{\mathcal{F}}(x) = \frac{1}{\Gamma(\alpha)}\int\limits_a^x(x-v)^{\alpha-1}{\mathcal{F}}(v){\mathrm{d}}v,\quad x > a, \end{equation*}$

and

$\begin{equation*} J_{b^-}^{\alpha}{\mathcal{F}}(x) = \frac{1}{\Gamma(\alpha)}\int\limits_x^b(v-x)^{\alpha-1}{\mathcal{F}}(v){\mathrm{d}}v,\quad x < b, \end{equation*}$

where

$\begin{align*} \Gamma(x) = \int_0^\infty e^{-v}v^{x-1}{\mathrm{d}}v, \end{align*}$

is the well known Gamma function.

Diaz et al. ^[8] introduced the notion of generalized $k$ -gamma function. The integral form of $\Gamma_{k}$ is given by:

$\begin{equation*} \Gamma_{k}(x) = \int\limits_{0}^{\infty}v^{x-1}e^{-\frac{v^{k}}{k}}{\mathrm{d}}v,\quad \Re(x) > 0. \end{equation*}$

Note that

$\Gamma_{k}(x) = k^{\frac{x}{k}-1}\Gamma\left(\frac{x}{k}\right).$

$k$ -Beta function is defined as:

$\begin{equation*} \beta_{k}(x,y) = \frac{1}{k}\int\limits_{0}^{1}v^{\frac{x}{k}-1}(1-v)^{\frac{x}{k}-1}{\mathrm{d}}v. \end{equation*}$

Obviously

$\beta_{k}(x,y) = \frac{1}{k}\beta\left(\frac{x}{k},\frac{y}{k}\right).$

Sarikaya et al. ^[18] extended the notion of Riemann-Liouville fractional integrals to $k$ -Riemann-Liouville fractional integrals and discussed some of its interesting properties.

To be more precise let ${\mathcal{F}}$ be piecewise continuous on $I^{*} = (0, \infty)$ and integrable on any finite subinterval of $I = [0, \infty]$ . Then for $v > 0$ , we consider $k$ -Riemann-Liouville fractional integral of ${\mathcal{F}}$ of order $\alpha$

$\begin{align*} _kJ_{a}^{\alpha}{\mathcal{F}}(x) = \frac{1}{k\Gamma_{k}(\alpha)}\int\limits_{a}^{x}(x-v)^{\frac{\alpha}{k}-1}{\mathcal{F}}(v){\mathrm{d}}v,\quad x > a,k > 0. \end{align*}$

It has been observed that $k$ -fractional integrals are significant generalizations of classical fractional integrals. For more details, see ^[18].

Ahmad et al. ^[1] defined fractional integral operators with an exponential kernel and obtained corresponding inequalities.

Definition 1.2. Let $\mathcal{F}\in [a, b]$ . The fractional left side integral $\, _{k}\mathcal{I}^{\alpha}_{a^{+}}\mathcal{F}$ and right side integral $\, _{k}\mathcal{I}^{\alpha}_{b^{-}}\mathcal{F}$ of order $\alpha\in(0, 1)$ are defined as follows:

$\mathcal{I}^{\alpha}_{a^{+}}\mathcal{F}(x) = \frac{1}{\alpha}\int\limits_{a}^{x}e^{-\frac{1-\alpha}{\alpha}(x-v)}\mathcal{F}(v){\mathrm{d}}v,\ \ \ \ x > a,$

and

$\mathcal{I}^{\alpha}_{b^{-}}\mathcal{F}(x) = \frac{1}{\alpha}\int\limits_{x}^{b}e^{-\frac{1-\alpha}{\alpha}(v-x)}\mathcal{F}(v){\mathrm{d}}v,\ \ \ \ x < b.$

Using the ideas of ^[1,18], we now introduce the notion of $k$ -fractional integral operators with an exponential kernel.

Definition 1.3. Let $\mathcal{F}\in L[a, b]$ . The $k$ -fractional left side integral $\, _{k}\mathcal{I}^{\alpha}_{a^{+}}\mathcal{F}$ and right side integral $\, _{k}\mathcal{I}^{\alpha}_{b^{-}}\mathcal{F}$ of order $\alpha\in(0, 1)$ for $k > 0$ are defined as follows

$_{k}\mathcal{I}^{\alpha}_{a^{+}}\mathcal{F}(x) = \frac{k}{\alpha}\int\limits_{a}^{x}e^{-\frac{k-\alpha}{\alpha}(x-v)}\mathcal{F}(v){\mathrm{d}}v,\ \ \ \ x > a,$

and

$_{k}\mathcal{I}^{\alpha}_{b^{-}}\mathcal{F}(x) = \frac{k}{\alpha}\int\limits_{x}^{b}e^{-\frac{k-\alpha}{\alpha}(v-x)}\mathcal{F}(v){\mathrm{d}}v,\ \ \ \ x < b.$

It is to be noted that by taking $k \to 1$ in Definition 1.3, we recapture Definition 1.2. Fractional analogues of integral inequalities have a great many applications in numerical quadrature, transform theory, probability, statistical problems etc. Therefore, a significant and rapid development in this field has been noticed, for details, see ^{[2,3,20,24,25]}. Sarikaya et al. ^[19] utilized the concepts of fractional integrals and obtained new fractional refinements of trapezium like inequalities. This article motivated many researchers and as a result several new fractional extensions of classical inequalities have been obtained in the literature, for example, see ^{[1,4,6,7,11,14,15,16,17,18,19,22,23]}. Recently Ahmad et al. ^[1] used fractional integral operators with an exponential kernel and obtained corresponding inequalities. Wu et al. ^[23] derived some new identities and bounds pertaining to fractional integrals with the exponential kernel.

The main motivation of this paper is to derive some new fractional refinements of trapezium like inequalities essentially using the new fractional integral operators with an exponential kernel to $k$ -fractional integral operators with an exponential kernel and the preinvexity property of the functions. In order to establish the significance of our main results, we offer some applications of our main results to means and $q$ -digamma functions. We hope that the ideas and techniques of this paper will inspire interested readers working in the field of inequalities.

Before we proceed further, we now recall some previously known concepts from convex analysis. We first, start with the definition of invex sets.

Definition 1.4 (^[10]). A set $K$ is said to be invex with respect to bifunction $\theta(., .)$ , if

$x+v\theta(y,x)\in K,\,\forall x,y\in K,\,v\in[0,1].$

The preinvexity of the functions is defined as:

Definition 1.5 (^[21]). A function $\mathcal{F}:K\to\mathbb{R}$ is said to be preinvex with respect to bifunction $\theta(., .)$ , if

${\mathcal{F}}(x+v\theta(y,x))\leq(1-v){\mathcal{F}}(x)+vf(y),\,\forall x,y\in K, v\in[0,1].$

In order to obtain some of the main results of the paper, we need the famous condition C, which was introduced by Mohan and Neogy ^[13]. This condition played a vital role in the development of several results involving preinvex functions.

Condition C. Let $\theta:K\times K\to\mathbb{R}^{n}$ . We say that the bifunction $\theta(., .)$ satisfies the condition C, if for any $x, y\in\mathbb{R}^{n}$

1. $\theta(x, x+v\theta(y, x)) = -v\theta(y, x)$ ,

2. $\theta(y, x+v\theta(y, x)) = (1-v)\theta(y, x)$ ,

for all $v\in[0, 1]$ .

Note that for any $x, y\in\mathbb{R}^{n}$ and $v_{1}, v_{2}\in[0, 1]$ and from the condition C, we have, see ^[12]

$\theta(x+v_{2}\theta(y,x),x+v_{1}\theta(y,x)) = (v_{2}-v_{1})\theta(y,x).$

2. Trapezium type inequalities

In this section, we derive some new fractional trapezium type inequalities involving the functions having preinvexity property. For the sake of simplicity, we set $\rho = \frac{k-\alpha}{\alpha}\theta(b, a)$ and $\rho_1 = \frac{1-\alpha}{\alpha}\theta(b, a)$ .

Theorem 2.1. Let $\mathcal{F}:[a, a+\theta(b, a)]\subseteq\mathbb{R}\to\mathbb{R}$ be a positive function with $\theta(b, a) > 0$ and $\mathcal{F}\in L[a, a+\theta(b, a)]$ . Suppose $\mathcal{F}$ is a preinvex function and $\theta(., .)$ satisfies condition C, then

$\begin{align} \mathcal{F}\left(\frac{2a+\theta(b,a)}{2}\right)\leq\frac{k-\alpha}{2k(1-e^{-\rho})}[\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{F}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{F}(a)]\leq\frac{\mathcal{F}(a)+\mathcal{F}(b)}{2}. \end{align}$

(2.1)

Proof. By preinvexity of $\mathcal{F}$ , we have for every $x, y\in [a, a+\theta(b, a)]$ with $\lambda = \frac{1}{2}$

$\begin{align*} 2\mathcal{F}\left(x+\frac{\theta(y,x)}{2}\right)\leq[\mathcal{F}(x)+\mathcal{F}(y)], \end{align*}$

with $x = a+v\theta(b, a), y = a+(1-v)\theta(b, a)$ and using the condition C, we have

$\begin{align} &2\mathcal{F}\left(a+v\theta(b,a)+\frac{\theta(a+(1-v)\theta(b,a),a+v\theta(b,a))}{2}\right)\\ & = 2\mathcal{F}\left(a+v\theta(b,a)+\frac{(1-2v)\theta(b,a)}{2}\right)\\ & = 2\mathcal{F}\left(\frac{2a+\theta(b,a)}{2}\right)\\ &\leq\mathcal{F}(a+v\theta(b,a))+\mathcal{F}(a+(1-v)\theta(b,a)). \end{align}$

(2.2)

Multiplying both sides of above inequality by $e^{-\rho v}$ and integrating with respect to $v$ over $[0, 1]$ , we have

$\begin{align*} \frac{2(1-e^{-\rho})}{\rho}\mathcal{F}\left(\frac{2a+\theta(b,a)}{2}\right)&\leq\int\limits^{1}_{0}e^{-\rho v}f(a+v\theta(b,a)){\mathrm{d}}v+\int\limits^{1}_{0}e^{-\rho v}\mathcal{F}(a+(1-v)\theta(b,a)){\mathrm{d}}v\\ & = \frac{1}{\theta(b,a)}\left[\int\limits^{a+\theta(b,a)}_{a}e^{-\rho(\frac{s-a}{\theta(b,a)})}\mathcal{F}(s){\mathrm{d}}s+\int\limits^{a+\theta(b,a)}_{a}e^{-\rho(\frac{a+\theta(b,a)-s}{\theta(b,a)})}\mathcal{F}(s){\mathrm{d}}s\right]\\ & = \frac{\alpha}{k\theta(b,a)}\left[\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{F}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{F}(a)\right]. \end{align*}$

As a result, we get

(2.3)

For the proof of second inequality, we note that $\mathcal{F}$ is a preinvex function, so we have

$\begin{align} \mathcal{F}(a+v\theta(b,a))+\mathcal{F}(a+(1-v)\theta(b,a))\leq\mathcal{F}(a)+\mathcal{F}(b). \end{align}$

(2.4)

Multiplying both sides by $e^{-\rho v}$ and integrating with respect to $v$ over $[0, 1]$ , we have

$\begin{align} \frac{\alpha}{k\theta(b,a)}\left[\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{F}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{F}(a)\right]\leq\frac{1-e^{-\rho}}{\rho}[\mathcal{F}(a)+\mathcal{F}(b)], \end{align}$

(2.5)

Combining (2.3) and (2.5) completes the proof.

Theorem 2.2. Let $\mathcal{F}:[a, a+\theta(b, a)]\subseteq\mathbb{R}\to\mathbb{R}$ be a positive and preinvex function with $\theta(b, a) > 0$ and $\mathcal{F}\in L[a, a+\theta(b, a)]$ . Let $W$ be a non-negative, integrable and symmetric with respect to $\frac{2a+\theta(b, a)}{2}$ , then using the condition C, we have

$\begin{align} &\mathcal{F}\left(\frac{2a+\theta(b,a)}{2}\right)[\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{W}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{W}(a)]\\ &\leq\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{(FW)}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{(FW)}(a)\\ &\leq\frac{\mathcal{F}(a)+\mathcal{F}(b)}{2}[\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{W}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{W}(a)]. \end{align}$

(2.6)

Proof. Since $\mathcal{F}$ is preinvex function on $L[a, a+\theta(b, a)]$ , so multiplying inequality (2.2) by

$\begin{align} e^{-\rho v}\mathcal{W}(a+v\theta(b,a)), \end{align}$

(2.7)

and then integrating with respect to $v$ over $[0, 1]$ , we get

$\begin{align*} &2\mathcal{F}\left(\frac{2a+\theta(b,a)}{2}\right)\int\limits_{0}^{1}e^{-\rho v}\mathcal{W}(a+v\theta(b,a)){\mathrm{d}}v\\ &\leq\int\limits^{1}_{0}e^{-\rho v}\mathcal{W}(a+v\theta(b,a))\mathcal{F}(a+v\theta(b,a)){\mathrm{d}}v +\int\limits^{1}_{0}e^{-\rho v}\mathcal{W}(a+v\theta(b,a))\mathcal{F}(a+(1-v)\theta(b,a)){\mathrm{d}}v\\ & = \int\limits^{1}_{0}e^{-\rho v}\mathcal{W}(a+v\theta(b,a))\mathcal{F}(a+v\theta(b,a)){\mathrm{d}}v+\int\limits^{1}_{0}e^{-\rho v}\mathcal{W}(a+(1-v)\theta(b,a))\mathcal{F}(a+(1-v)\theta(b,a)){\mathrm{d}}v\\ & = \frac{1}{\theta(b,a)}\left[\int\limits^{a+\theta(b,a)}_{a}e^{-\rho(\frac{s-a}{\theta(b,a)})}\mathcal{F}(s)\mathcal{W}(s){\mathrm{d}}s+\int\limits^{a+\theta(b,a)}_{a}e^{-\rho(\frac{a+\theta(b,a)-s}{\theta(b,a)})}\mathcal{F}(s)\mathcal{W}(s){\mathrm{d}}s\right]\\ & = \frac{\alpha}{k\theta{(b,a)}}[\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{(FW)}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{(FW)}(a)]. \end{align*}$

Thus

$\begin{align*} &2\mathcal{F}\left(\frac{2a+\theta(b,a)}{2}\right)\int\limits_{0}^{1}e^{-\rho v}\mathcal{W}(a+v\theta(b,a)){\mathrm{d}}v\\ &\leq\frac{\alpha}{k\theta{(b,a)}}[\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{FW}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{FW}(a)]. \end{align*}$

Since $\mathcal{W}$ is symmetric with respect to $\frac{2a+\theta(b, a)}{2},$ we have

$\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{W}(a+\theta(b,a)) = \,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{W}(a) = \frac{1}{2}[\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{W}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{W}(a)].$

Thus we get the left side of inequality (2.6).

For the proof of right side of inequality (2.6), we multiply (2.7) and (2.4) and then integrating the resulting inequality with respect to $v$ over $[0, 1],$ we get the required result.

Theorem 2.3. Let $\mathcal{F, W}:[a, a+\theta(b, a)]\subseteq\mathbb{R}\to\mathbb{R}$ be nonnegative and preinvex function on $L[a, a+\theta(b, a)]$ with $\theta(b, a) > 0$ . If $\theta(., .)$ satisfies condition C, then the following inequalities for the $k$ -fractional integrals with exponential kernel holds:

$\begin{align} &\frac{\alpha}{2k\theta(b,a)}[\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{(FW)}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{(FW)}(a)]\\ &\leq M(a,b)\frac{\rho^{2}-2\rho+4-e^{-\rho}(\rho^{2}+2\rho+4)}{2\rho^{3}}+N(a,b)\frac{\rho-2+e^{-\rho}(\rho+2)}{\rho^{3}}, \end{align}$

(2.8)

and

$\begin{align} &2\mathcal{F}\left(\frac{2a+\theta(b,a)}{2}\right)\mathcal{W}\left(\frac{2a+\theta(b,a)}{2}\right)\\ &\leq\frac{k-\alpha}{2k(1-e^{-\rho})}[\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{FW}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{FW}(a)]\\ &\leq M(a,b)\frac{\rho-2+e^{-\rho}(\rho+2)}{\rho^{2}(1-e^{-\rho})}+N(a,b)\frac{\rho^{2}-2\rho+4-e^{-\rho}(\rho^{2}+2\rho+4)}{2\rho^{2}(1-e^{-\rho})}, \end{align}$

(2.9)

where

$M(a,b) = \mathcal{F}(a)\mathcal{W}(a)+\mathcal{F}(b)\mathcal{W}(b),$

and

$N(a,b) = \mathcal{F}(a)\mathcal{W}(b)+\mathcal{F}(b)\mathcal{W}(a).$

Proof. Since $\mathcal{F}, \mathcal{W}$ are preinvex functions on $[a, a+\theta(b, a)]$ , then we have

$\begin{align*} \mathcal{F}(a+v\theta(b,a))\mathcal{W}(a+v\theta(b,a))\leq (1-v)^{2}\mathcal{F}(a)\mathcal{W}(a)+v^{2}\mathcal{F}(b)\mathcal{W}(b)+v(1-v)N(a,b), \end{align*}$

and

$\begin{align*} \mathcal{F}(a+(1-v)\theta(b,a))\mathcal{W}(a+(1-v)\theta(b,a))\leq v^{2}\mathcal{F}(a)\mathcal{W}(a)+(1-v)^{2}\mathcal{F}(b)\mathcal{W}(b)+v(1-v)N(a,b). \end{align*}$

Adding above inequalities, we have

$\begin{align*} \mathcal{F}(a+v\theta(b,a))\mathcal{W}(a+v\theta(b,a))&+\mathcal{F}(a+(1-v)\theta(b,a))\mathcal{W}(a+(1-v)\theta(b,a))\\ &\leq(2v^{2}-2v+1)M(a,b)+2v(1-v)N(a,b). \end{align*}$

Multiplying both sides of above inequality by $e^{-\rho v}$ and integrating with respect to $v$ over $[0, 1]$ , we have

$\begin{align*} &\int\limits_{0}^{1}e^{-\rho v}\mathcal{F}(a+v\theta(b,a))\mathcal{W}(a+v\theta(b,a)){\mathrm{d}}v\\ &\quad+\int\limits_{0}^{1}e^{-\rho v}\mathcal{F}(a+(1-v)\theta(b,a))\mathcal{W}(a+(1-v)\theta(b,a)){\mathrm{d}}v\\ &\leq M(a,b)\int\limits_{0}^{1}e^{-\rho v}(2v^{2}-2v+1){\mathrm{d}}v+M(a,b)\int\limits_{0}^{1}e^{-\rho v}2v(1-v){\mathrm{d}}v\\ & = M(a,b)\frac{\rho^{2}-2\rho+4-e^{-\rho}(\rho^{2}+2\rho+4)}{2\rho^{3}}+N(a,b)\frac{\rho-2+e^{-\rho}(\rho+2)}{\rho^{3}}. \end{align*}$

So,

$\begin{align*} &\frac{\alpha}{2k\theta(b,a)}[\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{(FW)}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{(FW)}(a)]\\ &\leq M(a,b)\frac{\rho^{2}-2\rho+4-e^{-\rho}(\rho^{2}+2\rho+4)}{2\rho^{3}}+N(a,b)\frac{\rho-2+e^{-\rho}(\rho+2)}{\rho^{3}}. \end{align*}$

For the proof of inequality (2.9), using the preinvexity of $\mathcal{F}, \mathcal{W}$ and condition C, we have

$\begin{align} &\mathcal{F}\left(\frac{2a+\theta(b,a)}{2}\right)\mathcal{W}\left(\frac{2a+\theta(b,a)}{2}\right)\\ & = \mathcal{F}\left(a+(1-v)\theta(b,a)+\frac{1}{2}\theta(a+v\theta(b,a),a+(1-v)\theta(b,a)\right)\\ &\times\mathcal{W}\left(a+(1-v)\theta(b,a)+\frac{1}{2}\theta(a+v\theta(b,a),a+(1-v)\theta(b,a)\right)\\ &\leq\left(\frac{\mathcal{F}(a+v\theta(b,a))+\mathcal{F}(a+(1-v)\theta(b,a))}{2}\right)\left(\frac{\mathcal{W}(a+v\theta(b,a))+\mathcal{W}(a+(1-v)\theta(b,a))}{2}\right)\\ &\leq\frac{\mathcal{F}(a+v\theta(b,a))\mathcal{W}(a+v\theta(b,a))}{4}+\frac{\mathcal{F}(a+(1-v)\theta(b,a))\mathcal{W}(a+(1-v)\theta(b,a))}{4}\\ &+\frac{v(1-v)}{2}M(a,b)+\frac{2v^{2}-2v+1}{4}N(a,b). \end{align}$

(2.10)

Multiplying both sides of inequality (2.10) by $e^{-\rho v}$ and integrating with respect to $v$ over $[0, 1]$ , we get

$\begin{align*} &\frac{1-e^{-\rho}}{\rho}\mathcal{F}\left(\frac{2a+\theta(b,a)}{2}\right)\mathcal{W}\left(\frac{2a+\theta(b,a)}{2}\right)\\ &\leq\int\limits_{0}^{1}e^{-\rho v}\frac{\mathcal{F}(a+v\theta(b,a))\mathcal{W}(a+v\theta(b,a))}{4}{\mathrm{d}}v\\ &+\int\limits_{0}^{1}e^{-\rho v}\frac{\mathcal{F}(a+(1-v)\theta(b,a))\mathcal{W}(a+(1-v)\theta(b,a))}{4}{\mathrm{d}}v\\ &+M(a,b)\int\limits_{0}^{1}e^{-\rho v}\frac{v(1-v)}{2}{\mathrm{d}}v+N(a,b)\int\limits_{0}^{1}e^{-\rho v}\frac{2v^{2}-2v+1}{4}{\mathrm{d}}v, \end{align*}$

which completes the proof.

3. Inequalities involving first and second order differentiable preinvex functions

Lemma 3.1. Assume that $\mathcal{F}:[a, a+\theta(b, a)]\subseteq\mathbb{R}\to\mathbb{R}$ is a differentiable function and $\mathcal{F^{\prime}}\in L[a, a+\theta(b, a)]$ . Then the following equality for the $k$ -fractional integrals with exponential kernel holds:

$\begin{align} R_{ab}& = \frac{\theta(b,a)}{2}\int\limits_{0}^{1}u\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v\\ &-\frac{\theta(b,a)}{2(1-e^{-\rho})}\left[\int\limits_{0}^{1}e^{-\rho v}\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v-\int\limits_{0}^{1}e^{-\rho(1-v)}\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right], \end{align}$

(3.1)

where

$\begin{align*} R_{ab} = \frac{k-\alpha}{2k(1-e^{-\rho})}[\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{F}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{F}(a)]-\mathcal{F}\left(\frac{2a+\theta(b,a)}{2}\right), \end{align*}$

and

$u = \left\{ \begin{array}{rcl} 1, \ & \mathit{{for}} & 0\leq v\leq\frac{1}{2}, \\\\ -1, & \mathit{{for}} & \frac{1}{2}\leq v\leq1. \end{array}\right.$

Proof. By simple calculations, we have

$\begin{align} &\int\limits_{0}^{1}e^{-\rho v}\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v\\ & = -\frac{1}{\theta(b,a)}\left[e^{-\rho v}\mathcal{F}(a+(1-v)\theta(b,a))\bigg|_{0}^{1}+\rho\int\limits_{0}^{1}e^{-\rho v}\mathcal{F}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]\\ & = -\frac{1}{\theta(b,a)}\left[\mathcal{F}(a+\theta(b,a))-e^{-\rho}\mathcal{F}(a)+\rho\int\limits_{0}^{1}e^{-\rho v}\mathcal{F}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]\\ & = \frac{1}{\theta(b,a)}\left[\mathcal{F}(a+\theta(b,a))-e^{-\rho}\mathcal{F}(a)-\frac{\rho}{\theta(b,a)}\int\limits_{a}^{a+\theta(b,a)}e^{-\rho \frac{a+\theta(b,a)-s}{\theta(b,a)}}\mathcal{F}(s){\mathrm{d}}s\right]\\ & = \frac{1}{\theta(b,a)}\left[\mathcal{F}(a+\theta(b,a))-e^{-\rho}\mathcal{F}(a)-\frac{k-\alpha}{\alpha}\int\limits_{a}^{a+\theta(b,a)}e^{-\frac{(k-\alpha)}{\alpha} (a+\theta(b,a)-s)}\mathcal{F}(s){\mathrm{d}}s\right]\\ & = \frac{1}{\theta(b,a)}\left[\mathcal{F}(a+\theta(b,a))-e^{-\rho}\mathcal{F}(a)-\frac{k-\alpha}{k}\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{F}(a+\theta(b,a))\right], \end{align}$

(3.2)

similarly

$\begin{align} &\int\limits_{0}^{1}e^{-\rho (1-v)}\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v\\ & = -\frac{1}{\theta(b,a)}\left[e^{-\rho (1-v)}\mathcal{F}(a+(1-v)\theta(b,a))\bigg|_{0}^{1}-\rho\int\limits_{0}^{1}e^{-\rho (1-v)}\mathcal{F}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]\\ & = -\frac{1}{\theta(b,a)}\left[\mathcal{F}(a)-e^{-\rho}\mathcal{F}(a+\theta(b,a))-\rho\int\limits_{0}^{1}e^{-\rho v}\mathcal{F}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]\\ & = \frac{1}{\theta(b,a)}\left[e^{-\rho}\mathcal{F}(a+\theta(b,a))-\mathcal{F}(a)+\frac{\rho}{\theta(b,a)}\int\limits_{a}^{a+\theta(b,a)}e^{-\rho \frac{s-a}{\theta(b,a)}}\mathcal{F}(s){\mathrm{d}}s\right]\\ & = \frac{1}{\theta(b,a)}\left[e^{-\rho}\mathcal{F}(a+\theta(b,a))-\mathcal{F}(a)+\frac{k-\alpha}{\alpha}\int\limits_{a}^{a+\theta(b,a)}e^{-\frac{(k-\alpha)}{\alpha} (s-a)}\mathcal{F}(s){\mathrm{d}}s\right]\\ & = \frac{1}{\theta(b,a)}\left[e^{-\rho}\mathcal{F}(a+\theta(b,a))-\mathcal{F}(a)+\frac{k-\alpha}{k}\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{F}(a)\right]. \end{align}$

(3.3)

Also note that

(3.4)

Substituting (3.2), (3.3) and (3.4) in (3.1) completes the proof.

Theorem 3.1. Assume that $\mathcal{F}:[a, a+\theta(b, a)]\subseteq\mathbb{R}\to\mathbb{R}$ is a differentiable function and $\mathcal{F^{\prime}}\in L[a, a+\theta(b, a)]$ and $|\mathcal{F^{\prime}}|$ is preinvex on $[a, a+\theta(b, a)]$ . Then the following inequality for the $k$ -fractional integrals with exponential kernel holds:

$\begin{align*} |R_{ab}|\leq\frac{\theta(b,a)}{2}\left(\frac{1}{2}-\frac{\tanh\left(\frac{\rho}{4}\right)}{\rho}\right)(|\mathcal{F}^{\prime}(a)|+|\mathcal{F}^{\prime}(b)|). \end{align*}$

Proof. Using Lemma 3.1, preinvexity of $|\mathcal{F^{\prime}}|$ and increasing property of exponential function, we have

$\begin{align*} |R_{ab}|& = \left|\frac{\theta(b,a)}{2}\int\limits_{0}^{1}u\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right.\\ &-\left.\frac{\theta(b,a)}{2(1-e^{-\rho})}\left[\int\limits_{0}^{1}e^{-\rho v}\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v-\int\limits_{0}^{1}e^{-\rho(1-v)}\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]\right|\\ &\leq\frac{\theta(b,a)}{2(1-e^{-\rho})}\left[\int\limits_{0}^{\frac{1}{2}}\left(1-e^{-\rho}-e^{-\rho v}+e^{-\rho(1-v)}\right)|\mathcal{F^{\prime}}(a+(1-v)\theta(b,a))|{\mathrm{d}}v\right.\\ &\left.-\int\limits_{\frac{1}{2}}^{1}\left(1-e^{-\rho}-e^{-\rho (1-v)}+e^{-\rho v}\right)|\mathcal{F^{\prime}}(a+(1-v)\theta(b,a))|{\mathrm{d}}v\right]\\ &\leq\frac{\theta(b,a)}{2(1-e^{-\rho})}\left[\int\limits_{0}^{\frac{1}{2}}\left(1-e^{-\rho}-e^{-\rho v}+e^{-\rho(1-v)}\right)(v|\mathcal{F^{\prime}}(a)|+(1-v)|\mathcal{F^{\prime}}(b)|){\mathrm{d}}v\right.\\ &\left.-\int\limits_{\frac{1}{2}}^{1}\left(1-e^{-\rho}-e^{-\rho (1-v)}+e^{-\rho v}\right)(v|\mathcal{F^{\prime}}(a)|+(1-v)|\mathcal{F}^{\prime}(b)|){\mathrm{d}}v\right] \\ & = \frac{\theta(b,a)}{2(1-e^{-\rho})}\left[\int\limits_{0}^{\frac{1}{2}}\left(1-e^{-\rho}-e^{-\rho v}+e^{-\rho(1-v)}\right)(v|\mathcal{F^{\prime}}(a)|+(1-v)|\mathcal{F^{\prime}}(b)|){\mathrm{d}}v\right.\\ &\left.-\int\limits_{0}^{\frac{1}{2}}\left(1-e^{-\rho}-e^{-\rho v}+e^{-\rho (1-v)}\right)((1-v)|\mathcal{F^{\prime}}(a)|+v|\mathcal{F}^{\prime}(b)|){\mathrm{d}}v\right]\\ & = \frac{\theta(b,a)}{2(1-e^{-\rho})}\int\limits_{0}^{\frac{1}{2}}\left(1-e^{-\rho}-e^{-\rho v}+e^{-\rho(1-v)}\right)(|\mathcal{F^{\prime}}(a)|+|\mathcal{F^{\prime}}(b)|){\mathrm{d}}v\\ & = \frac{\theta(b,a)}{2(1-e^{-\rho})}\left[\frac{1-e^{-\rho}}{2}-\frac{1}{\rho}(1-e^{\frac{-\rho}{2}})^{2}\right](|\mathcal{F^{\prime}}(a)|+|\mathcal{F^{\prime}}(b)|)\\ & = \frac{\theta(b,a)}{2}\left(\frac{1}{2}-\frac{\tanh\left(\frac{\rho}{4}\right)}{\rho}\right)(|\mathcal{F}^{\prime}(a)|+|\mathcal{F}^{\prime}(b)|), \end{align*}$

which completes the proof.

Lemma 3.2. Assume that $\mathcal{F}:[a, a+\theta(b, a)]\subseteq\mathbb{R}\to\mathbb{R}$ is a differentiable function and $\mathcal{F^{\prime}}\in L[a, a+\theta(b, a)]$ . Then the following equality for the fractional integrals with exponential kernel holds:

$\begin{align} L_{ab} = \frac{\theta(b,a)}{2(1-e^{-\rho})}\left[\int\limits_{0}^{1}e^{-\rho v}\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v-\int\limits_{0}^{1}e^{-\rho(1-v)}\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right], \end{align}$

(3.5)

where

$L_{ab} = \frac{\mathcal{F}(a)+\mathcal{F}(a+\theta(b,a)}{2}-\frac{k-\alpha}{2k(1-e^{-\rho})}[\,_{k}\mathcal{I}^{\alpha}_{(a)^{+}}\mathcal{F}(a+\theta(b,a))+\,_{k}\mathcal{I}^{\alpha}_{(a+\theta(b,a))^{-}}\mathcal{F}(a)].$

Proof. Using 3.2 and 3.3, we get the required result.

Theorem 3.2. Assume that $\mathcal{F}:[a, a+\theta(b, a)]\subseteq\mathbb{R}\to\mathbb{R}$ is a differentiable function and $\mathcal{F^{\prime}}\in L[a, a+\theta(b, a)]$ . If $|\mathcal{F^{\prime}}|$ is preinvex function on $[a, a+\theta(b, a)]$ , then the following inequality for the $k$ -fractional integrals with exponential kernel holds:

$\begin{align*} |L_{ab}|\leq\frac{\theta(b,a)}{2\rho}\tanh\left(\frac{\rho}{4}\right)(|\mathcal{F}^{\prime}(a)|+|\mathcal{F}^{\prime}(b)|). \end{align*}$

Proof. Using Lemma 3.2, preinvexity of $|\mathcal{F^{\prime}}|$ and increasing property of exponential function, we have

$\begin{align*} |L_{ab}|&\leq\frac{\theta(b,a)}{2}\int\limits_{0}^{1}\frac{|e^{-\rho v}-e^{-\rho(1-v)}|}{1-e^{-\rho}}|\mathcal{F^{\prime}}(a+v\theta(b,a))|{\mathrm{d}}v\\ &\leq\frac{\theta(b,a)}{2}\left[\int\limits_{0}^{1}\frac{|e^{-\rho v}-e^{-\rho(1-v)}|}{1-e^{-\rho}}v|\mathcal{F^{\prime}}(a)|{\mathrm{d}}v+\int\limits_{0}^{1}\frac{|e^{-\rho v}-e^{-\rho(1-v)}|}{1-e^{-\rho}}(1-v)|\mathcal{F^{\prime}}(b)|{\mathrm{d}}v\right] \\& = \frac{\theta(b,a)}{2}|\mathcal{F^{\prime}}(a)|\left[\int\limits_{0}^{\frac{1}{2}}\frac{e^{-\rho v}-e^{-\rho(1-v)}}{1-e^{-\rho}}v{\mathrm{d}}v+\int\limits^{1}_{\frac{1}{2}}\frac{e^{-\rho (1-v)}-e^{-\rho v}}{1-e^{-\rho}}v{\mathrm{d}}v\right]\\ &+\frac{\theta(b,a)}{2}|\mathcal{F^{\prime}}(b)|\left[\int\limits_{0}^{\frac{1}{2}}\frac{e^{-\rho v}-e^{-\rho(1-v)}}{1-e^{-\rho}}(1-v){\mathrm{d}}v+\int\limits^{1}_{\frac{1}{2}}\frac{e^{-\rho (1-v)}-e^{-\rho v}}{1-e^{-\rho}}(1-v){\mathrm{d}}v\right]\\ & = \frac{\theta(b,a)}{2\rho}\tanh\left(\frac{\rho}{4}\right)(|\mathcal{F}^{\prime}(a)|+|\mathcal{F}^{\prime}(b)|), \end{align*}$

which completes the proof.

Lemma 3.3. Assume that $\mathcal{F}:[a, a+\theta(b, a)]\subseteq\mathbb{R}\to\mathbb{R}$ is twice differentiable function and $\mathcal{F^{\prime\prime}}\in L[a, a+\theta(b, a)]$ . Then the following equality for the $k$ -fractional integrals with exponential kernel holds:

$\begin{align} L_{ab} = \frac{\theta^{2}(b,a)}{2\rho(1-e^{-\rho})}\int\limits_{0}^{1}\left(1+e^{-\rho}-e^{-\rho v}-e^{-\rho(1-v)}\right)\mathcal{F^{\prime\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v. \end{align}$

(3.6)

Proof. Using (3.5) and integration by parts, we have

$\begin{align} &\int\limits_{0}^{1}e^{-\rho v}\mathcal{F}^{\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\\ & = -\frac{1}{\rho}\left[e^{-\rho v}\mathcal{F}^{\prime}(a+(1-v)\theta(b,a))\bigg|_{0}^{1}+\theta(b,a)\int\limits_{0}^{1}e^{-\rho v}\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]\\ & = -\frac{1}{\rho}\left[e^{-\rho}\mathcal{F}^{\prime}(a)-\mathcal{F}^{\prime}(a+\theta(b,a))+\theta(b,a)\int\limits_{0}^{1}e^{-\rho v}\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right], \end{align}$

(3.7)

similarly

$\begin{align} &\int\limits_{0}^{1}e^{-\rho (1-v)}\mathcal{F}^{\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\\ & = \frac{1}{\rho}\left[e^{-\rho (1-v)}\mathcal{F}^{\prime}(a+(1-v)\theta(b,a))\bigg|_{0}^{1}+\theta(b,a)\int\limits_{0}^{1}e^{-\rho (1-v)}\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]\\ & = \frac{1}{\rho}\left[\mathcal{F}^{\prime}(a)-e^{-\rho}\mathcal{F}^{\prime}(a+\theta(b,a))+\theta(b,a)\int\limits_{0}^{1}e^{-\rho (1-v)}\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]. \end{align}$

(3.8)

Substituting (3.7) and (3.8) in (3.5), we have

$\begin{align*} L_{ab}& = \frac{\theta(b,a)}{2\rho(1-e^{-\rho})}\left[(1+e^{-\rho})(\mathcal{F}^{\prime}(a+\theta(b,a))-\mathcal{F}^{\prime}(a))\right.\\ &\left.-\theta(b,a)\int\limits_{0}^{1}(e^{-\rho v}+e^{-\rho (1-v)})\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]\\ & = \frac{\theta^{2}(b,a)}{2\rho(1-e^{-\rho})}\int\limits_{0}^{1}\left(1+e^{-\rho}-e^{-\rho v}-e^{-\rho(1-v)}\right)\mathcal{F^{\prime\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v, \end{align*}$

which completes the proof.

Theorem 3.3. Assume that $\mathcal{F}:[a, a+\theta(b, a)]\subseteq\mathbb{R}\to\mathbb{R}$ is twice differentiable function. If $\mathcal{F^{\prime\prime}}\in L[a, a+\theta(b, a)]$ and $|\mathcal{F^{\prime\prime}}|$ is preinvex on $[a, a+\theta(b, a)]$ , then the following inequality for the $k$ -fractional integrals with exponential kernel holds:

$\begin{align*} |L_{ab}|\leq\frac{\theta^{2}(b,a)}{2\rho(1-e^{-\rho})}\left(\frac{1+e^{-\rho}}{2}-\frac{1-e^{-\rho}}{\rho}\right)(|\mathcal{F}^{\prime\prime}(a)|+|\mathcal{F}^{\prime\prime}(b)|). \end{align*}$

Proof. It is to be noted that

$\begin{align} \int\limits_{0}^{1}\left(1+e^{-\rho}-e^{-\rho v}-e^{-\rho(1-v)}\right)v{\mathrm{d}}v = \frac{1+e^{-\rho}}{2}-\frac{1-e^{-\rho}}{\rho}, \end{align}$

(3.9)

and

$\begin{align} \int\limits_{0}^{1}\left(1+e^{-\rho}-e^{-\rho v}-e^{-\rho(1-v)}\right)(1-v){\mathrm{d}}v = \frac{1+e^{-\rho}}{2}-\frac{1-e^{-\rho}}{\rho}. \end{align}$

(3.10)

Using (3.6), (3.9), (3.10) and the preinvexity of $|\mathcal{F}^{\prime\prime}|,$ we have

$\begin{align*} L_{ab}& = \bigg|\frac{\theta^{2}(b,a)}{2\rho(1-e^{-\rho})}\int\limits_{0}^{1}\left(1+e^{-\rho}-e^{-\rho v}-e^{-\rho(1-v)}\right)\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\bigg|\\ &\leq\frac{\theta^{2}(b,a)}{2\rho(1-e^{-\rho})}\int\limits_{0}^{1}\left(1+e^{-\rho}-e^{-\rho v}-e^{-\rho(1-v)}\right)|\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a))|{\mathrm{d}}v\\ &\leq\frac{\theta^{2}(b,a)}{2\rho(1-e^{-\rho})}\int\limits_{0}^{1}\left(1+e^{-\rho}-e^{-\rho v}-e^{-\rho(1-v)}\right)(v|\mathcal{F}^{\prime\prime}(a)|+(1-v)|\mathcal{F}^{\prime\prime}(b)|){\mathrm{d}}v\\ & = \frac{\theta^{2}(b,a)}{2\rho(1-e^{-\rho})}\left(\frac{1+e^{-\rho}}{2}-\frac{1-e^{-\rho}}{\rho}\right)(|\mathcal{F}^{\prime\prime}(a)|+|\mathcal{F}^{\prime\prime}(b)|), \end{align*}$

the proof is complete.

Lemma 3.4. Assume that $\mathcal{F}:[a, a+\theta(b, a)]\subseteq\mathbb{R}\to\mathbb{R}$ is twice differentiable function and $\mathcal{F^{\prime\prime}}\in L[a, a+\theta(b, a)]$ . Then the following equality for the $k$ -fractional integrals with exponential kernel holds:

$\begin{align} R_{ab} = \frac{\theta^{2}(b,a)}{2}\int\limits_{0}^{1}h(v)\mathcal{F^{\prime\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v, \end{align}$

(3.11)

where

$h(v) = \left\{ \begin{array}{rcl} v-\frac{1+e^{-\rho}-e^{-\rho v}-e^{-\rho(1-v)}}{\rho(1-e^{-\rho})}, \ & \mathit{{for}} & 0\leq v\leq\frac{1}{2}, \\\\ (1-v)-\frac{1+e^{-\rho}-e^{-\rho v}-e^{-\rho(1-v)}}{\rho(1-e^{-\rho})}, & \mathit{{for}} & \frac{1}{2}\leq v\leq1. \end{array}\right.$

Proof. Using (3.1), we have

$\begin{align*} R_{ab}& = \frac{\theta(b,a)}{2}\int\limits_{0}^{1}u\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v\\ &-\frac{\theta(b,a)}{2(1-e^{-\rho})}\left[\int\limits_{0}^{1}e^{-\rho v}\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v-\int\limits_{0}^{1}e^{-\rho(1-v)}\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right], \end{align*}$

Thus

$\begin{align} &\frac{\theta(b,a)}{2}\int\limits_{0}^{1}u\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v\\ & = \frac{\theta(b,a)}{2}\left[\int\limits_{0}^{\frac{1}{2}}\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v-\int\limits_{\frac{1}{2}}^{1}\mathcal{F^{\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]\\ & = \frac{\theta(b,a)}{2}\left[v\mathcal{F}^{\prime}(a+(1-v)\theta(b,a))\bigg|_{0}^{\frac{1}{2}}+\theta(b,a)\int\limits_{0}^{\frac{1}{2}}v\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]\\ &-\frac{\theta(b,a)}{2}\left[v\mathcal{F}^{\prime}(a+(1-v)\theta(b,a))\bigg|_{\frac{1}{2}}^{1}+\theta(b,a)\int\limits_{\frac{1}{2}}^{1}v\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]\\ & = \frac{\theta(b,a)}{2}\left[\frac{1}{2}\mathcal{F}^{\prime}\left(\frac{2a+\theta(b,a)}{2}\right)+\theta(b,a)\int\limits_{0}^{\frac{1}{2}}v\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]\\ &-\frac{\theta(b,a)}{2}\left[\mathcal{F}^{\prime}(a)-\frac{1}{2}\mathcal{F}^{\prime}\left(\frac{2a+\theta(b,a)}{2}\right)+\theta(b,a)\int\limits_{\frac{1}{2}}^{1}v\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\right]\\ & = \frac{\theta(b,a)}{2}\left[\mathcal{F}^{\prime}\left(\frac{2a+\theta(b,a)}{2}\right)-\mathcal{F}^{\prime}(a)\right]\\ &+\frac{\theta^{2}(b,a)}{2}\int\limits_{0}^{\frac{1}{2}}v\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v-\frac{\theta^{2}(b,a)}{2}\int\limits_{\frac{1}{2}}^{1}v\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v \\ & = \frac{\theta^{2}(b,a)}{2}\int\limits_{\frac{1}{2}}^{1}\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\\ &+\frac{\theta^{2}(b,a)}{2}\int\limits_{0}^{\frac{1}{2}}v\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v-\frac{\theta^{2}(b,a)}{2}\int\limits_{\frac{1}{2}}^{1}v\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v\\ & = \frac{\theta^{2}(b,a)}{2}\int\limits_{0}^{\frac{1}{2}}v\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v+\frac{\theta^{2}(b,a)}{2}\int\limits_{\frac{1}{2}}^{1}(1-v)\mathcal{F}^{\prime\prime}(a+(1-v)\theta(b,a)){\mathrm{d}}v. \end{align}$

(3.12)

Substituting (3.7), (3.8) and (3.12) in (3.1), we get the required result.

Theorem 3.4. Assume that $\mathcal{F}:[a, a+\theta(b, a)]\subseteq\mathbb{R}\to\mathbb{R}$ is a twice differentiable function. If $\mathcal{F^{\prime\prime}}\in L[a, a+\theta(b, a)]$ and $|\mathcal{F^{\prime\prime}}|$ is preinvex on $[a, a+\theta(b, a)]$ , then the following inequality for the $k$ -fractional integrals with exponential kernel holds:

$\begin{align*} |R_{ab}|\leq\frac{\theta^{2}(b,a)}{2}\left(\frac{1}{8}+\frac{1+e^{-\rho}}{2\rho(1-e^{-\rho})}-\frac{1}{\rho^{2}}\right)(|\mathcal{F}^{\prime\prime}(a)|+|\mathcal{F}^{\prime\prime}(b)|). \end{align*}$

Proof. Using (3.11) and preinvexity of $|\mathcal{F^{\prime\prime}}|$ , we have

$\begin{align*} |R_{ab}|& = \bigg|\frac{\theta^{2}(b,a)}{2}\int\limits_{0}^{1}h(v)\mathcal{F^{\prime\prime}}(a+(1-v)\theta(b,a)){\mathrm{d}}v\bigg|\\ &\leq\frac{\theta^{2}(b,a)}{2}\int\limits_{0}^{1}h(v)|\mathcal{F^{\prime\prime}}(a+(1-v)\theta(b,a))|{\mathrm{d}}v\\ &\leq\frac{\theta^{2}(b,a)}{2}\int\limits_{0}^{\frac{1}{2}}\left(v-\frac{1+e^{-\rho}-e^{-\rho v}-e^{-\rho(1-v)}}{\rho(1-e^{-\rho})}\right)(v|\mathcal{F}^{\prime\prime}(a)|+(1-v)|\mathcal{F}^{\prime\prime}(b)|){\mathrm{d}}v\\ &+\frac{\theta^{2}(b,a)}{2}\int\limits_{\frac{1}{2}}^{1}\left(1-v-\frac{1+e^{-\rho}-e^{-\rho v}-e^{-\rho(1-v)}}{\rho(1-e^{-\rho})}\right)(v|\mathcal{F}^{\prime\prime}(a)|+(1-v)|\mathcal{F}^{\prime\prime}(b)|){\mathrm{d}}v\\ & = \frac{\theta^{2}(b,a)}{2}\left[\int\limits_{0}^{\frac{1}{2}}(v^{2}|\mathcal{F}^{\prime\prime}(a)|+v(1-v)|\mathcal{F}^{\prime\prime}(b)|){\mathrm{d}}v+\int\limits_{\frac{1}{2}}^{1}(v(1-v)|\mathcal{F}^{\prime\prime}(a)|+(1-v)^{2}|\mathcal{F}^{\prime\prime}(b)|){\mathrm{d}}v\right]\\ &+\frac{1}{\rho(1-e^{-\rho})}\int\limits_{0}^{1}\left(1+e^{-\rho}-e^{-\rho v}-e^{-\rho(1-v)}\right)(v|\mathcal{F}^{\prime\prime}(a)|+(1-v)|\mathcal{F}^{\prime\prime}(b)|){\mathrm{d}}v\\ & = \frac{\theta^{2}(b,a)}{2}\left(\frac{1}{8}+\frac{1+e^{-\rho}}{2\rho(1-e^{-\rho})}-\frac{1}{\rho^{2}}\right)(|\mathcal{F}^{\prime\prime}(a)|+|\mathcal{F}^{\prime\prime}(b)|). \end{align*}$

The proof is complete.

Remark 3.1. We would like to remark here that by taking $k\to1$ , new results can be obtained from our results.

Applications

We now discuss some applications of the results obtained in previous section. Before we proceed further let us recall the definition of arithmetic mean.

The arithmetic mean is defined as

$\mathcal{A}(a,b): = \frac{a+b}{2}, \quad a\neq b.$

Proposition 3.1. Suppose all the assumptions of Theorem 3.1 are satisfied, then

$\begin{align*} &\left|\alpha {\mathcal{A}}(a^2,b^2)+\frac{\alpha^2}{(1-\alpha)^2}[(a-b)(1-\alpha)+2\alpha]-{\mathcal{A}}^2(a,b)\right|\\ &\leq(b-a){\mathcal{A}}(a,b)\left(\frac{1}{2}-\frac{tanh\left(\frac{\rho_1}{4}\right)}{\rho_1}\right). \end{align*}$

Proof. The proof directly follows from Theorem 3.1 by setting $\theta(b, a) = b-a, k = 1$ and $\mathcal{F}(x) = x^2$ .

Proposition 3.2. Suppose all the assumptions of Theorem 3.2 are satisfied, then

$\begin{align*} &\left|(1-\alpha) {\mathcal{A}}(a^2,b^2)-\frac{\alpha^2}{(1-\alpha)^2}[(a-b)(1-\alpha)+2\alpha]\right|\\ &\leq(b-a){\mathcal{A}}(a,b)tanh\left(\frac{\rho_1}{4}\right). \end{align*}$

Proof. The proof directly follows from Theorem 3.2 by setting $\theta(b, a) = b-a, k = 1$ and $\mathcal{F}(x) = x^2$ .

Proposition 3.3. Suppose all the assumptions of Theorem 3.3 are satisfied, then

$\begin{align*} &\left|(1-\alpha) {\mathcal{A}}(a^2,b^2)-\frac{\alpha^2}{(1-\alpha)^2}[(a-b)(1-\alpha)+2\alpha]\right|\\ &\leq \frac{2(b-a)^2{\mathcal{A}}(a,b)}{\rho_1(1-e^{-\rho_1})}\left(\frac{1+e^{-\rho_1}}{2}-\frac{1-e^{-\rho_1}}{\rho_1}\right). \end{align*}$

Proof. The proof directly follows from Theorem 3.3 by setting $\theta(b, a) = b-a, k = 1$ and $\mathcal{F}(x) = x^2$ .

Proposition 3.4. Suppose all the assumptions of Theorem 3.4 are satisfied, then

$\begin{align*} &\left|\alpha {\mathcal{A}}(a^2,b^2)+\frac{\alpha^2}{(1-\alpha)^2}[(a-b)(1-\alpha)+2\alpha]-{\mathcal{A}}^2(a,b)\right|\\ &\leq{2(b-a)^2{\mathcal{A}}(a,b)}\left(\frac{1}{8}+\frac{1+e^{-\rho_1}}{2\rho_1(1-e^{-\rho_1})}-\frac{1}{\rho_1^{2}}\right). \end{align*}$

Proof. The proof directly follows from Theorem 3.4 by setting $\theta(b, a) = b-a, k = 1$ and $\mathcal{F}(x) = x^2$ . We now discuss applications to $q$ -digamma functions, which is defined as:

Suppose $0 < q < 1$ , the $q$ -digamma function $\chi_q(u)$ is given as

$\begin{align*} \chi_q(u)& = -\ln(1-q)+\ln(q)\sum\limits_{i = 0}^{\infty}\frac{q^{i+u}}{1-q^{i+u}}.\\ & = -\ln(1-q)+\ln(q)\sum\limits_{i = 0}^{\infty}\frac{q^{iu}}{1-q^{iu}}. \end{align*}$

For $q > 1, t > 0$ , then $q$ -digamma function $\chi_q$ can be given as

$\begin{align*} \chi_q(u)& = -\ln(q-1)+\ln(q)\left[u-\frac{1}{2}-\sum\limits_{i = 0}^{\infty}\frac{q^{-(i+u)}}{1-q^{-(i+u)}}\right].\\ & = -\ln(q-1)+\ln(q)\left[u-\frac{1}{2}-\sum\limits_{i = 0}^{\infty}\frac{q^{-iu}}{1-q^{-iu}}\right]. \end{align*}$

From the above definition, it is clear that $\chi'_q$ is completely monotone on $(0, \infty)$ for $q > 0$ . This implies that $\chi'_q$ is convex. For more details, see ^[5].

Proposition 3.5. Under the assumption of Theorem 2.1, the following inequality holds:

$\begin{align*} \chi_q\left(\frac{{a}+{b}}{2}\right)\leq\frac{1-\alpha}{2(1-e^{-\rho_1})}\left[\int_a^be^{-\frac{1-\alpha}{\alpha}(b-v)}\chi_q(v){\mathrm{d}}v +\int_a^be^{-\frac{1-\alpha}{\alpha}(v-a)}\chi_q(v){\mathrm{d}}v\right]\leq\frac{\chi_q(a)+\chi_q(b)}{2}. \end{align*}$

Proof. The proof is direct consequence of Theorem 2.1, by choosing $\theta(b, a) = b-a, k = 1$ and $\mathcal{F}(v)\rightarrow \chi_q(v)$ .

Proposition 3.6. Under the assumption of Theorem 3.1, the following inequality holds:

$\begin{align*} &\left|\frac{1-\alpha}{2(1-e^{-\rho_1})}\left[\int_a^be^{-\frac{1-\alpha}{\alpha}(b-v)}\chi_q(v){\mathrm{d}}v +\int_a^be^{-\frac{1-\alpha}{\alpha}(v-a)}\chi_q(v){\mathrm{d}}v\right]-\chi_q\left(\frac{a+b}{2}\right)\right|\\ &\leq\frac{b-a}{2}\left(\frac{1}{2}-\frac{\tanh\left(\frac{\rho_1}{4}\right)}{\rho_1}\right)(|\chi_q^{\prime}(a)|+|\chi_q^{\prime}(b)|). \end{align*}$

Proof. The proof is direct consequence of Theorem 3.1, by choosing $\theta(b, a) = b-a, k = 1$ and $\mathcal{F}(v)\rightarrow \chi_q(v)$ .

Proposition 3.7. Under the assumption of Theorem 3.1, the following inequality holds:

$\begin{align*} &\left|\frac{\chi_q(a)+\chi_q(b)}{2}-\frac{1-\alpha}{2(1-e^{-\rho_1})}\left[\int_a^be^{-\frac{1-\alpha}{\alpha}(b-v)}\chi_q(v){\mathrm{d}}v +\int_a^be^{-\frac{1-\alpha}{\alpha}(v-a)}\chi_q(v){\mathrm{d}}v\right]\right|\\ &\leq\frac{b-a}{2}\tanh\left(\frac{\rho_1}{4}\right)(|\chi_q^{\prime}(a)|+|\chi_q^{\prime}(b)|). \end{align*}$

Proof. The proof is direct consequence of Theorem 3.1, by choosing $\theta(b, a) = b-a, k = 1$ and $\mathcal{F}(v)\rightarrow \chi_q(v)$ .

4. Conclusions

In the article, we have extended the fractional integral operators with an exponential kernel to $k$ -fractional integral operators with an exponential kernel and derived several new trapezium type integral inequalities involving the new fractional integral operator essentially using the functions having preinvexity property. We have also discussed some interesting applications of our obtained results, which show the significance of our main results. It is also worth mentioning here that our obtained results are the generalizations of some previously known results and our ideas may lead to a lot of follow-up research.

Acknowledgment

The authors are thankful to the editor and anonymous reviewers for their valuable comments and suggestions. This research was funded by Dirección de Investigación from Pontificia Universidad Católica del Ecuador in the research project entitled, "Some integrals inequalities and generalized convexity" (Algunas desigualdades integrales para funciones con algún tipo de convexidad generalizada y aplicaciones).

References

[1]	A. W. Bergen, N. Caporaso, Cigarette smoking, J. Natl. Cancer Inst., 91 (1999), 1365–1376. https://doi.org/10.1093/jnci/91.16.1365 doi: 10.1093/jnci/91.16.1365
[2]	N. J. Wald, A. K. Hackshaw, Cigarette smoking: an epidemiological overview, Brit. Med. Bull., 52 (1996), 3–11. https://doi.org/10.1093/oxfordjournals.bmb.a011530 doi: 10.1093/oxfordjournals.bmb.a011530
[3]	B. Lloyd, K. Lucas, Smoking in adolescence: images and identities, London: Routledge, 1998.
[4]	S. Cohen, E. Lichtenstein, Perceived stress, quitting smoking, and smoking relapse, Health Psychol., 9 (1990), 466–478. https://doi.org/10.1037//0278-6133.9.4.466 doi: 10.1037//0278-6133.9.4.466
[5]	A. H. Mokdad, J. S. Marks, D. F. Stroup, J. L. Gerberding, Actual causes of death in the United States, JAMA, 291 (2004), 1238–1245. https://doi.org/10.1001/jama.291.10.1238 doi: 10.1001/jama.291.10.1238
[6]	A. Zeb, G. Zaman, S. Momani, Square-root dynamics of a giving up smoking model, Appl. Math. Model., 37 (2013), 5326–5334. https://doi.org/10.1016/j.apm.2012.10.005 doi: 10.1016/j.apm.2012.10.005
[7]	O. Sharomi, A. B. Gumel, Curtailing smoking dynamics: A mathematical modeling approach, Appl. Math. Comput., 195 (2008), 475–499. https://doi.org/10.1016/j.amc.2007.05.012 doi: 10.1016/j.amc.2007.05.012
[8]	Z. Alkhudhari, S. Al-Sheikh, S. Al-Tuwairqi, Stability analysis of a giving up smoking model, Int. J. Appl. Math. Res., 3 (2014), 168–177. http://doi.org/10.14419/ijamr.v3i2.2239 doi: 10.14419/ijamr.v3i2.2239
[9]	N. H. Shah, F. A. Thakkar, B. M. Yeolekar, Stability analysis of tuberculosis due to smoking, Int. J. Innov. Sci. Res. Tech., 3 (2018), 230–237.
[10]	Q. Din, M. Ozair, T. Hussain, U. Saeed, Qualitative behavior of asmoking model, Adv. Differ. Equ., 2016 (2016), 96. https://doi.org/10.1186/s13662-016-0830-6 doi: 10.1186/s13662-016-0830-6
[11]	A. M. Pulecio-Montoya, L. E. Lopez-Montenegro, L. M. Benavides, Analysis of a mathematical model of smoking, Contemp. Eng. Sci., 12 (2019), 117–129. https://doi.org/10.12988/ces.2019.9517 doi: 10.12988/ces.2019.9517
[12]	Z. Zhang, R. Wei, W. Xia, Dynamical analysis of a giving up smoking model with time delay, Adv. Differ. Equ., 2019 (2019), 505. https://doi.org/10.1186/s13662-019-2450-4 doi: 10.1186/s13662-019-2450-4
[13]	S. A. Khan, K. Shah, G. Zaman, F. Jarad, Existence theory and numerical solutions to smoking model under Caputo-Fabrizio fractional derivative, Chaos, 29 (2019), 013128. https://doi.org/10.1063/1.5079644 doi: 10.1063/1.5079644
[14]	S. Ucar, E. Ucar, N. ozdemir, Z. Hammouch, Mathematical analysis and numerical simulation for a smoking model with Atangana-Baleanu derivative, Chaos Soliton. Fract., 118 (2019), 300–306. https://doi.org/10.1016/j.chaos.2018.12.003 doi: 10.1016/j.chaos.2018.12.003
[15]	G. Rahman, R. P. Agarwal, Q. Din, Mathematical analysis of giving up smoking model via harmonic mean type incidence rate, Appl. Math. Comput., 354 (2019), 128–148. https://doi.org/10.1016/j.amc.2019.01.053 doi: 10.1016/j.amc.2019.01.053
[16]	C. Sun, J. Jia, Optimal control of a delayed smoking model with immigration, J. Biol. Dynam., 13 (2019), 447–460. https://doi.org/10.1080/17513758.2019.1629031 doi: 10.1080/17513758.2019.1629031
[17]	P. Veeresha, D. G. Prakasha, H. M. Baskonus, Solving smoking epidemic model of fractional order using a modified homotopy analysis transform method, Math. Sci., 13 (2019), 115–128. https://doi.org/10.1007/s40096-019-0284-6 doi: 10.1007/s40096-019-0284-6
[18]	A. M. S. Mahdy, N. H. Sweilam, M. Higazy, Approximate solution for solving nonlinear fractional order smoking model, Alex. Eng. J., 59 (2020), 739–752. https://doi.org/10.1016/j.aej.2020.01.049 doi: 10.1016/j.aej.2020.01.049
[19]	A. A. Alshareef, H. A. Batarfi, Stability analysis of chain, mild and passive smoking model, Amer. J. Comput. Math., 10 (2020), 31–42. https://doi.org/10.4236/ajcm.2020.101003 doi: 10.4236/ajcm.2020.101003
[20]	Z. Zhang, J. Zou, R. K. Upadhyay, A. Pratap, Stability and Hopf bifurcation analysis of a delayed tobacco smoking model containing snuffing class, Adv. Differ. Equ., 2020 (2020), 349. https://doi.org/10.1186/s13662-020-02808-5 doi: 10.1186/s13662-020-02808-5
[21]	A. Bernoussi, Global stability analysis of an SEIR epidemic model with relapse and general incidence rates, Appl. Sci., 21 (2019), 54–68.
[22]	M. Caputo, M. Fabrizio, A new definition of fractional derivative without singular kernel, Prog. Fract. Differ. Appl., 1 (2015), 73–85. https://doi.org/10.12785/pfda/010201 doi: 10.12785/pfda/010201
[23]	A. Atangana, D. Baleanu, New fractional derivatives with nonlocal and non-singular kernel: Theory and application to heat transfer model, Therm. Sci., 20 (2016), 763–769. https://doi.org/10.2298/TSCI160111018A doi: 10.2298/TSCI160111018A
[24]	D. Baleanu, A. Fernandez, On some new properties of fractional derivatives with Mittag-Leffler kernel, Commun. Nonlinear Sci. Numer. Simul., 59 (2018), 444–462. https://doi.org/10.1016/j.cnsns.2017.12.003 doi: 10.1016/j.cnsns.2017.12.003
[25]	T. Abdeljawad, Fractional operators with generalized Mittag-Leffler kernels and their iterated differintegrals, Chaos, 29 (2019), 023102. https://doi.org/10.1063/1.5085726 doi: 10.1063/1.5085726
[26]	T. Abdeljawad, D. Baleanu, On fractional derivatives with generalized Mittag-Leffler kernels, Adv. Differ. Equ., 2018 (2018), 468. https://doi.org/10.1186/s13662-018-1914-2 doi: 10.1186/s13662-018-1914-2
[27]	H. Khan, F. Jarad, T. Abdeljawad, A. Khan, A singular ABC-fractional differential equation with p-Laplacian operator, Chaos Soliton. Fract., 129 (2019), 56–61. https://doi.org/10.1016/j.chaos.2019.08.017 doi: 10.1016/j.chaos.2019.08.017
[28]	S. Rezapour, S. Etemad, H. Mohammadi, A mathematical analysis of a system of Caputo-Fabrizio fractional differential equations for the anthrax disease model in animals, Adv. Differ. Equ., 2020 (2020), 481. https://doi.org/10.1186/s13662-020-02937-x doi: 10.1186/s13662-020-02937-x
[29]	H. M. Alshehri, A. Khan, A fractional order Hepatitis C mathematical model with Mittag-Leffler kernel, J. Funct. Space., 2021 (2021), 2524027. https://doi.org/10.1155/2021/2524027 doi: 10.1155/2021/2524027
[30]	C. T. Deressa, S. Etemad, S. Rezapour, On a new four-dimensional model of memristor-based chaotic circuit in the context of nonsingular Atangana-Baleanu-Caputo operators, Adv. Differ. Equ., 2021 (2021), 444. https://doi.org/10.1186/s13662-021-03600-9 doi: 10.1186/s13662-021-03600-9
[31]	C. T. Deressa, S. Etemad, M. K. A. Kaabar, S. Rezapour, Qualitative analysis of a hyperchaotic Lorenz-Stenflo mathematical model via the Caputo fractional operator, J. Funct. Space., 2022 (2022), 4975104. https://doi.org/10.1155/2022/4975104 doi: 10.1155/2022/4975104
[32]	P. Kumar, V. S. Erturk, Environmental persistence influences infection dynamics for a butterfly pathogen via new generalised Caputo type fractional derivative, Chaos Soliton. Fract., 144 (2021), 110672. https://doi.org/10.1016/j.chaos.2021.110672 doi: 10.1016/j.chaos.2021.110672
[33]	A. Devi, A. Kumar, T. Abdeljawad, A. Khan, Stability analysis of solutions and existence theory of fractional Lagevin equation, Alex. Eng. J., 60 (2021), 3641–3647. https://doi.org/10.1016/j.aej.2021.02.011 doi: 10.1016/j.aej.2021.02.011
[34]	A. Pratap, R. Raja, R. P. Agarwal, J. Alzabut, M. Niezabitowski, E. Hincal, Further results on asymptotic and finite-time stability analysis of fractional-order time-delayed genetic regulatory networks, Neurocomputing, 475 (2022), 26–37. https://doi.org/10.1016/j.neucom.2021.11.088 doi: 10.1016/j.neucom.2021.11.088
[35]	H. Mohammadi, S. Kumar, S. Rezapour, S. Etemad, A theoretical study of the Caputo-Fabrizio fractional modeling for hearing loss due to Mumps virus with optimal control, Chaos Soliton. Fract., 144 (2021), 110668. https://doi.org/10.1016/j.chaos.2021.110668 doi: 10.1016/j.chaos.2021.110668
[36]	R. Begum, O. Tunc, H. Khan, H. Gulzar, A. Khan, A fractional order Zika virus model with Mittag-Leffler kernel, Chaos Soliton. Fract., 146 (2021), 110898. https://doi.org/10.1016/j.chaos.2021.110898 doi: 10.1016/j.chaos.2021.110898
[37]	A. Ali, Q. Iqbal, J. K. K. Asamoah, S. Islam, Mathematical modeling for the transmission potential of Zika virus with optimal control strategies, Eur. Phys. J. Plus, 137 (2022), 146. https://doi.org/10.1140/epjp/s13360-022-02368-5 doi: 10.1140/epjp/s13360-022-02368-5
[38]	P. Kumar, V. S. Erturk, H. Almusawa, Mathematical structure of mosaic disease using microbial biostimulants via Caputo and Atangana-Baleanu derivatives, Results Phys., 24 (2021), 104186. https://doi.org/10.1016/j.rinp.2021.104186 doi: 10.1016/j.rinp.2021.104186
[39]	R. Zarin, H. Khaliq, A. Khan, D. Khan, A. Akgul, U. W. Humphries, Deterministic and fractional modeling of a computer virus propagation, Results Phys., 33 (2022), 105130. https://doi.org/10.1016/j.rinp.2021.105130 doi: 10.1016/j.rinp.2021.105130
[40]	D. Baleanu, S. Etemad, S. Rezapour, A hybrid Caputo fractional modeling for thermostat with hybrid boundary value conditions, Bound. Value Probl., 2020 (2020), 64. https://doi.org/10.1186/s13661-020-01361-0 doi: 10.1186/s13661-020-01361-0
[41]	C. Thaiprayoon, W. Sudsutad, J. Alzabut, S. Etemad, S. Rezapour, On the qualitative analysis of the fractional boundary value problem describing thermostat control model via $\psi$ -Hilfer fractional operator, Adv. Differ. Equ., 2021 (2021), 201. https://doi.org/10.1186/s13662-021-03359-z doi: 10.1186/s13662-021-03359-z
[42]	J. Alzabut, G. M. Selvam, R. A. El-Nabulsi, D. Vignesh, M. E. Samei, Asymptotic stability of nonlinear discrete fractional pantograph equations with non-local initial conditions, Symmetry, 13 (2021), 473. https://doi.org/10.3390/sym13030473 doi: 10.3390/sym13030473
[43]	S. T. M. Thabet, S. Etemad, S. Rezapour, On a new structure of the pantograph inclusion problem in the Caputo conformable setting, Bound. Value Probl., 2020 (2020), 171. https://doi.org/10.1186/s13661-020-01468-4 doi: 10.1186/s13661-020-01468-4
[44]	P. Kumar, V. S. Erturk, A. Yusuf, K. S. Nisar, S. F. Abdelwahab, A study on canine distemper virus (CDV) and rabies epidemics in the red fox population via fractional derivatives, Results Phys., 25 (2021), 104281. https://doi.org/10.1016/j.rinp.2021.104281 doi: 10.1016/j.rinp.2021.104281
[45]	J. K. K. Asamoah, E. Okyere, E. Yankson, A. A. Opoku, A. Adom-Konadu, E. Acheampong, et al., Non-fractional and fractional mathematical analysis and simulations for Q fever, Chaos Soliton. Fract., 156 (2022), 111821. https://doi.org/10.1016/j.chaos.2022.111821 doi: 10.1016/j.chaos.2022.111821
[46]	H. Khan, C. Tunc, W. Chen, A. Khan, Existence theorems and Hyers-Ulam stability for a class of hybrid fractional differential equations with p-Laplacial operator, J. Appl. Anal. Comput., 8 (2018), 1211–1226. https://doi.org/10.11948/2018.1211 doi: 10.11948/2018.1211
[47]	A. Omame, U. K. Nwajeri, M. Abbas, C. P. Onyenegecha, A fractional order control model for Diabetes and COVID-19 co-dynamics with Mittag-Leffler function, Alex. Eng. J., 61 (2022), 7619–7635. https://doi.org/10.1016/j.aej.2022.01.012 doi: 10.1016/j.aej.2022.01.012
[48]	D. Baleanu, S. Etemad, H. Mohammadi, S. Rezapour, A novel modeling of boundary value problems on the glucose graph, Commun. Nonlinear Sci. Numer. Simul., 100 (2021), 105844. https://doi.org/10.1016/j.cnsns.2021.105844 doi: 10.1016/j.cnsns.2021.105844
[49]	S. Rezapour, B. Tellab, C. T. Deressa, S. Etemad, K. Nonlaopon, H-U-type stability and numerical solutions for a nonlinear model of the coupled systems of Navier BVPs via the generalized differential transform method, Fractal Fract., 5 (2021), 166. https://doi.org/10.3390/fractalfract5040166 doi: 10.3390/fractalfract5040166
[50]	E. Ucar, N. Özdemir, E. Altun, Fractional order model of immune cells influenced by cancer cells, Math. Model. Nat. Phenom., 14 (2019), 308. https://doi.org/10.1051/mmnp/2019002 doi: 10.1051/mmnp/2019002
[51]	E. Ucar, S. Ucar, F. Evirgen, N. Özdemir, A fractional SAIDR model in the frame of Atangana-Baleanu derivative, Fractal Fract., 50 (2021), 32. https://doi.org/10.3390/fractalfract5020032 doi: 10.3390/fractalfract5020032
[52]	S. Ucar, Existence and uniqueness results for a smoking model with determination and education in the frame of non-singular derivatives, Discrete Cont. Dyn. Syst. S, 14 (2021), 2571–2589. https://doi.org/10.3934/dcdss.2020178 doi: 10.3934/dcdss.2020178
[53]	S. Ucar, E. Ucar, N. Özdemir, Z. Hammouch, Mathematical analysis and numerical simulation for a smoking model with Atangana-Baleanu derivative, Chaos Soliton. Fract., 118 (2019), 300–306. https://doi.org/10.1016/j.chaos.2018.12.003 doi: 10.1016/j.chaos.2018.12.003
[54]	A. Khan, H. M. Alshehri, J. F. Gómez-Aguilar, Z. A. Khan, G. Fernández-Anaya, A predator-prey model involving variable-order fractional differential equations with Mittag-Leffler kernel, Adv. Differ. Equ., 2021 (2021), 183. https://doi.org/10.1186/s13662-021-03340-w doi: 10.1186/s13662-021-03340-w
[55]	H. M. Alshehri, A. Khan, A fractional order Hepatitis C mathematical model with Mittag-Leffler kernel, J. Funct. Space., 2021 (2021), 2524027. https://doi.org/10.1155/2021/2524027 doi: 10.1155/2021/2524027
[56]	P. Bedi, A. Kumar, A. Khan, Controllability of neutral impulsive fractional differential equations with Atangana-Baleanu-Caputo derivatives, Chaos Soliton. Fract., 150 (2021), 111153. https://doi.org/10.1016/j.chaos.2021.111153 doi: 10.1016/j.chaos.2021.111153
[57]	W. Chen, Time-space fabric underlying anomalous diffusion, Chaos Soliton. Fract., 28 (2006), 923–929. https://doi.org/10.1016/j.chaos.2005.08.199 doi: 10.1016/j.chaos.2005.08.199
[58]	R. Kanno, Representation of random walk in fractal space-time, Physica A, 248 (1998), 165–175. https://doi.org/10.1016/S0378-4371(97)00422-6 doi: 10.1016/S0378-4371(97)00422-6
[59]	W. Chen, H. G. Sun, X. Zhang, D. Korosak, Anomalous diffusion modeling by fractal and fractional derivatives, Comput. Math. Appl., 59 (2010), 1754–1758. https://doi.org/10.1016/j.camwa.2009.08.020 doi: 10.1016/j.camwa.2009.08.020
[60]	K. M. Owolabi, A. Atangana, A. Akgul, Modelling and analysis of fractal-fractional partial differential equations: Application to reaction-diffusion model, Alex. Eng. J., 59 (2020), 2477–2490. https://doi.org/10.1016/j.aej.2020.03.022 doi: 10.1016/j.aej.2020.03.022
[61]	Z. Ali, F. Rabiei, K. Shah, T. Khodadadi, Modeling and analysis of novel COVID-19 under fractal-fractional derivative with case study of malaysia, Fractals, 29 (2021), 2150020. https://doi.org/10.1142/S0218348X21500201 doi: 10.1142/S0218348X21500201
[62]	E. Bonyah, M. Yavuz, D. Baleanu, S. Kumar, A robust study on the listeriosis disease by adopting fractal-fractional operators, Alex. Eng. J., 61 (2022), 2016–2028. https://doi.org/10.1016/j.aej.2021.07.010 doi: 10.1016/j.aej.2021.07.010
[63]	M. Alqhtani, K. M. Saad, Numerical solutions of space-fractional diffusion equations via the exponential decay kernel, AIMS Mathematics, 7 (2022), 6535–6549. https://doi.org/10.3934/math.2022364 doi: 10.3934/math.2022364
[64]	M. Alqhtani, K. M. Saad, Fractal-fractional Michaelis–Menten enzymatic reaction model via different kernels, Fractal Fract., 6 (2021), 13. https://doi.org/10.3390/fractalfract6010013 doi: 10.3390/fractalfract6010013
[65]	K. M. Saad, J. F. Gomez-Aguilar, A. A. Almadiy, A fractional numerical study on a chronic hepatitis C virus infection model with immune response, Chaos Soliton. Fract., 139 (2020), 110062. https://doi.org/10.1016/j.chaos.2020.110062 doi: 10.1016/j.chaos.2020.110062
[66]	K. M. Saad, M. Alqhtani, Numerical simulation of the fractal-fractional reaction diffusion equations with general nonlinear, AIMS Mathematics, 6 (2021), 3788–3804. https://doi.org/10.3934/math.2021225 doi: 10.3934/math.2021225
[67]	A. Atangana, Fractal-fractional differentiation and integration: connecting fractal calculus and fractional calculus to predict complex system, Chaos Soliton. Fract., 102 (2017), 396–406. https://doi.org/10.1016/j.chaos.2017.04.027 doi: 10.1016/j.chaos.2017.04.027
[68]	M. Arfan, K. Shah, A. Ullah, Fractal-fractional mathematical model of four species comprising of prey-predation, Phys. Scr., 96 (2021), 124053. https://doi.org/10.1088/1402-4896/ac2f37 doi: 10.1088/1402-4896/ac2f37
[69]	M. Abdulwasaa, M. S. Abdo, K. Shah, T. A. Nofal, S. K. Panchal, S. V. Kawale, et al., Fractal-fractional mathematical modeling and forecasting of new cases and deaths of COVID-19 epidemic outbreaks in India, Results Phys., 20 (2021), 103702. https://doi.org/10.1016/j.rinp.2020.103702 doi: 10.1016/j.rinp.2020.103702
[70]	K. Shah, M. Arfan, I. Mahariq, A. Ahmadian, S. Salahshour, M. Ferrara, Fractal-fractional mathematical model addressing the situation of Corona virus in Pakistan, Results Phys., 19 (2020), 103560. https://doi.org/10.1016/j.rinp.2020.103560 doi: 10.1016/j.rinp.2020.103560
[71]	Z. A. Khan, M. ur Rahman, K. Shah, Study of a fractal-fractional smoking models with relapse and harmonic mean type incidence rate, J. Funct. Space., 2021 (2021) 6344079. https://doi.org/10.1155/2021/6344079 doi: 10.1155/2021/6344079
[72]	M. Arif, P. Kumam, W. Kumam, A. Akgul, T. Sutthibutpong, Analysis of newly developed fractal-fractional derivative with power law kernel for MHD couple stress fluid in channel embedded in a porous medium, Sci. Rep., 11 (2021), 20858. https://doi.org/10.1038/s41598-021-00163-3 doi: 10.1038/s41598-021-00163-3
[73]	H. Najafi, S. Etemad, N. Patanarapeelert, J. K. K. Asamoah, S. Rezapour, T. Sitthiwirattham, A study on dynamics of CD4 $^+$ T-cells under the effect of HIV-1 infection based on a mathematical fractal-fractional model via the Adams-Bashforth scheme and Newton polynomials. Mathematics, 10 (2022), 1366. https://doi.org/10.3390/math10091366 doi: 10.3390/math10091366
[74]	H. Khan, K. Alam, H. Gulzar, S. Etemad, S. Rezapour, A case study of fractal-fractional tuberculosis model in China: Existence and stability theories along with numerical simulations, Math. Comput. Simul., 198 (2022), 455–473. https://doi.org/10.1016/j.matcom.2022.03.009 doi: 10.1016/j.matcom.2022.03.009
[75]	A. Atangana, A. Akgul, K. M. Owolabi, Analysis of fractal fractional differential equations, Alex. Eng. J., 59 (2020), 1117–1134. https://doi.org/10.1016/j.aej.2020.01.005 doi: 10.1016/j.aej.2020.01.005
[76]	A. U. Awan, A. Sharif, K. A. Abro, M. Ozair, T. Hussain, Dynamical aspects of smoking model with cravings to smoke, Nonlinear Eng., 10 (2021), 91–108. http://doi.org/10.1515/nleng-2021-0008 doi: 10.1515/nleng-2021-0008

Reader Comments

Your name:*

Email:*
© 2022 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)