Fractional optimal control analysis of Covid-19 and dengue fever co-infection model with Atangana-Baleanu derivative

Asma Hanif; Azhar Iqbal Kashif Butt; Tariq Ismaeel; Asma Hanif; Azhar Iqbal Kashif Butt; Tariq Ismaeel

doi:10.3934/math.2024251

AIMS Mathematics

2024, Volume 9, Issue 3: 5171-5203. doi: 10.3934/math.2024251

Previous Article Next Article

Research article Special Issues

Fractional optimal control analysis of Covid-19 and dengue fever co-infection model with Atangana-Baleanu derivative

1.
Department of Mathematics, GC University, Lahore 54000, Pakistan
2.
Department of Mathematics and Statistics, College of Science, King Faisal University, Al-Ahsa 31982, Saudi Arabia

Received: 05 December 2023 Revised: 27 December 2023 Accepted: 08 January 2024 Published: 25 January 2024
MSC : 34H05, 49K15, 49K40

A co-infection with Covid-19 and dengue fever has had worse outcomes due to high mortality rates and longer stays either in isolation or at hospitals. This poses a great threat to a country's economy. To effectively deal with these threats, comprehensive approaches to prevent and control Covid-19/dengue fever co-infections are desperately needed. Thus, our focus is to formulate a new co-infection fractional model with the Atangana-Baleanu derivative to suggest effective and feasible approaches to restrict the spread of co-infection. In the first part of this paper, we present Covid-19 and dengue fever sub-models, as well as the co-infection model that is locally asymptotically stable when the respective reproduction numbers are less than unity. We establish the existence and uniqueness results for the solutions of the co-infection model. We extend the model to include a vaccination compartment for the Covid-19 vaccine to susceptible individuals and a treatment compartment to treat dengue-infected individuals as optimal control strategies for disease control. We outline the fundamental requirements for the fractional optimal control problem and illustrate the optimality system for the co-infection model using Pontraygin's principle. We implement the Toufik-Atangana approximating scheme to simulate the optimality system. The simulations show the effectiveness of the implemented strategy in determining optimal vaccination and treatment rates that decrease the cost functional to a minimum, thus significantly decreasing the number of infected humans and vectors. Additionally, we visualize a meaningful decrease in infection cases with an increase in the memory index. The findings of this study will provide reasonable disease control suggestions to regions facing Covid-19 and dengue fever co-infection.

Keywords:

Citation: Asma Hanif, Azhar Iqbal Kashif Butt, Tariq Ismaeel. Fractional optimal control analysis of Covid-19 and dengue fever co-infection model with Atangana-Baleanu derivative[J]. AIMS Mathematics, 2024, 9(3): 5171-5203. doi: 10.3934/math.2024251

Related Papers:

[1]	Sabila Ali, Shahid Mubeen, Rana Safdar Ali, Gauhar Rahman, Ahmed Morsy, Kottakkaran Sooppy Nisar, Sunil Dutt Purohit, M. Zakarya . Dynamical significance of generalized fractional integral inequalities via convexity. AIMS Mathematics, 2021, 6(9): 9705-9730. doi: 10.3934/math.2021565
[2]	Yousaf Khurshid, Muhammad Adil Khan, Yu-Ming Chu . Conformable integral version of Hermite-Hadamard-Fejér inequalities via η-convex functions. AIMS Mathematics, 2020, 5(5): 5106-5120. doi: 10.3934/math.2020328
[3]	Shuang-Shuang Zhou, Saima Rashid, Muhammad Aslam Noor, Khalida Inayat Noor, Farhat Safdar, Yu-Ming Chu . New Hermite-Hadamard type inequalities for exponentially convex functions and applications. AIMS Mathematics, 2020, 5(6): 6874-6901. doi: 10.3934/math.2020441
[4]	Maimoona Karim, Aliya Fahmi, Shahid Qaisar, Zafar Ullah, Ather Qayyum . New developments in fractional integral inequalities via convexity with applications. AIMS Mathematics, 2023, 8(7): 15950-15968. doi: 10.3934/math.2023814
[5]	Hari Mohan Srivastava, Soubhagya Kumar Sahoo, Pshtiwan Othman Mohammed, Bibhakar Kodamasingh, Kamsing Nonlaopon, Khadijah M. Abualnaja . Interval valued Hadamard-Fejér and Pachpatte Type inequalities pertaining to a new fractional integral operator with exponential kernel. AIMS Mathematics, 2022, 7(8): 15041-15063. doi: 10.3934/math.2022824
[6]	Muhammad Bilal Khan, Muhammad Aslam Noor, Thabet Abdeljawad, Bahaaeldin Abdalla, Ali Althobaiti . Some fuzzy-interval integral inequalities for harmonically convex fuzzy-interval-valued functions. AIMS Mathematics, 2022, 7(1): 349-370. doi: 10.3934/math.2022024
[7]	Yanping Yang, Muhammad Shoaib Saleem, Waqas Nazeer, Ahsan Fareed Shah . New Hermite-Hadamard inequalities in fuzzy-interval fractional calculus via exponentially convex fuzzy interval-valued function. AIMS Mathematics, 2021, 6(11): 12260-12278. doi: 10.3934/math.2021710
[8]	Sabir Hussain, Rida Khaliq, Sobia Rafeeq, Azhar Ali, Jongsuk Ro . Some fractional integral inequalities involving extended Mittag-Leffler function with applications. AIMS Mathematics, 2024, 9(12): 35599-35625. doi: 10.3934/math.20241689
[9]	Xiuzhi Yang, G. Farid, Waqas Nazeer, Muhammad Yussouf, Yu-Ming Chu, Chunfa Dong . Fractional generalized Hadamard and Fejér-Hadamard inequalities for m-convex functions. AIMS Mathematics, 2020, 5(6): 6325-6340. doi: 10.3934/math.2020407
[10]	Muhammad Bilal Khan, Pshtiwan Othman Mohammed, Muhammad Aslam Noor, Abdullah M. Alsharif, Khalida Inayat Noor . New fuzzy-interval inequalities in fuzzy-interval fractional calculus by means of fuzzy order relation. AIMS Mathematics, 2021, 6(10): 10964-10988. doi: 10.3934/math.2021637

Abstract

List of Abbreviations

The following abbreviations are used in this manuscript:

$\text{H-H} \quad \text{Hermite-Hadamard}\\ \text{H-H-M}\quad \text{Hermite-Hadamard-Mercer}\\ \text{H-H-F}\quad \text{Hermite-Hadamard-Fejér}$

1. Introduction and preliminaries

In science, convex functions have a long and distinguished history, and they have been the focus of study for almost a century. The rapid growth of convexity theory and applications of fractional calculus has kept the interest of a number of researchers on integral inequalities. Inequalities such as the H-H type, the Ostrowski-type, the H-H-M type, the Opial type, and other types, by using convex functions have been the focus of research for many years. The H-H inequality given in ^[1] has piqued the curiosity of most academics among all of these integral inequalities. Dragomir et al. ^[2] and Kirmaci et al. ^[3] presented some trapezoidal type inequalities and also some applications to special means. Following these articles, several mathematicians proposed new refinements of the Hermite-Hadamard inequality for various classes of convex functions and mappings such as quasi convex function ^[4], convex functions ^[5], $m$ -convex functions ^[6], $s$ -type convex functions of Raina type ^[7], $\sigma$ - $s$ -convex function ^[8] and harmonically convex functions ^[9]. Recently, this inequality was also investigated via different fractional integral operators, like Riemann-Liouville ^[10], $\psi$ -Riemann-Liouville ^[11], Proportional fractional ^[12,13], $k$ -Riemann-Liouville ^[14], Caputo-Fabrizio ^[15,16], generalized Atangana-Baleanu operator ^[17] to name a few.

It is important to emphasise that Leibniz and L'Hospital are credited with developing the idea of fractional calculus (1695). Other mathematicians, such as Riemann, Liouville, Letnikov, Erdéli, Grünwald, and Kober, have made significant inputs to the field of fractional calculus and its numerous applications. Many physical and engineering experts are interested in fractional calculus because of its behaviour and capacity to address a wide range of practical issues. Fractional calculus is currently concerned with the study of so-called fractional order integral and derivative functions over real and complex domains, as well as its applications. In many cases, fractional analysis requires the use of arithmetic from classical analysis to produce more accurate conclusions. Numerous mathematical models can be handled by differential equations of fractional order. Fractional mathematical models have more conclusive and precise results than classical mathematical models because they are particular examples of fractional order mathematical models. In classical analysis, integer orders do not serve as an adequate representation of nature. By using mathematical modelling, it is possible to identify the endemics' unique transmission dynamics and get insight into how infection impacts a new population. To enhance the actual phenomena to a higher degree of precision and accuracy, non-integer order fractional differential equations (FDEs) are applied. Additionally, ^{[18,19,20,21,22]} use and reference their utilization of fractional calculus. Other interesting results for fractional calclus can be found in ^[23,24,25]. However, fractional computation enables us to consider any number of orders and formulate far more measurable objectives. In recent years, mathematicians have become more and more interested in presenting well-known inequalities using a variety of novel theories of fractional integral operators. There are several different integral inequality results for fractional integrals. For generalizing significant and well-known integral inequalities, these operators are helpful. The Hermite-Hadamard integral inequality is a particular type of integral inequality. It is frequently used in the literature and outlines the necessary and sufficient conditions for a function to be convex. The Hermite-Hadamard inequalities were generalized by Sarikaya et al. ^[10] using Riemann-Liouville fractional integrals. Işcan ^[26] expanded Sarikaya et al. ^[10] 's findings to include Hermite-Hadamard-Fejer-type inequalities. By utilizing the product of two convex functions, Chen ^[27] produced fractional Hermite-Hadamard-type integral inequalities. Ögülmüs et al. ^[28] incorporated the Hermite-Hadamard and Jensen-Mercer inequalities to present Hermite-Hadamard-Mercer type inequalities for Riemann-Liouville fractional integrals. Motivated by the above articles, Butt et al. (see ^[29]), presented new versions of Jensen and Jensen-Mercer type inequalities in the fractal sense. New fractional versions of Hermite-Hadamard-Mercer and Pachpatte-Mercer type inclusions are established for convex ^[30] and harmonically convex functions ^[31] respectively. Latif et al. ^[32] established Hermite-Hadamard-Fejér type inequalities for convex harmonic and a positive symmetric increasing function. New refinements of Hermite-Mercer type inequalities are presented in ^[33], Mercer-Ostrowski type inequalities are presented in ^[34]. Further, the Hermite-Hadamard inequality is also generalized for convexity and quasi convexity ^[35] and differentiable convex functions ^[36]. For further information on other fractional-order integral inequalities, see the papers ^{[37,38,39,40,41]}.

Definition 1.1. (see ^[42]) Let $\mathscr{G}: \mathbb{X}\rightarrow \mathcal{R}$ be a function and $\mathbb{X}$ be a convex subset of a real vector space $\mathcal{R}$ . Then we say that the function $\mathscr{G}$ is convex if and only if the following condition:

$\begin{array}{l} \mathscr{G}\big( \Phi\mathfrak{r} +\left( 1-\Phi\right) \mathfrak{s} \big) \leq \Phi \mathscr{G}\left(\mathfrak{r}\right) +\left( 1-\Phi\right) \mathscr{G}\left( \mathfrak{s} \right), \end{array}$

holds true for all $\mathfrak{r}, \mathfrak{s}\in \mathbb{X}$ and $\Phi\in[0, 1].$

For further discussion, we first present the classical Hermite-Hadmard (H-H) inequality, which states that (see ^[1]):

If the function $\mathscr{G}:\mathbb{X}\subseteq{\mathcal{R}}\to{\mathcal{R}}$ is convex in $\mathbb{X}$ for $\mathfrak{r}, \mathfrak{s}\in \mathbb{X}$ and $\mathfrak{r} < \mathfrak{s},$ then

$\begin{array}{l} \mathscr{G}\biggl(\frac{\mathfrak{r} + \mathfrak{s}}{2}\biggl) \leq\frac{1}{\mathfrak{s} - \mathfrak{r}}\; \int_{\mathfrak{r}}^{\mathfrak{s}} \mathscr{G}(x) dx \leq \frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}. \end{array}$

(1.1)

Definition 1.2. (see ^[43]) Let there be a function $\mathscr{G}:[\mathfrak{r}, \mathfrak{s}]\to[0, \infty)$ and it is symmetric with respect to $\frac{\mathfrak{r}+\mathfrak{s}}{2}$ , if

$\begin{array}{l} \mathscr{G}\left(\mathfrak{r}+\mathfrak{s}-\mathrm{x}\right) = \mathscr{G}(\mathrm{x}). \end{array}$

In 1906, Fejér ^[44] preposed the following weighted variant of Hermite-Hadamard inequality famously known as Hermite-Hadamard-Fejér inequality, given as

Theorem 1.1. Let there be a convex function $\mathscr{G}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ with $\mathfrak{r} < \mathfrak{s}$ . If $\mathscr{D}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ be a convex symmetric and integrable function with respect to $\frac{\mathfrak{r}+\mathfrak{s}}{2}$ . Then

$\begin{array}{l} \mathscr{G}\biggl(\frac{\mathfrak{r} + \mathfrak{s}}{2}\biggl) \int_{\mathfrak{r}}^{\mathfrak{s}} \mathscr{D}(x) dx \leq\frac{1}{\mathfrak{s} - \mathfrak{r}}\; \int_{\mathfrak{r}}^{\mathfrak{s}} \mathscr{G}(x)\mathscr{D}(x) dx \leq \frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}\int_{\mathfrak{r}}^{\mathfrak{s}} \mathscr{D}(x) dx, \end{array}$

(1.2)

holds true.

Definition 1.3. (see ^[9]) Let $\mathscr{G}: \mathbb{X}\rightarrow \mathcal{R}$ be a function and $\mathbb{X}$ be a subset of a real vector space $\mathcal{R}$ . Then we say that the function $\mathscr{G}$ is harmonically convex if and only if the following condition

$\begin{array}{l} \mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{r} +\left( 1-\Phi\right) \mathfrak{s}} \right) \leq \Phi \mathscr{G}\left(\mathfrak{s}\right) +\left( 1-\Phi\right) \mathscr{G}\left( \mathfrak{r} \right), \end{array}$

holds true for all $\mathfrak{r}, \mathfrak{s}\in \mathbb{X}$ and $\Phi\in[0, 1].$

For further discussion, we first present the classical Hermite-Hadmard (H-H) inequality, which states that (see ^[9]):

If the function $\mathscr{G}:\mathbb{X}\subseteq{\mathcal{R}}\to{\mathcal{R}}$ is harmonically convex in $\mathbb{X}$ for $\mathfrak{r}, \mathfrak{s}\in \mathbb{X}$ and $\mathfrak{r} < \mathfrak{s},$ then

$\begin{array}{l} \mathscr{G}\biggl(\frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r} + \mathfrak{s}}\biggl) \leq\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s} - \mathfrak{r}}\; \int_{\mathfrak{r}}^{\mathfrak{s}} \frac{\mathscr{G}(x)}{x^{2}} dx \leq \frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}. \end{array}$

(1.3)

Definition 1.4. (see ^[45]) Let there be a function $\mathscr{G}:[\mathfrak{r}, \mathfrak{s}]\to[0, \infty)$ and it is harmonically symmetric with respect to $\frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}$ , if

$\begin{array}{l} \mathscr{G}\left({\Phi}\right) = \mathscr{G}\left(\frac{1}{\frac{1}{\mathfrak{r}}+\frac{1}{\mathfrak{s}}-\frac{1}{\Phi}}\right). \end{array}$

In the year 2014, Chen and Wu ^[46] proposed the following weighted variant of Hermite-Hadamard inequality for harmonically convex function, given as

Theorem 1.2. Let there be a convex function $\mathscr{G}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ with $\mathfrak{r} < \mathfrak{s}$ . If $\mathscr{D}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ be a convex symmetric and integrable function with respect to $\frac{\mathfrak{r}+\mathfrak{s}}{2}$ . Then,

$\begin{array}{l} \mathscr{G}\biggl(\frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r} + \mathfrak{s}}\biggl)\int_{\mathfrak{r}}^{\mathfrak{s}}\frac{\mathscr{D}(x)}{x^{2}} dx \leq\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s} - \mathfrak{r}}\; \int_{\mathfrak{r}}^{\mathfrak{s}} \frac{\mathscr{G}(x)\mathscr{D}(x)}{x^{2}} dx \leq \frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2} \int_{\mathfrak{r}}^{\mathfrak{s}}\frac{\mathscr{D}(x)}{x^{2}} dx, \end{array}$

(1.4)

holds true.

Definition 1.5. (see, for details, ^[10,47]; see also ^[48]) Let $\mathscr{G} \in \mathcal{L} \left\lbrack \mathfrak{r}, \ \mathfrak{s} \right\rbrack$. Then, the Riemann-Liouville fractional integrals of the order $\alpha > 0$ , are defined as follows:

$\begin{array}{l} {\mathtt{I}_{\mathfrak{r}^{+}}^{\alpha}{\mathscr{G}}\left( x \right)} = \frac{1}{\Gamma(\alpha)}\int_{\mathfrak{r}}^{x}\left( x - \mathfrak{m} \right)^{\alpha - 1} {\mathscr{G}\left( \mathfrak{m} \right)d\mathfrak{m}}\qquad \left(x > \mathfrak{r} \right), \end{array}$

and

$\begin{array}{l} {\mathtt{I}_{\mathfrak{s}^{-}}^{\alpha}{\mathscr{G}}\left( x \right)} = \frac{1}{\Gamma(\alpha)}\int_{x}^{\mathfrak{s}} \left( \mathfrak{m} - x \right)^{\alpha - 1}{\mathscr{G}\left( \mathfrak{m} \right)d\mathfrak{m}} \qquad (x < \mathfrak{s}), \end{array}$

respectively, where $\Gamma\left(\alpha \right) = \int_{0}^{\infty}{\Phi^{\alpha - 1}e^{- \Phi}}{d\Phi}$ is the Euler gamma function.

Definition 1.6. (see, ^[49] for details) Let ${\mathscr{G}}\in \mathcal{L}\left[\mathfrak{r}, \mathfrak{s}\right]$ . Then, the new left and right fractional integrals $\mathtt{I}_{\mathfrak{r}^{+}}^{\alpha}$ and $\mathtt{I}_{\mathfrak{s}^{-}}^{\alpha}$ of order $\alpha > 0$ are defined as

$\mathtt{I}_{\mathfrak{r}^{+}}^{\alpha}\mathscr{G}\left( x\right) : = \frac{1}{\alpha} \int_{\mathfrak{r}}^{x}e^{-\frac{1-\alpha}{\alpha}\left(x-\mathfrak{m}\right)}\; \mathscr{G}\left(\mathfrak{m}\right)d\mathfrak{m} \qquad (0 \leq \mathfrak{r} < x < \mathfrak{s}),$

and

$\mathtt{I}_{\mathfrak{s}^{-}}^{\alpha}\mathscr{G}\left(x\right) : = \frac{1}{\alpha} \int_{x}^{\mathfrak{s}}e^{-\frac{1-\alpha}{\alpha}\left(\mathfrak{m}-x\right)}\;\mathscr{G}\left(\mathfrak{m}\right)d\mathfrak{m} \qquad (0 \leq \mathfrak{r} < x < \mathfrak{s}),$

respectively.

It should be noted that

$\lim\limits_{\alpha\to 1} \mathtt{I}_{\mathfrak{r}^{+}}^{\alpha}\mathscr{G}\left( x\right) = \int_{\mathfrak{r}}^{x} \mathscr{G}\left(\mathfrak{m}\right)d\mathfrak{m} \ and \ \lim\limits_{\alpha\to 1} \mathtt{I}_{\mathfrak{s}^{-}}^{\alpha}\mathscr{G}\left(x\right) = \int_{x}^{\mathfrak{s}} \mathscr{G}\left(\mathfrak{m}\right)d\mathfrak{m}.$

Sarikaya et al. ^[37], in their article proved some interesting mid-point type Hermite-Hadamard inequalities. Here, we present one of his main results as follows:

Theorem 1.3. (see ^[37]) Let $\mathscr{G}:[\mathfrak{r}, \mathfrak{s}]\longrightarrow{\mathcal{R}}$ be a convex function with $0\leq\mathfrak{r}\leq\mathfrak{s}$ . If $\mathscr{G}\in\mathcal{L}[\mathfrak{r}, \mathfrak{s}]$ , then the following inequality for Riemann-Liouville fractional integral operator holds true:

$\mathscr{G}\left( \frac{\mathfrak{r}+\mathfrak{s}}{2}\right) \leq\frac{2^{\alpha -1}\Gamma(\alpha+1)}{(\mathfrak{s}-\mathfrak{r})^{\alpha}}\left[ \mathtt{I}_{(\frac {\mathfrak{r}+\mathfrak{s}}{2})^{+}}^{\alpha}\mathscr{G}(\mathfrak{s})+\mathtt{I}_{(\frac {\mathfrak{r}+\mathfrak{s}}{2})^{-}}^{\alpha}\mathscr{G}(\mathfrak{r})\right] \leq \frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}.$

The major goal of this paper is to establish Fejér type fractional inequalities using differintegrals of the $\left(\frac{\mathfrak{r}+\mathfrak{s}}{2}\right)$ type for both convex and harmonically convex functions via a novel fractional integral operator. In order to derive those inequalities, first we prove two new lemmas i.e., Lemmas 2.1 and 3.1 for convex and harmonic convex functions respectively.

In this study, we used a new fractional integral operator to achieve more generalized results. This is caused by the exponential function that makes up the kernel of this fractional operator. Our results differ from prior generalizations in that they do not lead to the aforementioned fractional integral inequalities. Numerous experts have suggested utilizing different fractional integral operators to extend the Hermite-Hadamard and Fejér type inequalities, however, none of their findings exhibit an exponential property. This study generated interest in using an exponential function as the kernel to create more generalized fractional inequalities. Furthermore, the application of symmetric and harmonically symmetric functions to the main results gives the study of inequalities a new path. For other generalization regarding exponential kernel interested reader can see e.g., on distributed-order fractional derivative in ^[50]. There are many research gaps to be filled for integral inequalities involving fractional calculus for different types of convex functions, despite the fact that there exist many different forms of research on the growth of fractional integral inequalities. As a result, the main purpose of this research is to find new Hermite-Hadamard and Fejér type inequalities for positive symmetric functions using fractional integral operators.

Our present investigation is structured as follows. In Sections 2 and 3, we discuss two additional characteristics of the relevant fractional operator before proving some enhanced versions (mid-point types) of the Fejér and Hermite-Hadamard type inequalities for convex and harmonically convex functions respectively. Some applications are also taken into consideration in Section 4 to determine whether the predetermined results are appropriate. Section 5 explores a brief conclusion and possible areas for additional research that is related to the findings in this paper are discussed in Section 6.

2. Improved Fejér type results for convex functions

In this section for simplicity, we denote $\rho_{c} = \frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r}).$ If $\alpha \to 1$ , then $\rho_{c} = \frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\to 0.$

Lemma 2.1. Let $\mathscr{D}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ be a symmetric convex function with respect to $\frac{\mathfrak{r}+\mathfrak{s}}{2}$ . Then for $\alpha > 0$ , the following equality holds true:

$\mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}\mathscr{D}(\mathfrak{s}) = \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}\mathscr{D}(\mathfrak{r}) = \frac{1}{2} \left[ \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}\mathscr{D}(\mathfrak{r})+\mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}\mathscr{D}(\mathfrak{s})\right].$

Proof. Since $\mathscr{D} :[\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ is integrable and symmetric to $\frac{\mathfrak{r}+\mathfrak{s}}{2}$ we have $\mathscr{D}(\mathfrak{r}+\mathfrak{s}-\mathrm{x}) = \mathscr{D}(\mathrm{x})$ . Also, Setting $\Phi = \mathfrak{r}+\mathfrak{s}-\mathrm{x}$ and $d\Phi = -d\mathrm{x}$ , we have

$\begin{align*} \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}\mathscr{D}(\mathfrak{s}) & = \frac{1}{\alpha}\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2}}^{\mathfrak{s}} e^{-\frac{1-\alpha}{\alpha}\left( \mathfrak{s}-\Phi\right) } \mathscr{D}(\Phi) d\Phi\\ & = - \frac{1}{\alpha}\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2}}^{\mathfrak{r}} e^{-\frac{1-\alpha}{\alpha}\left( \mathfrak{s}-(\mathfrak{r}+\mathfrak{s}-\mathrm{x})\right) } \mathscr{D}(\mathfrak{r}+\mathfrak{s}-\mathrm{x}) d\mathrm{x} \nonumber\\ & = \frac{1}{\alpha}\int_{\mathfrak{r}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}} e^{-\frac{1-\alpha}{\alpha}(\mathrm{x}-\mathfrak{r})} \mathscr{D}(\mathfrak{r}+\mathfrak{s}-\mathrm{x}) d\mathrm{x} \nonumber\\ & = \frac{1}{\alpha}\int_{\mathfrak{r}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}} e^{-\frac{1-\alpha}{\alpha}(\mathrm{x}-\mathfrak{r})} \mathscr{D}(\mathrm{x}) d\mathrm{x} \nonumber\\ & = \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}\mathscr{D}(\mathfrak{r}). \end{align*}$

$\implies \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}\mathscr{D}(\mathfrak{s}) = \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}\mathscr{D}(\mathfrak{r}) = \frac{1}{2} \left[ \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}\mathscr{D}(\mathfrak{r})+\mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}\mathscr{D}(\mathfrak{s})\right].$

This led us to the desired equality.

First, we prove both the first and second kind Fejér type inequalities in a different approach. Then, we also prove the Hermite-Hadamard inequality using symmetric convex functions.

Theorem 2.1. Let there be a convex function $\mathscr{G}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ with $\mathfrak{r} < \mathfrak{s}$ . If $\mathscr{D}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ be a convex symmetric and integrable function with respect to $\frac{\mathfrak{r}+\mathfrak{s}}{2}$ . Then for $\alpha > 0$ , the following inequality holds true:

$\begin{array}{l} \mathscr{G}\left( \frac{\mathfrak{r}+\mathfrak{s}}{2}\right)\left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}\mathscr{D}(\mathfrak{r})+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}\mathscr{D}(\mathfrak{s})\right]\leq \left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}(\mathscr{G}\mathscr{D})(\mathfrak{r})+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}(\mathscr{G}\mathscr{D})(\mathfrak{s})\right]. \end{array}$

Proof. Using the convexity of $\mathscr{G}$ on $[\mathfrak{r}, \mathfrak{s}]$ , we have

$\begin{array}{l} 2\mathscr{G}\left( \frac{\mathfrak{r}+\mathfrak{s}}{2}\right)\leq \mathscr{G}\left( \Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}\right)+\mathscr{G}\left( \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right). \end{array}$

(2.1)

Upon multiplication of both sides of the inequality (2.1) by $e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}\mathscr{D}\left(\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right)$ and then integrating the resultant over $\left[0, \frac{1}{2}\right]$ , we obtain

$\begin{array}{l} &2\mathscr{G}\left( \frac{\mathfrak{r}+\mathfrak{s}}{2}\right)\int_{0}^{\frac{1}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}\mathscr{D}\left( \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right)d\Phi\\ &\leq \int_{0}^{\frac{1}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}\mathscr{G}\left( \Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}\right)\mathscr{D}\left( \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right)d\Phi \\ &+\int_{0}^{\frac{1}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}\mathscr{G}\left( \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right)\mathscr{D}\left( \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right)d\Phi. \end{array}$

(2.2)

Since $\mathscr{D}$ is symmetric with respect to $\frac{\mathfrak{r}+\mathfrak{s}}{2}$ , we have $\mathscr{D}(\mathrm{x}) = \mathscr{D}(\mathfrak{r}+\mathfrak{s}-\mathrm{x})$ .

Moreover, setting $\mathrm{x} = \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}$ and $d\mathrm{x} = (\mathfrak{s}-\mathfrak{r})d\Phi$ in (2.2), we have

$\begin{align*} &2\mathscr{G}\left( \frac{\mathfrak{r}+\mathfrak{s}}{2}\right)\frac{1}{\mathfrak{s}-\mathfrak{r}}\int_{\mathfrak{r}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathrm{x}-\mathfrak{r})}\mathscr{D}(\mathrm{x})d\mathrm{x} = 2\mathscr{G}\left( \frac{\mathfrak{r}+\mathfrak{s}}{2}\right)\frac{\alpha}{\mathfrak{s}-\mathfrak{r}}\left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}\mathscr{D}(\mathfrak{r})\right] \nonumber\\ & \leq \frac{1}{\mathfrak{s}-\mathfrak{r}}\int_{\mathfrak{r}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathrm{x}-\mathfrak{r})}\mathscr{G}(\mathfrak{r}+\mathfrak{s}-\mathrm{x})\mathscr{D}(\mathrm{x})d\mathrm{x} + \frac{1}{\mathfrak{s}-\mathfrak{r}}\int_{\mathfrak{r}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathrm{x}-\mathfrak{r})}\mathscr{G}(\mathrm{x})\mathscr{D}(\mathrm{x})d\mathrm{x} \nonumber\\ & = \frac{1}{\mathfrak{s}-\mathfrak{r}}\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2}}^{\mathfrak{s}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathrm{x})}\mathscr{G}(\mathrm{x})\mathscr{D}(\mathfrak{r}+\mathfrak{s}-\mathrm{x})d\mathrm{x} + \frac{1}{\mathfrak{s}-\mathfrak{r}}\int_{\mathfrak{r}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathrm{x}-\mathfrak{r})}\mathscr{G}(\mathrm{x})\mathscr{D}(\mathrm{x})d\mathrm{x} \nonumber\\ & = \frac{\alpha}{\mathfrak{s}-\mathfrak{r}}\left[ \frac{1}{\alpha}\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2}}^{\mathfrak{s}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathrm{x})}\mathscr{G}(\mathrm{x})\mathscr{D}(\mathrm{x})d\mathrm{x} + \frac{1}{\alpha}\int_{\mathfrak{r}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathrm{x}-\mathfrak{r})}\mathscr{G}(\mathrm{x})\mathscr{D}(\mathrm{x})d\mathrm{x}\right] . \end{align*}$

It follows from the above developments and Lemma 2.1 that,

$\begin{array}{l} \frac{\alpha}{\mathfrak{s}-\mathfrak{r}}\mathscr{G}\left( \frac{\mathfrak{r}+\mathfrak{s}}{2}\right)\left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}\mathscr{D}(\mathfrak{r})+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}\mathscr{D}(\mathfrak{s})\right] \leq \frac{\alpha}{\mathfrak{s}-\mathfrak{r}}\left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}(\mathscr{G}\mathscr{D})(\mathfrak{r})+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}(\mathscr{G}\mathscr{D})(\mathfrak{s})\right]. \end{array}$

This concludes the proof of the required result.

Example 2.1. Let $\mathscr{G}\left(\mathfrak{m}\right) = e^{\mathfrak{m}}$ , $\mathfrak{m}\in\left[1, 9\right]$ and $\mathscr{D}\left(\mathfrak{m}\right) = \left(5-\mathfrak{m}\right) ^{2}$ , is non-negative symmetric about $\mathfrak{m} = 5$ . Let $0 < \alpha < 1$ , then

$\begin{align*} & \mathscr{G}\left( \frac{\mathfrak{r}+\mathfrak{s}}{2}\right) \left[ \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}^{\alpha}\mathscr{D} (\mathfrak{r})+\mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}^{\alpha} \mathscr{D}(\mathfrak{s})\right] \\ & = e^{5}\left[ \mathtt{I}_{5-}^{\alpha}\mathscr{D} (\mathit{1})+\mathtt{I}_{5+}^{\alpha}\mathscr{D}(\mathit{9})\right] \\ & = e^{5}\left[ \frac{1}{\alpha}\int_{1}^{\mathit{5}} e^{-\frac{1-\alpha}{\alpha}\left( \mathfrak{m}-1\right) }\;\left( 5-\mathfrak{m}\right) ^{2}d\mathfrak{m}+\frac{1}{\alpha}\int_{\mathit{5}} ^{9}e^{-\frac{1-\alpha}{\alpha}\left( 9-\mathfrak{m}\right) }\;\left( 5-\mathfrak{m}\right) ^{2}d\mathfrak{m}\right] . \end{align*}$

and

$\begin{align*} & \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}^{\alpha}(\mathscr{G} \mathscr{D})(\mathfrak{r})+\mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+} ^{\alpha}(\mathscr{G}\mathscr{D})(\mathfrak{s})\\ & = \mathtt{I}_{5-}^{\alpha}(\mathscr{G}\mathscr{D})(\mathit{1})+\mathtt{I} _{5+}^{\alpha}(\mathscr{G}\mathscr{D})(\mathit{9})\\ & = \frac{1}{\alpha}\int_{1}^{\mathit{5}}e^{-\frac{1-\alpha}{\alpha}\left( \mathfrak{m}-1\right) }\;e^{\mathfrak{m}}\left( 5-\mathfrak{m}\right) ^{2}d\mathfrak{m}+ \frac{1}{\alpha}\int_{\mathit{5}}^{9}e^{-\frac{1-\alpha}{\alpha}\left( 9-\mathfrak{m}\right) }\;e^{\mathfrak{m}}\left( 5-\mathfrak{m}\right) ^{2}d\mathfrak{m}. \end{align*}$

The graphical representation of Theorem 2.1 is shown in the graph below (see ) for $0 < \alpha < 1$ :

Figure 1. The graphical representation of Theorem 2.1 for

$0 < \alpha < 1$ .

DownLoad: Full-Size Img PowerPoint

Theorem 2.2. Let there be a convex function $\mathscr{G}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ with $\mathfrak{r} < \mathfrak{s}$ . If $\mathscr{D}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ be a convex symmetric and integrable function with respect to $\frac{\mathfrak{r}+\mathfrak{s}}{2}$ . Then for $\alpha > 0$ , the following inequality holds true:

$\begin{array}{l} \left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}(\mathscr{G}\mathscr{D})(\mathfrak{r})+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}(\mathscr{G}\mathscr{D})(\mathfrak{s})\right] \leq \frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}\left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}(\mathscr{D})(\mathfrak{r})+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}(\mathscr{D})(\mathfrak{s})\right]. \end{array}$

Proof. Since $\mathscr{G}$ is convex function, we have

$\begin{array}{l} \mathscr{G}\left( \Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}\right)+\mathscr{G}\left( \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right) \leq \mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s}). \end{array}$

(2.3)

Multiplying both side of the above inequality (2.3) by $e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}\mathscr{D}\left(\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right)$ and upon integration of the obtained result over $\left[0, \frac{1}{2}\right]$ , one has

$\begin{align*} &\int_{0}^{\frac{1}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}\mathscr{G}\left( \Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}\right) \mathscr{D}\left( \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right)d\Phi\\ &+\int_{0}^{\frac{1}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}\mathscr{G}\left( \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right)\mathscr{D}\left( \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right) d\Phi \\ &\leq [\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})]\int_{0}^{\frac{1}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}\mathscr{D}\left( \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right) d\Phi. \end{align*}$

It follows

$\begin{array}{l} \frac{\alpha}{\mathfrak{s}-\mathfrak{r}}\left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}(\mathscr{G}\mathscr{D})(\mathfrak{r})+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}(\mathscr{G}\mathscr{D})(\mathfrak{s})\right] \leq \alpha\frac{[\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})]}{\mathfrak{s}-\mathfrak{r}}\left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}(\mathscr{D})(\mathfrak{r})\right]. \end{array}$

Furthermore, using the Lemma 2.1, we obtain

$\begin{array}{l} \frac{\alpha}{\mathfrak{s}-\mathfrak{r}}\left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}(\mathscr{G}\mathscr{D})(\mathfrak{r})+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}(\mathscr{G}\mathscr{D})(\mathfrak{s})\right] \leq \frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}\frac{\alpha}{\mathfrak{s}-\mathfrak{r}}\left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}(\mathscr{D})(\mathfrak{r})+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}(\mathscr{D})(\mathfrak{s})\right]. \end{array}$

This concludes the proof of the desired result.

Example 2.2. Let $\mathscr{G}\left(\mathfrak{m}\right) = e^{\mathfrak{m}}$ , $\mathfrak{m}\in\left[1, 9\right]$ and $\mathscr{D}\left(\mathfrak{m}\right) = \left(5-\mathfrak{m}\right) ^{2}$ , is non-negative symmetric about $\mathfrak{m} = 5$ . Let $0 < \alpha < 1$ , then

$\begin{align*} & \frac{e+e^{9}}{2}\left[ \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2} -}^{\alpha}(\mathscr{D})(\mathfrak{r})+\mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s} }{2}+}^{\alpha}(\mathscr{D})(\mathfrak{s})\right] \\ & = \frac{e+e^{9}}{2}\left[ \mathtt{I}_{5-}^{\alpha}\mathscr{D} (\mathit{1})+\mathtt{I}_{5+}^{\alpha}\mathscr{D}(\mathit{9})\right] \\ & = \frac{e+e^{9}}{2}\left[ \frac{1}{\alpha}\int_{1}^{\mathit{5}} e^{-\frac{1-\alpha}{\alpha}\left( \mathfrak{m}-1\right) }\;\left( 5-\mathfrak{m}\right) ^{2}d\mathfrak{m}+\frac{1}{\alpha}\int_{\mathit{5}} ^{9}e^{-\frac{1-\alpha}{\alpha}\left( 9-\mathfrak{m}\right) }\;\left( 5-\mathfrak{m}\right) ^{2}d\mathfrak{m}\right] . \end{align*}$

And

The graphical representation of Theorem 2.2 is shown in the graph below (see ) for $0 < \alpha < 1$ :

Figure 2. The graphical representation of Theorem 2.2 for

$0 < \alpha < 1$ .

DownLoad: Full-Size Img PowerPoint

Theorem 2.3. Let there be a convex function $\mathscr{G}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ with $\mathfrak{r} < \mathfrak{s}$ . Then for $\alpha > 0$ , the following fractional integral inequality holds true:

$\begin{array}{l} \mathscr{G}\left( \frac{\mathfrak{r}+\mathfrak{s}}{2}\right) \leq\frac{1-\alpha}{2(1-e^{-\frac{\rho_{c}}{2}})}\left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}\mathscr{G}(\mathfrak{r})+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}\mathscr{G}(\mathfrak{s})\right] \leq \frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}. \end{array}$

(2.4)

Proof. By the hypothesis of convexity, we have

$\begin{array}{l} \mathscr{G}\left( \frac{\mathfrak{r}+\mathfrak{s}}{2}\right) & = \mathscr{G}\left( \frac{\Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}+\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}{2}\right)\\ &\leq \frac{\mathscr{G}\left( \Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}\right)+\mathscr{G}\left( \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right) }{2}. \end{array}$

(2.5)

Upon multiplication of both sides of the inequality (2.5) by $2e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}$ and then integrating the obtained result over $\left[0, \frac{1}{2}\right]$ , we have

$\begin{array}{l} &2\mathscr{G}\left( \frac{\mathfrak{r}+\mathfrak{s}}{2}\right)\int_{0}^{\frac{1}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}d\Phi\\ & \leq \int_{0}^{\frac{1}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}\mathscr{G}\left( \Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}\right)d\Phi+\int_{0}^{\frac{1}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}\mathscr{G}\left( \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right) d\Phi. \end{array}$

(2.6)

Furthermore, let $\mathfrak{m} = \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\implies d\Phi = \frac{d\mathfrak{m}}{\mathfrak{s}-\mathfrak{r}}$ . Then inequality (2.6) gives

$\begin{array}{l} &\frac{2(1-e^{-\frac{\rho_{c}}{2}})}{\rho}\mathscr{G}\left( \frac{\mathfrak{r}+\mathfrak{s}}{2}\right)\\ &\leq \left[ \frac{1}{\mathfrak{s}-\mathfrak{r}}\int_{\mathfrak{r}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\frac{\mathfrak{m}-\mathfrak{r}}{\mathfrak{s}-\mathfrak{r}}}\mathscr{G}\left( \mathfrak{r}+\mathfrak{s}-\mathfrak{m}\right)d\mathfrak{m} +\frac{1}{\mathfrak{s}-\mathfrak{r}}\int_{\mathfrak{r}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\frac{\mathfrak{m}-\mathfrak{r}}{\mathfrak{s}-\mathfrak{r}}}\mathscr{G}\left( \mathfrak{m}\right)d\mathfrak{m}\right] \\ & = \frac{\alpha}{\mathfrak{s}-\mathfrak{r}}\left[\frac{1}{\alpha}\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2}}^{\mathfrak{s}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{m})}\mathscr{G}\left( \mathfrak{m}\right)d\mathfrak{m} +\frac{1}{\alpha}\int_{\mathfrak{r}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{m}-\mathfrak{r})}\mathscr{G}\left( \mathfrak{m}\right)d\mathfrak{m} \right] \\ & = \frac{\alpha}{\mathfrak{s}-\mathfrak{r}}\left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}\mathscr{G}(\mathfrak{r})+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}\mathscr{G}(\mathfrak{s})\right]. \end{array}$

(2.7)

This concludes the proof of the first part of the inequality (2.4). To prve the next part of inequality, under the given hypothesis, we have

(2.8)

Upon multiplication of both sides of the inequality (2.8) by $e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}$ and integrating over $\left[0, \frac{1}{2}\right]$ , we have

$\begin{align*} &\int_{0}^{\frac{1}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}\mathscr{G}\left( \Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}\right)d\Phi +\int_{0}^{\frac{1}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}\mathscr{G}\left( \Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}\right) d\Phi\\ &\leq [\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})]\int_{0}^{\frac{1}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{r})\Phi}d\Phi. \end{align*}$

Consequently, we obtain

$\begin{array}{l} \frac{\alpha}{\mathfrak{s}-\mathfrak{r}}\left[ \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}\mathscr{G}(\mathfrak{r})+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}\mathscr{G}(\mathfrak{s})\right] \leq \frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}\frac{2(1-e^{-\frac{\rho_{c}}{2}})}{\rho_{c}}. \end{array}$

(2.9)

Consequently, it follows from the above developments (2.7) and (2.9) that

This concludes the proof of the required result.

Remark 2.1. If one chooses $\alpha\to 1$ i.e., $\frac{\rho_{c}}{2}\to 0$ , then

$\lim\limits_{\alpha\to 1}\frac{1-\alpha}{2(1-e^{-\frac{\rho_{c}}{2}})} = \frac{1}{\mathfrak{s}-\mathfrak{r}}$

and hence Theorem 2.3 retrieves the classical Hermite-Hadamard inequality (1.1).

Example 2.3. Let $\mathscr{G}\left(\mathfrak{m}\right) = e^{\mathfrak{m}}$ , $\mathfrak{m}\in\left[1, 9\right]$ and $0 < \alpha < 1$ , then

$\frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2} = \frac{e+e^{9}}{2},$

$\mathscr{G}\left( \frac{\mathfrak{r}+\mathfrak{s}}{2}\right) = e^{{5}}$

and

$\begin{align*} & \frac{1-\alpha}{2(1-e^{-\frac{\rho_{c}}{2}})}\left[ \mathtt{I} _{\frac{\mathfrak{r}+\mathfrak{s}}{2}-}^{\alpha}\mathscr{G}(\mathfrak{r} )+\mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}+}^{\alpha}\mathscr{G} (\mathfrak{s})\right] \\ & = \frac{1-\alpha}{2\left( 1-e^{-\frac{4\left( 1-\alpha\right) }{\alpha} }\right) }\left[ \mathtt{I}_{5-}^{\alpha}(\mathscr{G})(\mathit{1} )+\mathtt{I}_{5+}^{\alpha}(\mathscr{G})(\mathit{9})\right] \\ & = \frac{1-\alpha}{2\left( 1-e^{-\frac{4\left( 1-\alpha\right) }{\alpha} }\right) }\left[ \frac{1}{\alpha}\int_{1}^{\mathit{5}}e^{-\frac{1-\alpha }{\alpha}\left( \mathfrak{m}-1\right) }\;e^{\mathfrak{m}}d\mathfrak{m}+\frac{1}{\alpha} \int_{\mathit{5}}^{9}e^{-\frac{1-\alpha}{\alpha}\left( 9-\mathfrak{m}\right) }\;e^{\mathfrak{m}}d\mathfrak{m}\right] . \end{align*}$

The graphical representation of Theorem 2.3 is shown in the graph below (see ) for $0 < \alpha < 1$ :

Figure 3. The graphical representation of Theorem 2.3 for

$0 < \alpha < 1$ .

DownLoad: Full-Size Img PowerPoint

3. Improved inclusions via harmonically convex functions

The family of Lebesgue measurable functions is represented here by $\mathcal{L} \left\lbrack \mathfrak{r}, \mathfrak{s} \right\rbrack $. In this section, for brevity we use, $\rho_{h} = \frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}$ wherever needed. If $\alpha \to 1$ , then $\rho_{h} = \frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\to 0.$

Lemma 3.1. Let $\mathscr{D}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ be a harmonically symmetric and integrable function with respect to $\frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}$ . Then for $\alpha > 0$ , the following equality holds true:

$\begin{array}{l} \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right) = \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right) = \frac{1}{2}\left[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right)\right], \end{array}$

where $\mathtt{K}(\mathrm{x}) = \frac{1}{\mathrm{x}}, \; \; \mathrm{x}\in\; \; \left[\frac{1}{\mathfrak{s}}, \frac{1}{\mathfrak{r}}\right]$ .

Proof. Let $\mathscr{D}$ be a harmonically symmetric function with respect to $\frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}$ . Then using the harmonically symmetric property of $\mathscr{D}$ , given as $\mathscr{D}\left(\frac{1}{\Phi}\right) = \mathscr{D}\left(\frac{1}{\frac{1}{\mathfrak{r}}+\frac{1}{\mathfrak{s}}-\Phi}\right)$ for $\alpha > 0$ .

$\begin{align*} \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right) & = \frac{1}{\alpha}\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}}^{\frac{1}{\mathfrak{r}}}\mathrm{e}^{-\frac{1-\alpha}{\alpha}\left(\frac{1}{\mathfrak{r}}-\Phi \right) }\mathscr{D}\left( \frac{1}{\Phi}\right) d\Phi \nonumber\\ & = \frac{1}{\alpha}\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}}^{\frac{1}{\mathfrak{r}}}\mathrm{e}^{-\frac{1-\alpha}{\alpha}\left(\frac{1}{\mathfrak{r}}-\Phi \right) }\mathscr{D}\left( \frac{1}{\frac{1}{\mathfrak{r}}+\frac{1}{\mathfrak{s}}-\Phi}\right) d\Phi \nonumber\\ & = -\frac{1}{\alpha}\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}}^{\frac{1}{\mathfrak{s}}}\mathrm{e}^{-\frac{1-\alpha}{\alpha}\left(\mathrm{x}-\frac{1}{\mathfrak{s}}\right) }\mathscr{D}\left( \frac{1}{\mathrm{x}}\right) d\mathrm{x} \nonumber\\ & = \frac{1}{\alpha}\int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}}\mathrm{e}^{-\frac{1-\alpha}{\alpha}\left(\mathrm{x}-\frac{1}{\mathfrak{s}}\right) }\mathscr{D}\left( \frac{1}{\mathrm{x}}\right) d\mathrm{x} \nonumber\\ & = \mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right). \end{align*}$

Consequently, it follows from the above developments that

where $\mathtt{K}(\mathrm{x}) = \frac{1}{\mathrm{x}}, \; \; \mathrm{x}\in\; \; \left[\frac{1}{\mathfrak{s}}, \frac{1}{\mathfrak{r}}\right].$

Now, we use the above result to produce new Hadamard-Fejér type inequalities of both first and second kind for harmonically convex functions.

Let us begin with the Hadamard-Fejér type inequality of the first kind.

Theorem 3.1. Let there be a harmonically convex function $\mathscr{G}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ with $\mathfrak{r} < \mathfrak{s}$ . If $\mathscr{D}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ be a harmonically symmetric and integrable function with respect to $\frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}$ . Then for $\alpha > 0$ , the following inequality holds true:

$\begin{align*} \mathscr{G}\left( \frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}\right)\bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right)\bigg]\\ \leq \bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right)\bigg]. \end{align*}$

Proof. Since $\mathscr{G}$ is harmonically convex function on $[\mathfrak{r}, \mathfrak{s}]$ , we have

$\mathscr{G}\left( \frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}\right) \leq \frac{\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}}\right)+\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)}{2}.$

Multiplying both side by $2e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}\mathscr{D}\left(\frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)$ and then integrating over $[0, \frac{1}{2}]$ , we find

$\begin{align*} &2\mathscr{G}\left( \frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}\right) \int_{0}^{\frac{1}{2}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}\mathscr{D}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)d\Phi \nonumber\\ &\leq \int_{0}^{\frac{1}{2}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}\mathscr{D}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}}\right)d\Phi\nonumber\\ &+ \int_{0}^{\frac{1}{2}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}\mathscr{D}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)d\Phi. \end{align*}$

Since, $\mathscr{D}$ is harmonically symmetric with respect to $\frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}$ i.e

$\mathscr{D}\left( \frac{1}{\mathrm{x}}\right) = \mathscr{D}\left( \frac{1}{\frac{1}{\mathfrak{r}}+\frac{1}{\mathfrak{s}}-\mathrm{x}}\right).$

Also, setting $\mathrm{x} = \frac{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}$ $\implies d\Phi = \frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}d\mathrm{x}$ the above developments proceed as follows:

$\alpha \frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\mathscr{G}\left( \frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}\right)\frac{1}{\alpha}\int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\left( \mathrm{x}-\frac{1}{\mathfrak{s}}\right) } \mathscr{D}\left( \frac{1}{\mathrm{x}}\right) d\mathrm{x}\nonumber\\ \leq \alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}} \bigg[\frac{1}{\alpha}\int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\left( \mathrm{x}-\frac{1}{\mathfrak{s}}\right) } \mathscr{G}\left( \frac{1}{\frac{1}{\mathfrak{r}}+\frac{1}{\mathfrak{s}}-\mathrm{x}}\right)\mathscr{D}\left( \frac{1}{\mathrm{x}}\right) d\mathrm{x} +\frac{1}{\alpha} \int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\left( \mathrm{x}-\frac{1}{\mathfrak{s}}\right) } \mathscr{G}\left( \frac{1}{\mathrm{x}}\right)\mathscr{D}\left( \frac{1}{\mathrm{x}}\right) d\mathrm{x}\bigg] \nonumber\\ = \alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\bigg[\frac{1}{\alpha}\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}}^{\frac{1}{\mathfrak{r}}} e^{-\frac{1-\alpha}{\alpha}\left( \frac{1}{\mathfrak{r}}-\mathrm{x}\right) } \mathscr{G}\left(\frac{1}{\mathrm{x}} \right)\mathscr{D}\left(\frac{1}{\frac{1}{\mathfrak{r}}+\frac{1}{\mathfrak{s}}-\mathrm{x}} \right) d\mathrm{x} +\frac{1}{\alpha}\int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\left( \mathrm{x}-\frac{1}{\mathfrak{s}}\right) } \mathscr{G}\left( \frac{1}{\mathrm{x}}\right)\mathscr{D}\left( \frac{1}{\mathrm{x}}\right) d\mathrm{x}\bigg] \nonumber\\ = \alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\bigg[\frac{1}{\alpha}\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}}^{\frac{1}{\mathfrak{r}}} e^{-\frac{1-\alpha}{\alpha}\left( \frac{1}{\mathfrak{r}}-\mathrm{x}\right) } \mathscr{G}\left(\frac{1}{\mathrm{x}} \right)\mathscr{D}\left(\frac{1}{\mathrm{x}} \right) d\mathrm{x} +\frac{1}{\alpha}\int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\left( \mathrm{x}-\frac{1}{\mathfrak{s}}\right) } \mathscr{G}\left( \frac{1}{\mathrm{x}}\right)\mathscr{D}\left( \frac{1}{\mathrm{x}}\right) d\mathrm{x}\bigg] \nonumber\\ = \alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right)\bigg].$

From the above developments and Lemma 3.1, we have

$\begin{align*} &\alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\mathscr{G}\left( \frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}\right)\bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right)\bigg]\\ &\leq \alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right)\bigg]. \end{align*}$

This concludes the proof of the required result.

Example 3.1. Let $\mathscr{G}\left(\mathfrak{m}\right) = \mathfrak{m}^{2}$ , $\mathfrak{m}\in\left[1, 4\right]$ , $\mathscr{D}\left(\mathfrak{m}\right) = \left(\frac{5\mathfrak{m}-8}{8\mathfrak{m}}\right) ^{2}$ and $0 < \alpha < 1$ , then

$\begin{align*} & \mathscr{G}\left( \frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s} }\right) \left[ \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r} \mathfrak{s}}-}^{\alpha}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s} }\right) +\mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}} +}^{\alpha}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right) \right] \\ & = \frac{64}{25}\left[ \frac{1}{\alpha}\int_{\frac{1}{4}}^{\frac{\mathit{5} }{8}}e^{-\frac{1-\alpha}{\alpha}\left( \mathfrak{m}-\frac{\mathit{1}} {4}\right) }\;\left( \frac{5-8\mathfrak{m}}{8}\right) ^{2}d\mathfrak{m}+\frac{1}{\alpha }\int_{\frac{\mathit{5}}{8}}^{1}e^{-\frac{1-\alpha}{\alpha}\left( 1-\mathfrak{m}\right) }\;\left( \frac{5-8\mathfrak{m}}{8}\right) ^{2}d\mathfrak{m}\right] \end{align*}$

and

$\begin{align*} & \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}^{\alpha }\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right) +\mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}^{\alpha }\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right) \\ & = \mathtt{I}_{\frac{\mathit{5}}{8}-}^{\alpha}\mathscr{G}\mathscr{D} \circ\mathtt{K}\left( \frac{1}{\mathit{4}}\right) +\mathtt{I}_{\frac {\mathit{5}}{8}+}^{\alpha}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( 1\right) \\ &{ = \frac{1}{\alpha}\int_{\frac{1}{4}}^{\frac{\mathit{5}}{8}}e^{-\frac {1-\alpha}{\alpha}\left( \mathfrak{m}-\frac{\mathit{1}}{4}\right) }\;\left(\frac{1}{\mathfrak{m}}\right)^{2}\left( \frac{5-8\mathfrak{m}}{8}\right) ^{2}d\mathfrak{m}+\frac{1}{\alpha}\int_{\frac{\mathit{5} }{8}}^{1}e^{-\frac{1-\alpha}{\alpha}\left( 1-\mathfrak{m}\right) }\;\left(\frac{1}{\mathfrak{m}}\right)^{2}\left( \frac{5-8\mathfrak{m}}{8}\right) ^{2}d\mathfrak{m}.} \end{align*}$

The graphical representation of Theorem 3.1 is shown in the graph below (see ) for $0 < \alpha < 1$ :

Figure 4. The graphical representation of Theorem 3.1 for

$0 < \alpha < 1$ .

DownLoad: Full-Size Img PowerPoint

Now, we will establish the Fejér type inequality of the second kind.

Theorem 3.2. Let there be a harmonically convex function $\mathscr{G}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ with $\mathfrak{r} < \mathfrak{s}$ . If $\mathscr{D}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ be a harmonically symmetric and integrable function with respect to $\frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}$ . Then for $\alpha > 0$ , the following inequality holds true:

$\begin{align*} & \bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right)\bigg]\\ &\leq\frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}\bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right)\bigg]. \end{align*}$

Proof. Since $\mathscr{G}$ is harmonically convex function

$\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}}\right)+\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right) \leq\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s}).$

Multiplying both side by $e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}\mathscr{D}\left(\frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)$ and then integrating the resultant over $\left[0, \frac{1}{2}\right]$ , we find

$\begin{array}{l} &\int_{0}^{\frac{1}{2}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}}\right)\mathscr{D}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)d\Phi\\ &+\int_{0}^{\frac{1}{2}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)\mathscr{D}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)d\Phi \\ &\leq [\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})]\int_{0}^{\frac{1}{2}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}\mathscr{D}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)d\Phi. \end{array}$

(3.1)

Setting $\mathrm{x} = \frac{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}$ and $\mathscr{D}\left(\frac{1}{\mathrm{x}}\right) = \mathscr{D}\left(\frac{1}{\frac{1}{\mathfrak{r}}+\frac{1}{\mathfrak{s}}-\mathrm{x}} \right)$ in (3.1), we have

$\begin{array}{l} &\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}} \bigg[\int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\left( \mathrm{x}-\frac{1}{\mathfrak{s}}\right) } \mathscr{G}\left( \frac{1}{\frac{1}{\mathfrak{r}}+\frac{1}{\mathfrak{s}}-\mathrm{x}}\right)\mathscr{D}\left( \frac{1}{\mathrm{x}}\right) d\mathrm{x} \int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\left( \mathrm{x}-\frac{1}{\mathfrak{s}}\right) } \mathscr{G}\left( \frac{1}{\mathrm{x}}\right)\mathscr{D}\left( \frac{1}{\mathrm{x}}\right) d\mathrm{x}\bigg] \\ & = \alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\bigg[\frac{1}{\alpha}\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}}^{\frac{1}{\mathfrak{r}}} e^{-\frac{1-\alpha}{\alpha}\left( \frac{1}{\mathfrak{r}}-\mathrm{x}\right) } \mathscr{G}\left(\frac{1}{\mathrm{x}} \right)\mathscr{D}\left(\frac{1}{\mathrm{x}} \right) d\mathrm{x} +\frac{1}{\alpha}\int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\left( \mathrm{x}-\frac{1}{\mathfrak{s}}\right) } \mathscr{G}\left( \frac{1}{\mathrm{x}}\right)\mathscr{D}\left( \frac{1}{\mathrm{x}}\right) d\mathrm{x}\bigg] \\ & = \alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right)\bigg]. \end{array}$

(3.2)

Also,

$\begin{array}{l} &[\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})]\int_{0}^{\frac{1}{2}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}\mathscr{D}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)d\Phi\\ & = \alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}[\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})]\frac{1}{\alpha}\int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\left( \mathrm{x}-\frac{1}{\mathfrak{s}}\right) } \mathscr{D}\left(\frac{1}{\mathrm{x}} \right) d\mathrm{x} \\ & = \alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}[{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}]\bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right)\bigg]. \end{array}$

(3.3)

From the above developments (3.2), (3.3) and Lemma 3.1, we have

$\begin{align*} &\alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right)\bigg]\\ &\leq\alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}\bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right)\bigg]. \end{align*}$

This concludes the proof of the required result.

Example 3.2. Let $\mathscr{G}\left(\mathfrak{m}\right) = \mathfrak{m}^{2}$ , $\mathfrak{m}\in\left[1, 4\right]$ , $\mathscr{D}\left(\mathfrak{m}\right) = \left(\frac{5\mathfrak{m}-8}{8\mathfrak{m}}\right) ^{2}$ and $0 < \alpha < 1$ , then

$\begin{align*} & \frac{17}{2}\left[ \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r} \mathfrak{s}}-}^{\alpha}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s} }\right) +\mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}} +}^{\alpha}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right) \right] \\ & = \frac{17}{2}\left[ \frac{1}{\alpha}\int_{\frac{1}{4}}^{\frac{\mathit{5} }{8}}e^{-\frac{1-\alpha}{\alpha}\left( \mathfrak{m}-\frac{1}{4}\right) }\;\left( \frac{5-8\mathfrak{m}}{8}\right) ^{2}d\mathfrak{m}+\frac{1}{\alpha}\int _{\frac{\mathit{5}}{8}}^{1}e^{-\frac{1-\alpha}{\alpha}\left( 1-\mathfrak{m} \right) }\;\left( \frac{5-8\mathfrak{m}}{8}\right) ^{2}d\mathfrak{m}\right] \end{align*}$

and

$\begin{align*} & \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}^{\alpha }\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{s}}\right) +\mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}^{\alpha }\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( \frac{1}{\mathfrak{r}}\right) \\ & = \mathtt{I}_{\frac{\mathit{5}}{8}-}^{\alpha}\mathscr{G}\mathscr{D} \circ\mathtt{K}\left( \frac{1}{\mathit{4}}\right) +\mathtt{I}_{\frac {\mathit{5}}{8}+}^{\alpha}\mathscr{G}\mathscr{D}\circ\mathtt{K}\left( 1\right) \\ & = \frac{1}{\alpha}\int_{\frac{1}{4}}^{\frac{\mathit{5}}{8}}e^{-\frac {1-\alpha}{\alpha}\left( \mathfrak{m}-\frac{1}{4}\right) }\;\left( \frac{5-8\mathfrak{m}}{8}\right) ^{4}d\mathfrak{m}+\frac{1}{\alpha}\int_{\frac{\mathit{5}}{8} }^{1}e^{-\frac{1-\alpha}{\alpha}\left( 1-\mathfrak{m}\right) }\;\left( \frac{5-8\mathfrak{m}}{8}\right) ^{4}d\mathfrak{m}. \end{align*}$

The graphical representation of Theorem 3.2 is shown in the graph below (see ) for $0 < \alpha < 1$ :

Figure 5. The graphical representation of Theorem 3.2 for

$0 < \alpha < 1$ .

DownLoad: Full-Size Img PowerPoint

Theorem 3.3. Let there be a harmonically convex function $\mathscr{G}: [\mathfrak{r}, \mathfrak{s}]\subseteq \mathcal{R}\setminus\{0\}\to\mathcal{R}$ with $\mathfrak{r} < \mathfrak{s}$ . Then for $\alpha > 0$ ,

$\begin{array}{l} \mathscr{G}\left( \frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}\right)\leq\frac{(1-\alpha)}{2\left( 1-e^{\frac{\rho_{h}}{2}}\right) }\bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{G}\circ\mathscr{K}\left( \frac{1}{\mathfrak{r}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{G}\circ\mathscr{K}\left( \frac{1}{\mathfrak{s}}\right)\bigg]\leq\frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}, \end{array}$

(3.4)

holds true.

Proof. Since $\mathscr{G}$ is harmonically convex function on $[\mathfrak{r}, \mathfrak{s}]$ , we have

$\begin{array}{l} \mathscr{G}\left( \frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}\right) = \mathscr{G}\left( \frac{2\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}}\right)\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right) }{\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}}\right)+\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)}\right) \leq \frac{\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}}\right)+\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)}{2}{.} \end{array}$

(3.5)

Multiplying both side of the inequality (3.5) by $2e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}$ and integrating over $[0, \frac{1}{2}]$ , we obtain

$\begin{align*} 2\mathscr{G}\left( \frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}\right) \int_{0}^{1/2}e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}d\Phi& \leq \int_{0}^{1/2}e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi} \mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}}\right) d\Phi\\ &+\int_{0}^{1/2}e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi} \mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right) d\Phi. \end{align*}$

Let $\mathfrak{m} = \frac{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}$ , then $d\mathfrak{m} = \frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}d\Phi$

$\begin{array}{l} &\frac{2(1-e^{-\frac{\rho_{h}}{2}})}{\rho_{h}}\mathscr{G}\left( \frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}\right)\\ &\leq \int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\left(\mathfrak{m}-\frac{1}{\mathfrak{s}} \right) } \left( \frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\right) \mathscr{G}\left( \frac{1}{\frac{1}{\mathfrak{r}}+\frac{1}{\mathfrak{s}}-\mathfrak{m}}\right) d\mathfrak{m}+ \int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\left(\mathfrak{m}-\frac{1}{\mathfrak{s}} \right) } \left( \frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\right) \mathscr{G}\left( \frac{1}{\mathfrak{m}}\right) d\mathfrak{m} \\ & = \frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\bigg[\int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\left(\mathfrak{m}-\frac{1}{\mathfrak{s}} \right) } \mathscr{G}\left( \frac{1}{\frac{1}{\mathfrak{r}}+\frac{1}{\mathfrak{s}}-\mathfrak{m}}\right) d\mathfrak{m} +\int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\left(\mathfrak{m}-\frac{1}{\mathfrak{s}} \right) } \mathscr{G}\left( \frac{1}{\mathfrak{m}}\right) d\mathfrak{m}\bigg] \\ & = \frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\bigg[\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}}^{\frac{1}{\mathfrak{r}}} e^{-\frac{1-\alpha}{\alpha}\left(\frac{1}{\mathfrak{r}} -\mathfrak{m}\right) } \mathscr{G}\left( \frac{1}{\mathfrak{m}}\right) d\mathfrak{m} +\int_{\frac{1}{\mathfrak{s}}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}} e^{-\frac{1-\alpha}{\alpha}\left(\mathfrak{m}-\frac{1}{\mathfrak{s}} \right) } \mathscr{G}\left( \frac{1}{\mathfrak{m}}\right) d\mathfrak{m}\bigg] \\ & = \alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{G}\circ\mathscr{K}\left( \frac{1}{\mathfrak{r}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{G}\circ\mathscr{K}\left( \frac{1}{\mathfrak{s}}\right)\bigg]. \end{array}$

(3.6)

This gives us the first part of the inequality (3.4). Now, for the next part, we use the hypotheses of harmonically convex function i.e.

$\begin{array}{l} \mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}}\right)+\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right) \leq\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s}). \end{array}$

(3.7)

Multiplying both side of the above inequality (3.7) by $e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}$ and then integrating the resultant over $[0, 1]$ , we obtain

$\begin{align*} &\int_{0}^{\frac{1}{2}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}}\right)d\Phi +\int_{0}^{\frac{1}{2}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)d\Phi\nonumber\\ &\leq [\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})]\int_{0}^{\frac{1}{2}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}d\Phi \nonumber\\ & = \frac{2(1-e^{-\frac{\rho_{h}}{2}})}{\rho_{h}}\frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}. \end{align*}$

Consequently from the first inequality (3.6), we have

$\begin{array}{l} &\int_{0}^{\frac{1}{2}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{r}+(1-\Phi)\mathfrak{s}}\right)d\Phi +\int_{0}^{\frac{1}{2}} e^{-\frac{1-\alpha}{\alpha}\frac{\mathfrak{s}-\mathfrak{r}}{\mathfrak{r}\mathfrak{s}}\Phi}\mathscr{G}\left( \frac{\mathfrak{r}\mathfrak{s}}{\Phi\mathfrak{s}+(1-\Phi)\mathfrak{r}}\right)d\Phi\\ & = \alpha\frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}\bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{G}\circ\mathscr{K}\left( \frac{1}{\mathfrak{r}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{G}\circ\mathscr{K}\left( \frac{1}{\mathfrak{s}}\right)\bigg]\leq\frac{2(1-e^{-\frac{\rho_{h}}{2}})}{\rho_{h}}\frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}. \end{array}$

(3.8)

From the above developments (3.6) and (3.8), it follows

$\mathscr{G}\left( \frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}\right)\leq\frac{(1-\alpha)}{2\left( 1-e^{\frac{\rho_{h}}{2}}\right) }\bigg[\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}\mathscr{G}\circ\mathscr{K}\left( \frac{1}{\mathfrak{r}}\right)+\mathtt{I}^\alpha_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}\mathscr{G}\circ\mathscr{K}\left( \frac{1}{\mathfrak{s}}\right)\bigg]\leq\frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2}.$

This concludes the proof of the required result.

Remark 3.1. If one chooses $\alpha\to 1$ i.e., $\frac{\rho_{h}}{2}\to 0$ , then

$\lim\limits_{\alpha\to 1}\frac{1-\alpha}{2(1-e^{-\frac{\rho_{h}}{2}})} = \frac{\mathfrak{r}\mathfrak{s}}{\mathfrak{s}-\mathfrak{r}}$

and hence Theorem 2.3 retrieves the classical Hermite-Hadamard inequality (1.3) for harmonically convex function.

Example 3.3. Let $\mathscr{G}\left(\mathfrak{m}\right) = \mathfrak{m}^{2}$ , $\mathfrak{m}\in\left[1, 4\right]$ , $\mathtt{K}\left(\mathfrak{m}\right) = \frac{1}{\mathfrak{m}}$ and $0 < \alpha < 1$ , then

$\mathscr{G}\left( \frac{2\mathfrak{r}\mathfrak{s}}{\mathfrak{r}+\mathfrak{s}}\right) = \frac{64}{25},$

$\begin{align*} & \frac{(1-\alpha)}{2\left( 1-e^{\frac{\rho_{h}}{2}}\right) }\left[ \mathtt{I}_{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}+}^{\alpha }\mathscr{G}\circ\mathscr{K}\left( \frac{1}{\mathfrak{r}}\right) +\mathtt{I} _{\frac{\mathfrak{r}+\mathfrak{s}}{2\mathfrak{r}\mathfrak{s}}-}^{\alpha} \mathscr{G}\circ\mathscr{K}\left( \frac{1}{\mathfrak{s}}\right) \right] \\ & = \frac{(1-\alpha)}{2\left( 1-e^{\frac{3\left( 1-\alpha\right) }{8\alpha }}\right) }\left[ \mathtt{I}_{\frac{\mathit{5}}{8}+}^{\alpha}\mathscr{G} \circ\mathscr{K}\left( \frac{1}{\mathfrak{r}}\right) +\mathtt{I} _{\frac{\mathit{5}}{8}-}^{\alpha}\mathscr{G}\circ\mathscr{K}\left( \frac {1}{\mathfrak{s}}\right) \right] \\ & = \frac{(1-\alpha)}{2\left( 1-e^{\frac{3\left( 1-\alpha\right) }{8\alpha }}\right) }\left[ \frac{1}{\alpha}\int_{\frac{1}{4}}^{\frac{\mathit{5}}{8} }e^{-\frac{1-\alpha}{\alpha}\left( \mathfrak{m}-\frac{1}{4}\right) }\;\left( \frac{1}{\mathfrak{m}}\right) ^{2}d\mathfrak{m}+\frac{1}{\alpha}\int_{\frac{\mathit{5} }{8}}^{1}e^{-\frac{1-\alpha}{\alpha}\left( 1-\mathfrak{m}\right) }\;\left( \frac{1}{\mathfrak{m}}\right) ^{2}d\mathfrak{m}\right] \end{align*}$

and

$\frac{\mathscr{G}(\mathfrak{r})+\mathscr{G}(\mathfrak{s})}{2} = \frac{17}{2}.$

The graphical representation of Theorem 3.3 is shown in the graph below (see ) for $0 < \alpha < 1$ :

Figure 6. The graphical representation of Theorem 3.3 for

$0 < \alpha < 1$ .

DownLoad: Full-Size Img PowerPoint

4. Applications

Example 4.1. Let $\mathbb{C}^{n}$ be the set of $n \times n$ complex matrices, $\mathbb{M}_{n}$ denote the algebra of $n \times n$ complex matrices, and $\mathbb{M}_{n}^{+}$ denote the strictly positive matrices in $\mathbb{M}_{n}$ . That is, for any nonzero $u\in \mathbb{C}^{n}$ , $\mathrm{A} \in \mathbb{M}^{+}_{n}$ if $\big < \mathrm{A}u, u\big > > 0$ .

Sababheh ^[51], proved that $\mathscr{G}(\kappa) = \parallel\mathrm{A}^{\kappa}\mathrm{X}\mathrm{B}^{1-\kappa}+\mathrm{A}^{1-\kappa}\mathrm{X}\mathrm{B}^{\kappa}\parallel$ , $\mathrm{A}, \mathrm{B}\in \mathbb{M}^{+}_{n}, \mathrm{X}\in\mathbb{M}_{n}$ is convex for all $\kappa\in[0, 1]$ .

Then, by using Theorem 2.3, we have

$\begin{align*} &\parallel\mathrm{A}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}}\mathrm{X}\mathrm{B}^{1-\left(\frac{\mathfrak{r}+\mathfrak{s}}{2}\right)}+\mathrm{A}^{1-\left(\frac{\mathfrak{r}+\mathfrak{s}}{2}\right)}\mathrm{X}\mathrm{B}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}}\parallel\\ &\leq \frac{1-\alpha}{2(1-e^{-\frac{\rho_{c}}{2}})}\bigg[\;{\mathtt{I}}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}^{+}}^\alpha \parallel\mathrm{A}^{\mathfrak{s}}\mathrm{X}\mathrm{B}^{1-\mathfrak{s}}+\mathrm{A}^{1-\mathfrak{s}}\mathrm{X}\mathrm{B}^{\mathfrak{s}}\parallel + \;{\mathtt{I}}_{\frac{\mathfrak{r}+\mathfrak{s}}{2}^{-}} ^\alpha \parallel\mathrm{A}^{\mathfrak{r}}\mathrm{X}\mathrm{B}^{1-\mathfrak{r}}+\mathrm{A}^{1-\mathfrak{r}}\mathrm{X}\mathrm{B}^{\mathfrak{r}}\parallel\bigg]\\ &\leq \frac{\parallel\mathrm{A}^{\mathfrak{r}}\mathrm{X}\mathrm{B}^{1-\mathfrak{r}}+\mathrm{A}^{1-\mathfrak{r}}\mathrm{X}\mathrm{B}^{\mathfrak{r}}\parallel+\parallel\mathrm{A}^{\mathfrak{s}}\mathrm{X}\mathrm{B}^{1-\mathfrak{s}}+\mathrm{A}^{1-\mathfrak{s}}\mathrm{X}\mathrm{B}^{\mathfrak{s}}\parallel}{2}. \end{align*}$

Example 4.2. The $q$ -digamma(psi) function $\psi_\Phi$ given as (see ^[52]):

$\begin{align*} \psi_\Phi(\zeta)& = -\ln \left( 1-\Phi\right) +\ \ln \Phi \sum\limits_{k = 0}^{\infty}\frac{\Phi^{k+\zeta}}{1-\Phi^{k+\zeta}} \\& \ \ \ \ = -\ln \left( 1-\Phi\right) +\ \ln \Phi \sum\limits_{k = 1}^{\infty}\frac{\Phi^{k\zeta}}{1-\Phi^{k\zeta}}. \end{align*}$

For $\Phi > 1$ and $\zeta > 0,$ $\Phi$ -digamma function $\psi_\Phi$ can be given as:

$\begin{align*} \psi_\Phi(\zeta)& = -\ln \left( \Phi-1\right) +\ \ln \Phi\bigg[\zeta-\frac{1}{2} - \sum\limits_{k = 0}^{\infty}\frac{\Phi^{-(k+\zeta)}}{1-\Phi^{-(k+\zeta)}}\bigg] \\& \ \ \ \ = -\ln \left( \Phi-1\right) +\ \ln \Phi\bigg[\zeta-\frac{1}{2} - \sum\limits_{k = 1}^{\infty}\frac{\Phi^{-k\zeta}}{1-\Phi^{-k\zeta}}\bigg]. \end{align*}$

If we set $\mathscr{G}(\zeta) = \psi^{\prime}_{\Phi}(\zeta)$ in Theorem 2.3, then we have the following inequality.

$\begin{array}{l} \psi^{\prime}_{\Phi}\left(\frac{\mathfrak{r}+\mathfrak{s}}{2}\right) \leq\frac{1-\alpha}{2\alpha(1-e^{-\frac{\Phi_{c}}{2}})}\left[\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2}}^{\mathfrak{s}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{m})}\psi^{\prime}_{\Phi}\left( \mathfrak{m}\right)d\mathfrak{m} +\int_{\mathfrak{r}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{m}-\mathfrak{r})}\psi^{\prime}_{\Phi}\left( \mathfrak{m}\right)d\mathfrak{m} \right] \leq \frac{\psi^{\prime}_{\Phi}(\mathfrak{r})+\psi^{\prime}_{\Phi}(\mathfrak{s})}{2}. \end{array}$

Modified Bessel functions

Example 4.3. Let the function $\mathscr{J}_{\rho}:\mathcal{R}\to [1, \infty)$ be defined ^[52] as

$\mathscr{J}_{\rho}(\mathfrak{m}) = 2^{\rho}\Gamma(\rho+1)\mathfrak{m}^{-\rho}\mathrm{I}_{\rho}(\mathfrak{m}), \ \ \mathfrak{m}\in{\mathcal{R}}.$

Here, we consider the modified Bessel function of first kind given in

$\mathscr{J}_{\rho}(\mathfrak{m}) = \sum\limits_{n = 0}^{\infty}\frac{\left(\frac{\mathfrak{m}}{2}\right)^{\rho+2n}}{n!\Gamma(\rho+n+1)}.$

The first and second order derivative are given as

$\mathscr{J}^{\prime}_{\rho}(\mathfrak{m}) = \frac{\mathfrak{m}}{2(\rho+1)}\mathscr{J}_{\rho+1}(\mathfrak{m}).$

$\mathscr{J}^{\prime\prime}_{\rho}(\mathfrak{m}) = \frac{1}{4(\rho+1)}\left[\frac{u^{2}}{(\rho+1)}\mathscr{J}_{\rho+2}(\mathfrak{m})+2\mathscr{J}_{\rho+1}(\mathfrak{m})\right].$

If we use, $\mathscr{G}(\mathfrak{m}) = \mathscr{J}^{\prime}_{\rho}(\mathfrak{m})$ and the above functions in Theorem 2.3, we have

$\begin{array}{l} \frac{\mathfrak{r}+\mathfrak{s}}{2}\mathscr{J}_{\rho+1}\left(\frac{\mathfrak{r}+\mathfrak{s}}{2}\right)\\ \leq \frac{1-\alpha}{2\alpha(1-e^{-\frac{\rho_{c}}{2}})} \left[\int_{\frac{\mathfrak{r}+\mathfrak{s}}{2}}^{\mathfrak{s}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{s}-\mathfrak{m})}\mathfrak{m}\mathscr{J}_{\rho+1}\left(\mathfrak{m}\right)d\mathfrak{m} +\int_{\mathfrak{r}}^{\frac{\mathfrak{r}+\mathfrak{s}}{2}}e^{-\frac{1-\alpha}{\alpha}(\mathfrak{m}-\mathfrak{r})}\mathfrak{m}\mathscr{J}_{\rho+1}\left(\mathfrak{m}\right)d\mathfrak{m}\right]\\ \leq \frac{\mathfrak{r}\mathscr{J}_{\rho+1}\left(\mathfrak{r}\right)+\mathfrak{s}\mathscr{J}_{\rho+1}\left(\mathfrak{s}\right)}{2} . \end{array}$

5. Concluding remarks

The use of fractional calculus for finding various integral inequalities via convex functions has skyrocketed in recent years. This paper addresses a novel sort of Fejér type integral inequalities. In order to generalize some H-H-F (Hermite-Hadamard-Fejér) type inequalities, a new fractional integral operator with exponential kernel is employed. New midpoint type inequalities for both convex and harmonically convex functions are studied. Applications related to matrices, q-digamma and modifed Bessel functions are presented as well.

6. Future scopes

We will use our theories and methods to create new inequalities for future research by combining these new weighted generalized fractional integral operators with Chebyshev, Simpson, Jensen-Mercer Markov, Bullen, Newton, and Minkowski type inequalities. Quantum calculus, fuzzy interval-valued analysis, and interval-valued analysis can all be used to establish these kinds of inequalities. The idea of Digamma functions and other special functions will be integrated with this kind of inequality as the major focus. We also aim to find other novel inequalities using finite products of functions. In future, we will employ the concept of $\mathit{cr}$ -order defined by Bhunia and Samanta ^[53] to present different inequalities for $\mathit{cr}$ -convexity and $\mathit{cr}$ -harmonically convexity ^[54].

Acknowledgments

This research received funding support from the NSRF via the Program Management Unit for Human Resources & Institutional Development, Research and Innovation, (grant number B05F650018).

Conflict of interest

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

References

[1]	World Health Organization, Epidemiological bulletin WHO health emergencies programme WHO regional office for south-east Asia, 8 Eds., 2023. Available from: https://iris.who.int/handle/10665/372269
[2]	S. Naseer, S. Khalid, S. Parveen, K. Abbass, H. Song, M. V. Achim, COVID-19 outbreak: Impact on global economy, Front. Public Health, 10 (2023), 1009393. https://doi.org/10.3389/fpubh.2022.1009393 doi: 10.3389/fpubh.2022.1009393
[3]	Y. Shang, H. Li, R. Zhang, Effects of Pandemic outbreak on economies: Evidence from business history context, Front. Public Health, 9 (2021), 632043. https://doi.org/10.3389/fpubh.2021.632043 doi: 10.3389/fpubh.2021.632043
[4]	S. Murthy, C. D. Gomersall, R. A. Fowler, Care for critically ill patients with Covid-19, JAMA, 323 (2020), 1499–1500. https://doi.org/10.1001/jama.2020.3633 doi: 10.1001/jama.2020.3633
[5]	H. S. Rodrigues, M. T. T. Monteiro, D. F. M. Torres, Sensitivity analysis in a dengue epidemiological model, Conf. Pap. Math., 2013 (2013), 721406. https://doi.org/10.1155/2013/721406 doi: 10.1155/2013/721406
[6]	World Health Organization, Dengue guidelines, for diagnosis, treatment, prevention and control, 2009. Available from: https://www.who.int/publications-detail-redirect/9789241547871
[7]	M. Verduyn, N. Allou, V. Gazaille, M. Andre, T. Desroche, M. C. Jaffar et al., Co-infection of dengue and Covid-19: A case report, PLoS Negl. Trop. Dis., 14 (2020), e0008476. https://doi.org/10.1371/journal.pntd.0008476 doi: 10.1371/journal.pntd.0008476
[8]	S. Masyeni, M. S. Santoso, P. D. Widyaningsih, D. G. W. Asmara, F. Nainu, H. Harapan, et al., Serological cross-reaction and coinfection of dengue and COVID-19 in Asia: Experience from Indonesia, Int. J. Infect. Dis., 102 (2021), 152–154. https://doi.org/10.1016j.ijid.2020.10.043
[9]	L. Lansbury, B. Lim, V. Baskaran, W. S. Lim, Co-infections in people with Covid-19: A systematic review and meta-analysis, J. infect., 81 (2020), 266–275. https://doi.org/10.1016/j.jinf.2020.05.046 doi: 10.1016/j.jinf.2020.05.046
[10]	A. A. Malibari, F. Al-Husayni, A. Jabri, A. Al-Amri, M. Alharbi, A patient with Dengue Fever and COVID-19: coinfection or not? Cureus, 12 (2020), e11955. https://doi.org/10.7759/cureus.11955 doi: 10.7759/cureus.11955
[11]	A. Saddique, M. S. Rana, M. M. Alam, A. Ikram, M. Usman, M. Salman, et al., Emergence of co-infection of COVID-19 and dengue: A serious public health threat, J. Infect., 81 (2020), e16–e18. https://doi.org/10.1016/j.jinf.2020.08.009 doi: 10.1016/j.jinf.2020.08.009
[12]	M. I. H. Hussein, A. A. D. Albashir, O. A. M. A. Elawad, A. Homeida, Malaria and COVID-19: Unmasking their ties, Malaria J., 19 (2020), 457. https://doi.org/10.1186/s12936-020-03541-w doi: 10.1186/s12936-020-03541-w
[13]	K. Sverdrup, S. J. Kimmerle, P. Berg, Computational investigation of the stability and dissolution of nanobubbles, Appl. Math. Model., 49 (2017), 199–219. https://doi.org/10.1016/j.apm.2017.05.006 doi: 10.1016/j.apm.2017.05.006
[14]	M. Farman, A. Akgül, M. T. Tekin, M. M. Akram, A. Ahmad, E. E. Mahmoud, et al., Fractal fractional-order derivative for HIV/AIDS model with Mittag-Leffler kernel, Alex. Eng. J., 61 (2022), 10965–10980. https://doi.org/10.1016/j.aej.2022.04.030 doi: 10.1016/j.aej.2022.04.030
[15]	S. Ahmad, A. Ullah, A. Akgül, D. Baleanu, Analysis of the fractional tumour-immune-vitamins model with Mittag-Leffler kernel, Results Phys., 19 (2020), 103559. https://doi.org/10.1016/j.rinp.2020.103559 doi: 10.1016/j.rinp.2020.103559
[16]	Y. Sabbar, A. Din, D. Kiouach, Influence of fractal-fractional differentiation and independent quadratic Levy jumps on the dynamics of a general epidemic model with vaccination strategy, Chaos Solitons Fractals, 171 (2023), 113434. https://doi.org/10.1016/j.chaos.2023.113434 doi: 10.1016/j.chaos.2023.113434
[17]	A. I. K. Butt, Atangana-Baleanu fractional dynamics of predictive whooping Cough model with optimal control analysis, Symmetry, 15 (2023), 1773. https://doi.org/10.3390/sym15091773 doi: 10.3390/sym15091773
[18]	B. Fatima, M. Yavuz, M Ur Rahman, A. Althobaiti, S. Althobaiti, Predictive modeling and control strategies for the transmission of middle east respiratory syndrome coronavirus, Math. Comput. Appl., 28 (2023), 98. https://doi.org/10.3390/mca28050098 doi: 10.3390/mca28050098
[19]	A. Hanif, A. I. K. Butt, Atangana-Baleanu fractional dynamics of dengue fever with optimal control strategies, AIMS Mathematics, 8 (2023), 15499–15535. https://doi.org/10.3934/math.2023791 doi: 10.3934/math.2023791
[20]	A. Omame, D. Okuonghae, A co-infection model for oncogenic Human papillomavirus and tuberculosis with optimal control and cost-effectiveness analysis, Optim. Control. Appl. Methods, 42 (2021), 1081–1101. https://doi.org/10.1002/oca.2717 doi: 10.1002/oca.2717
[21]	Y. Guo, T. Li, Fractional-order modeling and optimal control of a new online game addiction model based on real data, Commun. Nonlinear Sci. Numer. Simul., 121 (2023), 107221. https://doi.org/10.1016/j.cnsns.2023.107221 doi: 10.1016/j.cnsns.2023.107221
[22]	D. Baleanu, F. A. Ghassabzade, J. J. Nieto, A. Jajarmi, On a new and generalized fractional model for a real cholera outbreak, Alex. Eng. J., 61 (2022), 9175–9186. https://doi.org/10.1016/j.aej.2022.02.054 doi: 10.1016/j.aej.2022.02.054
[23]	B. Razia, T. Osman, H. Khan, G. Haseena, A. Khan, A fractional order Zika virus model with Mittag-Leffler kernel, Chaos Solitons Fractals, 146 (2021), 110898. https://doi.org/10.1016/j.chaos.2021.110898 doi: 10.1016/j.chaos.2021.110898
[24]	A. Atangana, I. Koca, Chaos in a simple nonlinear system with Atangana-Baleanu derivative with fractional order, Chaos Solitons Fractals, 89 (2016), 447–454. https://doi.org/10.1016/j.chaos.2016.02.012 doi: 10.1016/j.chaos.2016.02.012
[25]	A. I. K. Butt, W. Ahmad, M. Rafiq, N. Ahmad, M. Imran, Optimally analyzed fractional Coronavirus model with Atangana–Baleanu derivative, Results Phys., 53 (2023), 106929. https://doi.org/10.1016/j.rinp.2023.106929 doi: 10.1016/j.rinp.2023.106929
[26]	M. Caputo, M. Fabrizio, A new definition of fractional derivative without singular kernel, Progr. Fract. Differ. Appl., 1 (2015), 73–85. http://doi.org/10.12785/pfda/010201 doi: 10.12785/pfda/010201
[27]	D. Baleanu, A. Jajarmi, H. Mohammadi, S. Rezapour, A new study on the mathematical modelling of human liver with Caputo-Fabrizio fractional derivative, Chaos Solitons Fractals, 134 (2020), 109705. https://doi.org/10.1016/j.chaos.2020.109705 doi: 10.1016/j.chaos.2020.109705
[28]	T. Khan, R. Ullah, G. Zaman, J. Alzabut, A mathematical model for the dynamics of SARS-CoV-2 virus using the Caputo-Fabrizio operator, Math. Biosci. Eng., 18 (2021), 6095–6116. https://doi.org/10.3934/mbe.2021305 doi: 10.3934/mbe.2021305
[29]	A. I. K. Butt, M. Imran, S. Batool, M. Al Nuwairan, Theoretical analysis of a COVID-19 CF-fractional model to optimally control the spread of pandemic, Symmetry, 15 (2023), 380. https://doi.org/10.3390/sym15020380 doi: 10.3390/sym15020380
[30]	A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Elsevier Science, 2006.
[31]	A. Atangana, D. Balneau, New fractional derivative with nonlocal and nonsingular kernel: Theory and application to heat transfer model, Thermal Sci., 20 (2016), 763–769. http://doi.org/10.2298/TSCI160111018A doi: 10.2298/TSCI160111018A
[32]	M. A. Al Shorbagy, M. Ur Rahman, Y. Karaca, A computational analysis fractional complex-order values by ABC operator and Mittag-Leffler Kernel modeling, Fractals, 31 (2023), 2340164. https://doi.org/10.1142/S0218348X23401643 doi: 10.1142/S0218348X23401643
[33]	P. Liu, Mati ur Rahman, A. Din, Fractal fractional based transmission dynamics of COVID-19 epidemic model, Comput. Methods Biomec. Biomed. Eng., 25 (2022), 1852–1869. https://doi.org/10.1080/10255842.2022.2040489 doi: 10.1080/10255842.2022.2040489
[34]	A. Hanif, A. I. K. Butt, W. Ahmad, Numerical approach to solve Caputo‐Fabrizio‐fractional model of corona pandemic with optimal control design and analysis, Math. Methods Appl. Sci., 46 (2023), 9751–9782. https://doi.org/10.1002/mma.9085 doi: 10.1002/mma.9085
[35]	A. Hanif, A. I. K. Butt, S. Ahmad, R. U. Din, M. Inc, A new fuzzy fractional order model of transmission of Covid-19 with quarantine class, Eur. Phys. J. Plus, 136 (2021), 1179. https://doi.org/10.1140/epjp/s13360-021-02178-1 doi: 10.1140/epjp/s13360-021-02178-1
[36]	A. I. K. Butt, M. Rafiq, W. Ahmad, N. Ahmad, Implementation of computationally efficient numerical approach to analyze a Covid-19 pandemic model, Alex. Eng. J., 69 (2023), 341–362. https://doi.org/10.1016/j.aej.2023.01.052 doi: 10.1016/j.aej.2023.01.052
[37]	M. L. Diagne, H. Rwezaura, S. Y. Tchoumi, J. M. Tchuenche, A mathematical model of Covid-19 with vaccination and treatment, Comput. Math. Methods Med., 2021 (2021), 1250129. https://doi.org/10.1155/2021/1250129 doi: 10.1155/2021/1250129
[38]	A. I. K. Butt, W. Ahmad, M. Rafiq, N. Ahmad, M. Imran, Computationally efficient optimal control analysis for the mathematical model of Coronavirus pandemic, Expert Syst. Appl., 234 (2023), 121094. https://doi.org/10.1016/j.eswa.2023.121094 doi: 10.1016/j.eswa.2023.121094
[39]	K. Diethelm, A fractional calculus based model for the simulation of an outbreak of dengue fever, Nonlinear Dynam., 71 (2013), 613–619. https://doi.org/10.1007/s11071-012-0475-2 doi: 10.1007/s11071-012-0475-2
[40]	Y. Chu, S. Sultana, S. Rashid, M. S. Alharthi, Dynamical analysis of the stochastic COVID-19 model using piecewise differential equation technique, Comput. Model. Eng. Sci., 137 (2023), 2427–2464. https://doi.org/10.32604/cmes.2023.028771 doi: 10.32604/cmes.2023.028771
[41]	Y. Chu, S. Rashid, A. O. Akdemir, A. Khalid, D. Baleanu, B. R. Al-Sinan, et al., Predictive dynamical modeling and stability of the equilibria in a discrete fractional difference COVID-19 epidemic model, Results Phys., 49 (2023), 106467. https://doi.org/10.1016/j.rinp.2023.106467 doi: 10.1016/j.rinp.2023.106467
[42]	T. Li, Y. Guo, Modeling and optimal control of mutated COVID-19 (Delta strain) with imperfect vaccination, Chaos Solitons Fractals, 156 (2022), 111825. https://doi.org/10.1016/j.chaos.2022.111825 doi: 10.1016/j.chaos.2022.111825
[43]	A. I. K. Butt, M. Imran, B. A. McKinney, S. Batool, H. Aftab, Mathematical and stability analysis of dengue-malaria co-infection with disease control strategies, Mathematics, 11 (2023), 4600. https://doi.org/10.3390/math11224600 doi: 10.3390/math11224600
[44]	E. Mtisi, H. Rwezaura, J. M. Tchuenche, A mathematical analysis of malaria and tuberculosis co-dynamics, Discrete Contin. Dyn. Syst. Ser. B, 12 (2009), 827–864. https://doi.org/10.3934/dcdsb.2009.12.827 doi: 10.3934/dcdsb.2009.12.827
[45]	A. Omame, M. Abbas, C. P. Onyenegecha, A fractional-order model for Covid-19 and tuberculosis co-infection using Atangana-Baleanu derivative, Chaos Solitons Fractals, 153 (2001), 111486. https://doi.org/10.1016/j.chaos.2021.111486 doi: 10.1016/j.chaos.2021.111486
[46]	A. Omame, H. Rwezaura, M. L. Diagne, S. C. Inyama, J. M. Tchuenche, COVID-19 and dengue co-infection in Brazil: Optimal control and cost-effectiveness analysis, Eur. Phys. J. Plus, 136 (2021), 1090. https://doi.org/10.1140/epjp/s13360-021-02030-6 doi: 10.1140/epjp/s13360-021-02030-6
[47]	A. Atangana, D. Baleanu, New fractional derivatives with nonlocal and non-singular kernel: Theory and applications to heat transfer model, Therm Sci., 20 (2016), 763–769. http://dx.doi.org/10.2298/TSCI160111018A doi: 10.2298/TSCI160111018A
[48]	J. K. K. Asamoah, Fractal fractional model and numerical scheme based on Newton polynomial for Q fever disease under Atangana-Baleanu derivative, Results Phys., 34 (2022), 105189. https://doi.org/10.1016/j.rinp.2022.105189 doi: 10.1016/j.rinp.2022.105189
[49]	G. Birkhoff, G. C. Rota, Ordinary differential equations, 4 Eds., John Wiley & Sons, 1991.
[50]	C. M. A. Pinto, J. A. T. Machado, Fractional model for malaria transmission under control strategies, Comput. Math. Appl., 66 (2013), 908–916. https://doi.org/10.1016/j.camwa.2012.11.017 doi: 10.1016/j.camwa.2012.11.017
[51]	E. Mtisi, H. Rwezaura, J.M. Tchuenche, A mathematical analysis of malaria and tuberculosis co-dynamics, Discrete Contin. Dyn. Syst. Ser. B, 12 (2009), 827–864. https://doi.org/10.3934/dcdsb.2009.12.827 doi: 10.3934/dcdsb.2009.12.827
[52]	M. Toufik, A. Atangana, New numerical approximation of fractional derivative with non-local and non-singular kernel: Application to chaotic models, Eur. Phys. J. Plus, 132 (2017), 444. https://doi.org/10.1140/epjp/i2017-11717-0 doi: 10.1140/epjp/i2017-11717-0
[53]	R. Kamocki, Pontryagin maximum principle for fractional ordinary optimal control problems, Math. Methods Appl. Sci., 37 (2014), 1668–1686. https://doi.org/10.1002/mma.2928 doi: 10.1002/mma.2928
[54]	S. Lenhart, J. T. Workman, Optimal control applied to biological models, CRC Press, 2007.
[55]	H. M. Ali, F. L. Pereira, S. M. A. Gama, A new approach to the Pontryagin maximum principle for nonlinear fractional optimal control problems, Math. Methods Appl. Sci., 39 (2016), 3640–3649. https://doi.org/10.1002/mma.3811 doi: 10.1002/mma.3811

This article has been cited by:

1.	Soubhagya Kumar Sahoo, Artion Kashuri, Munirah Aljuaid, Soumyarani Mishra, Manuel De La Sen, On Ostrowski–Mercer’s Type Fractional Inequalities for Convex Functions and Applications, 2023, 7, 2504-3110, 215, 10.3390/fractalfract7030215
2.	Artion Kashuri, Soubhagya Kumar Sahoo, Pshtiwan Othman Mohammed, Eman Al-Sarairah, Y. S. Hamed, Some New Hermite-Hadamard Type Inequalities Pertaining to Fractional Integrals with an Exponential Kernel for Subadditive Functions, 2023, 15, 2073-8994, 748, 10.3390/sym15030748
3.	Artion Kashuri, Soubhagya Kumar Sahoo, Munirah Aljuaid, Muhammad Tariq, Manuel De La Sen, Some New Hermite–Hadamard Type Inequalities Pertaining to Generalized Multiplicative Fractional Integrals, 2023, 15, 2073-8994, 868, 10.3390/sym15040868
4.	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, 2023, 8, 2473-6988, 13785, 10.3934/math.2023700
5.	Muhammad Tariq, Sotiris K. Ntouyas, Asif Ali Shaikh, A Comprehensive Review on the Fejér-Type Inequality Pertaining to Fractional Integral Operators, 2023, 12, 2075-1680, 719, 10.3390/axioms12070719
6.	Çetin Yıldız, Gauhar Rahman, Luminiţa-Ioana Cotîrlă, On Further Inequalities for Convex Functions via Generalized Weighted-Type Fractional Operators, 2023, 7, 2504-3110, 513, 10.3390/fractalfract7070513
7.	Yu Peng, Serap Özcan, Tingsong Du, Symmetrical Hermite–Hadamard type inequalities stemming from multiplicative fractional integrals, 2024, 183, 09600779, 114960, 10.1016/j.chaos.2024.114960
8.	XIAOHUA ZHANG, YUNXIU ZHOU, TINGSONG DU, PROPERTIES AND 2α̃-FRACTAL WEIGHTED PARAMETRIC INEQUALITIES FOR THE FRACTAL (m,h)-PREINVEX MAPPINGS, 2023, 31, 0218-348X, 10.1142/S0218348X23501347
9.	Muhammad Bilal Khan, Ali Althobaiti, Cheng-Chi Lee, Mohamed S. Soliman, Chun-Ta Li, Some New Properties of Convex Fuzzy-Number-Valued Mappings on Coordinates Using Up and Down Fuzzy Relations and Related Inequalities, 2023, 11, 2227-7390, 2851, 10.3390/math11132851

Reader Comments

Your name:*

Email:*
© 2024 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)