Investigation of fractional-order pantograph delay differential equations using Sumudu decomposition method

Asrar Saleh Alsulami; Mariam Al-Mazmumy; Maryam Ahmed Alyami; Mona Alsulami; Asrar Saleh Alsulami; Mariam Al-Mazmumy; Maryam Ahmed Alyami; Mona Alsulami

doi:10.3934/math.20241702

AIMS Mathematics

2024, Volume 9, Issue 12: 35910-35930. doi: 10.3934/math.20241702

Previous Article Next Article

Research article Special Issues

Investigation of fractional-order pantograph delay differential equations using Sumudu decomposition method

Department of Mathematics and Statistics, Faculty of Science, University of Jeddah, Jeddah 23218, Saudi Arabia

This paper combines the Sumudu transform with the Adomian decomposition method to address Caputo-type fractional-order pantograph delay differential equations. It features numerical evaluations that confirm the effectiveness of the proposed methods. The study introduces a powerful computational technique for solving these equations, providing results that establish its efficiency and relevance through comparisons with existing methods. The findings underscore both the efficiency and accuracy of the proposed algorithm.

Keywords:

Citation: Asrar Saleh Alsulami, Mariam Al-Mazmumy, Maryam Ahmed Alyami, Mona Alsulami. Investigation of fractional-order pantograph delay differential equations using Sumudu decomposition method[J]. AIMS Mathematics, 2024, 9(12): 35910-35930. doi: 10.3934/math.20241702

Related Papers:

[1]	Maysaa Al-Qurashi, Saima Rashid, Fahd Jarad, Madeeha Tahir, Abdullah M. Alsharif . New computations for the two-mode version of the fractional Zakharov-Kuznetsov model in plasma fluid by means of the Shehu decomposition method. AIMS Mathematics, 2022, 7(2): 2044-2060. doi: 10.3934/math.2022117
[2]	Saima Rashid, Sobia Sultana, Bushra Kanwal, Fahd Jarad, Aasma Khalid . Fuzzy fractional estimates of Swift-Hohenberg model obtained using the Atangana-Baleanu fractional derivative operator. AIMS Mathematics, 2022, 7(9): 16067-16101. doi: 10.3934/math.2022880
[3]	Rasool Shah, Abd-Allah Hyder, Naveed Iqbal, Thongchai Botmart . Fractional view evaluation system of Schrödinger-KdV equation by a comparative analysis. AIMS Mathematics, 2022, 7(11): 19846-19864. doi: 10.3934/math.20221087
[4]	Wedad Albalawi, Muhammad Imran Liaqat, Fahim Ud Din, Kottakkaran Sooppy Nisar, Abdel-Haleem Abdel-Aty . Well-posedness and Ulam-Hyers stability results of solutions to pantograph fractional stochastic differential equations in the sense of conformable derivatives. AIMS Mathematics, 2024, 9(5): 12375-12398. doi: 10.3934/math.2024605
[5]	Yousef Jawarneh, Humaira Yasmin, Abdul Hamid Ganie, M. Mossa Al-Sawalha, Amjid Ali . Unification of Adomian decomposition method and ZZ transformation for exploring the dynamics of fractional Kersten-Krasil'shchik coupled KdV-mKdV systems. AIMS Mathematics, 2024, 9(1): 371-390. doi: 10.3934/math.2024021
[6]	Muhammad Imran Liaqat, Sina Etemad, Shahram Rezapour, Choonkil Park . A novel analytical Aboodh residual power series method for solving linear and nonlinear time-fractional partial differential equations with variable coefficients. AIMS Mathematics, 2022, 7(9): 16917-16948. doi: 10.3934/math.2022929
[7]	Mariam Al-Mazmumy, Mona Alsulami, Norah AL-Yazidi . Algorithms and applications of the new modified decomposition method to solve initial-boundary value problems for fractional partial differential equations. AIMS Mathematics, 2025, 10(3): 5806-5829. doi: 10.3934/math.2025267
[8]	Yousef Jawarneh, Humaira Yasmin, M. Mossa Al-Sawalha, Rasool Shah, Asfandyar Khan . Numerical analysis of fractional heat transfer and porous media equations within Caputo-Fabrizio operator. AIMS Mathematics, 2023, 8(11): 26543-26560. doi: 10.3934/math.20231356
[9]	Muammer Ayata, Ozan Özkan . A new application of conformable Laplace decomposition method for fractional Newell-Whitehead-Segel equation. AIMS Mathematics, 2020, 5(6): 7402-7412. doi: 10.3934/math.2020474
[10]	M. Mossa Al-Sawalha, Azzh Saad Alshehry, Kamsing Nonlaopon, Rasool Shah, Osama Y. Ababneh . Approximate analytical solution of time-fractional vibration equation via reliable numerical algorithm. AIMS Mathematics, 2022, 7(11): 19739-19757. doi: 10.3934/math.20221082

Abstract

1. Introduction

Fractional calculus is a branch of mathematics that generalizes integration and differentiation to arbitrary orders, leading to the study of fractional differential equations (FDEs), which have gained prominence in recent times. Its origins trace back to the early conjectures of Leibniz in 1695 and in 1730. The literature offers numerous definitions of fractional derivatives, such as Riemann-Liouville, Caputo, Grunwald-Letnikov, Weyl, Marchaud, and Prabhakar, among others (see ^[1,2,3] for a historical overview). Fractional calculus finds practical applications across various fields, including economics ^[4], chaotic systems ^[5], dynamics of viscoelastic materials ^[6], electrochemical processes ^[7], traffic models in fluid dynamics ^[8], and optics ^[9]. Fractional derivatives are particularly useful for describing the historical profiles and genetic characteristics of diverse processes. The significance of FDEs has prompted numerous investigations into the existence and uniqueness of valid solutions ^[2,10]. Traditional mathematical modeling often assumes that a system's current state is sufficient for description; however, this assumption can be inadequate. In many cases, especially for dynamical systems exhibiting delays common in fields like economics, control systems, thermoelasticity, hydraulic networks, and biology, it becomes necessary to incorporate past states for its accurate modeling. These time-delay systems can be represented by delay differential equations (DDEs), a specific class of differential equations that relate the derivative of a function at a given time to its values at earlier times. DDEs encompass hereditary systems, systems with dead time, and equations with different arguments. They are particularly useful for modeling time-dependent processes, as they integrate past information into ordinary differential equations. Fractional delay differential equations (FDDEs), which incorporate delay terms within the framework of fractional calculus, offer greater flexibility in modeling than standard DDEs. Recently, FDDEs have attracted significant attention from researchers, as even minor delays can substantially impact system behavior. These equations find applications in diverse scientific fields, including control theory, chemistry, economics, finance, and vibration theory ^[11,12]. For mathematical insights into FDDEs, we recommend ^[13] for information on existence and uniqueness and ^[14,15] for findings related to the stability of solutions.

Pantograph equations are a specific type of DDEs that incorporate proportional delay terms and are widely used to model various real-life phenomena. The term "pantograph" originates from a device first used in constructing electric locomotives in 1851. The mathematical modeling of this device was established in 1971 ^[16]. These equations frequently arise in diverse scientific fields, including population dynamics, electrodynamics, dynamical systems, number theory, and quantum mechanics. Both analytical and numerical techniques have been employed historically and currently to derive closed-form and approximate solutions for fractional pantograph equations. For instance, multi-pantograph systems were addressed using the spectral Tau method in ^[17], while the Genocchi operational matrix in the collocation method was utilized in ^[18]. Additionally, the classical operational matrix method ^[19] and Montez-Legendre polynomials ^[20] have been applied to computationally tackle pantograph equations with fractional-order derivatives.

In the mid-1980s, George Adomian Introduced a valuable technique known as the Adomian decomposition method (ADM) ^[21], which effectively provides both closed-form and approximate solutions for linear and nonlinear functional equations. The ADM has gained increasing attention for its convergence properties and the stability of its solutions. Numerous studies have explored these aspects, emphasizing the method's robustness without reiterating the same findings (see ^{[22,23,24,25]}). The ADM has become widely recognized as a powerful tool for solving various FDEs. Applications include the analysis of fractional Bernoulli's equations ^[26] and the study of heat transfer processes with fractional components ^[27]. To enhance its accuracy and efficiency, the ADM has undergone several modifications designed to enhance convergence speed and reduce computational time. These advancements have led to significant progress, showing faster convergence of series solutions compared to the standard ADM. The modified ADM has proven to be computationally efficient across various models, making it invaluable for researchers in applied science. Recently, Masood et al. ^[28] proposed effective modifications to the ADM in addressing initial-boundary value problems for diffusion equations with fractional-order derivatives. Additionally, another efficient modification was applied by ^[29] to solve nonlinear FDEs. Furthermore, researchers have integrated the ADM with various integral transformations, such as Laplace, Sumudu, Natural, and Elzaki, which have proven to be robust for solving different FDE models ^[30,31,32].

In this context, Watugala ^[33] introduced the Sumudu integral transform as a modification of the well-known Laplace transform. This new approach offers several motivational advantages, including:

● It preserves units, facilitating problem-solving without the need for frequency domain conversion, which is beneficial in physical sciences where dimensions are crucial.

● The Sumudu integral transform is a linear operator, and conserves linear functions; this means that the units and dimensions remain unchanged.

● The Sumudu transform (in contrast with the Laplace transform) accurately handles initial conditions with singularities in differential equations, typical in practical engineering models.

● Calculating the inverse using the Sumudu transform avoids the complexities of contour integration, offering a simpler solution approach.

Furthermore, the authors in ^[34,35] introduced several fundamental properties of the Sumudu transformation, which are instrumental in constructing and solving mathematical models. Additionally, the authors in ^[36,37] discussed the combination of the Sumudu transform and the ADM to tackle various models of FDEs.

However, this study aims to utilize the coupling of the ADM and the Sumudu transform, referred to as the Sumudu decomposition method (SDM), to find approximate closed-form solutions for FPDDEs. The proposed method effectively addresses the complexity of these problems while maintaining ease of use. Error comparisons demonstrate that the devised SDM is exceptionally precise.

The structure of the manuscript is as follows: Section 2 provides definitions related to fractional calculus and the Sumudu transform. Section 3 discusses the solution method, while Section 4 applies the SDM to various numerical problems, determining approximate solutions for different FPDDEs. Finally, Section 5 presents the conclusions of the study.

2. Main concepts and theorems

This section presents key definitions and properties of fractional calculus. Additionally, it discusses fundamental features of the Sumudu transform (ST).

2.1. Some basic definitions and properties of the fractional calculus

Definition 2.1. Let $y(x) \in C([a, b])$ and $a < x < b$ . The Riemann-Liouville fractional integral operator of the fractional-order $\alpha > 0$ is given by ^[1,2,10]

$\begin{equation} I^{\alpha}_{a}y(x) = \frac{1}{\Gamma(\alpha)}\displaystyle{\int}_{a}^{x} y(t)(x - t)^{\alpha - 1} dt. \end{equation}$

(2.1)

For $\alpha = 0$ , this reduces to the identity operator, i.e., $I^{0}_{a}y(x) = y(x)$ . Moreover, when $\alpha \in \mathbb{N}$ , $I^{\alpha}_{a}y(x)$ coincides with the classical integral operator.

Definition 2.2. Let $y(x) \in C([a, b])$ . The Caputo fractional derivative of the fractional-order $\alpha > 0$ is defined as ^[1,2,10]

$\begin{equation} D^{\alpha}_{*}y(x) = \begin{cases} I^{m - \alpha} y^{(m)}(x), & m-1 < \alpha < m, \\ \frac{d^{m}}{dx^{m}} y(x), & \alpha = m, \end{cases} \end{equation}$

(2.2)

where $m = \lceil \alpha \rceil$ and $m \in \mathbb{N}$ .

Additionally, several crucial properties of the fractional differential and integral operators based on the previously mentioned definitions are outlined as follows ^[1,2,10]:

(1)

$\begin{equation} I^{\alpha}_{a}I^{\beta}_{a}y(x) = I^{\alpha + \beta}_{a}y(x), \quad \alpha, \beta \geq 0; \end{equation}$

(2.3)

(2)

$\begin{equation} I^{\alpha}_{a}D^{\alpha}_{*}y(x) = y(x) - \sum\limits_{k = 0}^{m-1} y^{(k)}(0^+) \frac{x^k}{\Gamma(k+1)}, \text{ where } m-1 < \alpha\leq m ; \end{equation}$

(2.4)

(3)

$\begin{equation} I^{\alpha}_{a} x^n = \frac{\Gamma(n+1)}{\Gamma(\alpha + n + 1)} x^{n+\alpha}, \text{ where } x > 0 , \ m-1 < \alpha\leq m, \ n > -1 ; \end{equation}$

(2.5)

(4)

$\begin{equation} D^{\alpha}_{*} C = 0, \text{ where } C \text{ is a real constant;} \end{equation}$

(2.6)

(5)

$\begin{equation} D^{\alpha}_{*}x^{n} = \begin{cases} 0, & n < m-1, \\ \frac{\Gamma(n+1)}{\Gamma(n+1-\alpha)} x^{n-\alpha}, & n \geq m-1. \end{cases} \end{equation}$

(2.7)

Moreover, since this study focuses on solving Caputo-type FPDDEs, the Caputo fractional derivative is chosen. This derivative allows for the specification of additional initial conditions, ensuring a unique and precise solution for the FDEs. It is widely recognized for its capacity to define initial conditions based on the integer-order derivatives of the functions involved in the models under consideration. Additionally, it is worth noting that numerous studies have thoroughly explored the geometric interpretation of fractional derivatives ^[38].

2.2. Some basic definitions and properties of the Sumudu transform

Definition 2.3. The ST over the following set ^[33,34]:

$\begin{eqnarray} A = \{y(x)\; |\; \exists \; M, \nu_{1},\nu_{2} > 0, \; \vert y(x)\vert < Me^{\vert x \vert/\nu_{i}}, \; \text{if}\; \; x \in (-1)^{i} \times{[0,\infty)} \}, \end{eqnarray}$

(2.8)

of functions is defined by

$\begin{eqnarray} \mathbb{S}[y(x)] = \int_{0}^{\infty}y(ux)e^{-x}dx = G(u), \; \; \; u\in(-\nu_{1},\nu_{2}), \end{eqnarray}$

(2.9)

where $u$ is a ST parameter, real or complex, and independent of $x$ . In addition, the formula for inverting the ST is expressed as follows:

$\begin{eqnarray} \mathbb{S}^{-1}[G(u)] = \frac{1}{2 \pi i }\int_{\lambda-i\infty}^{\lambda+i\infty} \frac{e^{ux}}{u} G\left(\frac{1}{u}\right) du = y(x). \end{eqnarray}$

(2.10)

Definition 2.4. The ST of the Caputo derivative is expressed as ^[33,34]

$\begin{eqnarray} \mathbb{S}[D^{\alpha}y(x)] = u^{-\alpha}\mathbb{S}[y(x)]-\sum\limits_{k = 0}^{m-1} u^{-\alpha+k}y^{(k)}(0),\; \; m-1 < \alpha\leq m. \end{eqnarray}$

(2.11)

Furthermore, to illustrate how the adopted integral transform operates, Table 1 presents the Sumudu transform of several elementary functions ^[33,34]. Additionally, the authors in ^[34] demonstrate the equivalence of the Laplace and Sumudu transformations. However, one may prove more convenient for computations than the other, particularly when addressing FDEs. Consequently, this study selects the Sumudu transform for solving the FPDDEs.

Table 1. ST of some functions.

$y(x)$	$\mathbb{S}\{y(x)\}=Y(u)$
1	1
$e^{ax}$	$\frac{1}{1-au}$
$\frac{x^{n-1}}{\Gamma(n)}$	$u^{n-1}$ for $n > 0$
$\cos ax$	$\frac{1}{1+a^{2}u^{2}}$
$\sin ax$	$\frac{au}{1+a^{2}u^{2}}$

| Show Table

DownLoad: CSV

3. Treatment of fractional pantograph delay differential equations via Sumudu decomposition method

To give the outline of the deployed SDM on the class of FPDDEs, the current section makes consideration of the generalized initial-value problem for FDDE with delay terms $q_1$ and $q_2$ as follows:

$\begin{equation} D^{\alpha}_{*}y(x) = F(x^{\alpha},y(x),y(q_{1}x),y(q_{2}x))+g(x), \end{equation}$

(3.1)

$\begin{equation} y^{i}(0) = y_{i}, \ \ i = 1,\dots, m-1, \end{equation}$

(3.2)

where $x\in I = [0, T]$ , $m-1 < \alpha\leq m$ , $0 < q_{1}, q_{2} < 1$ , $F$ is an arbitrary continuous nonlinear function, $g(x)$ is the given source function, $y_{i}\in \mathbb{R},$ while $D^{\alpha}_{*}$ is the Caputo fractional derivative of order $\alpha.$

To solve the problem expressed in (3.1) and (3.2) via the deployed method, SDM starts by operating the ST on both sides of the governing equation to obtain the following:

$\begin{equation} \mathbb{S}[D^{\alpha}_{*}y(x)] = \mathbb{S}[F(x^{\alpha},y(x),y(q_{1}x),y(q_{2}x))]+\mathbb{S}[g(x)]. \end{equation}$

(3.3)

Further, on using (2.11) together with the prescribed initial conditions, one obtains

$\begin{equation} u^{-\alpha}\mathbb{S}[y(x)]-\sum\limits_{i = 1}^{m-1}u^{i-\alpha}y^{(i)}(0) = \mathbb{S}[F(x^{\alpha},y(x),y(q_{1}x),y(q_{2}x))]+\mathbb{S}[g(x)], \end{equation}$

(3.4)

or alternatively

$\begin{equation} \mathbb{S}[y(x)]-\sum\limits_{i = 1}^{m-1}u^{i}y_{i} = \frac{1}{u^{-\alpha}}\mathbb{S}[F(x^{\alpha},y(x),y(q_{1}x),y(q_{2}x))]+\frac{1}{u^{-\alpha}}\mathbb{S}[g(x)]. \end{equation}$

(3.5)

Next, upon applying the inverse ST on both sides of the latter equation, one then obtains

$\begin{equation} y(x)-\mathbb{S}^{-1}\bigg[\sum\limits_{i = 1}^{m-1}u^{i}y_{i}\bigg] = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[F(x^{\alpha},y(x),y(q_{1}x),y(q_{2}x))\bigg]\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[g(x)\bigg]\bigg]. \end{equation}$

(3.6)

Further, the standard ADM procedure proceeds by expressing the solution $y(x)$ by the following infinite series:

$\begin{equation} y(x) = \sum\limits_{n = 0}^{\infty}y_{n}(x), \end{equation}$

(3.7)

while the nonlinear term $F(x^{\alpha}, y(x), y(q_{1}x), y(q_{2}x))$ is expressed through the summation series of Adomian polynomials $A_{n}$ as follows:

$\begin{equation} F(x^{\alpha},y(x),y(q_{1}x),y(q_{2}x)) = \sum\limits_{n = 0}^{\infty}A_{n}, \end{equation}$

(3.8)

where $A_{n}$ are the Adomian polynomials, unequivocally expressed as follows:

$\begin{equation} A_{n} = \frac{1}{n!}\frac{d^{n}}{d\lambda^{n}}\bigg[N(\sum\limits_{i = 0}^{\infty}\lambda^{i}y_{i})\bigg]_{\lambda = 0}, \ \ n = 0,1,2,\dots \end{equation}$

(3.9)

Therefore, with the above-decomposed series, Eq (3.6) is thus re-expressed as follows:

$\begin{equation} \sum\limits_{n = 0}^{\infty}y_{n}(x) = \mathbb{S}^{-1}\bigg[\sum\limits_{i = 1}^{m-1}u^{i}y_{i}\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[\sum\limits_{n = 0}^{\infty}A_{n}\bigg]\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[g(x)\bigg]\bigg], \end{equation}$

(3.10)

upon which the formal recursive relation is acquired from the standard ADM procedure as follows:

$\begin{equation} \begin{cases} y_{0}(x) = \mathbb{S}^{-1}\bigg[\sum\nolimits_{i = 1}^{m-1}u^{i}y_{i}\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\big[g(x)\big]\bigg] = p(x),\\ y_{n}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\big[A_{n-1}\big]\bigg], \ n\geq 1. \end{cases} \end{equation}$

(3.11)

Notably, it is crucial to note that the standard ADM involves defining the zeroth component $y_{0}(x)$ based on the function $p(x)$ derived from the source term and the prescribed initial conditions. However, in this context, Wazwaz ^[39] proposed an extension of the standard ADM by further decomposing the zeroth component $y_{0}(x)$ into two components as follows:

$\begin{equation*} y_{0}(x) = p_{0}(x)+p_{1}(x), \end{equation*}$

where the first part is associated with the initial term, and the second part is attached to the subsequent iterate. Thus, under this assumption, the new modified recursive algorithm from (3.11) takes the following form:

$\begin{equation} \begin{cases} y_{0}(x) = p_{0}(x),\\ y_{1}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\big[A_{0}\big]\bigg]+p_{1}(x),\\ y_{n}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\big[A_{n-1}\big]\bigg], \ n\geq 2. \end{cases} \end{equation}$

(3.12)

Certainly, the closed-form solution determined in Eq (3.12) (of course after taking the net sum of the components, $y(x) = \sum_{n = 0}^{\infty}y_{n}(x)$ ) converges more rapidly than the standard ADM, as shown in ^[39]. In addition, the modified scheme by Wazwaz heavily relies on the selection of the right $p_{0}(x)$ and $p_{1}(x)$ functions; this dependency will be demonstrated in the next application section. In Algorithm 1, we systematically outline our method for solving FPDDSs.

Algorithm 1. Solution procedure for FPDDEs using SDM
Input: Define the FPDDE with given fractional derivative, nonlinear terms, delay terms, and initial conditions.
Step 1: Apply the Sumudu transform (ST) to the given equation, incorporating the initial conditions.
Step 2: Use the inverse ST to reformulate the problem, expressing the solution in terms of a series and Adomian polynomials.
Step 3: Decompose the solution and nonlinear terms into infinite series. Represent the nonlinear terms using Adomian polynomials.
Step 4: Derive recursive relations for the solution components using inverse ST.
Step 5: For Wazwaz's modified approach, decompose the initial component into two parts for improved convergence.
Step 6: Iteratively compute the solution components until convergence is achieved.
Output: Construct the approximate solution as the sum of all computed components.

4. Applications

This section presents a numerical demonstration of the derived Sumudu decomposition method (SDM) on several initial-value problems (IVPs) for fractional pantograph delay differential equations (FPDDEs). Additionally, graphical illustrations will be provided to compare the obtained approximate solutions with available exact solutions. Furthermore, comparison tables will be included to evaluate the effectiveness of the deployed SDM against other methods found in the literature. All calculations are performed, and the plots are generated by using Maple 22 with 8 digits precision.

Example 4.1. Consider the IVP for FPDDE as follows ^[40,41,42]:

$\begin{equation} \begin{cases} D^{1/2}_{*}y(x) = -y(x)+y(\frac{x}{2})+\frac{7}{8}x^{3}+\frac{16}{5\Gamma(\frac{1}{2})}x^{\frac{5}{2}},\\ y(0) = 0. \end{cases} \end{equation}$

(4.1)

To solve (4.1) by SDM, we apply ST on both sides of the equation, which when using (2.11) gives

$\begin{equation} \mathbb{S}[y(x)] = \frac{1}{u^{-1/2}}\mathbb{S}[-y(x)]+\frac{1}{u^{-1/2}}\mathbb{S}\bigg[y(\frac{x}{2})\bigg]+\frac{1}{u^{-1/2}}\mathbb{S}\bigg[\frac{7}{8}x^{3}+\frac{16}{5\Gamma(\frac{1}{2})}x^{\frac{5}{2}}\bigg]. \end{equation}$

(4.2)

Now, applying the inverse ST on both sides of the above equation, one gets

$\begin{equation} y(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}[-y(x)]\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[y(\frac{x}{2})\bigg]\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[\frac{7}{8}x^{3}+\frac{16}{5\Gamma(\frac{1}{2})}x^{\frac{5}{2}}\bigg]\bigg], \end{equation}$

(4.3)

or equally expressed after deploying the standard ADM procedure on the latter equation as follows:

$\begin{equation} \sum\limits_{n = 0}^{\infty}y_{n}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[-\sum\limits_{n = 0}^{\infty}y_{n}(x)\bigg]\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[\sum\limits_{n = 0}^{\infty}y_{n}(\frac{x}{2})\bigg]\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[\frac{7}{8}x^{3}+\frac{16}{5\Gamma(\frac{1}{2})}x^{\frac{5}{2}}\bigg]\bigg]. \end{equation}$

(4.4)

Moreover, the formal recursive relationship for the governing fractional IVP is obtained as follows:

$\begin{equation} \begin{cases} y_{0}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[\frac{7}{8}x^{3}+\frac{16}{5\Gamma(\frac{1}{2})}x^{\frac{5}{2}}\bigg]\bigg] = x^{3}+\frac{4}{5\sqrt{\pi}}x^{\frac{7}{2}},\\ y_{n}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[-y_{n-1}(x)\bigg]\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[y_{n-1}(\frac{x}{2})\bigg]\bigg], \ n\geq 1. \end{cases} \end{equation}$

(4.5)

In addition, upon using the modified recursive algorithm by Wazwaz in (3.12), the governing model via the modified SDM yields the following recursive algorithm as follows:

$\begin{equation} \begin{cases} &y_{0}(x) = x^{3},\\ &y_{1}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[-x^{3}\bigg]\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[\frac{x^{3}}{8}\bigg]+\frac{4}{5\sqrt{\pi}}x^{\frac{7}{2}} \\ &\qquad \ \ = -\frac{4}{5\sqrt{\pi}}x^{\frac{7}{2}}+\frac{4}{5\sqrt{\pi}}x^{\frac{7}{2}} \\ &\qquad \ \ = 0, \\ & y_{n}(x) = 0, \ \ n\ge 2, \end{cases} \end{equation}$

(4.6)

upon which the exact solution is arrived at after taking the net sum of the above solution components as follows:

$\begin{equation} y(x) = x^{3}. \end{equation}$

(4.7)

In this context, the authors in ^[40,41,42] employed various methods to tackle the same model in (4.1), using the regenerative kernel method (RKM), the piecewise Picard iteration method (PPIM), and the generalized Legendre polynomial configuration method (GLPCM), respectively. Each of these authors obtained approximate solutions by calculating nine terms, while our algorithm achieves an exact solution using only two terms. This advantage leads to more effective convergence to the exact solution compared to numerical methods, significantly reducing computational complexity. Additionally, the study compares the resulting errors from these methods with those produced by the devised SDM, as detailed in Table 2.

Table 2. Comparison of absolute errors at various

$x_i$ values for the fractional IVP (4.1).

$x$	RKM ^[40]	PPIM ^[41]	GLPCM ^[42]	Present method (SDM)
0.1	$1.475\times 10^{-7}$	$1.303\times 10^{-13}$	$6.328\times 10^{-15}$	0
0.2	$4.492\times 10^{-6}$	$3.731\times 10^{-11}$	$6.212\times 10^{-15}$	0
0.3	$2.361 \times 10^{-5}$	$1.171 \times 10^{-9}$	$5.648 \times 10^{-15}$	0
0.4	$8.239 \times 10^{-5}$	$1.351 \times 10^{-8}$	$2.914 \times 10^{-15}$	0
0.5	$2.107 \times 10^{-4}$	$2.107 \times 10^{-4}$	$5.065 \times 10^{-15}$	0
0.6	$4.403 \times 10^{-4}$	$4.403 \times 10^{-4}$	$4.219 \times 10^{-15}$	0
0.7	$8.037 \times 10^{-4}$	$8.037 \times 10^{-4}$	$4.163 \times 10^{-15}$	0
0.8	$1.333 \times 10^{-3}$	$1.333 \times 10^{-3}$	$3.220 \times 10^{-15}$	0
0.9	$2.064 \times 10^{-3}$	$2.064 \times 10^{-3}$	$2.665 \times 10^{-15}$	0
1	$1.303 \times 10^{-3}$	$1.303 \times 10^{-3}$	$2.331 \times 10^{-15}$	0

| Show Table

DownLoad: CSV

Example 4.2. Consider the nonlinear IVP for FPDDE as follows ^{[42,43,44,45]}:

$\begin{equation} \begin{cases} D^{\alpha}_{*}y(x) = 1-2y^{2}\bigg(\frac{x}{2}\bigg), \; \; 0 < \alpha\leq 1,\\ y(0) = 0. \end{cases} \end{equation}$

(4.8)

To solve (4.8) by SDM, we apply ST on both sides of the equation to accordingly obtain

$\begin{equation} \mathbb{S}[y(x)] = \frac{1}{u^{-\alpha}}\mathbb{S}[1]-\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[2y^{2}\bigg(\frac{x}{2}\bigg)\bigg]. \end{equation}$

(4.9)

Next, upon applying the inverse ST on the latter equation, one obtains

$\begin{equation} y(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}[1]\bigg]-\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[2y^{2}\bigg(\frac{x}{2}\bigg)\bigg]\bigg], \end{equation}$

(4.10)

such that when the ADM is deployed reveals

$\begin{equation} \sum\limits_{n = 0}^{\infty}y_{n}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}[1]\bigg]-\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[2\sum\limits_{n = 0}^{\infty}A_{n}\bigg(\frac{x}{2}\bigg)\bigg]\bigg], \end{equation}$

(4.11)

where $A_{n}$ are the Adomian polynomials corresponding to the nonlinear terms $y^{2}(\frac{x}{2})$ defined as follows:

$\begin{equation} \begin{split} A_0& = y_{0}\! \left(\frac{x}{2}\right)^{2},\\ A_1& = 2 y_{0}\! \left(\frac{x}{2}\right) y_{1}\! \left(\frac{x}{2}\right),\\ A_2& = y_{1}\! \left(\frac{x}{2}\right)^{2}+2 y_{0}\! \left(\frac{x}{2}\right) y_{2}\! \left(\frac{x}{2}\right),\\ & \vdots \end{split} \end{equation}$

(4.12)

Moreover, the formal recursive relationship for the governing fractional IVP is obtained as follows:

$\begin{equation} \begin{cases} y_{0}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}[1]\bigg],\\ y_{n}(x) = -\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[2A_{n-1}\bigg(\frac{x}{2}\bigg)\bigg]\bigg], \ n\geq 1. \end{cases} \end{equation}$

(4.13)

Accordingly, the recursive scheme is expressed more plainly as follows:

$\begin{equation} \begin{split} y_{0}(x)& = \frac{1}{\Gamma(\alpha+1)}x^{\alpha},\\ y_{1}(x)& = -\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[2A_{0}\bigg(\frac{x}{2}\bigg)\bigg]\bigg] \\ & = -\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[2\frac{x^{2\alpha}}{4^{\alpha}\Gamma^{2}(\alpha+1)}\bigg]\bigg] \\ & = -\frac{2\Gamma(2\alpha+1)}{4^{\alpha}\Gamma^{2}(\alpha+1)\Gamma(3\alpha+1)}x^{3\alpha},\\ y_{2}(x)& = \frac{2^{3}\Gamma(2\alpha+1)\Gamma(4\alpha+1)}{4^{3\alpha}\Gamma^{3}(\alpha+1)\Gamma(3\alpha+1)\Gamma(5\alpha+1)}x^{5\alpha},\\ & \vdots \end{split} \end{equation}$

(4.14)

Hence, on taking the net sum of the latter solution components, one acquires

$\begin{equation} y(x) = \frac{x^{\alpha}}{\Gamma(\alpha+1)}-\frac{2\Gamma(2\alpha+1)}{4^{\alpha}\Gamma^{2}(\alpha+1)\Gamma(3\alpha+1)}x^{3\alpha}+\frac{2^{3}\Gamma(2\alpha+1)\Gamma(4\alpha+1)}{4^{3\alpha}\Gamma^{3}(\alpha+1)\Gamma(3\alpha+1)\Gamma(5\alpha+1)}x^{5\alpha}+\dots \end{equation}$

(4.15)

Remarkably, when the fractional-order $\alpha$ takes unity in the above series solution, that is, $\alpha = 1$ , one obtains the following compacted closed-form solution:

$\begin{equation} y(x) = x-\frac{1}{3!}x^{3}+\dots = \sum\limits_{n = 0}^{\infty}\frac{(-1)^{n}x^{2n+1}}{(2n+1)!} = \sin(x). \end{equation}$

(4.16)

Similarly, the authors in ^{[42,43,44,45]} addressed the governing model in (4.8) using the Chebyshev wavelet method (CWM), the Jacobi spectral configuration method (JSCM), the rescaled spectral configuration method (RSCM), and the generalized Legendre polynomial configuration method (GLPCM), respectively. These authors obtained approximate solutions by calculating nine terms, while our algorithm provides an exact solution for the model at $\alpha = 1$ . For error analysis, refer to Table 3, which compares the proposed SDM with other competing approaches. Additionally, Figure 1 presents a graphical depiction of the approximate SDM solution and the absolute errors for the governing nonlinear initial-value problem (IVP) for FPDDE (4.8).

Table 3. Comparison of absolute errors at various

$x_i$ values for the fractional IVP (4.8).

$x$	JSCM ^[43]	CWM ^[44]	RSCM ^[45]	GLPCM ^[42]	Present method (SDM)
0.125	$2.05 \times 10^{-12}$	$1.71 \times 10^{-12}$	$1.81 \times 10^{-13}$	$1.61 \times 10^{-13}$	0
0.250	$1.61 \times 10^{-12}$	$2.47 \times 10^{-12}$	$1.34 \times 10^{-13}$	$3.14 \times 10^{-14}$	0
0.375	$3.48 \times 10^{-11}$	$9.36 \times 10^{-11}$	$1.50 \times 10^{-13}$	$5.29 \times 10^{-14}$	0
0.500	$5.31 \times 10^{-11}$	$1.78 \times 10^{-11}$	$1.61 \times 10^{-13}$	$2.38 \times 10^{-14}$	0
0.625	$2.07 \times 10^{-11}$	$1.88 \times 10^{-11}$	$1.69 \times 10^{-13}$	$5.35 \times 10^{-14}$	0
0.750	$2.47 \times 10^{-11}$	$3.04 \times 10^{-12}$	$1.30 \times 10^{-13}$	$6.05 \times 10^{-14}$	0
0.875	$1.59 \times 10^{-11}$	$1.64 \times 10^{-12}$	$1.38 \times 10^{-13}$	$1.61 \times 10^{-14}$	0
1	$3.85 \times 10^{-10}$	$4.51 \times 10^{-10}$	$8.98 \times 10^{-13}$	$3.11 \times 10^{-15}$	0

| Show Table

DownLoad: CSV

Figure 1. SDM solutions and error curves for the fractional IVP (4.8).

DownLoad: Full-Size Img PowerPoint

Example 4.3. Consider the nonhomogeneous nonlinear IVP for FPDDE as follows ^[42,43,46]:

$\begin{equation} \begin{cases} D^{1/2}_{*}y(x) = y(\frac{x}{3})+3y^{2}(x)+g(x),\\ y(0) = 0, \end{cases} \end{equation}$

(4.17)

where

$g(x) = \frac{8}{3\sqrt{\pi}}x^{{\frac{3}{2}}}-\frac{2}{\sqrt{\pi}}x^{{\frac{1}{2}}}-x^{4}+2x^{3}-\frac{10}{9}x^{2}+\frac{1}{3}x.$

Thus, without much delay, the formal recursive relationship for the governing fractional model is obtained as follows:

$\begin{equation} \begin{cases} y_{0}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[g(x)\bigg]\bigg] = x^{2}-x-\frac{256}{315\sqrt{\pi}}x^{\frac{9}{2}}+\frac{64}{35\sqrt{\pi}}x^{\frac{7}{2}}-\frac{32}{27\sqrt{\pi}}x^{\frac{5}{2}}+\frac{4}{9\sqrt{\pi}}x^{\frac{3}{2}} = p(x),\\ y_{n}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[y_{n-1}(\frac{1}{3}x)\bigg]\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[3A_{n-1}\bigg]\bigg], \ n\geq 1, \end{cases} \end{equation}$

(4.18)

where the Adomian polynomials $A_{n}$ in the above scheme for the nonlinear term $y^{2}(x)$ , take the following forms:

$\begin{equation} \begin{split} A_0& = y_{0}\! \left(x\right)^{2},\\ A_1& = 2 y_{0}\! \left(x\right) y_{1}\! \left(x\right),\\ A_2& = y_{1}\! \left(x\right)^{2}+2 y_{0}\! \left(x\right) y_{2}\! \left(x\right),\\ & \vdots \end{split} \end{equation}$

(4.19)

Accordingly, when considering the Wazwaz's modification for ADM (3.12) by choosing the functions $p_0(x)$ and $p_1(x)$ as follows:

$\begin{equation} \begin{split} p_0(x)& = x^{2}-x,\\ p_{1}(x)& = -\frac{256}{315\sqrt{\pi}}x^{\frac{9}{2}}+\frac{64}{35\sqrt{\pi}}x^{\frac{7}{2}}-\frac{32}{27\sqrt{\pi}}x^{\frac{5}{2}}+\frac{4}{9\sqrt{\pi}}x^{\frac{3}{2}}, \end{split} \end{equation}$

(4.20)

the modified SDM algorithm then yields the following scheme:

$\begin{equation} \begin{cases} y_{0}(x) = p_{0}(x) = x^{2}-x,\\ y_{1}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[y_{0}(x)\bigg]\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[{3A_{0}(x)}\bigg]\bigg]+p_{1}(x) \\ \qquad = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[x^{2}-x\bigg]\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-1/2}}\mathbb{S}\bigg[{3(x^{2}-x)^{2}}\bigg]\bigg] -\frac{256}{315\sqrt{\pi}}x^{\frac{9}{2}}+\frac{64}{35\sqrt{\pi}}x^{\frac{7}{2}}-\frac{32}{27\sqrt{\pi}}x^{\frac{5}{2}}+\frac{4}{9\sqrt{\pi}}x^{\frac{3}{2}} \\ \qquad = \frac{256}{315\sqrt{\pi}}x^{\frac{9}{2}}-\frac{64}{35\sqrt{\pi}}x^{\frac{7}{2}}+\frac{32}{27\sqrt{\pi}}x^{\frac{5}{2}}-\frac{4}{9\sqrt{\pi}}x^{\frac{3}{2}}-\frac{256}{315\sqrt{\pi}}x^{\frac{9}{2}}+\frac{64}{35\sqrt{\pi}}x^{\frac{7}{2}}-\frac{32}{27\sqrt{\pi}}x^{\frac{5}{2}}+\frac{4}{9\sqrt{\pi}}x^{\frac{3}{2}} \\ \qquad = 0,\\ y_{n}(x) = 0, \ n\ge 2, \end{cases} \end{equation}$

(4.21)

which then sums to yield the resulting exact solution as follows:

$\begin{equation} y(x) = x^{2}-x. \end{equation}$

(4.22)

Consequently, the references ^[42,43,46] also addressed the governing model in (4.17) using the Jacobi spectral configuration method (JSCM), the Jacobi spectral Galerkin method (JSGM), and the generalized Legendre polynomial spectral method (GLPSM), respectively. These authors obtained approximate solutions by calculating varying numbers of terms, while our algorithm provides an exact solution using only two terms. This advantage leads to more effective convergence to the exact solution compared to numerical methods, significantly reducing computational complexity. For a comparative analysis of the aforementioned methods with respect to different error norms, see Table 4.

Table 4. Comparison of maximum errors for different methods for the fractional IVP (4.17).

Maximum Error	$n = 6$	$n = 8$	$n = 10$
JSCM ^[43] ( $L^2$ )	$3.6822 \times 10^{-6}$	$1.20720 \times 10^{-8}$	$1.3049 \times 10^{-11}$
JSCM ^[43] ( $L^\infty$ )	$1.0832 \times 10^{-5}$	$2.7460 \times 10^{-8}$	$3.9471 \times 10^{-11}$
JSGM ^[46] ( $L^2$ )	$2.3632 \times 10^{-6}$	$1.0032 \times 10^{-8}$	$1.1006\times 10^{-11}$
JSGM ^[46] ( $L^\infty$ )	$1.2851 \times 10^{-5}$	$1.4350 \times 10^{-8}$	$2.9356 \times 10^{-11}$
GLPCM ^[42]	$4.5103 \times 10^{-16}$	$4.7184 \times 10^{-16}$	$5.3429 \times 10^{-16}$

| Show Table

DownLoad: CSV

Example 4.4. Consider the IVP for FPDDE of the following form:

$\begin{equation} \begin{cases} D^{\alpha}_{*}y(x) = 1+2y(\frac{x}{2})-y(x),\\ y(0) = 0, \end{cases} \end{equation}$

(4.23)

where the fractional-order $\alpha$ assumes a general value over the interval $0 < \alpha \le 1.$

Subsequently, upon applying the ST on the fractional IVP in via (2.11), one obtains

$\begin{equation} \mathbb{S}[y(x)] = \frac{1}{u^{-\alpha}}+\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[2y(\frac{x}{2})\bigg]-\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[y(x)\bigg]. \end{equation}$

(4.24)

Now, on applying the inverse ST on the latter equation, one obtains

$\begin{equation} y(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[2y(\frac{x}{2})\bigg]\bigg]-\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[y(x)\bigg]\bigg], \end{equation}$

(4.25)

such that when the ADM is deployed, one obtains

$\begin{equation} \sum\limits_{n = 0}^{\infty}y_{n}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[2\sum\limits_{n = 0}^{\infty}y_{n}(\frac{x}{2})\bigg]\bigg]-\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[\sum\limits_{n = 0}^{\infty}y_{n}(x)\bigg]\bigg]. \end{equation}$

(4.26)

Accordingly, the formal recursive relation for the governing fractional IVP is obtained as follows:

$\begin{equation} \begin{cases} y_{0}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\bigg],\\ y_{n}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[2y_{n-1}(\frac{x}{2})\bigg]\bigg]-\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[y_{n-1}(x)\bigg]\bigg], \ n\geq 1, \end{cases} \end{equation}$

(4.27)

or more openly from the above recursive relation as follows:

$\begin{equation} \begin{split} y_{0}(x)& = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\bigg] = \frac{x^{\alpha}}{\Gamma(\alpha+1)},\\ y_{1}(x)& = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[2y_{0}(\frac{x}{2})\bigg]\bigg]-\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[y_{0}(x)\bigg]\bigg] = \frac{x^{2\alpha}}{\Gamma(2\alpha+1)}\bigg(\frac{1}{2^{\alpha-1}}-1\bigg),\\ y_{2}(x)& = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[2y_{1}(\frac{x}{2})\bigg]\bigg]-\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[y_{1}(x)\bigg]\bigg] \\ & = \frac{x^{3\alpha}}{\Gamma(3\alpha+1)}\bigg(\frac{1}{2^{\alpha-1}}-1\bigg)\bigg(\frac{1}{2^{\alpha-1}}-1\bigg),\\ y_{3}(x)& = \frac{x^{4\alpha}}{\Gamma(4\alpha+1)}\bigg(\frac{1}{2^{\alpha-1}}-1\bigg)\bigg(\frac{1}{2^{\alpha-1}}-1\bigg)\bigg(\frac{1}{2^{\alpha-1}}-1\bigg),\\ & \vdots\\ \end{split} \end{equation}$

(4.28)

Hence, on taking the net sum of the above solution components, one obtains

$\begin{equation} y(x) = \frac{x^{\alpha}}{\Gamma(\alpha+1)}+\sum\limits_{i = 1}^{\infty}\bigg(\prod\limits_{l = 1}^{i} \frac{1}{2^{l\alpha-1}}-1\bigg)\frac{x^{(i+1)\alpha}}{\Gamma((i+1)\alpha+1)}. \end{equation}$

(4.29)

Additionally, when the fractional-order $\alpha$ assumes full integer-order, that is, when $\alpha = 1$ , the second term in the above series solution vanishes, which then leads to the acquisition of the exact solution of the integer-order model as follows:

$\begin{equation} y(x) = x. \end{equation}$

(4.30)

Moreover, see for the pictorial representation of the approximate SDM solution, together with the resulting absolute error for various values of the fractional-order $\alpha$ for the governing IVP in (4.23).

Figure 2. SDM solution and absolute error of the fractional IVP (4.23) for various values of

$\alpha$ .

DownLoad: Full-Size Img PowerPoint

Example 4.5. Consider the inhomogeneous IVP for FPDDE ^[18,47]

$\begin{equation} \begin{cases} D^{5/2}_{*}y(x) = -y(x)-y(x-0.5)+g(x),\; \; x \in [0,1],\\ y(0) = y'(0) = y''(0) = 0, \end{cases} \end{equation}$

(4.31)

where

$g(x) = \frac{\Gamma(4)}{\Gamma(3/2)}x^{1/2}+x^{3}+(x-0.5)^{3}.$

Accordingly, upon applying the ST and its subsequent inverse transform on (4.31) as previously highlighted, one thus obtains

$\begin{equation} y(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-5/2}}\mathbb{S}[-y(x)]\bigg]-\mathbb{S}^{-1}\bigg[\frac{1}{u^{-5/2}}\mathbb{S}\big[y(x-0.5)\big]\bigg]+\mathbb{S}^{-1}\bigg[\frac{1}{u^{-5/2}}\mathbb{S}\big[g(x)\big]\bigg], \end{equation}$

(4.32)

such that the standard ADM reveals the resulting recurrent scheme as follows:

$\begin{equation} \begin{cases} y_{0}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-5/2}}\mathbb{S}[g(x)]\bigg] = p(x),\\ y_{n}(x) = -\mathbb{S}^{-1}\bigg[\frac{1}{u^{-5/2}}\mathbb{S}\bigg[y_{n-1}(x)\bigg]\bigg]-\mathbb{S}^{-1}\bigg[\frac{1}{u^{-5/2}}\mathbb{S}\bigg[y_{n-1}(x-0.5)\bigg]\bigg], \ n\geq 1. \end{cases} \end{equation}$

(4.33)

In the same fashion, when considering the Wazwaz's modification for ADM (3.12) by choosing the functions $p_0(x)$ and $p_1(x)$ as follows:

$\begin{equation} \begin{split} p_{0}(x)& = x^3,\\ p_{1}(x)& = \frac{\Gamma(4)}{\Gamma(13/2)}x^{11/2}+\mathbb{S}^{-1}\bigg[u^{5/2}(u-0.5)^{3}\Gamma(4)\bigg], \end{split} \end{equation}$

(4.34)

the modified SDM scheme is then expressed as follows:

$\begin{equation} \begin{split} y_{0}(x)& = p_{0}(x) = x^3,\\ y_{1}(x)& = -\mathbb{S}^{-1}\bigg[\frac{1}{u^{-5/2}}\mathbb{S}\bigg[y_{0}(x)\bigg]\bigg]-\mathbb{S}^{-1}\bigg[\frac{1}{u^{-5/2}}\mathbb{S}\bigg[y_{0}(x-0.5)\bigg]\bigg]+p_{1}(x),\\ & = -\mathbb{S}^{-1}\bigg[\frac{1}{u^{-5/2}}\mathbb{S}\bigg[x^{3}\bigg]\bigg]-\mathbb{S}^{-1}\bigg[\frac{1}{u^{-5/2}}\mathbb{S}\bigg[(x-0.5)^{3}\bigg]\bigg]+\frac{\Gamma(4)}{\Gamma(13/2)}x^{11/2}\\ &\quad+\mathbb{S}^{-1}\bigg[u^{5/2}(u-0.5)^{3}\Gamma(4)\bigg],\\ & = -\frac{\Gamma(4)}{\Gamma(13/2)}x^{11/2}-\mathbb{S}^{-1}\bigg[u^{5/2}(u-0.5)^{3}\Gamma(4)\bigg]+\frac{\Gamma(4)}{\Gamma(13/2)}x^{11/2}\\ &\quad+\mathbb{S}^{-1}\bigg[u^{5/2}(u-0.5)^{3}\Gamma(4)\bigg],\\ & = 0, \\ y_{n}(x)& = 0, \ n\geq 2, \end{split} \end{equation}$

(4.35)

such that the closed-form solution of the governing IVP is attained by summing the above solution components as follows:

$\begin{equation} y(x) = x^{3}. \end{equation}$

(4.36)

Moreover, upon simulating the derived SDM scheme for the examining model, together with some contending computational methods, then provides the absolute error contrast between the deployed SDM (with only a few terms, $N = 2$ ), the results reported in ^[47], where $N = 22$ using different Laguerre's parameters $\beta$ through the application of the Laguerre-Gauss collocation scheme (L-GCS), and the approximate solution in ^[18] using the operational matrix via the Genocchi polynomials collocation method (GPCM), where $N = 4$ . Certainly, the results show that the proposed SDM method provides an almost exact solution over the contending approaches.

Table 5. Comparison of the absolute errors at various

$x_i$ values for the fractional IVP (4.31).

$x$	L-GCS ^[47] $N = 22$			GPCM ^[18] $N = 4$	Present method (SDM) $N = 2$
$x$	$\beta = 2$	$\beta = 3$	$\beta = 5$	GPCM ^[18] $N = 4$	Present method (SDM) $N = 2$
0.1	$1.030 \times 10^{-4}$	$1.019 \times 10^{-5}$	$6.273 \times 10^{-6}$	$6.17040 \times 10^{-9}$	0
0.2	$6.510 \times 10^{-4}$	$6.051 \times 10^{-5}$	$3.892 \times 10^{-5}$	$4.93630 \times 10^{-8}$	0
0.3	$1.740 \times 10^{-3}$	$1.495 \times 10^{-4}$	$1.023 \times 10^{-4}$	$1.66600 \times 10^{-7}$	0
0.4	$3.283 \times 10^{-3}$	$2.559 \times 10^{-4}$	$1.901 \times 10^{-4}$	$3.94910 \times 10^{-7}$	0
0.5	$5.138 \times 10^{-3}$	$3.546 \times 10^{-4}$	$2.944 \times 10^{-4}$	$7.71300 \times 10^{-7}$	0
0.6	$7.175 \times 10^{-3}$	$4.261 \times 10^{-4}$	$4.088 \times 10^{-4}$	$1.33280 \times 10^{-6}$	0
0.7	$9.303 \times 10^{-3}$	$4.592 \times 10^{-4}$	$5.306 \times 10^{-4}$	$2.11640 \times 10^{-6}$	0
0.8	$1.147 \times 10^{-2}$	$4.510 \times 10^{-4}$	$6.597 \times 10^{-4}$	$3.15920 \times 10^{-6}$	0
0.9	$1.367 \times 10^{-2}$	$4.055 \times 10^{-4}$	$7.977 \times 10^{-4}$	$4.49820 \times 10^{-6}$	0
1.0	$1.589 \times 10^{-2}$	$3.311 \times 10^{-4}$	$9.468 \times 10^{-4}$	$6.17040 \times 10^{-6}$	0

| Show Table

DownLoad: CSV

Example 4.6. Consider the coupled system of nonhomogeneous IVP for FPDDE as follows ^[48]:

$\begin{equation} \begin{cases} D^{\alpha}_{*}y_{1}(x) = e^{-x^{\alpha}}-e^{\frac{x^{\alpha}}{2}}+y_{1}(\frac{x}{2})+y_{1}(x)-y_{2}(x),\; \; y_{1}(0) = 1,\\ D^{\alpha}_{*}y_{2}(x) = e^{x^{\alpha}}+e^{\frac{x^{\alpha}}{2}}-y_{2}(\frac{x}{2})-y_{1}(x)-y_{2}(x), \ \ \ \ y_{2}(0) = 1. \end{cases} \end{equation}$

(4.37)

Accordingly, the deployed SDM posses the recursive scheme for the governing coupled fractional model as follows:

$\begin{equation} \begin{cases} y_{1,0}(x) = 1+ \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}[e^{-x^{\alpha}}-e^{\frac{x^{\alpha}}{2}}]\bigg],\\ y_{2,0}(x) = 1+ \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}[e^{x^{\alpha}}+e^{\frac{x^{\alpha}}{2}}]\bigg], \end{cases} \end{equation}$

(4.38)

and $\forall$ $n\geq 1$ as follows:

$\begin{equation} \begin{cases} y_{1,n}(x) = \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}[y_{1,n-1}(\frac{x}{2})+y_{1,n-1}(x)-y_{2,n-1}(x)]\bigg],\\ y_{2,n}(x) = -\mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}[y_{2,n-1}(\frac{x}{2})+y_{1,n-1}(x)+y_{2,n-1}(x)]\bigg]. \end{cases} \end{equation}$

(4.39)

Certainly, the current coupled system of fractional IVP extends the corresponding integer-order model examined in ^[49]. Thus, upon considering the integer-order version in ^[49], where the fractional-order takes unity, that is, when $\alpha = 1$ , the model admits the following exact solution $y_{1}(x) = e^{x},$ $y_{2}(x) = e^{-x}$ . What is more, upon simplifying the derived scheme further, one replaces the exponential terms in (4.38) with the first seven terms of their Taylor series expansions to eventually obtain as follows:

$\begin{equation} \begin{cases} y_{1,0}(x) = 1+ \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[\frac{3x^{\alpha}}{2}-\frac{3x^{2\alpha}}{8}+\frac{3x^{3\alpha}}{16}-\frac{5x^{4\alpha}}{128}+\frac{11x^{5\alpha}}{1280}-\frac{7x^{6\alpha}}{5120}\bigg]\bigg],\\ y_{2,0}(x) = 1- \mathbb{S}^{-1}\bigg[\frac{1}{u^{-\alpha}}\mathbb{S}\bigg[2+\frac{x^{\alpha}}{2}+\frac{5x^{2\alpha}}{8}+\frac{7x^{3\alpha}}{48}+\frac{17x^{4\alpha}}{384}+\frac{31x^{5\alpha}}{3840}+\frac{13x^{6\alpha}}{9216}\bigg]\bigg], \end{cases} \end{equation}$

(4.40)

that is,

$\begin{align*} y_{1,0}(x) & = 1 + \frac{3x^{2\alpha} \Gamma(1 + \alpha)}{2\Gamma(1 + 2\alpha)} - \frac{3x^{3\alpha} \Gamma(1 + \alpha)}{8\Gamma(1 + 3\alpha)}+ \frac{3x^{4\alpha} \Gamma(1 + 3\alpha)}{16\Gamma(1 + 4\alpha)} - \frac{5x^{4\alpha} \Gamma(1 + \alpha)}{128\Gamma(1 + 128\alpha)} \\&\quad+ \frac{11x^{6\alpha} \Gamma(1 + 5\alpha)}{1280\Gamma(1 + 6\alpha)} - \frac{7x^{7\alpha} \Gamma(1 + 6\alpha)}{5120\Gamma(1 + 7\alpha)}, \\ y_{2,0}(x) & = 1 -\frac{2x^\alpha}{\Gamma(1 + \alpha)} - \frac{x^{2\alpha} \Gamma(1 + \alpha)}{2\Gamma(1 + 2\alpha)} - \frac{5x^{3\alpha} \Gamma(1 + 2\alpha)}{8\Gamma(1 + 3\alpha)} - \frac{7x^{4\alpha} \Gamma(1 + 3\alpha)}{48\Gamma(1 + 4\alpha)} \\&\quad- \frac{17x^{5\alpha} \Gamma(1 + 4\alpha)}{384\Gamma(1 + 5\alpha)} - \frac{31x^{6\alpha} \Gamma(1 + 5\alpha)}{3840\Gamma(1 + 6\alpha)} - \frac{13x^{7\alpha} \Gamma(1 + 6\alpha)}{9216\Gamma(1 + 7\alpha)},\\ y_{1,1}(x)& = \frac{x^\alpha}{\Gamma(1 + \alpha)} - \frac{2x^{2\alpha}}{\Gamma(1 + 2\alpha)} + \frac{2^\alpha x^{2\alpha}}{\Gamma(1 + 2\alpha)} + \dots - \frac{x^{8\alpha} \Gamma(1 + 6\alpha)}{23040\Gamma(1 + 8\alpha)} + \frac{x^{8\alpha} \Gamma(1 + 6\alpha)}{72^{7\alpha+10} \times 5\Gamma(1 + 8\alpha)}, \\ y_{2,1}(x)& = \frac{x^\alpha}{\Gamma(1 + \alpha)} - \frac{2^\alpha x^{2\alpha}}{\Gamma(1 + 2\alpha)} - \frac{x^{2\alpha} \Gamma(1 + \alpha)}{2\Gamma(1 + 2\alpha)} - \dots - \frac{x^{8\alpha} \Gamma(1 + 6\alpha)}{360\Gamma(1 + 8\alpha)} - \frac{x^{8\alpha} \Gamma(1 + 6\alpha)}{132^{7\alpha+10} \times 9\Gamma(1 + 8\alpha)},\\ \vdots \end{align*}$

In addition, depicts the pictorial view of the obtained SDM solution for the fractional-order system (4.37). Furthermore, with the consideration of the corresponding integer-order case, that is, when $\alpha = 1,$ the first 6 terms of the solution components by SDM sum to yield the following series solution:

$\begin{equation} \begin{cases} y_{1}(x) = 1+x+\frac{1}{2} x^{2}+\frac{1}{6} x^{3}+\frac{1}{24} x^{4}+\frac{1}{120} x^{5}+\dots\approx e^{x},\\ y_{2}(x) = 1-x+\frac{1}{2} x^{2}-\frac{1}{6} x^{3}+\frac{1}{24} x^{4}-\frac{1}{120} x^{5}+\dots\approx e^{-x}. \end{cases} \end{equation}$

(4.41)

Figure 3. Exact and SDM solutions of the coupled fractional IVP (4.37) for various values

$\alpha$ .

DownLoad: Full-Size Img PowerPoint

5. Conclusions

The current manuscript introduces a novel computational approach for solving the class of FPDDEs by combining the Sumudu transform with the reliable ADM. This approach enables the derivation of approximate or exact solutions through straightforward iterative calculations that are easily implementable in computer programs. The numerical results confirm the precision of the proposed method. Furthermore, the study successfully achieves its goal by obtaining numerical solutions in the form of rapidly converging series with easily calculable components. To validate the accuracy and efficiency of the algorithms adopted via the SDM, we present several important and interesting problems, comparing the obtained solutions both numerically and graphically. The results demonstrate that our method is highly effective, straightforward, and yields approximate solutions that closely match the exact solutions. Additionally, the paper recommends applying this technique to more complex neutral pantograph equations involving fractional derivatives on the right-hand side, utilizing various fractional operators.

Author contributions

M. Al-Mazmumy: Conceptualization; A. S. Alsulami: Formal analysis, Writing–original draft; A. S. Alsulami, M. Al-Mazmumy, M. A. Alyami and M. Alsulami: Methodology, Investigation, Writing–review & editing; A. S. Alsulami and M. Al-Mazmumy: Software. All authors have read and agreed to the published version of the manuscript.

Use of Generative-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.

References

[1]	K. S. Miller, B. Ross, An introduction to the fractional calculus and fractional differential equations, New York: Wiley, 1993.
[2]	I. Podlubny, Fractional differential equations, Academic Press, 1999.
[3]	J. T. Machado, V. Kiryakova, F. Mainardi, Recent history of fractional calculus, Commun. Nonlinear Sci. Numer. Simul., 16 (2011), 1140–1153. https://doi.org/10.1016/j.cnsns.2010.05.027 doi: 10.1016/j.cnsns.2010.05.027
[4]	R. T. Baillie, Long memory processes and fractional integration in econometrics, J. Econom., 73 (1996), 5–59. https://doi.org/10.1016/0304-4076(95)01732-1 doi: 10.1016/0304-4076(95)01732-1
[5]	G. C. Wu, Z. G. Deng, D. Baleanu, D. Q. Zeng, New variable-order fractional chaotic systems for fast image encryption, Chaos, 29 (2019), 083103. https://doi.org/10.1063/1.5096645 doi: 10.1063/1.5096645
[6]	A. Bonfanti, J. L. Kaplan, G. Charras, A. Kabla, Fractional viscoelastic models for power-law materials, Soft Matter, 16 (2020), 6002–6020. https://doi.org/10.1039/D0SM00354A doi: 10.1039/D0SM00354A
[7]	J. F. Gomez-Aguilar, J. E. Escalante-Martinez, C. Calderon-Ramon, L. J. Morales-Mendoza, M. Benavidez-Cruz, M. Gonzalez-Lee, Equivalent circuits applied in electrochemical impedance spectroscopy and fractional derivatives with and without singular kernel, Adv. Math. Phys., 2016 (2016), 9720181. https://doi.org/10.1155/2016/9720181 doi: 10.1155/2016/9720181
[8]	N. Singh, K. Kumar, P. Goswami, H. Jafari, Analytical method to solve the local fractional vehicular traffic flow model, Math. Methods Appl. Sci., 45 (2022), 3983–4001. https://doi.org/10.1002/mma.8027 doi: 10.1002/mma.8027
[9]	H. Bulut, T. A. Sulaiman, H. M. Baskonus, H. Rezazadeh, M. Eslami, M. Mirzazadeh, Optical solitons and other solutions to the conformable space-time fractional Fokas-Lenells equation, Optik, 172 (2018), 20–27. https://doi.org/10.1016/j.ijleo.2018.06.108 doi: 10.1016/j.ijleo.2018.06.108
[10]	A. A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and applications of fractional differential equations, Elsevier, 2006.
[11]	I. T. Huseynov, N. I. Mahmudov, A class of Langevin time-delay differential equations with general fractional orders and their applications to vibration theory, J. King Saud Univ. Sci., 33 (2021), 101596. https://doi.org/10.1016/j.jksus.2021.101596 doi: 10.1016/j.jksus.2021.101596
[12]	S. Mashayekhi, S. Sedaghat, Fractional model of stem cell population dynamics, Chaos Solitons Fract., 146 (2021), 110919. https://doi.org/10.1016/j.chaos.2021.110919 doi: 10.1016/j.chaos.2021.110919
[13]	I. T. Huseynov, N. I. Mahmudov, Delayed analogue of three-parameter Mittag-Leffler functions and their applications to Caputo-type fractional time delay differential equations, Math. Methods Appl. Sci., 47 (2024), 11019–11043. https://doi.org/10.1002/mma.6761 doi: 10.1002/mma.6761
[14]	I. T. Huseynov, N. I. Mahmudov, Analysis of positive fractional-order neutral time-delay systems, J. Franklin Inst., 359 (2022), 294–330. https://doi.org/10.1016/j.jfranklin.2021.07.001 doi: 10.1016/j.jfranklin.2021.07.001
[15]	M. Bohner, O. Tunc, C. Tunc, Qualitative analysis of Caputo fractional integro-differential equations with constant delays, Comput. Appl. Math., 40 (2021), 214. https://doi.org/10.1007/s40314-021-01595-3 doi: 10.1007/s40314-021-01595-3
[16]	J. R. Ockendon, A. B. Tayler, The dynamics of a current collection system for an electric locomotive, Proc. R. Soc. Lond. A, 322 (1971), 447–468. https://doi.org/10.1098/rspa.1971.0078 doi: 10.1098/rspa.1971.0078
[17]	S. S. Ezz-Eldien, On solving systems of multi-pantograph equations via spectral tau method, Appl. Math. Comput., 321 (2018), 63–73. https://doi.org/10.1016/j.amc.2017.10.014 doi: 10.1016/j.amc.2017.10.014
[18]	A. Isah, C. Phang, P. Phang, Collocation method based on Genocchi operational matrix for solving generalized fractional pantograph equations, Int. J. Differ. Equ., 2017 (2017), 2097317. https://doi.org/10.1155/2017/2097317 doi: 10.1155/2017/2097317
[19]	L. Shi, X. H. Ding, Z. Chen, Q. Ma, A new class of operational matrices method for solving fractional neutral pantograph differential equations, Adv. Differ. Equ., 2018 (2018), 1–17. https://doi.org/10.1186/s13662-018-1536-8 doi: 10.1186/s13662-018-1536-8
[20]	M. T. Kajani, Numerical solution of fractional pantograph equations via Müntz-Legendre polynomials, Math. Sci., 18 (2024), 387–395. https://doi.org/10.1007/s40096-022-00507-8 doi: 10.1007/s40096-022-00507-8
[21]	G. Adomian, Nonlinear stochastic system theory and applications to physics, Dordrecht: Kluwer Academic Publishers, 1989.
[22]	A. Abdelrazec, D. Pelinovsky, Convergence of the Adomian decomposition method for initial-value problems, Numer. Methods Partial Differ. Equ., 27 (2011), 749–766. https://doi.org/10.1002/num.20549 doi: 10.1002/num.20549
[23]	V. Daftardar-Gejji, H. Jafari, Adomian decomposition: a tool for solving a system of fractional differential equations, J. Math. Anal. Appl., 301 (2005), 508–518. https://doi.org/10.1016/j.jmaa.2004.07.039 doi: 10.1016/j.jmaa.2004.07.039
[24]	S. O. Adesanya, Linear stability analysis of a plane-Poiseuille hydromagnetic flow using Adomian decomposition method, U.P.B. Sci. Bull. Series A, 75 (2013), 99–106.
[25]	H. Thabet, S. Kendre, New modification of Adomian decomposition method for solving a system of nonlinear fractional partial differential equations, Int. J. Adv. Appl. Math. Mech., 6 (2019), 1–13.
[26]	M. D. Johansyah, A. K. Supriatna, E. Rusyaman, J. Saputra, Bernoulli fractional differential equation solution using Adomian decomposition method, IOP Conf. Ser. Mater. Sci. Eng., 1115 (2021), 012015. https://doi.org/10.1088/1757-899X/1115/1/012015 doi: 10.1088/1757-899X/1115/1/012015
[27]	R. I. Nuruddeen, J. F. Gómez-Aguilar, J. R. Razo-Hernández, Fractionalizing, coupling and methods for the coupled system of two-dimensional heat diffusion models, AIMS Math., 8 (2023), 11180–11201. https://doi.org/10.3934/math.2023566 doi: 10.3934/math.2023566
[28]	S. Masood, Hajira, H. Khan, R. Shah, S. Mustafa, Q. Khan, et al., A new modified technique of Adomian decomposition method for fractional diffusion equations with initial-boundary conditions, J. Funct. Spaces, 2022 (2022), 6890517. https://doi.org/10.1155/2022/6890517 doi: 10.1155/2022/6890517
[29]	M. Al-Mazmumy, M. A. Alyami, M. Alsulami, A. S. Alsulami, Efficient modified Adomian decomposition method for solving nonlinear fractional differential equations, Int. J. Anal. Appl., 22 (2024), 76. https://doi.org/10.28924/2291-8639-22-2024-76 doi: 10.28924/2291-8639-22-2024-76
[30]	Y. F. Zhang, L. Z. Wang, Application of Laplace Adomian decomposition method for fractional Fokker-Planck equation and time fractional coupled Boussinesq-Burger equations, Eng. Comput., 2024. https://doi.org/10.1108/ec-06-2023-0275
[31]	A. Khalouta, A. Kadem, Fractional natural decomposition method for solving a certain class of nonlinear time-fractional wave-like equations with variable coefficients, Acta Univ. Sapientiae Math., 11 (2019), 99–116. https://doi.org/10.2478/ausm-2019-0009 doi: 10.2478/ausm-2019-0009
[32]	O. E. Ige, R. A. Oderinu, T. M. Elzaki, Numerical simulation of the nonlinear coupled Jaulent-Miodek equation by Elzaki transform-Adomian polynomial method, Adv. Math. Sci. J., 9 (2020), 10335–10355. https://doi.org/10.37418/amsj.9.12.25 doi: 10.37418/amsj.9.12.25
[33]	G. K. Watugala, Sumudu transform: a new integral transform to solve differential equations and control engineering problems, Int. J. Math. Ed. Sci. Technol., 24 (1993), 35–43. https://doi.org/10.1080/0020739930240105 doi: 10.1080/0020739930240105
[34]	F. B. M. Belgacem, A. A. Karaballi, Sumudu transform fundamental properties investigations and applications, Int. J. Stoch. Anal., 2006 (2006), 091083. https://doi.org/10.1155/JAMSA/2006/91083 doi: 10.1155/JAMSA/2006/91083
[35]	M. A. Asiru, Further properties of the Sumudu transform and its applications, Int. J. Math. Ed. Sci. Technol., 33 (2002), 441–449. https://doi.org/10.1080/002073902760047940 doi: 10.1080/002073902760047940
[36]	S. A. Ahmed, M. Elbadri, Solution of Newell-Whitehead-Segal equation of fractional order by using Sumudu decomposition method, Math. Statist., 8 (2020), 631–636. https://doi.org/10.13189/ms.2020.080602 doi: 10.13189/ms.2020.080602
[37]	T. Patel, R. Meher, Adomian decomposition Sumudu transform method for solving a solid and porous fin with temperature-dependent internal heat generation, SpringerPlus, 5 (2016), 1–18. https://doi.org/10.1186/s40064-016-2106-8 doi: 10.1186/s40064-016-2106-8
[38]	M. Caputo, Linear models of dissipation whose $Q$ is almost frequency independent–Ⅱ, Geophys. J. Int., 13 (1967), 529–539. https://doi.org/10.1111/j.1365-246X.1967.tb02303.x doi: 10.1111/j.1365-246X.1967.tb02303.x
[39]	A. M. Wazwaz, A reliable technique for solving linear and nonlinear Schrödinger equations by Adomian decomposition method, Bull. Inst. Math. Acad. Sin., 29 (2001), 125–134.
[40]	M. Ghasemi, M. Fardi, R. K. Ghaziani, Numerical solution of nonlinear delay differential equations of fractional order in reproducing kernel Hilbert space, Appl. Math. Comput., 268 (2015), 815–831. https://doi.org/10.1016/j.amc.2015.06.012 doi: 10.1016/j.amc.2015.06.012
[41]	Z. Chen, Q. Q. Gou, Piecewise Picard iteration method for solving nonlinear fractional differential equation with proportional delays, Appl. Math. Comput., 348 (2019), 465–478. https://doi.org/10.1016/j.amc.2018.10.058 doi: 10.1016/j.amc.2018.10.058
[42]	X. Y. Cui, Y. Q. Feng, J. Jiang, Generalized Legendre polynomial configuration method for solving numerical solutions of fractional pantograph delay differential equations, J. Adv. Appl. Comput. Math., 10 (2023), 88–101. https://doi.org/10.15377/2409-5761.2023.10.9 doi: 10.15377/2409-5761.2023.10.9
[43]	C. Q. Yang, J. H. Hou, X. G. Lv, Jacobi spectral collocation method for solving fractional pantograph delay differential equations, Eng. Comput., 38 (2022), 1985–1994. https://doi.org/10.1007/s00366-020-01193-7 doi: 10.1007/s00366-020-01193-7
[44]	U. Saeed, M. ur Rehman, M. A. Iqbal, Modified Chebyshev wavelet methods for fractional delay-type equations, Appl. Math. Comput., 264 (2015), 431–442. https://doi.org/10.1016/j.amc.2015.04.113 doi: 10.1016/j.amc.2015.04.113
[45]	N. A. Elkot, E. H. Doha, I. G. Ameen, A. S. Hendy, M. A. Zaky, A re-scaling spectral collocation method for the nonlinear fractional pantograph delay differential equations with non-smooth solutions, Commun. Nonlinear Sci. Numer. Simul., 118 (2023), 107017. https://doi.org/10.1016/j.cnsns.2022.107017 doi: 10.1016/j.cnsns.2022.107017
[46]	C. Q. Yang, X. G. Lv, Generalized Jacobi spectral Galerkin method for fractional pantograph differential equation, Math. Methods Appl. Sci., 44 (2021), 153–165. https://doi.org/10.1002/mma.6718 doi: 10.1002/mma.6718
[47]	A. H. Bhrawy, A. A. Al-Zahrani, Y. A. Alhamed, D. Baleanu, A new generalized Laguerre-Gauss collocation scheme for numerical solution of generalized fractional pantograph equations, Rom. J. Phys., 59 (2014), 646–657.
[48]	A. El-Ajou, M. N. Oqielat, Z. Al-Zhour, S. Momani, Analytical numerical solutions of the fractional multi-pantograph system: two attractive methods and comparisons, Results Phys., 14 (2019), 102500. https://doi.org/10.1016/j.rinp.2019.102500 doi: 10.1016/j.rinp.2019.102500
[49]	S. Widatalla, M. A. Koroma, Approximation algorithm for a system of pantograph equations, J. Appl. Math., 2012 (2012), 714681. https://doi.org/10.1155/2012/714681 doi: 10.1155/2012/714681

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)

通讯作者: 陈斌, bchen63@163.com

1.
沈阳化工大学材料科学与工程学院沈阳 110142

AIMS Mathematics

1.8 3.4

Metrics

Article views(585) PDF downloads(26) Cited by(0)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

AIMS Mathematics

Investigation of fractional-order pantograph delay differential equations using Sumudu decomposition method

Related Papers:

Abstract

1. Introduction

2. Main concepts and theorems

2.1. Some basic definitions and properties of the fractional calculus

2.2. Some basic definitions and properties of the Sumudu transform

3. Treatment of fractional pantograph delay differential equations via Sumudu decomposition method

4. Applications

5. Conclusions

Author contributions

Use of Generative-AI tools declaration

Conflict of interest

References

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Catalog

Abstract

1. Introduction

2. Main concepts and theorems

2.1. Some basic definitions and properties of the fractional calculus

2.2. Some basic definitions and properties of the Sumudu transform

3. Treatment of fractional pantograph delay differential equations via Sumudu decomposition method

4. Applications

5. Conclusions

Author contributions

Use of Generative-AI tools declaration

Conflict of interest

References

AIMS Mathematics

Investigation of fractional-order pantograph delay differential equations using Sumudu decomposition method

Related Papers:

Abstract

1. Introduction

2. Main concepts and theorems

2.1. Some basic definitions and properties of the fractional calculus

2.2. Some basic definitions and properties of the Sumudu transform

3. Treatment of fractional pantograph delay differential equations via Sumudu decomposition method

4. Applications

5. Conclusions

Author contributions

Use of Generative-AI tools declaration

Conflict of interest

References

Reader Comments

通讯作者: 陈斌, bchen63@163.com

Metrics

Figures and Tables

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog

Abstract

1. Introduction

2. Main concepts and theorems

2.1. Some basic definitions and properties of the fractional calculus

2.2. Some basic definitions and properties of the Sumudu transform

3. Treatment of fractional pantograph delay differential equations via Sumudu decomposition method

4. Applications

5. Conclusions

Author contributions

Use of Generative-AI tools declaration

Conflict of interest

References