An analysis of the isoparametric bilinear finite volume element method by applying the Simpson rule to quadrilateral meshes

1.   Introduction

Fractional calculus (FC) has long been devoted to studying integrals and derivatives of non-integer order. In recent years, researchers have been growing interest in finding the most suitable fractional derivative (FD) for modeling real-world problems. While classical calculus equations are sometimes inadequate for capturing complex systems, fractional calculus methods have proven to be highly effective in modeling phenomena observed in various disciplines such as physics, chemistry, electricity, and mechanics, thanks to their ability to handle long-memory processes ^{[1,2,3,4,5,6,7,8]}.

Traditionally, the focus of FC was primarily on Riemann-Liouville (RL) and Caputo fractional derivatives. However, as we entered the new millennium, researchers began introducing new types of fractional operators (FO), expanding the field of fractional calculus ^{[9,10,11,12,13,14,15]}. It is noteworthy that the FOs proposed in this work are distinct instances of fractional derivatives/integrals compared to those described in ^[16,17,18]. Nevertheless, they share one crucial characteristic with other FOs, which is nonlocality.

In ^[19], the authors introduced a local FD, known as the conformable derivative (CD). The conceptualization of these local FD eventually led to the rediscovery of nonlocal FD previously described in ^[14]. We provide an overview of the fundamental principles behind the CD and propose a derivative consistent with both the left and right versions. Furthermore, we find that the nonlocal FD version proposed in ^[20,21] can be derived from ^[13].

In any form of FC or calculus with derivatives, it is essential for a function's order zero to be equal to the function itself. However, the CD lacks this crucial property, which can be seen as a deficiency. To address this, the authors in ^[22,23] redefined the CD to ensure that it yields the function itself when the local FD has an order of zero. Subsequently, Sadek ^[24] introduced the FD version of this redefined CD called cotangent fractional derivative which features the that they achieve a semi-group property, the kernel operator is the exponential of the cotangent function, and is a generalization of the Riemann-Liouville FD and Caputo FD.

The existence and uniqueness of solutions represent significant qualitative properties of fractional differential problems. Several mathematicians have addressed the existence and uniqueness of solutions for fractional differential problems involving various types of FD and boundary/initial conditions (see ^{[25,26,27,28,29,30,31,32,33,34,35]}).

From the works of ^[36,37], we introduce an innovative FD called the Hilfer cotangent fractional derivative. Building upon this newly defined derivative, we investigate the solutions' existence and uniqueness for a distinct class of nonlinear fractional differential problems that feature nonlocal initial conditions. The Hilfer cotangent fractional differential equation is formulated as follows:

In this equation, $\mathcal{D}_{c^{+}}^{G, q, r}(\cdot)$ represents the Hilfer cotangent fractional derivative (HCFD), and $\mathcal{I}_{c^{+}}^{1-\sigma, r}(\cdot)$ denotes the cotangent fractional integral (CFI) with an order of $1-\sigma > 0$. Here, $c_i$ belongs to the set of real numbers, $f: L \times \mathbb{R} \rightarrow \mathbb{R}$ is a continuous function, and $\tau_i$ are elements of the interval $L$ satisfying $c < \tau_1 < \cdots < \tau_m < d$, $i = 1, \ldots, m$. The Hilfer fractional model, along with other fractional calculus models, provides a powerful mathematical framework for describing and understanding complex systems and processes that do not conform to traditional integer-order calculus. Its applications are diverse and span various scientific and engineering disciplines. Notably, the existence and uniqueness of solutions for Eq (1.1) have not been previously discussed using this new fractional definition.

The Hilfer fractional model, named after the mathematician Rudolf Gorenflo Hilfer, is a mathematical framework used to describe systems and processes that exhibit non-integer order fractional behavior. Unlike traditional integer-order calculus, fractional calculus, including the Hilfer fractional model, deals with derivatives and integrals of non-integer orders, which allows it to capture more complex dynamics in various domains. Here are some applications of the Hilfer fractional model:

● Anomalous diffusion: The Hilfer fractional model is often used to describe anomalous diffusion, where particles or information spread in a non-standard way. This can be applied to fields like physics (e.g., diffusion in porous media), biology (e.g., cell migration), and finance (e.g., modeling price movements in financial markets).

● Viscoelasticity: The Hilfer fractional model is used to describe the behavior of viscoelastic materials, which exhibit both viscous (liquid-like) and elastic (solid-like) properties. This is relevant in materials science and engineering, especially for modeling polymers and complex fluids.

● Fractional differential equations (FDEs): FDEs based on the Hilfer fractional model are used to describe a wide range of dynamic systems, such as electrical circuits, mechanical systems, and chemical reactions. These equations can capture non-standard responses and memory effects in these systems.

● Control systems: The Hilfer fractional model has applications in control theory for describing and controlling complex systems with fractional dynamics. It's used in areas like robotics, automation, and mechatronics to design controllers for systems with non-integer order behavior.

● Signal processing: Fractional calculus, including the Hilfer fractional model, is applied to signal processing for tasks like noise reduction, image processing, and data compression, where it can capture long-range dependencies in signals.

● Biomedical engineering: The Hilfer fractional model is used in modeling biological and physiological systems, such as modeling drug delivery in the human body or describing the dynamics of neurons in the brain.

● Geophysics: In geophysics, the Hilfer fractional model can be applied to model seismic phenomena, groundwater flow, and other geological processes that exhibit complex, non-integer order behavior.

● Environmental science: Environmental processes that exhibit anomalous behavior, such as pollutant transport in groundwater, can be described using the Hilfer fractional model.

● Economics: Fractional calculus models, including the Hilfer fractional model, have been applied to financial time series analysis and modeling to capture long-term dependencies and improve predictions in financial markets.

● Renewable energy: The modeling of energy storage systems, battery behavior, and renewable energy generation can benefit from the Hilfer fractional model to account for the complex and non-standard behavior of these systems.

The remainder of the manuscript is structured as follows. In Section 2, we provide a review of essential definitions and theoretical results necessary for our subsequent analysis. In Section 3, we introduce our proposed derivatives, namely the Hilfer cotangent fractional derivatives, along with their preliminary properties. Additionally, we explore the connection between an initial value problem and a cotangent Volterra integral equation, establishing uniqueness and the existence of the solution using Banach's fixed point theorem and Kransnoselskii's fixed point theorem. Furthermore, we present two examples to illustrate our findings. Finally, we conclude the paper in Section 4.

2.   Preliminaries of cotangent fractional derivative

In this section, we present preliminary information, definitions, and results related to cotangent fractional derivative, which is fundamental for the rest of the paper. Let us consider finite and infinite intervals $-\infty < c < d < \infty$ within the real numbers $\mathbb{R}_{+}$. We denote the space of continuous functions $x$ defined on $[c, d]$ as $\mathcal{C}[c, d]$. The norm for this space is defined as in ^[4]:

Additionally, we define $\mathcal{A C}^n[c, d]$, the space of functions $x$ that are $n$ times absolutely continuous differentiable, as:

The weighted space $\mathcal{C}_\sigma[c, d]$ for functions $f$ on $(c, d]$ is defined as:

where the corresponding norm:

Similarly, the weighted space $\mathcal{C}_\sigma^n[c, d]$ for functions $x$ on $(c, d]$ is defined as:

where the corresponding norm:

It is worth noting that $\mathcal{C}_\sigma^0[c, d] = \mathcal{C}_\sigma[c, d]$ when $n = 0$.

Definition 2.1. ^[4] Suppose that $x \in L^1([c, d], \mathbb{R})$. Then the RL integral of order $G > 0$ of the function $x$ is defined by

with $\Gamma(\cdot)$ is the gamma function.

Definition 2.2. ^[4] Let $x\in \mathcal{C}([c, d])$. The RL fractional derivative of order $G$ such that $n-1 < G < n, n \in \mathbb{N}$ of the function $x$ is defined by

Definition 2.3. ^[4] Let $x \in \mathcal{C}^n([c, d])$ and $p > 0$ such that $n-1 < G < n, n \in \mathbb{N}$. Then the Caputo fractional derivative of order $G$ of the function $x$ is defined by

where $x^n(\mu) = \frac{d^nx(\mu)}{d\mu^n}$.

Definition 2.4. ^[24] Let $G\in \mathbb{C}, \operatorname{Re}(G) > 0$ and $r \in(0, 1]$. The CFI of order $G$ of the function $x$ is defined by

Remark 2.5. ^[24] Let

and

Definition 2.6. ^[24] The left cotangent fractional derivative of order $G\in \mathbb{C}$ ($\operatorname{Re}(G) > 0$) and $r \in(0, 1]$ of a function $x$ is defined by

where $n = [G]+1$.

Definition 2.7. ^[24] Let $r \in(0, 1]$. Then the left-sided cotangent fractional derivative in the sense of Caputo of order $G$ of the function $x$ is defined by

where $n = [G]+1$.

Remark 2.8. For $r = 1$, Definitions 2.4–2.7 coincide the Definitions 2.1–2.3, respectively.

The cotangent fractional derivative and integral exhibit the following significant properties:

Proposition 2.9. ^[24] Let $r \in \; (0, 1]$ and $\delta, G \in \mathbb{C}$ such that $\operatorname{Re}(\delta) > 0$ and $\operatorname{Re}(G) \geq 0$. Then we have

Theorem 2.10. ^[24] Let $r \in(0, 1], \operatorname{Re}(G) > 0$, $\operatorname{Re}(q) > 0$ and $x\in \mathcal{C}([c, d], \mathbb{R})$, then we have

Theorem 2.11. ^[24] Let $x\in L^1([c, d])$. Suppose that $r \in(0, 1]$ and $0 \leq m < [\operatorname{Re}(G)]+1$. Then

Corollary 2.12. ^[24] Let $0 < \operatorname{Re}(q) < \operatorname{Re}(G)$ and $m-1 < \operatorname{Re}(q) \leq m$. Then we get

Theorem 2.13. ^[24] Let $x\in L^1([c, d]), \operatorname{Re}(G) > 0$, $n = [\operatorname{Re}(G)]+1$ and $r \in(0, 1]$. Then

Lemma 2.14. ^[24] Let $m \in \mathbb{Z}_{+}$, $G > 0$ and $r \in(0, 1]$. Then

In particular, if $m = 1$, we get

Theorem 2.15. ^[24] Let $x\in L_1(c, d)$ and $\left(I_{c^+}^{G, r}x\right)(\ell) \in A C^n[c, d]$, $\operatorname{Re}(G) > 0, n = -[-\operatorname{Re}(G)]$. Then

3.   Main results

In this section, we present the HCFD and explore its properties. Furthermore, we establish the equivalence between the given problem (1.1) and a cotangent fractional integral equation of Volterra type. By utilizing fixed point theorems, we then demonstrate the existence and uniqueness of solutions for Eq (1.1).

Definition 3.1. Let $n \in \mathbb{N}$ such that $n-1 < G < n, r \in(0, 1]$ and $0 \leq q \leq 1$. The left/right HCFD of order $G$ and type $q$ of a function $x$ is defined by

where $\mathcal{I}$ is the CFI defined in Eq (2.1).

In particular, if $n = 1$, Eq (3.1) is equivalent with

Thus, throughout this work, we discuss the case where $0 < G < 1, n = 1, 0\leq q \leq 1$ and $\sigma = G+q-G q$.

This definition can be used in many ways, including: A robust scheme for variable-order time-fractional diffusion-type equations ^[38], a fractal-fractional-order modified predator-prey mathematical model with immigrations ^[39], a pseudo-operational collocation method for variable-order time-space fractional KdV-Burgers-Kuramoto equation ^[40], Leffler-Galerkin method to solve Riccati differential equation of fractional order ^[41], stability of fractional linear systems ^[42], and controllability, observability of fractional linear systems ^[43].

Remark 3.2. It is important to mention that:

(1) If $q = 0$, we get Definition 2.6.

(2) If $q = 1$, we get Definition 2.7.

(3) The parameter $\sigma$ satisfies

Proposition 3.3. We have

Proof. From Eq (3.2) and Definition 2.6, we get

□

Let us examine the weighted function spaces defined on $(c, d]$ as follows:

and

since

Lemma 3.4. Suppose $r \in(0, 1], 0 < G < 1$, $0 \leq \sigma < 1$ and $x\in \mathcal{C}_\sigma[c, d]$, then

Proof. If we consider $x \in \mathcal{C}[c, d]$, it indicates that $x$ belongs to $\mathcal{C}_\sigma[c, d]$ and $(\ell-c)^\sigma$ belongs to $\mathcal{C}[c, d]$. Consequently, there exists a positive constant $M$ such that

and

It follows from Proposition 2.9, that

this implies that as $\ell$ approaches to $a$ from the right, the right-hand side tends to zero. □

Lemma 3.5. Let $0 < G < 1, 0 \leq q \leq 1, r \in(0, 1]$, $\sigma = G+q-G q$ and $x\in \mathcal{C}_{1-\sigma}^\sigma[c, d]$, then

and

Proof. By utilizing Theorem 2.10 and Proposition 3.3, we can proceed to

Moreover, considering Theorem 2.10 and Eq (3.2), it is evident that

□

Lemma 3.6. Suppose $x \in L^1(c, d)$ such that $\mathcal{D}_{c^{+}}^{q(1-G), p} x$ exists in $L^1(c, d)$. Then

Proof. By utilizing Definition 2.6 and Eq (3.2), we can deduce that

□

Lemma 3.7. Let $0 < G < 1$ and $0 \leq \sigma < 1$, $r \in(0, 1]$. Let $x\in \mathcal{C}_\sigma[c, d]$ and $\mathcal{I}_{c^{+}}^{1-p, r}x\in \mathcal{C}_\sigma^1[c, d]$, then

for all $\ell\in(c, d]$.

Proof. The proof follows a similar approach to the ones presented in ^[24]. □

Lemma 3.8. Let $0 < G < 1, r \in(0, 1], 0 \leq q \leq 1$ and $\sigma = G+q-G q$. Let $x \in \mathcal{C}_{1-\sigma}[c, d]$ and $\mathcal{D}_{c^{+}}^{G, q, r} x$ then $\mathcal{D}_{c^{+}}^{G, q, r} \mathcal{I}_{c^{+}}^{G, r}x$ exists in $(c, d)$ and

Proof. By employing Lemmas 3.4, 3.6, and 3.7, we obtain the following result

□

Lemma 3.9. Let $r \in(0, 1], 0 < G < 1, 0 \leq q \leq 1$, $0 < \sigma < 1$, $x\in \mathcal{C}_{1-\sigma}[c, d]$ and $\mathcal{I}_{c^{+}}^{1-\sigma, r}x$, then

Proof. From Definition 3.1 and Lemma 3.7 we have

□

3.1.   Cotangent Volterra integral equation (CVIE)

The subsequent lemma establishes the equivalence between the proposed Eq (1.1) and a CVIE.

Lemma 3.10. Consider $0 < G < 1$, $0 \leq q \leq 1$, and $\sigma = G+q-Gq$. Let $f:L \times \mathbb{R} \rightarrow \mathbb{R}$ such that $f\in \mathcal{C}_{1-\sigma}[L, \mathbb{R}]$ for any $x \in \mathcal{C}_{1-\sigma}[L, \mathbb{R}]$. If $x \in \mathcal{C}_{1-\sigma}^\sigma[L, \mathbb{R}]$, then $x$ satisfies Eq (1.1) if and only if it satisfies the CVIE given by:

where

Proof. Let $x \in \mathcal{C}_{1-\sigma}^\sigma[L, \mathbb{R}]$ be a solution of (1.1). We aim to demonstrate that $x$ is also a solution of (3.3). Based on Lemma 3.9, we can conclude that

By substituting $\ell = \tau_i$ and multiplying both sides of Eq (3.5) by $c_i$, we obtain

which implies that

From the initial condition $\mathcal{I}_{c^{+}}^{1-\sigma, r} x(c) = \sum_{i = 1}^m c_i x\left(\tau_i\right)$, we get

Hence, by substituting Eq (3.7) into Eq (3.5), we can conclude that $x(\ell)$ satisfies Eq (3.3).

Conversely, let us assume that $x \in \mathcal{C}_{1-\sigma}^\sigma$ satisfies Eq (3.3). Our goal is to show that $x$ also satisfies Eq (1.1).

By applying $\mathcal{D}_{c^{+}}^{\sigma, r}$ of Eq (3.3), and considering Proposition 2.9, Theorem 2.11, and Definition 3.1, we obtain

Since $\mathcal{D}_{c^{+}}^{G, q, r} x \in \mathcal{C}_{1-\sigma}[L, \mathbb{R}]$, by the definition of $\mathcal{C}_{1-\sigma}^\sigma[L, \mathbb{R}]$ Eq (3.8) implies that

Considering $f \in \mathcal{C}_{1-\sigma}[L, \mathbb{R}]$ and using Theorem 2.13, we observe that $\mathcal{I}_{c^{+}}^{1-q(1-G), r} f \in \mathcal{C}_{1-\sigma, r}[L, \mathbb{R}]$. Consequently, we can conclude that $\mathcal{I}_{c^{+}}^{1-q(1-G), r} f \in \mathcal{C}_{1-\sigma}^1[L, \mathbb{R}]$ based on the definition of $\mathcal{C}_\sigma^n[L, \mathbb{R}]$.

By applying $\mathcal{I}_{c^{+}}^{q(1-G), r}$ to both sides of Eq (3.8), considering Proposition 2.9, Lemma 3.7, and Definition 3.1, we obtain

Hence, its remains to show that if $x \in \mathcal{C}_{1-\sigma}^\sigma[L, \mathbb{R}]$ satisfies (3.3), it also satisfies the initial condition. So, by applying $\mathcal{I}_{c^{+}}^{1-\sigma, r}$ to both sides of Eq (3.3) and using Proposition 2.9, Theorem 2.10 and Corollary 2.12, we obtain

Taking the limit as $\ell \rightarrow c^{+}$in Eq (3.10) and the fact that $1-q < 1-G(1-r)$ get

Substituting $t = \tau_i$ and multiplying through by $c_i$ in Eq (3.3),

which implies that

Thus

So, in view of (3.11) and (3.14), we have

□

Remark 3.11. The introduced Hilfer cotangent fractional derivative (as defined in Definition 3.1) consolidates the previously existing Riemann-Liouville, cotangent, and Hilfer fractional derivatives.

3.2.   Uniqueness

In this subsection, we will present a comprehensive proof of the uniqueness of solutions for the given Eq (1.1) by employing the principles of the Banach contraction. To accomplish this, we rely on the following assumptions:

$\left(H_1\right)$ Let $f: L \times \mathbb{R} \rightarrow \mathbb{R}$ such that $f\in \mathcal{C}_{1-\sigma}^{q(1-G)}[L, \mathbb{R}]$ for any $x \in \mathcal{C}_{1-\sigma}^\sigma[L, \mathbb{R}]$.

$\left(H_2\right)$ There exists a constant $K > 0$ such that

for any $z, \bar{z} \in \mathbb{R}$ and $\ell\in L$.

$\left(H_3\right)$ Suppose that $K \psi < 1,$ where

and $\mathcal{B}(\sigma, G)$ is the Beta function ^[4]

Theorem 3.12. Let $0 < G < 1$, $0 \leq q \leq 1$, and $\sigma = G+q-G q$. Assuming that the conditions $\left(H_1\right)$–$\left(H_3\right)$ hold. Then the Eq (1.1) possesses a unique solution within the space $\mathcal{C}_{1-\sigma}^\sigma[L, \mathbb{R}]$.

Proof. Let the operator $T: \mathcal{C}_{1-\sigma}[L, \mathbb{R}] \rightarrow \mathcal{C}_{1-\sigma}[L, \mathbb{R}]$ by

Consequently, the operator $T$ is unambiguously defined. For arbitrary $x_1, x_2 \in \mathcal{C}_{1-\sigma}[L, \mathbb{R}]$ and $\ell\in L$, the following expression holds:

Since $\left|e^{-\cot(r\frac{\pi}{2})\ell}\right| < 1$, we get

Therefore,

Therefore, by considering (3.16), it can be deduced that $T$ behaves as a contraction map. As a result of the Banach contraction principle, it can be firmly stated that Eq (1.1) possesses a unique solution. □

3.3.   Existence

We utilize the Kransnoselskii's fixed point theorem ^[44] to establish the existence of solutions for the given Eq (1.1).

$\left(H_4\right)$ Suppose that $K \Delta < 1,$ where

Theorem 3.13. Assume that the conditions $\left(H_1\right), \left(H_2\right)$, and $\left(H_4\right)$ are fulfilled, let $0 < G < 1, 0 \leq q \leq 1$, and $\sigma = G+q-Gq$. Then problem (1.1) possesses at least one solution within the space $\mathcal{C}_{1-\sigma}^\sigma[L, \mathbb{R}]$.

Proof. In Step 1, we have $\|\eta\|_{\mathcal{C}_{1-\sigma}[L, \mathbb{R}]} = \sup_{\ell \in L}\left|(\ell-c)^{1-\sigma} \eta(\ell)\right|$ and choose $\kappa \geq M\|\eta\|_{\mathcal{C}_{1-\sigma}[L, \mathbb{R}]}$, where

Let us $\mathbf{B}_\kappa = \left\{x \in \mathbb{C}[L, \mathbb{R}]:\|x\|_{\mathcal{C}_{1-\sigma}[L, \mathbb{R}]} \leq \kappa\right\}$ and the operators $\mathbb{T}_1$ and $\mathbb{T}_2$ on $\mathbf{B}_k$ by

for each $\ell\in[c, d]$, the following inequality holds true. Now, considering any $x$ and $y$ belonging to the set $\mathbf{B}_k$, we have the following:

This implies that $\mathbb{T}_1 x+\mathbb{T}_2 y \in \mathbf{B}_\kappa$.

In Step 2, we demonstrate that the operator $\mathbb{T}_2$ is a contraction.

Now, let $x, y \in \mathcal{C}_{1-\sigma}[L, \mathbb{R}]$ and $\ell\in L$, then

Consequently, it can be deduced from $\left(H_4\right)$ that $\mathbb{T}_2$ is a contraction.

In Step 3, we establish the continuity and compactness of the operator $\mathbb{T}_1$.

Evidently, the operator $\mathbb{T}_1$ is continuous since the function $f$ is continuous. Therefore, for any $x \in \mathcal{C}_{1-\sigma}[L, \mathbb{R}]$, we obtain:

This shows that the operator $\mathbb{T}_1$ is uniformly bounded on $\mathbf{B}_k$. Thus, it remains to show that $\mathbb{T}_1$ is compact. Let $\sup_{(\ell, x) \in L \times \mathbf{B}_{\mathbb{R}}}|f(\ell, x(\ell))| = \delta < \infty$ and for any $c < \tau_1 < \tau_2 < d$,

By applying the Arzelà-Ascoli theorem ^[45], we can conclude that the operator $\mathbb{T}_1$ is compact on $\mathbf{B}\kappa$. Consequently, Eq (1.1) possesses at least one solution. □

3.4.   Examples

Example 3.14. Let us consider a fractional differential equation that incorporates the Hilfer cotangent fractional derivative in the following manner:

By comparing (1.1) with (3.26), we get $G = \frac{2}{3}, q = \frac{1}{2}, r = 1, \sigma = \frac{5}{6}, c = 0, d = 2, c_1 = 2$ since $m = 1, \tau_1 = \frac{2}{5} \in L$ and

Consequently, the function $f$ is continuous, and for any $u, v \in \mathbb{R}_{+}$ and $\ell\in L$, the inequality $|f(\ell, u)-f(\ell, v)| \leq \frac{1}{25}|u-v|$ holds. Hence, we can affirm that both conditions $\left(H_1\right)$ and $\left(H_3\right)$ are satisfied with $K = \frac{1}{25}$. Through straightforward calculations, we obtain

and

These values imply that $K \psi = 0.1345251830 < 1$. As a result, all the requirements of Theorem 3.12 are fulfilled, leading to the conclusion that problem (1.1) possesses a unique solution on $L$.

Similarly, we determine that

and

Since all the conditions of Theorem 3.13 are met, we can infer that problem (1.1) has at least one solution on $L$.

Example 3.15. Let's examine the Hilfer cotangent fractional differential equation given by the following expression:

By following the same procedure as demonstrated in Example 3.14, we obtain the values

and

Consequently, we have

Based on Theorem 3.12, it can be concluded that problem (1.1) possesses a unique solution on $L$. Moreover, we have

which implies, according to Theorem 3.13, that problem (1.1) has at least one solution on $L$.

4.   Conclusions

In the present study, we introduced the Hilfer cotangent fractional derivatives within the framework of fractional calculus. Leveraging well-established theorems from fixed point theory, we were able to establish the existence and uniqueness of solutions for a particular class of fractional initial value problems incorporating the Hilfer cotangent fractional derivative. To illustrate the efficacy of our findings, we provided illustrative examples. Notably, the Hilfer cotangent fractional derivative encompasses three parameters, offering greater flexibility. This expanded parameter space becomes particularly valuable when considering stability and other qualitative characteristics of differential equations involving fractional derivatives.

Use of AI tools declaration

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

Conflict of interest

The authors declare that they have no conflicts of interest.