In this paper, we aim to establish new inequalities of Hermite–Hadamard (H.H) type for harmonically convex functions using proportional Caputo-Hybrid (P.C.H) fractional operators. Parameterized by α, these operators offer a unique flexibility: setting α=1 recovers the classical inequalities for harmonically convex functions, while setting α=0 yields inequalities for differentiable harmonically convex functions. This framework allows us to unify classical and fractional cases within a single operator. To validate the theoretical results, we provide several illustrative examples supported by graphical representations, marking the first use of such visualizations for inequalities derived via P.C.H operators. Additionally, we demonstrate practical applications of the results by deriving new fractional-order recurrence relations for the modified Bessel function of type-1, which are useful in mathematical modeling, engineering, and physics. The findings contribute to the growing body of research in fractional inequalities and harmonic convexity, paving the way for further exploration of generalized convexities and higher-order fractional operators.
Read full abstract