Adapted block hybrid method for the numerical solution of Duffing equations and related problems

1.   Introduction

Convexities and preinvexities significantly contribute to advancing the optimization theory in both pure and applied mathematics. The increasing demand for fractional operators stems from the extensions and generalizations of fractional inequalities. Modifications applied to fractional operators enhance their kernels, and multi-dimensional series of special functions play a pivotal role in deriving new kernel versions, thereby improving the precision and applicability of fractional operators. Over the past two decades, extensive research has been conducted on extending generalized fractional operators, utilizing multi-index special functions as kernels. Researchers have diligently explored the use of these generalized operators to illustrate inequalities, garnering considerable attention in the field of fractional order analysis ^[1,2,3,4,5]. The interconnection between control optimization and the theory of inequalities is closely linked to the presence of (non)convex functions within their structures. This relationship has been thoroughly examined by researchers ^[6,7,8]. The classical theory of convexities facilitates numerous generalizations and extensions. One such essential generalization is the concept of an invex function, introduced by Hanson in ^[9]. This concept serves as a foundational extension of convex functions.

Furthermore, the fundamental characteristics of preinvex mappings and their applications in the realms of optimization and variational inequalities have been extensively explored. Noor ^[10,11], along with the works of Weir and Mond ^[12], has addressed these preinvex mappings, solving a multitude of problems related to equilibria. The application of preinvex mappings has proven to be instrumental in addressing challenges within optimization and variational inequalities, contributing to the resolution of complex equilibrium-related issues.

Diverse strategies are employed to tackle issues related to community specifications, decision-making, and optimizations, with their significance growing daily. The conventional crisp sets, which provide membership values of $0$ or $1$, pose a challenge when an object exhibits a property to varying degrees. To address this limitation, Zadeh introduced the concept of the fuzzy set theory, designed to handle ambiguous or uncertain information by allowing membership values in the range $[0, 1]$. The degree of membership within this interval reflects the level of assurance or conviction associated with the corresponding output value.

The most commonly concept of the fuzzy number theory is based on the characteristics function, but in fuzzy domain, the characteristics function is known as membership function that depends on the unit interval $[0, 1]$. This means that each subset of $\mathbb{R}$, and any real number, are exponential cases of the fuzzy number theory. In decision-making processes entailing unclear or imprecise information, fuzzy interval-valued functions prove invaluable. They enable the representation and analysis of hazy preferences or ambiguous results. Pioneering works by Moore in 1966, Kulish and Miranker ^[13] delved into the idea of interval analysis, uncovering its applications in arithmetic. Furthermore, a multitude of studies in ^[14] have explored fuzzy-set theory and system development across various fields, making noteworthy contributions to solving a wide array of problems in applied and pure mathematics. These contributions extend to disciplines such as computer science, operational research, management sciences, control optimization theory, decision sciences, and artificial intelligence.

So far, numerous researchers have explored the diverse properties and applications of variational inequalities within a set of functions exhibiting generalized convexity ^[15,16]. Additionally, in the realm of studying fractional inequalities using fuzzy interval-valued functions (FIV-functions), which involve fractional integrals with the integrand as an FIV-function, several noteworthy works have been addressed by researchers. For instance, Flores et al. ^[17] extended the Gronwall inequality to encompass interval-valued functions (Ⅳ-functions). Furthermore, Chalco et al. ^[18] applied generalized Hukuhara derivatives to establish Ostrowski-type inequalities for FIV-functions and delved into their numerical applications. Zhao et al. ^[19,20] directed their attention to Ⅳ-functions for the modification of well-known inequalities such as Jensen's, Chebyshev, and Hermite-Hadamard type inequalities. Following this line of research, Zhang et al. ^[21] derived Jensen's inequalities through the utilization of fuzzy set-valued functions.

The primary motivation of this study stems from the essential and widely applicable nature of two classes of functions characterized by convexity and preinvexity properties within the realm of optimization. In the context of fuzzy optimization, we aim to articulate optimality conditions under such functions and explore the derivation of various variational inequalities of fuzzy type. FIV-functions emerge as a pivotal fuzzy notion, offering a means to establish innovative applied fractional fuzzy inequalities through their combination with fractional operators. The incorporation of FIV-functions into these inequalities provides a unique avenue for investigating uncertainties in prediction scenarios, departing from the semantic norms of standard non-fuzzy inequalities. We delve into the utilization of FIV-functions to formulate fuzzy integral inequalities, where uncertainties manifest in the integrand rather than the measures of integration. To achieve this, we draw on a specific group of $h$-convex and $h$-preinvex FIV-functions of the Godunova-Levin type. This strategic utilization facilitates the derivation of notable fuzzy-based integral inequalities. In essence, the research endeavors to bridge the conceptual gap between convexity, preinvexity, and fuzzy optimization, offering a novel perspective on uncertainties in prediction situations through the lens of FIV-functions and fractional fuzzy inequalities.

In the following section, certain essential concepts are revisited to serve as a reminder. Section 3 explores Prabhakar fuzzy-based Hermite-Hadamard (H-H) type inequalities utilizing $h$-Godunova-Levin FIV-functions endowed with convexity. The convexity property is formulated explicitly within this section. In Section 4, the focus shifts to Prabhakar fuzzy-based trapezoidal-type inequalities, this time incorporating $h$-Godunova-Levin FIV-functions with the preinvexity property. The study culminates in Section 5, presenting concluding remarks and insights drawn from the findings.

2.   Preliminaries

In the present section, we focus on some basic definitions and properties, which are helpful to understand the contents of the paper. Throughout the paper, $K \subseteq \mathbb{R}$ is a convex set.

Definition 2.1. ^[22] The function $\varkappa:K\rightarrow \mathbb{R}$ is said to be convex if

where $t\in[0, 1]$ and $u, v \in K$ are arbitrary.

By considering the convex function defined above, the Hermite-Hadamard (H-H) type inequality

was discussed many remarkable results in ^[23], where $u, v \in K \subseteq \mathbb{R}$ and $u < v$.

Definition 2.2. ^[24] The preinvex function $y: K \rightarrow \mathbb{R}$, where $u, v \in K$ (Here, $K \subset \mathbb{R}$ is an invex set w.r.t. the bi-function $\xi$ on $K \times K$) and $\lambda \in [0, 1]$, is defined as

Definition 2.3. ^[25] Consider the functions $h:(0, 1)\rightarrow \mathbb{R}$ and $y:K \rightarrow \mathbb{R}^{\geq 0}$. Then, the $h$-Godunova-Levin function $y$ satisfies

for all $u, v\in K$ and $t\in(0, 1)$.

Definition 2.4. ^[25] Consider the functions $h:(0, 1)\rightarrow \mathbb{R}$ and $y:K \rightarrow \mathbb{R}$. Then the $h$-Godunova-Levin preinvex function $y$ w.r.t. $\xi$ satisfies

so that $u, v \in K, t\in (0, 1)$ are arbitrary.

Definition 2.5. ^[26] Let $y\in L^{1}[u, v]$. The Riemann-Liouville left and right integrals $\mathbb{I}_{u^{+}}^{b} y$ and $\mathbb{I}_{v^{-}}^{b} y$ of order $b > 0$ are defined as follows

and

By considering $F \subseteq \mathbb{R}$ as the fuzzy interval, let $\hbar_1, \hbar_2 \in F$. The partial order relation $\preccurlyeq$ is given by

where $[\hbar]^\ell$ is called the $\ell$-level set of $\hbar$ defined by

Proposition 2.6. ^[27] For $\hbar_1, \hbar_2 \in F$, $c\in \mathbb{R}$ and $\ell \in [0, 1]$, we have

Let $\hbar_1 = \hbar_2 \widetilde{+} H$, where $H \in F$. In this case, the Hukuhara difference ($H$-difference) of $\hbar_1$ and $\hbar_2$ is represented by the symbol $\hbar_1 \widetilde{-} \hbar_2$ and we have

The left fuzzy-based Riemann-Liouville integral of the function $y$ is defined as

where

and

We denote the family of all bounded closed fuzzy intervals of $\mathbb{R}$ by $K_{bc}$.

Definition 2.7. ^[28] For all $\ell \in [0, 1]$, the FIV-function $y:K\subseteq \mathbb{R} \rightarrow K_{bc}$ is formulated by $y_{\ell}(x) = \big[y_{*}(x, \ell), y^{*}(x, \ell)\big], \, \forall\, x\in K$ in which the left and right functions $y^{*}(x, \ell), y_{*}(x, \ell): K \rightarrow \mathbb{R}$ are named as the upper and lower functions of $y$.

Note that the FIV-function $y:K\subseteq \mathbb{R} \rightarrow K_{bc}$ is continuous at $x\in K$ if for all $\ell \in [0, 1]$, both $y_{*}(x, \ell)$ and $y^{*}(x, \ell)$ are continuous at $x \in K$.

Definition 2.8. ^[29] Assume that $K_{bc}$ is the aforesaid class of intervals and $H\in K_{bc}$. Then

If $\Lambda_{*} = \Lambda^{*}$, then we say that $\Lambda$ is degenerate.

In this research, we consider all the fuzzy intervals as the non-degenerate intervals. If $\Lambda_{*}\geq0$, then $[\Lambda_{*}, \Lambda^{*}]$ is named as a positive interval, and all of them can be specified by the symbol $K_{bc}^{+}$ so that

Let $\sigma\in\mathbb{R}$. Then, we define the scalar multiplication as ^[29]

Also, the Minkowski algebraic difference and addition actions $\Lambda_1-\Lambda_2$ and $\Lambda_1 + \Lambda_2$, where $\Lambda_1, \Lambda_2\in K_{bc}$, are defined as

Moreover, the inclusion "$\subseteq$" means that $\Lambda_1\subseteq \Lambda_2$ if and only if $[\Lambda_{1*}, \Lambda_1^{*}]\subseteq[\Lambda_{2*}, \Lambda_2^{*}]$ iff $\Lambda_{1*}\preccurlyeq \Lambda_{2*}, \Lambda_1^{*}\preccurlyeq \Lambda_2^{*}$ ^[29].

The fuzzy order relation $\preccurlyeq_{1}$ is defined on $K_{bc}$ as

for all $[\Lambda_{1*}, \Lambda_1^{*}], [\Lambda_{2*}, \Lambda_2^{*}]\in K_{bc}$ ^[30].

Definition 2.9. ^[31] The FIV-function $y:[u, v]\rightarrow K_{bc}$ is convex on $[u, v]$ if

so that $u, v\in[u, v], t\in[0, 1]$, $y(v)\succcurlyeq0$ for all $v\in[u, v]$.

Definition 2.10. ^[32] For $u, v \in K_{bc}$, $\lambda \in [0, 1]$, if $\xi:[\alpha, \beta]\rightarrow \mathbb{R}$ is the bi-function, then the invex set $K_{bc}\subset \mathbb{R}$ is defined as follows:

Definition 2.11. ^[32] The FIV-function $y:[\alpha, \beta]\rightarrow K_{bc}$ is preinvex w.r.t. $\xi$ on invex interval $[\alpha, \beta]$ if

so that $u, v\in[\alpha, \beta], t\in[0, 1]$, and $y(u)\succcurlyeq0, \forall\, v\in[\alpha, \beta]$.

Definition 2.12. ^[29] Let $P$ be a partition on $[\alpha, \beta]$ as $P = \{u = x_{1} < x_{2} < x_{3} < x_{4} < x_{5} < \dots < x_{m} = v\}$. The subintervals including the point $p$, having the maximum length, are termed as the mesh of partition $P$ given by

The collection of all possible partitions for $[u, v]$ is denoted by $P\{\delta, [\alpha, \beta]\}$ and for $p\in P(\delta, [\alpha, \beta])$, we have $\text{Mesh}(p) < \delta$. For every point $b_{j}$ belonging to the subinterval $[x_{j-1}, x_{j}]$, where $1 \leq j\leq m$, we define the Riemann sum of the function $y$ corresponding to $p \in P\{\delta, [\alpha, \beta]\}$ as

so that $y:[\alpha, \beta]\rightarrow K_{bc}$.

Definition 2.13. ^[33] Assume $y:[\alpha, \beta]\rightarrow K_{bc}$. We say that $y$ is Riemann-integrable ($\mathcal{R}$-integrable) on $[\alpha, \beta]$ if

where $B\in\mathbb{R}$ and $\delta > 0$. For each choice of $b_{j}\in[y_{j-1}, y_{j}]$, $1\leq j \leq m$, $B$ is called the $\mathcal{R}$-integral of $y$ w.r.t. $[\alpha, \beta]$, and is denoted by $B = (\mathcal{R}) \int_{u}^{v} y(t)dt$.

Definition 2.14. The integral form of gamma function is defined as given below

for $\mathbb{R}(k) > 0$

Definition 2.15. The Pochammer's symbol is defined as follows

For $k \in \mathbb{N}$ and $y \in \mathbb{C}:$

where $\Gamma$ is the gamma function.

Theorem 2.16. ^[34] Let $y:[\alpha, \beta]\subseteq \mathbb{R} \rightarrow K_{bc}$ be an FIV-function such that $y(x) = [y_{*}, y^{*}]$. Then $y$ is $\mathcal{R}$-integrable on $[\alpha, \beta]$ iff $y_{*}$ and $y^{*}$ are $\mathcal{R}$-integrable on $[\alpha, \beta]$, and

Definition 2.17. ^[35] The Prabhakar fractional linear operator on a space $L$ of functions by integral and employs fractional operator integration $I^{\mu}:L\rightarrow L$ to prove results for $a < x$, $\alpha, \beta, \rho, \lambda \in \mathbb{C}$ and $\Re(\alpha) > 0, \Re(\rho) > 0$ as follows:

and operator for $a > x$

for any real number $\alpha\geq0$ where the function

is an entire function of order $(\Re (\alpha))^{-1}$.

Definition 2.18. The left- and right-sided Prabhakar fuzzy fractional integral operators pertaining left, right end point functions, are defined as given below

where

And

Remark 2.19. Prabhakar fuzzy fractional integral operators are represented with the difference $\xi(x, a)$ in short notation as follows:

3.   On the (H-H)-type inequalities using the Prabhakar FIV-functions endowed with $h$-Godunova-Levin convexity

In the present section, we introduce a novel category of Prabhakar FIV-functions characterized by the convexity property aligned with the $h$-Godunova-Levin type. Through the application of these Prabkakar FIV-functions, we obtain (H-H)-type inequalities and delve into their applications, illustrated through various examples and cases.

Definition 3.1. The FIV-function $y:[\alpha, \beta]\rightarrow K_{bc}$ is said to be the Godunova-Levin FIV-function if

is satisfied for all $u, v \in [u, v], t \in (0, 1)$.

Definition 3.2. Consider $h:(0, 1)\rightarrow \mathbb{R}$. The FIV-function $y:[\alpha, \beta]\rightarrow K_{bc}$ is named as a convex FIV-function of the $h$-Godunova-Levin type if

is satisfied for all $u, v \in [u, v], t \in (0, 1)$.

Definition 3.3. Consider $h:(0, 1)\rightarrow \mathbb{R}$. The FIV-function $y:[\alpha, \beta]\rightarrow K_{bc}$ is named as a preinvex FIV-function of the $h$-Godunova-Levin type w.r.t. $\xi$ if

is satisfied for all $u, v \in [u, v], t \in (0, 1)$.

Now, we are going to prove the (H-H)-type inequality for convex FIV-functions of the $h$-Godunova-Levin type.

Theorem 3.4. Consider $h:(0, 1)\rightarrow \mathbb{R}$ and $y\in L_{1}[u, v]$ such that $h(w) \neq 0$ and $y:[u, v] \rightarrow K_{bc}$ is a convex FIV-function of the $h$-Godunova-Levin type with $0 < u < v$. Then the fuzzy fractional integral inequalities

hold.

Proof. In view of the $h$-Godunova-Levin convexity for the FIV-function $y$ on the closed interval $[u, v]$, for each $x, y \in [u, v]$, we have

Put $x = w u+(1-w)v, y = (1-w)u+w v$ and $\delta = \frac{1}{2}$. Hence

For each $\ell\in[0, 1]$ and for the lower FIV-function $y_*$, we can write

Multiplying $w^{\beta-1} E^{\rho}_{\alpha, \beta} \lambda w^{\alpha}$ by Eq (3.2) and integrating w.r.t. $w$ on $[0, 1]$, we obtain

Now, by substituting $r = w u+(1-w)v$ and $t = (1-w)u+w v$, it becomes

and so,

Similarly, we continue the calculation for upper FIV-function $y^{*}$, and we obtain

From inequalities (3.5) and (3.6), the inequality

is derived; that is,

Now, we again consider the $h$-Godunova-Levin convexity for the FIV-function $y$ to prove the second part of the inequality (3.1). We have

and

By adding (3.8) and (3.9), we get

For each $\ell\in[0, 1]$ and for the lower FIV-function $y_*$, we have

By multiplying $w^{\beta-1} E^{\rho}_{\alpha, \beta} \lambda w^{\alpha}$ by both sides of (3.10) and integrating on $[0, 1]$ w.r.t. $w$, we get

After simplification, we obtain

Similarly, we continue our computation for the upper FIV-function $y^{*}$, and we obtain

Combining inequalities (3.11) and (3.12), the inequality

is derived; that is,

From (3.7) and (3.13), we get

which is our desired result. □

Some special cases of the conclusion of the above theorem can be stated in the framework of several examples.

Example 3.5. In Theorem 3.4, put $h(w) = w^{s}$. Then the (H-H)-type inequality

is derived for the FIV-function with convexity of the $s$-Godunova-Levin type.

Example 3.6. In Theorem 3.4, put $h(w) = 1$. Then

Example 3.7. In Theorem 3.4, put $h(w) = \dfrac{1}{w}$. Then the (H-H)-type inequality

is derived for the convex FIV-function.

Example 3.8. In Theorem 3.4, put $h(w) = w$. Then the (H-H)-type inequality

is derived for the FIV-function with convexity of the Godunova-Levin type.

Example 3.9. (1) By choosing the value $h(w) = \dfrac{1}{w^{s}}$ in Theorem 3.4, the (H-H)-type inequality

is derived for the FIV-function with $s$-convexity.

(2) If we take the value $\lambda = 0, h(w) = \dfrac{1}{w^{s}}$ in Theorem 3.4, then

is derived for the FIV-function with $s$-convexity.

(3) If we take the value $\lambda = 0, h(w) = \dfrac{1}{w^{s}}, \beta = 1$ in Theorem 3.4, then

4.   On the trapezoidal-type inequalities using Prabhakar FIV-functions characterized by the preinvexity property of the $h$-Godunova-Levin type

In this context, we formulate a lemma specific to Prabhakar FIV-functions featuring $h$-Godunova-Levin preinvexity, providing valuable support in establishing our key findings. Indeed, within this section, we deduce fuzzy-based trapezoidal-type inequalities and explore their applications through illustrated examples and cases.

Lemma 4.1. Let $y:K = [u, u+\xi(v, u)]\rightarrow \mathbb{R}$ be a differentiable function, and $K$ be an invex set w.r.t. $\xi:K\times K\rightarrow \mathbb{R}$ so that $\xi(v, u) > 0$ for $u, v\in K$. Then

Proof. Consider the integral

For every $\ell\in[0, 1]$, we have

Now we solve the integrals of the above equation by considering it as $I_{1}$ and $I_{2}$. Simply, for $I_{1}$, we get

and so,

Using the similar inference to that of $I_{2}$, we have

Consequently,

Similarly, we continue our computations for the upper FIV-function $y^{'*}$, and we get

By combining Eqs (4.1) and (4.2), then we have

That is,

This completes the proof. □

By Lemma 4.1, we can state the following theorem.

Theorem 4.2. Consider the function $y:K = [u, u+\xi(v, u)]\rightarrow (0, \infty)$. Assume that the differentiable FIV-function $| y'|$ has the preinvexity property of the $h$-Godunova-Levin type on $K$. Then

Proof. Using Lemma 4.1, we obtain that for every $\ell \in [0, 1]$

The $h$-Godunova-Levin preinvexity property of the FIV-function $\big| y^{'}_{*}(u+w\xi(v, u), \ell)\big|$ gives

Similarly, we repeat the computations to solve it this time for the upper FIV-function $y^{'*}$. We have

From inequalities (4.4) and (4.5), we get

as required. □

Some special cases of the conclusion of the above theorem can be stated in the framework of several examples.

Example 4.3. In Theorem 4.2, put $\xi(v, u) = v-u$. Then

Theorem 4.4. Consider the function $y:K = [u, u+\xi(v, u)]\rightarrow (0, \infty)$. Assume that the differentiable FIV-function $| y'|^{q}$ has the preinvexity property of the $h$-Godunova-Levin type on $K$ so that $q = \frac{p}{p-1}$ with $p > 1$. Then

Proof. Using Lemma 4.1, and Hölder's integral inequality for every $\ell\in[0, 1]$, we obtain

where $\frac{1}{p}+\frac{1}{q} = 1$.

The $h$-Godunova-Levin preinvexity property of the FIV-function $\big| y^{'}_{*}(u+w\xi(v, u)), \ell\big|^{q}$ gives

Similarly, we solve it for the upper FIV-function $\big| y^{'*}\big|^{q}$, and we get

By combining (4.7) and (4.8), we have

Using (4.9) in (4.6), we get the required result. □

Theorem 4.5. The H-H inequality related to the hypotheses of Theorem 4.4 is as follows

Proof. Using Lemma 4.1, and thanks to the power-mean inequality, the last inequality becomes for every $\ell\in[0, 1]$

The $h$-Godunova-Levin preinvexity of the FIV-function $\big| y^{'}(u+w\xi(v, u))\big|^{q}$, we can obtain

Similarly, for the upper FIV-function $| y^{'*}|^{q}$, we get

By combining (4.12) and (4.13), we have

On the other hand, using the equation

and (4.14) and (4.12), we get the required result. □

5.   Conclusions

In this research, we introduced a novel classification of convex FIV-functions and preinvex FIV-functions, employing an arbitrary auxiliary $h$-function of the Godunova-Levin type. Utilizing this category, we derived updated versions of Hermite-Hadamard (H-H) and trapezoidal-type inequalities by implementation of newly described fuzzy fractional operators (Prabhakar fuzzy-based operators). The application of these findings was demonstrated through various examples, highlighting diverse cases. We conclude that extending Prabhakar (H-H) and Prabhakar trapezoidal-type fractional inequalities to encompass different types of fuzzy-based convexities and preinvexities under Prabhakar FIV-functions is both feasible and applicable. Looking ahead, we plan to explore this concept further within the domain of generalized convex FIV-functions, with a focus on applications in fuzzy-interval-valued non-linear programming. We anticipate that this exploration will contribute valuable insights to the broader field of fuzzy-based optimization theory and provide a foundation for other researchers in fulfilling their roles within this realm. Moreover, we conclude that such type of fractional operators, extensions, and generalizations of inequalities and its refinments could be discussed for fuzzy interval valued functions. The validated results are expected to open new avenues for future research in this evolving area. These outcomes indicate a promising route for further exploration in the sphere of integral inequalities.

Author contributions

Iqra Nayab: Conceptualization, Formal Analysis, Software, Software, Writing-Original Draft; Shahid Mubeen: Conceptualization, Formal Analysis, Writing-Review & Editing; Rana Safdar Ali: Conceptualization, Formal Analysis, Investigation, Methodology, Writing-Original Draft; Faisal Zahoor: Formal Analysis, Investigation, Software; Muath Awadalla: Investigation, Writing-Original Draft; Abd Elmotaleb A. M. A. Elamin: Investigation, Methodology, Writing-Review & Editing. All authors read and approved the final manuscript.

Use of AI tools declaration

The authors declare that they have not used Artificial Intelligence (AI) tools in the creation of this article.

Acknowledgments

This work was supported by the Deanship of Scientific Research, Vice Presidency for Graduate Studies and Scientific Research, King Faisal University, Saudi Arabia [Grant No. A166]. This study is supported via funding from Prince Sattam bin Abdulaziz University, project number (PSAU/2024/R/1445).

Conflict of interest

The authors declare no conflicts of interest.