Some New Harmonically Convex Function Type Generalized Fractional Integral Inequalities

In this article, we established a new version of generalized fractional Hadamard and Fejér–Hadamard type integral inequalities. A fractional integral operator (FIO) with a non-singular function (multi-index Bessel function) as its kernel and monotone increasing functions is utilized to obtain the new version of such fractional inequalities. Our derived results are a generalized form of several proven inequalities already existing in the literature. The proven inequalities are useful for studying the stability and control of corresponding fractional dynamic equations.


Introduction
In the present era, fractional integral operators involving inequalities are widely derived by [1][2][3][4]. These fractional integral operators of any arbitrary real or complex order involve a different type of kernel. The field of fractional calculus has gained considerable importance among mathematicians and scientists due to its wide applications in sciences, engineering, and many other fields [5][6][7][8][9]. Hadamard and Fejér-Hadamard type inequalities have been discussed for many functions using different fractional operators with different kernels. Abbas and Farid [10] proposed the Hadamard and Fejér-Hadamard type integral inequalities for harmonically convex functions using the two-sided generalized fractional integral operator. Farid et al. [11,12] discussed these results in generalized form with an extended generalized Mittag-Leffler function. Hadamard and Fejér-Hadamard type inequalities are widely studied by the researchers [12][13][14][15][16][17][18][19]. The objective of this paper is to derive Hadamard, Fejér-Hadamard, and some other related type inequalities for the harmonically convex function via a generalized fractional operator with a nonsingular function as its kernel, which involves a multi-index Bessel function. For a recent related weighted fractional generalized approach, we refer to [20].
Hermite-Hadamard inequality and Fejér-Hadamard inequality are given by Theorem 1 ([21-23]). The inequality derived on the interval I = [u, v] ⊆ R called Hermite Hadamard inequality is given by where u, v ∈ I, with u = v and ρ : I → R is a convex function.
Theorem 2 ([21,24,25]). The Fejér-Hadamard inequality is defined for a convex function ρ : I → R and for a function µ : I → R, which is non-negative, integrable, and symmetric about u+v 2 , defined by where u, v ∈ I, with u = v.
Definition 2 ( [21,22]). Let I be an interval of nonzero real numbers. Then a function ρ : I → R is said to be harmonically convex if holds for all u, v ∈ I and t ∈ [0, 1].

Main Results
In this section, we present Hadamard, and Fejér-Hadamard type inequalities for harmonically convex functions by employing the new generalized fractional integral operators with a multi-index Bessel function as its kernel. We also establish a new version of inequalities by expressing the generalized fractional integral operator as the sum of two fractional integrals. where Proof. If θ is harmonically convex on [a, b], for every x, y ∈ [a, b], the following inequality holds and By multiplying by ( Solving the integrals involved in right side of inequality (13) by making substitution in first integral and 1 To obtain the second part of the inequality, the harmonic convexity of θ, we have the following relation Solving the integrals involved in the left side of inequality (16) by making substitution in first integral and in the second integral, we obtain .

Corollary 1.
If ψ(x) = 1 x in Theorem 3 then the following inequality holds Now, we derive the following Lemma before giving the next result.
, and ψ is differentiable and strictly increasing. If θ is a harmonically convex function on [a,b] and satisfies θ 1 , we have where Putting u = 1 ψ(a) in Equation (20), we have By the addition of T (21) on both sides, we have the required result.
] is a positive function, ψ is a differentiable and strictly increasing function and η is nonnegative and integrable and satisfies η 1 , then the following inequality holds where Proof. By using the harmonic convexity of θ, we have By multiplying by( in Equation (23) and then integrating over the closed interval [0, 1], we have By making a substitution of in the first integral and 1 in second integrals occurring at right side and in the integral occurring at left side of inequality (25) and using η 1 Now, we take By ) in Equation (27) and then integrating over [0, 1], we get dt.
Solving the integrals involved in left side of inequality (28) by making substitution in the first integral and 1 in the second integral and in the integral on the right side of the inequality and using η 1 Combining (26) and (29), we have the required result.
] is a positive and harmonically convex function and ψ is differentiable and strictly increasing, then the following inequality holds for the operators defined in Definition 6 ), where ).
on both sides and then integrating over Solving the integrals involved in the right side of inequality (32) by making a sub- in the first integral and the second integral as well as in the integral occurring at the left side of inequality (32), we have To obtain the second part of inequality, the harmonic convexity of θ gives the following relation Combining (33) and (36), we have the result.
Remark 2. 1. If ψ(x) = x, m = 1, ς = 0 and τ j is replaced by τ j − 1, it reduces to the result produced by Mehmet et al. [31] θ 2ab a + b ≤ Γ(τ j + 1) , and ψ is differentiable and strictly increasing. If θ is a harmonically convex function on [a,b] and satisfies θ 1 , we have Proof. Consider − ψ(x) and using θ 1 By the addition of T in Equation (39) on both sides, we have the required result.
] is a positive function, ψ is a differentiable, strictly increasing function and η is nonnegative and integrable and satisfies η 1 , then the following inequality holds for the operators defined in Definition 6. where Proof. By the harmonic convexity of θ, we have 2θ 2ψ(a)ψ(b) in the Equation (41) and then integrating over the closed interval [ 1 2 , 1], we have By substituting in the first integral and 1 in the second integrals occurring at the right side and in the integral occurring at left side of inequality (42), we have Now, we take ) in Equation (44) ≤ (θ(ψ(a)) + θ(ψ(b))) Solving the integrals involved in left side of inequality (45) by making substitution in the first integral and 1 in the second integral and in the integral on the right side of the inequality and using the above lemma and the condition 1 , we have .

Conclusion Remarks
In this article, we established Hadamard and Fejér-Hadamard type inequalities via a new generation of the generalized fractional integral operators (8) and (9) with a nonsingular function (multi-index Bessel function) as its kernel for harmonically convex functions. It is concluded that many classical inequalities cited in the literature can be easily derived by employing certain conditions on generalized fractional integral operators (8) and (9). We believe that our formulated inequalities will be useful to investigate the stability of certain fractional controlled systems.