<em>S</em>-function associated with fractional derivative and double Dirichlet average

Jitendra Daiya; Dinesh Kumar; Jitendra Daiya; Dinesh Kumar

doi:10.3934/math.2020094

AIMS Mathematics

2020, Volume 5, Issue 2: 1372-1382. doi: 10.3934/math.2020094

Previous Article Next Article

Research article

S-function associated with fractional derivative and double Dirichlet average

Jitendra Daiya ¹,
Dinesh Kumar ^{2
,
,}

1.
Department of Mathematics & Statistics, Jai Narain Vyas University, Jodhpur-342005 (Raj.), India
2.
Department of Applied Sciences, College of Agriculture, Sumerpur-Pali, Agriculture University of Jodhpur, Jodhpur-342304 (Raj.), India

Received: 30 September 2019 Accepted: 08 January 2020 Published: 20 January 2020
MSC : Primary: 26A33, 33C99; Secondary: 33E12, 33E99

The object of this article is to investigate the double Dirichlet averages of S-functions. Representations of such relations are obtained in terms of fractional derivative. Some interesting special cases are also stated.

Keywords:

Citation: Jitendra Daiya, Dinesh Kumar. S-function associated with fractional derivative and double Dirichlet average[J]. AIMS Mathematics, 2020, 5(2): 1372-1382. doi: 10.3934/math.2020094

Related Papers:

[1]	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
[2]	Musawa Yahya Almusawa, Hassan Almusawa . Numerical analysis of the fractional nonlinear waves of fifth-order KdV and Kawahara equations under Caputo operator. AIMS Mathematics, 2024, 9(11): 31898-31925. doi: 10.3934/math.20241533
[3]	Aslı Alkan, Halil Anaç . A new study on the Newell-Whitehead-Segel equation with Caputo-Fabrizio fractional derivative. AIMS Mathematics, 2024, 9(10): 27979-27997. doi: 10.3934/math.20241358
[4]	Humaira Yasmin, Aljawhara H. Almuqrin . Analytical study of time-fractional heat, diffusion, and Burger's equations using Aboodh residual power series and transform iterative methodologies. AIMS Mathematics, 2024, 9(6): 16721-16752. doi: 10.3934/math.2024811
[5]	Mariam Sultana, Muhammad Waqar, Ali Hasan Ali, Alina Alb Lupaş, F. Ghanim, Zaid Ameen Abduljabbar . Numerical investigation of systems of fractional partial differential equations by new transform iterative technique. AIMS Mathematics, 2024, 9(10): 26649-26670. doi: 10.3934/math.20241296
[6]	Humaira Yasmin, Aljawhara H. Almuqrin . Efficient solutions for time fractional Sawada-Kotera, Ito, and Kaup-Kupershmidt equations using an analytical technique. AIMS Mathematics, 2024, 9(8): 20441-20466. doi: 10.3934/math.2024994
[7]	M. Mossa Al-Sawalha, Khalil Hadi Hakami, Mohammad Alqudah, Qasem M. Tawhari, Hussain Gissy . Novel Laplace-integrated least square methods for solving the fractional nonlinear damped Burgers' equation. AIMS Mathematics, 2025, 10(3): 7099-7126. doi: 10.3934/math.2025324
[8]	Meshari Alesemi . Innovative approaches of a time-fractional system of Boussinesq equations within a Mohand transform. AIMS Mathematics, 2024, 9(10): 29269-29295. doi: 10.3934/math.20241419
[9]	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
[10]	Maysaa Al Qurashi, Saima Rashid, Sobia Sultana, Fahd Jarad, Abdullah M. Alsharif . Fractional-order partial differential equations describing propagation of shallow water waves depending on power and Mittag-Leffler memory. AIMS Mathematics, 2022, 7(7): 12587-12619. doi: 10.3934/math.2022697

Abstract

1. Introduction

A novel mathematical tool for characterizing non-local structures is fractional calculus (FC). Mathematical explanations of many physical problems using fractional derivatives have proved successful in recent generations when applied to situations close to reality. Many authors, including Hadamard, Riemann-Liouville, Coimbra, Grunwald-Letnikov, Riesz, Weyl, Liouville Caputo, Atangana-Baleanu, and Caputo-Fabrizio, have offered crucial definitions of fractional operators ^[1,2,3,4]. The underlying principle of these traditional differential equations is their reliance on integer-order derivatives, which give the order an integer numerical value indicating the number of times that a function is differentiated. As opposed to fractional partial differential equations (FPDEs), these concepts are expanded by the addition of fractional derivatives. This kind of model is necessary for the description of delayed or dependent responses, non-local interactions, and anomalous diffusion. These relationships not only have the power to explain the complex phenomena of physics, biology, finance, and engineering, but they also can predict special events. The solution of FPDEs is based on a series of special numerical methods and analytical techniques, all adapted to face the scale-free and non-integer properties of fractional derivatives. Often, the systems required for modeling real-world phenomena follow interrelated processes that can be described using systems of partial differential equations ^[1,2,3,4,9]. Such systems have multiple differential equations, with each one describing how the value of a specific physical quantity or any of the interacting systemic components changes with time as they progress.

Computational models often include partial differential equations (PDEs), which are important in applications such as fluid flow, electromagnetics, population dynamics, and quantized mechanics. By studying the behavior of the components, their interactions, and their relationships with each other, one can understand how patterns and dynamics are formed and how the system stabilizes. Studying PDE systems with coupled terms is extremely complex and requires advanced mathematics, such as numerical simulations, perturbation methods, and symmetry analysis, to discern a solution ^{[10,11,12,13,14]}. The Hermite colocation method ^[10], the optimal homotopy asymptotic technique ^[11], the Adomian decomposition method ^[12], the homotopy perturbation transform method ^[13], the Pade approximation and homotopy-Pade technique ^[14], the invariant subspace method ^[20], the q-homotopy analysis transform method ^[21], the homotopy analysis Sumudu transform method ^[22], and the Sumudu transform series expansion method ^[23] are some of the sophisticated approaches developed for finding exact solutions to nonlinear FPDE models ^[24,25,26]. If perturbation methods are not used, the homotopy analysis method breaks a problem into an endless series of linear problems. This method employs the concept of homotopy from topology to derive a convergent series solution ^[27,28]. An approach to homotopy analysis proposed by Liao ^[29], and the Laplace transform ^[30] are combined in the Laplace homotopy perturbation method.

A gradient of chemical molecules guides the movement of cells, a process known as chemotaxis, which is essential for cell population self-organization and developmental biology in general. In 1970, Lee Segel and Evelyn Keller presented the first mathematical model of chemotaxis. To further understand how the mould aggregation process works in the chemical-attraction-based cellular slime, they used parabolic approaches ^[31]. Here, we take a look at the fractional-order system of a KS model that goes like this:

$\begin{equation} \begin{split} &D_{{\varphi}}^{p} {\beta_{1}}({\psi}, {\varphi})-a \frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2}+\frac{\partial}{\partial \psi}\Big({\beta_{1}}(\psi,\varphi) \frac{\partial\varpi({\beta_{2}})}{\partial \psi}\Big) = 0,\\ &D_{{\varphi}}^{p} {\beta_{2}}({\psi}, {\varphi})-b \frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}-c {\beta_{1}}({\psi}, {\varphi})+d {\beta_{2}}({\psi}, {\varphi}) = 0, \ \ \text { where }\ \ 0 < p \leq 1, \end{split} \end{equation}$

(1.1)

having IC's:

$\begin{equation} \begin{split} {\beta_{1}}({\psi}, 0) = {\beta_{1}}_0(\psi),\\ {\beta_{2}}({\psi}, 0) = {\beta_{2}}_0(\psi). \end{split} \end{equation}$

(1.2)

The concentration of amoebae are indicated by the unknown term ${\beta_{1}}(\psi, \varphi)$ , while the chemical substance of concentration is expressed by ${\beta_{2}}(\psi, \varphi)$ ; $\frac{\partial}{\partial \psi}\Big({\beta_{1}}(\psi, \varphi) \frac{\partial\varpi({\beta_{2}})}{\partial \psi}\Big)$ ; stands for the chemotactic word, indicating that the chemicals are attractive to and sensitive to the cells. The sensitivity function is denoted by $\varpi({\beta_{2}})$ , and $a, b, c,$ and $d$ are positive constants. The parameter $0 < p \leq 1$ represents the order of the fractional derivative. Much recent research has focused on the KS model. For example, to solve the KS model, Atangana used a combination of methods, including a modified homotopy perturbation, the homotopy decomposition, and the Laplace transform approach ^[32,33,34]. Zayernouri established a fractional class of implicit Adams-Moulton and explicit Adams-Bashforth methods in ^[35] and so on ^[36,37,38].

According to ^[39], the residual power series method (RPSM) was developed in 2013 by a Jordanian mathematician named Omar Abu Arqub. The RPSM is a semi-analytical approach that uses Taylor's series to integrate the residual error function. It finds convergence series solutions for differential equations. In 2013, RPSM was first used to resolve fuzzy differential equations. A new RPSM method was created by Arqub et al. ^[40] to quickly get power series solutions for ordinary differential equations (DEs). A new and attractive RPSM approach for fractional DEs problems was developed by Arqub et al. ^[41]. A novel iterative technique to estimate fractional KdV-burgers equations was presented by El-Ajou et al. ^[42] utilizing RPSM. A unique method was developed by Xu et al. ^[43] for solving Boussinesq DEs with fractional power series. Zhang et al. ^[44] stated that a trustworthy numerical approach was developed. More readings on RPSM may be found in ^[45,46,47].

To resolve fractional-order differential equations (FODEs), the research team used two separate approaches. One approach to solving the updated equation is to project it into the space generated by the Aboodh transform. Next, the original equation may be solved by using the inverse Aboodh transform ^[48]. This novel methodology combines the Sumudu transform with the homotopy perturbation method. Without discretization, linearization, or perturbation, this novel approach may solve PDEs as power series expansions, irrespective of their linearity or nonlinearity. There is a significant reduction in the computations needed to find the coefficients compared to RPSM, which requires several repetitions of calculating distinct fractional derivatives throughout the solution phases. The proposed approach has the potential to provide an accurate and closed-form approximation solution.

The Aboodh transform iterative technique (ATIM) is a significant mathematical achievement for fractional partial differential equations. Complexity and convergence issues may develop when using traditional techniques to solve partial differential equations with fractional derivatives. Keeping a steady computational economy while continually improving approximations allows our new strategy to improve accuracy continuously, avoiding these limits. Due to this discovery, we can tackle difficult problems in applied mathematics, engineering, and physics, which enhances our capacity to identify and understand complex systems governed by fractional partial differential equations ^[49,50,51].

The two most basic approaches to solving fractional differential equations are the Aboodh transform iterative technique (ATIM) and the Aboodh residual power series method (ARPSM) ^{[49,50,51,52,53]}, respectively. These techniques not only provide numerical solutions to PDEs that do not need discretization or linearization but also make the symbolic terms in analytical solutions instantly and visible. The primary objective of this study is to compare and contrast the performance of ARPSM and ATIM in solving the Keller-Segel (KS) model. It is worth mentioning that several linear and nonlinear fractional differential problems have been solved using these two approaches.

2. Basic concepts

Definition 2.1. ^[54] Let us assume that the function ${\beta_{1}}({{\psi}}, {{\varphi}})$ is piecewise continuous with exponential order. The Aboodh transform (AT) is defined as follows, assuming $\tau\geq0$ for ${\beta_{1}}({{\psi}}, {{\varphi}})$ ,

$A[{\beta_{1}}({{\psi}}, {{\varphi}})] = \Psi({{\psi}},\xi) = \frac{1}{\xi}\int_{0}^{\infty}{\beta_{1}}({{\psi}}, {{\varphi}})e^{-{{\varphi}}\xi}d{{\varphi}},\ \ {r}_{1}\leq\xi\leq{r}_{2}.$

The Aboodh inverse transform (AIT) is specifically described as follows:

$A^{-1}[\Psi({{\psi}},\xi)] = {\beta_{1}}({{\psi}}, {{\varphi}}) = \frac{1}{2\pi i}\int_{u-i\infty}^{u+i\infty}\Psi({{\psi}}, {{\varphi}})\xi e^{{{\varphi}}\xi}d{{\varphi}},$

where ${{\psi}} = ({{\psi}}_{1}, {{\psi}}_{2}, \cdots, {{\psi}}_{p})\in \mathbb{R}$ and $p\in\mathbb{N}$ .

Lemma 2.2. Let ^[55,56] ${\beta_{1}}_{1}({{\psi}}, \varphi)$ and ${\beta_{1}}_{2}({{\psi, \varphi}})$ are two functions. It is assumed that they are piecewise continuous on $[0, \infty[$ and exponentially ordered. Let $A[{\beta_{1}}_{1}({{\psi}}, {{\varphi}})] = \Psi_{1}({{\psi}}, {{\varphi}}), A[{\beta_{1}}_{2}({{\psi}}, {{\varphi}})] = \Psi_{2}({{\psi}}, {{\varphi}})$ and $\chi_{1}, \chi_{2}$ are constants. Thus, the following characteristics are true:

(1) $A[\chi_{1}{\beta_{1}}_{1}({{\psi}}, {{\varphi}})+\chi_{2}{\beta_{1}}_{2}({{\psi}}, {{\varphi}})] = \chi_{1}\Psi_{1}({{\psi}}, \xi)+\chi_{2}\Psi_{2}({{\psi}}, {{\varphi}})$ ,

(2) $A^{-1}[\chi_{1}\Psi_{1}({{\psi}}, {{\varphi}})+\chi_{2}\Psi_{2}({{\psi}}, {{\varphi}})] = \chi_{1}{\beta_{1}}_{1}({{\psi}}, \xi)+\chi_{2}{\beta_{1}}_{2}({{\psi}}, {{\varphi}})$ ,

(3) $A[J_{{{\varphi}}}^{{p}}{\beta_{1}}({{\psi}}, {{\varphi}})] = \frac{\Psi({{\psi}}, \xi)}{\xi^{{p}}}$ ,

(4) $A[D_{{{\varphi}}}^{{p}}{\beta_{1}}({{\psi}}, {{\varphi}})] = \xi^{{p}}\Psi({{\psi}}, \xi)-\sum_{K = 0}^{r-1}\frac{{\beta_{1}}^{K}({{\psi}}, 0)}{\xi^{K-{p}+2}}, r-1 < {p}\leq r, \ r\in \mathbb{N}$ .

Definition 2.3. ^[57] In terms of order $p$ , the Caputo defines the fractional derivative of the function ${\beta_{1}}({{\psi}}, {{\varphi}})$ as:

$D_{{{\varphi}}}^{{p}}{\beta_{1}}({{\psi}}, {{\varphi}}) = J_{{{\varphi}}}^{m-{p}}{\beta_{1}}^{(m)}({{\psi}}, {{\varphi}}), \ r\geq 0, \ m-1 < {p}\leq m,$

where ${{\psi}} = ({{\psi}}_{1}, {{\psi}}_{2}, \cdots, {{\psi}}_{p})\in \mathbb{R}^{p}$ and $m, p\in R, J_{{{\varphi}}}^{m-{p}}$ is the R-L integral of ${\beta_{1}}({{\psi}}, {{\varphi}})$ .

Definition 2.4. ^[58] Following is the structure of the power series notation:

$\sum\limits_{r = 0}^{\infty}\hbar_{r}({{\psi}})({{\varphi}}-{{\varphi}}_{0})^{r{p}} = \hbar_{0}({{\varphi}}-{{\varphi}}_{0})^{0}+\hbar_{1}({{\varphi}}-{{\varphi}}_{0})^{{p}}+\hbar_{2}({{\varphi}}-{{\varphi}}_{0})^{2{p}} +\cdots,$

where ${{\psi}} = ({{\psi}}_{1}, {{\psi}}_{2}, \cdots, {{\psi}}_{p})\in \mathbb{R}^{p}$ and $p\in \mathbb{N}$ . The series concerning ${{\varphi}}_{0}$ is referred to as a multiple fractional power series (MFPS), where the series coefficients are $\hbar_{r}({{\psi}})'s$ and ${{\varphi}}$ is variable.

Lemma 2.5. Let us suppose that the exponential order function is ${\beta_{1}}({{\psi}}, {{\varphi}})$ . In this case, the AT is defined as: $A[{\beta_{1}}({{\psi}}, {{\varphi}})] = \Psi({{\psi}}, \xi)$ . Hence,

$\begin{equation} A[D_{{{\varphi}}}^{r{p}}{\beta_{1}}({{\psi}}, {{\varphi}})] = \xi^{r{p}}\Psi({{\psi}},\xi)-\sum\limits_{j = 0}^{r-1}\xi^{{p}(r-j)-2}D_{{{\varphi}}}^{j{p}}{\beta_{1}}({{\psi}},0),\ \ 0 < {p}\leq 1, \end{equation}$

(2.1)

where ${{\psi}} = ({{\psi}}_{1}, {{\psi}}_{2}, \cdots, {{\psi}}_{p})\in \mathbb{R}^{p}$ and $p\in \mathbb{N}$ and $D_{{{\varphi}}}^{r{p}} = D_{{{\varphi}}}^{{p}}.D_{{{\varphi}}}^{{p}}.\cdots.D_{{{\varphi}}}^{{p}}(r-times)$ .

Proof. By induction, we are able to illustrate Eq (2.5). When $r = 1$ is used in Eq (2.5), the following results are obtained:

$A[D_{{{\varphi}}}^{2{p}}{\beta_{1}}({{\psi}}, {{\varphi}})] = \xi^{2{p}}\Psi({{\psi}},\xi)-\xi^{2{p}-2}{\beta_{1}}({{\psi}},0)-\xi^{{p}-2}D_{{{\varphi}}}^{{p}}{\beta_{1}}({{\psi}},0).$

Equation (2.5) is true for $r = 1$ , according to Lemma 2.2, part (4). After substituting $r = 2$ in Eq (2.5), we get:

$\begin{equation} A[D_{r}^{2{p}}{\beta_{1}}({{\psi}}, {{\varphi}})] = \xi^{2{p}}\Psi({{\psi}},\xi)-\xi^{2{p}-2}{\beta_{1}}({{\psi}},0)-\xi^{{p}-2}D_{{{\varphi}}}^{{p}}{\beta_{1}}({{\psi}},0). \end{equation}$

(2.2)

Equation (2.2) L.H.S. enables us to determine

$\begin{equation} L.H.S = A[D_{{{\varphi}}}^{2{p}}{\beta_{1}}({{\psi}}, {{\varphi}})]. \end{equation}$

(2.3)

The following way is used to express Eq (2.3):

$\begin{equation} L.H.S = A[D_{{{\varphi}}}^{{p}}{\beta_{1}}({{\psi}}, {{\varphi}})]. \end{equation}$

(2.4)

Assume

$\begin{equation} z({{\psi}}, {{\varphi}}) = D_{{{\varphi}}}^{{p}}{\beta_{1}}({{\psi}}, {{\varphi}}). \end{equation}$

(2.5)

Thus, Eq (2.4) becomes

$\begin{equation} L.H.S = A[D_{{{\varphi}}}^{{p}}z({{\psi}}, {{\varphi}})]. \end{equation}$

(2.6)

Implementing the Caputo derivative led to a modification in Eq (2.6).

$\begin{equation} L.H.S = A[J^{1-{p}}z^{'}({{\psi}}, {{\varphi}})]. \end{equation}$

(2.7)

Equation (2.7) provides the R-L integral for AT, which allows us to deduce the following:

$\begin{equation} L.H.S = \frac{A[z^{'}({{\psi}}, {{\varphi}})]}{\xi^{1-{p}}}. \end{equation}$

(2.8)

Equation (2.8) is changed into the following form by using the differential characteristic of the AT:

$\begin{equation} L.H.S = \xi^{{p}}Z({{\psi}},\xi)-\frac{z({{\psi}},0)}{\xi^{2-{p}}}, \end{equation}$

(2.9)

from Eq (2.5), we obtain:

$Z({{\psi}},\xi) = \xi^{{p}}\Psi({{\psi}},\xi)-\frac{{\beta_{1}}({{\psi}},0)}{\xi^{2-{p}}},$

where $A[z({{\psi}}, {{\varphi}})] = Z({{\psi}}, \xi)$ . Therefore, Eq (2.9) is transformed to

$\begin{equation} L.H.S = \xi^{2{p}}\Psi({{\psi}},\xi)-\frac{{\beta_{1}}({{\psi}},0)}{\xi^{2-2{p}}}-\frac{D_{{{\varphi}}}^{{p}}{\beta_{1}}({{\psi}},0)}{\xi^{2-{p}}}, \end{equation}$

(2.10)

when $r = K$ . Equations (2.5) and (2.10) are compatible. For $r = K$ , let's assume that Eq (2.5) holds. Therefore, we substitute $r = K$ into Eq (2.5):

$\begin{equation} A[D_{{{\varphi}}}^{K{p}}{\beta_{1}}({{\psi}}, {{\varphi}})] = \xi^{K{p}}\Psi({{\psi}},\xi)-\sum\limits_{j = 0}^{K-1}\xi^{{p}(K-j)-2}D_{{{\varphi}}}^{j{p}}D_{{{\varphi}}}^{j{p}}{\beta_{1}}({{\psi}},0),\ 0 < {p}\leq1. \end{equation}$

(2.11)

Next, we will show how to solve Eq (2.5) for $r = K+1$ . Based on Eq (2.5), we may express

$\begin{equation} A[D_{{{\varphi}}}^{(K+1){p}}{\beta_{1}}({{\psi}}, {{\varphi}})] = \xi^{(K+1){p}}\Psi({{\psi}},\xi)-\sum\limits_{j = 0}^{K}\xi^{{p}((K+1)-j)-2}D_{{{\varphi}}}^{j{p}}{\beta_{1}}({{\psi}},0). \end{equation}$

(2.12)

After examining the left side of Eq (2.12), we get

$\begin{equation} L.H.S = A[D_{{{\varphi}}}^{K{p}}(D_{{{\varphi}}}^{K{p}})], \end{equation}$

(2.13)

let

$D_{{{\varphi}}}^{K{p}} = g({{\psi}}, {{\varphi}}),$

by Eq (2.13), we drive

$\begin{equation} L.H.S = A[D_{{{\varphi}}}^{{p}}g({{\psi}}, {{\varphi}})]. \end{equation}$

(2.14)

Equation (2.14) is modified to provide the following result by using the R-L integral and Caputo derivative:

$\begin{equation} L.H.S = \xi^{{p}}A[D_{{{\varphi}}}^{K{p}}{\beta_{1}}({{\psi}}, {{\varphi}})]-\frac{g({{\psi}},0)}{\xi^{2-{p}}}. \end{equation}$

(2.15)

Equation (2.15) is derived from Eq (2.11),

$\begin{equation} L.H.S = \xi^{r{p}}\Psi({{\psi}},\xi)-\sum\limits_{j = 0}^{r-1}\xi^{{p}(r-j)-2}D_{{{\varphi}}}^{j{p}}{\beta_{1}}({{\psi}},0). \end{equation}$

(2.16)

In addition, the outcome that follows is obtained from Eq (2.16):

$L.H.S = A[D_{{{\varphi}}}^{r{p}}{\beta_{1}}({{\psi}},0)].$

Thus, for $r = K+1$ , the Eq (2.5) is valid. Equation (2.5) is valid for all positive integers according to the mathematical induction method. □

Here, we find another novel way of looking to MFTS, or multiple fractional Taylor's series. The ARPSM, which will be discussed in more depth later on, will benefit from this formula.

Lemma 2.6. Assume that ${\beta_{1}}({{\psi}}, {{\varphi}})$ represents the exponential order function. The expression $A[{\beta_{1}}({{\psi}}, {{\varphi}})] = \Psi({{\psi}}, \xi)$ is the AT of ${\beta_{1}}({{\psi}}, {{\varphi}})$ . The AT MFTS notation looks like this:

$\begin{equation} \Psi({{\psi}},\xi) = \sum\limits_{r = 0}^{\infty}\frac{\hbar_{r}({{\psi}})}{\xi^{r{p}+2}}, \xi > 0, \end{equation}$

(2.17)

where, ${{\psi}} = (s_{1}, {{\psi}}_{2}, \cdots, {{\psi}}_{p})\in\mathbb{R}^{p}, \ p\in \mathbb{N}$ .

Proof. Let us investigate Taylor's series' fractional order expression:

$\begin{equation} {\beta_{1}}({{\psi}}, {{\varphi}}) = \hbar_{0}({{\psi}})+\hbar_{1}({{\psi}})\frac{{{\varphi}}^{{p}}}{\Gamma[{p}+1]}++\hbar_{2}({{\psi}})\frac{{{\varphi}}^{2{p}}}{\Gamma[2{p}+1]}+\cdots. \end{equation}$

(2.18)

The following equality is obtained by applying the AT to Eq (2.18):

$A\left[{\beta_{1}}({{\psi}}, {{\varphi}})\right] = A\left[\hbar_{0}({{\psi}})\right]+A\left[\hbar_{1}({{\psi}})\frac{{{\varphi}}^{{p}}}{\Gamma[{p}+1]}\right]+A\left[\hbar_{1}({{\psi}}) \frac{{{\varphi}}^{2{p}}}{\Gamma[2{p}+1]}\right]+\cdots$

This is accomplished by using the AT's characteristics.

$A\left[{\beta_{1}}({{\psi}}, {{\varphi}})\right] = \hbar_{0}({{\psi}})\frac{1}{\xi^{2}}+\hbar_{1}({{\psi}})\frac{\Gamma[{p}+1]}{\Gamma[{p}+1]}\frac{1}{\xi^{{p}+2}}+\hbar_{2}({{\psi}}) \frac{\Gamma[2{p}+1]}{\Gamma[2{p}+1]}\frac{1}{\xi^{2{p}+2}}\cdots$

A distinct variant of Taylor's series in the AT is therefore obtained. □

Lemma 2.7. As stated in the new form of Taylor's series 2.17, the MFPS may be represented as $A[{\beta_{1}}(\psi, \varphi)] = \Psi(\psi, \xi)$ .

$\begin{equation} \hbar_{0}({{\psi}}) = \lim\limits_{\xi\rightarrow \infty}\xi^2\Psi({{\psi}},\xi) = {\beta_{1}}({{\psi}},0). \end{equation}$

(2.19)

Proof. This can be determined from the revised version of Taylor's series:

$\begin{equation} \hbar_{0}({{\psi}}) = \xi^2\Psi({{\psi}},\xi)-\frac{\hbar_{1}({{\psi}})}{\xi^{{p}}}-\frac{\hbar_{2}({{\psi}})}{\xi^{2{p}}}-\cdots \end{equation}$

(2.20)

As shown in Eq (2.20), the necessary solution may be obtained by evaluating $\lim_{\xi\rightarrow \infty}$ into Eq (2.19) and doing a quick computation. □

Theorem 2.8. The function $A[{\beta_{1}}({{\psi}}, {{\varphi}})] = \Psi({{\psi}}, \xi)$ may be expressed in MFPS form as follows:

$\Psi({{\psi}},\xi) = \sum\limits_{0}^{\infty}\frac{\hbar_{r}({{\psi}})}{\xi^{r{p}+2}}, \ \xi > 0,$

where ${{\psi}} = ({{\psi}}_{1}, {{\psi}}_{2}, \cdots, {{\psi}}_{p})\in\mathbb{R}^p$ and $p\in\mathbb{N}$ . Then we have

$\hbar_{r}({{\psi}}) = D_{r}^{r{p}}{\beta_{1}}({{\psi}},0),$

where, $D_{{{\varphi}}}^{r{p}} = D_{{{\varphi}}}^{{p}}.D_{{{\varphi}}}^{{p}}.\cdots.D_{{{\varphi}}}^{{p}}(r-times)$ .

Proof. The new Taylor's series is as follows:

$\begin{equation} \hbar_{1}({{\psi}}) = \xi^{{p}+2}\Psi({{\psi}},\xi)-\xi^{{p}}\hbar_{0}({{\psi}})-\frac{\hbar_{2}({{\psi}})}{\xi^{{p}}}-\frac{\hbar_{3}({{\psi}})}{\xi^{2{p}}}- \cdots \end{equation}$

(2.21)

$\lim_{\xi\rightarrow \infty}$ , is applied to (2.21), we get

$\hbar_{1}({{\psi}}) = \lim\limits_{\xi\rightarrow \infty}(\xi^{{p}+2}\Psi({{\psi}},\xi)-\xi^{{p}}\hbar_{0}({{\psi}}))-\lim\limits_{\xi\rightarrow \infty}\frac{\hbar_{2}({{\psi}})}{\xi^{{p}}} -\lim\limits_{\xi\rightarrow \infty}\frac{\hbar_{3}({{\psi}})}{\xi^{2{p}}}- \cdots$

After calculating the limit, we have the following equality:

$\begin{equation} \hbar_{1}({{\psi}}) = \lim\limits_{\xi\rightarrow \infty}(\xi^{{p}+2}\Psi({{\psi}},\xi)-\xi^{{p}}\hbar_{0}({{\psi}})). \end{equation}$

(2.22)

The result of inserting Lemma 2.5 into Eq (2.22) is as follows:

$\begin{equation} \hbar_{1}({{\psi}}) = \lim\limits_{\xi\rightarrow \infty}(\xi^2A[D_{{{\varphi}}}^{{p}}{\beta_{1}}({{\psi}}, {{\varphi}})](\xi)). \end{equation}$

(2.23)

Furthermore, it is transformed into by using Lemma 2.6 to Eq (2.23),

$\hbar_{1}({{\psi}}) = D_{{{\varphi}}}^{{p}}{\beta_{1}}({{\psi}},0).$

Again, applying limit $\xi\rightarrow \infty$ and using the new form of Taylor's series, we obtain:

$\hbar_{2}({{\psi}}) = \xi^{2{p}+2}\Psi({{\psi}},\xi)-\xi^{2{p}}\hbar_{0}({{\psi}})-\xi^{{p}}\hbar_{1}({{\psi}})-\frac{\hbar_{3}({{\psi}})}{\xi^{{p}}}- \cdots$

We get the result from Lemma 2.6.

$\begin{equation} \hbar_2({{\psi}}) = \lim\limits_{\xi\rightarrow \infty}\xi^2(\xi^{2{p}}\Psi({{\psi}},\xi)-\xi^{2{p}-2}\hbar_{0}({{\psi}})-\xi^{{p}-2}\hbar_{1}({{\psi}})). \end{equation}$

(2.24)

Using Lemmas 2.5 and 2.7, we convert Eq (2.24) into

$\hbar_{2}({{\psi}}) = D_{{{\varphi}}}^{2{p}}{\beta_{1}}({{\psi}},0),$

when the new Taylor's series is put through the same process, the following results are obtained:

$\hbar_{3}({{\psi}}) = \lim\limits_{\xi\rightarrow \infty}\xi^2(A[D_{{{\varphi}}}^{2{p}}{\beta_{1}}({{\psi}},{p})](\xi)),$

Lemma 2.7 is used to derive the final equation:

$\hbar_{3}({{\psi}}) = D_{{{\varphi}}}^{3{p}}{\beta_{1}}({{\psi}},0),$

in general

$\hbar_{r}({{\psi}}) = D_{{{\varphi}}}^{r{p}}{\beta_{1}}({{\psi}},0).$

Consequently, proof ends here. □

The principles regulating the convergence of Taylor's series in its new form are explained and proven in the following theorem.

Theorem 2.9. Presented in Lemma 2.6, the formula for multiple fractional Taylor's series may be represented in the following new form: $A [{\beta_{1}}({{\psi}}, {{\varphi}})] = \Psi({{\psi}}, \xi)$ . When $|\xi^{a}A[D_{{{\varphi}}}^{(K+1){p}}{\beta_{1}}({{\psi}}, {{\varphi}})]|\leq T$ , for all $0 < \xi \leq s$ and $0 < {p} \leq 1$ , the following inequality satisfies the residual $R_K({{\psi}}, \xi)$ of the new MFTS:

$|R_{K}({{\psi}},\xi)|\leq \frac{T}{\xi^{(K = 1){p}+2}}, \ 0 < \xi\leq s.$

Proof. Let $A[D_{{{\varphi}}}^{r{p}}{\beta_{1}}({{\psi}}, {{\varphi}})](\xi)$ is defined on $0 < \xi \leq s$ for $r = 0, 1, 2, \cdots, K + 1$ . Let us assume that $|{\xi^{2}A[{D_{{psi}}, \varphi}]^{K+1}}{\beta_{1}}^{1}]| \leq T, \ on\ 0 < \xi \leq s$ . Determine the following relation using the new Taylor's series:

$\begin{equation} R_K({{\psi}},\xi) = \Psi({{\psi}},\xi)-\sum\limits_{r = 0}^K\frac{\hbar_r({{\psi}})}{\xi^{r{p}+2}}. \end{equation}$

(2.25)

Equation (2.25) is converted using Theorem 2.8,

$\begin{equation} R_K({{\psi}},\xi) = \Psi({{\psi}},\xi)-\sum\limits_{r = 0}^{K}\frac{D_{{{\varphi}}}^{r{p}}{\beta_{1}}({{\psi}},0)}{\xi^{r{p}+2}}. \end{equation}$

(2.26)

To solve Eq (2.26), multiply $\xi^{(K+1)a+2}$ on both sides,

$\begin{equation} \xi^{(K+1){p}+2}R_K({{\psi}},\xi) = \xi^{2}(\xi^{(K+1){p}}\Psi({{\psi}},\xi)-\sum\limits_{r = 0}^{K}\xi^{(K+1-r){p}-2}D_{{{\varphi}}}^{r{p}}{\beta_{1}}({{\psi}},0)). \end{equation}$

(2.27)

Lemma 2.5 applied to Eq (2.27) yields

$\begin{equation} \xi^{(K+1){p}+2}R_K({{\psi}},\xi) = \xi^{2}A[D_{{{\varphi}}}^{(K+1){p}}{\beta_{1}}({{\psi}}, {{\varphi}})]. \end{equation}$

(2.28)

Taking absolute of Eq (2.28), we get

$\begin{equation} |\xi^{(K+1){p}+2}R_K({{\psi}},\xi)| = |\xi^{2}A[D_{{{\varphi}}}^{(K+1){p}}{\beta_{1}}({{\psi}}, {{\varphi}})]|. \end{equation}$

(2.29)

Applying the criteria listed in Eq (2.29) yields the following result:

$\begin{equation} \frac{-T}{\xi^{(K+1){p}+2}}\leq R_K({{\psi}},\xi)\leq \frac{T}{\xi^{(K+1){p}+2}}. \end{equation}$

(2.30)

We use Eq (2.30) to get the necessary result,

$|R_K({{\psi}},\xi)|\leq \frac{T}{\xi^{(K+1){p}+2}}.$

Thus, a new series convergence criteria is developed. □

3. Proposed methodologies

3.1. Aboodh residual power series method (ARPSM)

The ARPSM rules served as the foundation for our general model solution, which we describe below.

Step 1. The general equation may be simplified to obtain:

$\begin{equation} D_{{{\varphi}}}^{{p}}{\beta_{1}}({{\psi}}, {{\varphi}})+\vartheta({{\psi}})N({\beta_{1}})-\delta({{\psi}},{\beta_{1}}) = 0. \end{equation}$

(3.1)

Step 2. The two sides of Eq (3.1) are evaluated using the AT in order to get

$\begin{equation} A[D_{{{\varphi}}}^{{p}}{\beta_{1}}({{\psi}}, {{\varphi}})+\vartheta({{\psi}})N({\beta_{1}})-\delta({{\psi}},{\beta_{1}})] = 0, \end{equation}$

(3.2)

transformation of Eq (3.2) by using Lemma 2.5. Thus,

$\begin{equation} \Psi({{\psi}},{s}) = \sum\limits_{j = 0}^{q-1}\frac{D_{{{\varphi}}}^{j}{\beta_{1}}({{\psi}},0)}{{s}^{q{p}+2}}-\frac{\vartheta({{\psi}})Y({s})}{{s}^{q{p}}}+\frac{F({{\psi}},{s})}{{s}^{q{p}}}, \end{equation}$

(3.3)

where, $A[\delta({{\psi}}, {\beta_{1}})] = F({{\psi}}, {s}), A[N({\beta_{1}})] = Y({s})$ .

Step 3. Examine the form that the solution to Eq (3.3) takes:

$\Psi({{\psi}}, {s}) = \sum\limits_{r = 0}^{\infty}\frac{\hbar_r({{\psi}})}{{s}^{r{p}+2}}, \ {s} > 0,$

Step 4. To proceed, follow these steps:

$\hbar_0({{\psi}}) = \lim\limits_{{s}\rightarrow \infty}{s}^2 \Psi({{\psi}}, {s}) = {\beta_{1}}({{\psi}}, 0).$

The following outcome by using Theorem 2.9:

$\hbar_{1}({{\psi}}) = D_{{{\varphi}}}^{{p}}{\beta_{1}}({{\psi}}, 0),$

$\hbar_{2}({{\psi}}) = D_{{{\varphi}}}^{2{p}}{\beta_{1}}({{\psi}}, 0),$

$\vdots$

$\hbar_{w}({{\psi}}) = D_{{{\varphi}}}^{w{p}}{\beta_{1}}({{\psi}}, 0).$

Step 5. The $\Psi({{\psi}}, {s})$ series that has been $K^{th}$ truncated may be found using the formula below:

$\Psi_{K}({{\psi}}, {s}) = \sum\limits_{r = 0}^{K}\frac{\hbar_r({{\psi}})}{{s}^{r{p}+2}}, \ {s} > 0,$

$\Psi_{K}({{\psi}}, {s}) = \frac{\hbar_{0}({{\psi}})}{{s}^{2}}+\frac{\hbar_{1}({{\psi}})}{{s}^{{p}+2}}+\cdots+\frac{\hbar_{w}({{\psi}})}{{s}^{w{p}+2}}+\sum\limits_{r = w+1}^{K} \frac{\hbar_r({{\psi}})}{{s}^{r{p}+2}}.$

Step 6. Remember that, in order to derive the following, you must take into consideration both the $K^{th}$ -truncated Aboodh residual function and the Aboodh residual function (ARF) from (3.3) separately:

$ARes({{\psi}}, {s}) = \Psi({{\psi}},{s})-\sum\limits_{j = 0}^{q-1}\frac{D_{{{\varphi}}}^{j}{\beta_{1}}({{\psi}},0)}{{s}^{j{p}+2}}+\frac{\vartheta({{\psi}})Y({s})}{{s}^{j{p}}}-\frac{F({{\psi}},{s})}{{s}^{j{p}}},$

and

$\begin{equation} ARes_{K}({{\psi}}, {s}) = \Psi_{K}({{\psi}},{s})-\sum\limits_{j = 0}^{q-1}\frac{D_{{{\varphi}}}^{j}{\beta_{1}}({{\psi}},0)}{{s}^{j{p}+2}}+\frac{\vartheta({{\psi}})Y({s})}{{s}^{j{p}}}-\frac{F({{\psi}},{s})}{{s}^{j{p}}}. \end{equation}$

(3.4)

Step 7. Put $\Psi_{K}({{\psi}}, {s})$ into Eq (3.4) rather than use its expansion form,

$\begin{equation} \begin{split} ARes_{K}({{\psi}}, {s}) = &\Big(\frac{\hbar_{0}({{\psi}})}{{s}^{2}}+\frac{\hbar_{1}({{\psi}})}{{s}^{{p}+2}}+\cdots+\frac{\hbar_{w}({{\psi}})}{{s}^{w{p}+2}}+\sum\limits_{r = w+1}^{K} \frac{\hbar_r({{\psi}})}{{s}^{r{p}+2}}\Big)\\&-\sum\limits_{j = 0}^{q-1}\frac{D_{{{\varphi}}}^{j}{\beta_{1}}({{\psi}},0)}{{s}^{j{p}+2}}+\frac{\vartheta({{\psi}})Y({s})}{{s}^{j{p}}} -\frac{F({{\psi}},{s})}{{s}^{j{p}}}. \end{split} \end{equation}$

(3.5)

Step 8. Equation (3.5) may be solved by multiplying both sides by ${s}^{K{p}+2}$ ,

$\begin{equation} \begin{split} {s}^{K{p}+2}ARes_{K}({{\psi}}, {s}) = &{s}^{K{p}+2}\Big(\frac{\hbar_{0}({{\psi}})}{{s}^{2}}+\frac{\hbar_{1}({{\psi}})}{{s}^{{p}+2}}+\cdots+\frac{\hbar_{w}({{\psi}})}{{s}^{w{p}+2}}+\sum\limits_{r = w+1}^{K} \frac{\hbar_r({{\psi}})}{{s}^{r{p}+2}}\\&-\sum\limits_{j = 0}^{q-1}\frac{D_{{{\varphi}}}^{j}{\beta_{1}}({{\psi}},0)}{{s}^{j{p}+2}}+\frac{\vartheta({{\psi}})Y({s})}{{s}^{j{p}}} -\frac{F({{\psi}},{s})}{{s}^{j{p}}}\Big). \end{split} \end{equation}$

(3.6)

Step 9. After taking $\lim_{{s}\rightarrow \infty}$ , we calculate the solution to Eq (3.6), which is:

$\begin{equation*} \begin{split} \lim\limits_{{s}\rightarrow \infty}{s}^{K{p}+2}ARes_{K}({{\psi}}, {s}) = &\lim\limits_{{s}\rightarrow \infty}{s}^{K{p}+2}\Big(\frac{\hbar_{0}({{\psi}})}{{s}^{2}}+\frac{\hbar_{1}({{\psi}})}{{s}^{{p}+2}}+\cdots+\frac{\hbar_{w}({{\psi}})}{{s}^{w{p}+2}}+\sum\limits_{r = w+1}^{K} \frac{\hbar_r({{\psi}})}{{s}^{r{p}+2}}\\&-\sum\limits_{j = 0}^{q-1}\frac{D_{{{\varphi}}}^{j}{\beta_{1}}({{\psi}},0)}{{s}^{j{p}+2}}+\frac{\vartheta({{\psi}})Y({s})}{{s}^{j{p}}} -\frac{F({{\psi}},{s})}{{s}^{j{p}}}\Big). \end{split} \end{equation*}$

Step 10. Solving the above equation will provide the value of $\hbar_K({{\psi}})$ ,

$\lim\limits_{{s}\rightarrow \infty}({s}^{K{p}+2}ARes_{K}({{\psi}}, {s})) = 0,$

where $K = w+1, w+2, \cdots.$

Step 11. Using a $K$ -truncated series of $\Psi({{\psi}}, {s})$ , replace the values of $\hbar_K({{\psi}})$ to get the $K$ -approximate solution of Eq (3.3).

Step 12. Solve $\Psi_K({{\psi}}, {s})$ using the AIT to get the $K$ -approximate solution ${\beta_{1}}_K({{\psi}}, {{\varphi}})$ .

3.2. Aboodh transform iteration method (ATIM)

Consider the following PDE of space and time fractional order:

$\begin{equation} D_{{{\varphi}}}^{{p}}{{{\beta_{1}}}}({{\psi}},{{\varphi}}) = \Phi\Big({{{\beta_{1}}}}({{\psi}},{{\varphi}}),D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}({{\psi}},{{\varphi}}) ,D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}({{\psi}},{{\varphi}}),D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}({{\psi}},{{\varphi}})\Big),\ 0 < {p},\varphi\leq1, \end{equation}$

(3.7)

having the IC's

$\begin{equation} {{{\beta_{1}}}}^{(k)}({{\psi}},0) = h_{k},\ k = 0,1,2,\cdots,m-1, \end{equation}$

(3.8)

the function ${{{\beta_{1}}}}({{\psi}}, {{\varphi}})$ is unknown, while $\Phi\Big({{{\beta_{1}}}}({{\psi}}, {{\varphi}}), D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}({{\psi}}, {{\varphi}}), D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}({{\psi}}, {{\varphi}}), D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}({{\psi}}, {{\varphi}})\Big)$ may be a nonlinear operator or linear of ${{{\beta_{1}}}}({{\psi}}, {{\varphi}}), D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}({{\psi}}, {{\varphi}}), D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}({{\psi}}, {{\varphi}})$ and $D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}({{\psi}}, {{\varphi}})$ . Applying the AT to both sides of Eq (3.7) yields the following equation; for convenience, we will denote ${{{\beta_{1}}}}({{\psi}}, {{\varphi}})$ using the symbol ${{{\beta_{1}}}}$ ,

$\begin{equation} {A} [{{{\beta_{1}}}}({{\psi}},{{\varphi}})] = \frac{1}{{s}^{{p}}}\Big(\sum\limits_{k = 0}^{m-1}\frac{{{{\beta_{1}}}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}+{A} \Big[\Phi\Big({{{\beta_{1}}}}({{\psi}},{{\varphi}}),D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}({{\psi}},{{\varphi}}),D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}({{\psi}},{{\varphi}}),D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}({{\psi}},{{\varphi}})\Big)\Big]\Big), \end{equation}$

(3.9)

as a result of using the AIT to solve this problem,

$\begin{equation} {{{\beta_{1}}}}({{\psi}},{{\varphi}}) = {A}^{-1}\Big[\frac{1}{{s}^{{p}}}\Big(\sum\limits_{k = 0}^{m-1}\frac{{{{\beta_{1}}}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}+{A} \Big[\Phi\Big({{{\beta_{1}}}}({{\psi}},{{\varphi}}),D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}({{\psi}},{{\varphi}}),D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}({{\psi}},{{\varphi}}),D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}({{\psi}},{{\varphi}})\Big)\Big]\Big)\Big]. \end{equation}$

(3.10)

The solution obtained by using the iterative Aboodh transform method is represented as an infinite series,

$\begin{equation} {{{\beta_{1}}}}({{\psi}},{{\varphi}}) = \sum\limits_{i = 0}^{\infty}{{{\beta_{1}}}}_{i}. \end{equation}$

(3.11)

Since $\Phi\Big({{{\beta_{1}}}}, D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}, D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}, D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}\Big)$ is either a nonlinear or linear operator, which can be decomposed as follows:

$\begin{equation} \begin{split} &\Phi \Big({{{\beta_{1}}}},D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}},D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}},D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}\Big)\\ = &\Phi \Big({{{\beta_{1}}}}_{0},D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}_{0},D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}_{0},D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}_{0}\Big)\\&+ \sum\limits_{i = 0}^{\infty} \left( \Phi \Big(\sum\limits_{k = 0}^{i}\Big({{{\beta_{1}}}}_{k},D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}_{k},D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}_{k},D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}_{k}\Big)\Big) -\Phi \Big(\sum\limits_{k = 1}^{i-1}\Big({{{\beta_{1}}}}_{k},D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}_{k},D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}_{k},D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}_{k}\Big)\Big)\right). \end{split} \end{equation}$

(3.12)

Equation (3.10) is changed to the following equation by substituting the values of (3.12) and (3.11).

$\begin{equation} \begin{split} \sum\limits_{i = 0}^{\infty}{{{\beta_{1}}}}_i({{\psi}},{{\varphi}}) = &{A}^{-1}\Big[ \frac{1}{{s}^{{p}}}\Big( \sum\limits_{k = 0}^{m-1}\frac{{{{\beta_{1}}}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}+{A}[\Phi ({{{\beta_{1}}}}_{0},D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}_{0},D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}_{0},D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}_{0})]\Big)\Big]\\&+{A}^{-1}\Big[ \frac{1}{{s}^{{p}}}\Big({A}\Big[ \sum\limits_{i = 0}^{\infty}\Big(\Phi \sum\limits_{k = 0}^{i}({{{\beta_{1}}}}_{k},D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}_{k},D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}_{k},D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}_{k})\Big)\Big]\Big)\Big]\\&-{A}^{-1}\Big[ \frac{1}{{s}^{{p}}}\Big({A}\Big[\Big(\Phi \sum\limits_{k = 1}^{i-1}({{{\beta_{1}}}}_{k},D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}_{k},D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}_{k},D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}_{k})\Big)\Big]\Big)\Big], \end{split} \end{equation}$

(3.13)

$\begin{equation} \begin{split} &{{{\beta_{1}}}}_0({{\psi}},{{\varphi}}) = {A}^{-1}\Big[ \frac{1}{{s}^{{p}}}\Big( \sum\limits_{k = 0}^{m-1}\frac{{{{\beta_{1}}}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}\Big)\Big],\\ &{{{\beta_{1}}}}_1({{\psi}},{{\varphi}}) = {A}^{-1}\Big[ \frac{1}{{s}^{{p}}}\Big( {A}[\Phi ({{{\beta_{1}}}}_{0},D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}_{0},D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}_{0},D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}_{0})]\Big)\Big],\\ \vdots\\ {{{\beta_{1}}}}_{m+1}({{\psi}},{{\varphi}}) = &{A}^{-1}\Big[ \frac{1}{{s}^{{p}}}\Big({A}\Big[ \sum\limits_{i = 0}^{\infty}\Big(\Phi \sum\limits_{k = 0}^{i}({{{\beta_{1}}}}_{k},D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}_{k},D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}_{k},D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}_{k})\Big)\Big]\Big)\Big]\\&-{A}^{-1}\Big[ \frac{1}{{s}^{{p}}}\Big({A}\Big[\Big(\Phi \sum\limits_{k = 1}^{i-1}({{{\beta_{1}}}}_{k},D_{{{\psi}}}^{\varphi}{{{\beta_{1}}}}_{k},D_{{{\psi}}}^{2\varphi}{{{\beta_{1}}}}_{k},D_{{{\psi}}}^{3\varphi}{{{\beta_{1}}}}_{k})\Big)\Big]\Big)\Big],\ m = 1,2,\cdots. \end{split} \end{equation}$

(3.14)

The m-term of Eq (3.7) may be analytically approximated using the following expression:

$\begin{equation} {{{\beta_{1}}}}({{\psi}},{{\varphi}}) = \sum\limits_{i = 0}^{m-1}{{{\beta_{1}}}}_{i}. \end{equation}$

(3.15)

4. Problem 1

4.1. Problem 1 with ARPSM

Examine the time-fractional KS model with sensitivity term $\varpi({\beta_{2}}) = 1$ , as shown in []. Then, $\frac{\partial}{\partial \psi}\Big({\beta_{1}}(\psi, \varphi) \frac{\partial\varpi({\beta_{2}})}{\partial \psi}\Big) = 0$ ,

$\begin{equation} \begin{split} &D_{{\varphi}}^{p} {\beta_{1}}({\psi}, {\varphi})-a\frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2} = 0,\\ &D_{{\varphi}}^{p} {\beta_{2}}({\psi}, {\varphi})-b\frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}-c{\beta_{1}}({\psi}, {\varphi})+d{\beta_{2}}({\psi}, {\varphi}) = 0, \ \ \text { where }\ \ 0 < p \leq 1, \end{split} \end{equation}$

(4.1)

having IC's:

$\begin{equation} \begin{split} {\beta_{1}}({\psi}, 0) = {l_1} e^{-\psi^2},\\ {\beta_{2}}({\psi}, 0) = {l_2} e^{-\psi^2}. \end{split} \end{equation}$

(4.2)

Using Eq (4.2), AT is applied to Eq (4.1) in order to get

$\begin{equation} \begin{split} &{\beta_{1}}({\psi}, s)-\frac{{l_1} e^{-\psi^2}}{s^2}-\frac{a}{s^p}\Big[\frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2}\Big] = 0,\\ &{\beta_{2}}({\psi}, s)-\frac{{l_2} e^{-\psi^2}}{s^2}-\frac{b}{s^p}\Big[\frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}\Big]-\frac{c}{s^p}\Big[{\beta_{1}}({\psi}, {\varphi})\Big]+\frac{d}{s^p}\Big[{\beta_{2}}({\psi}, {\varphi})\Big] = 0. \end{split} \end{equation}$

(4.3)

The $k^{th}$ -truncated term series are

$\begin{equation} \begin{split} &{\beta_{1}}({{\psi}},s) = \frac{{l_1} e^{-\psi^2}}{s^2}+\sum\limits_{r = 1}^{k}\frac{f_{r}({{\psi}},s)}{s^{rp+1}},\\ &{\beta_{2}}({{\psi}},s) = \frac{{l_2} e^{-\psi^2}}{s^2}+\sum\limits_{r = 1}^{k}\frac{j_{r}({{\psi}},s)}{s^{rp+1}}, \ \ r = 1,2,3,4\cdots \end{split} \end{equation}$

(4.4)

Aboodh residual functions (ARFs) are

$\begin{equation} \begin{split} &{A}_{{{\varphi}}}Res({{\psi}},s) = {\beta_{1}}({\psi}, s)-\frac{{l_1} e^{-\psi^2}}{s^2}-\frac{a}{s^p}\Big[\frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2}\Big] = 0,\\ &{A}_{{{\varphi}}}Res({{\psi}},s) = {\beta_{2}}({\psi}, s)-\frac{{l_2} e^{-\psi^2}}{s^2}-\frac{b}{s^p}\Big[\frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}\Big]-\frac{c}{s^p}\Big[{\beta_{1}}({\psi}, {\varphi})\Big]+\frac{d}{s^p}\Big[{\beta_{2}}({\psi}, {\varphi})\Big] = 0, \end{split} \end{equation}$

(4.5)

and the ${k}^{th}$ -LRFs as:

$\begin{equation} \begin{split} &{A}_{{{\varphi}}}Res_{k}({{\psi}},s) = {\beta_{1}}_{k}({\psi}, s)-\frac{{l_1} e^{-\psi^2}}{s^2}-\frac{a}{s^p}\Big[\frac{\partial ^2{\beta_{1}}_{k}(\psi,\varphi)}{\partial \psi^2}\Big] = 0,\\ &{A}_{{{\varphi}}}Res_{k}({{\psi}},s) = {\beta_{2}}_{k}({\psi}, s)-\frac{{l_2} e^{-\psi^2}}{s^2}-\frac{b}{s^p}\Big[\frac{\partial ^2{\beta_{2}}_{k}({\psi}, {\varphi})}{\partial \psi^2}\Big]-\frac{c}{s^p}\Big[{\beta_{1}}_{k}({\psi}, {\varphi})\Big]+\frac{d}{s^p}\Big[{\beta_{2}}_{k}({\psi}, {\varphi})\Big] = 0, \end{split} \end{equation}$

(4.6)

To determine $f_{r}({{\psi}}, s)$ and $j_{r}({{\psi}}, s)$ , for $r = 1, 2, 3, ...$ . Then, we iteratively solve $\lim_{s\rightarrow \infty}(s^{r{p}+1})$ by multiplying the resulting equation by $s^{r{p}+1}$ , substituting the $r^{th}$ -Aboodh residual function Eq (4.6) for the $r^{th}$ -truncated series Eq (4.4). ${A}_{{{\varphi}}}Res_{{\beta_{1}}, r}({{\psi}}, s)) = 0$ and ${A}_{{{\varphi}}}Res_{{\beta_{2}}, r}({{\psi}}, s)) = 0$ , and $r = 1, 2, 3, \cdots$ . Putting $a = 0.5, b = 3, c = 1$ and $d = 0.8$ and taking the values of $l_1 = 160$ and $l_2 = 120$ , we find the first few terms as:

$\begin{equation} \begin{split} & f_{1}({{\psi}},s) = e^{-\psi^2} \left(320 {\beta_{2}}^2-160\right),\\ & j_{1}({{\psi}},s) = e^{-\psi^2} \left(1440 \psi^2-656\right), \end{split} \end{equation}$

(4.7)

$\begin{equation} \begin{split} & f_{2}({{\psi}},s) = e^{-\psi^2} \left(640 \psi^4-1920 \psi^2+480\right),\\ & j_{2}({{\psi}},s) = e^{-\psi^2} \left(17280 \psi^4-51904 \psi^2+12941\right), \end{split} \end{equation}$

(4.8)

and so on.

Putting $f_{r}({{\psi}}, s)$ , for $r = 1, 2, 3, \cdots,$ in Eq (4.4), we get

$\begin{equation} \begin{split} &{\beta_{1}}({\psi},s) = \frac{e^{-\psi^2} \left(320 \psi^2-160\right)}{s^{p +1}}+\frac{e^{-\psi^2} \left(640 \psi^4-1920 \psi^2+480\right)}{s^{2 p +1}}+\frac{160 e^{-\psi^2}}{s^2}+\cdots,\\ &{\beta_{2}}({\psi},s) = \frac{e^{-\psi^2} \left(1440 \psi^2-656\right)}{s^{p +1}}+\frac{e^{-\psi^2} \left(17280 \psi^4-51904 \psi^2+12941\right)}{s^{2 p +1}}+\frac{120 e^{-\psi^2}}{s^2}+\cdots. \end{split} \end{equation}$

(4.9)

The AIT may be used to get

$\begin{equation} \begin{split} &{\beta_{1}}({\psi},{\varphi}) = \frac{e^{-\psi^2}\varphi^{p} \left(320 \psi^2-160\right)}{\Gamma(p+1)}+\frac{e^{-\psi^2}\varphi^{2p} \left(640 \psi^4-1920 \psi^2+480\right)}{\Gamma(2p+1)}+160 e^{-\psi^2}+\cdots,\\ &{\beta_{2}}({\psi},s) = \frac{e^{-\psi^2}\varphi^p \left(1440 \psi^2-656\right)}{\Gamma(p+1)}+\frac{e^{-\psi^2}\varphi^{2p} \left(17280 \psi^4-51904 \psi^2+12941\right)}{\Gamma(2p+1)}+120 e^{-\psi^2}+\cdots. \end{split} \end{equation}$

(4.10)

4.2. Problem 1 with ATIM

$\begin{equation} \begin{split} &D_{{\varphi}}^{p} {\beta_{1}}({\psi}, {\varphi}) = a\frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2},\\ &D_{{\varphi}}^{p} {\beta_{2}}({\psi}, {\varphi}) = b\frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}+c{\beta_{1}}({\psi}, {\varphi})-d{\beta_{2}}({\psi}, {\varphi}), \ \ \text { where }\ \ 0 < p \leq 1, \end{split} \end{equation}$

(4.11)

having IC's:

$\begin{equation} \begin{split} {\beta_{1}}({\psi}, 0) = {l_1} e^{-\psi^2},\\ {\beta_{2}}({\psi}, 0) = {l_2} e^{-\psi^2}, \end{split} \end{equation}$

(4.12)

By using the AT on each side of Eq (4.11), we are able to get the following result:

$\begin{equation} \begin{split} &{A} [D_{{\varphi}}^{p} {\beta_{1}}({\psi}, {\varphi})] = \frac{1}{{s}^{{p}}}\Big(\sum\limits_{k = 0}^{m-1}\frac{{\beta_{1}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}+{A} \Big[a\frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2}\Big]\Big),\\ &{A} [D_{{\varphi}}^{p} {\beta_{2}}({\psi}, {\varphi})] = \frac{1}{{s}^{{p}}}\Big(\sum\limits_{k = 0}^{m-1}\frac{{\beta_{2}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}+{A} \Big[b\frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}+c{\beta_{1}}({\psi}, {\varphi})-d{\beta_{2}}({\psi}, {\varphi})\Big]\Big), \end{split} \end{equation}$

(4.13)

using the AIT on each side of 4.13, we get the following result:

$\begin{equation} \begin{split} &{\beta_{1}}({\psi}, {\varphi}) = {A}^{-1}\Big[\frac{1}{{s}^{{p}}}\Big(\sum\limits_{k = 0}^{m-1}\frac{{\beta_{1}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}+{A} \Big[a\frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2}\Big]\Big)\Big],\\ &{\beta_{2}}({\psi}, {\varphi}) = {A}^{-1}\Big[\frac{1}{{s}^{{p}}}\Big(\sum\limits_{k = 0}^{m-1}\frac{{\beta_{2}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}+{A} \Big[b\frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}+c{\beta_{1}}({\psi}, {\varphi})-d{\beta_{2}}({\psi}, {\varphi})\Big]\Big)\Big]. \end{split} \end{equation}$

(4.14)

The equation that is produced as a consequence of applying the AT in an iterative manner is as follows:

$\begin{equation*} \begin{split} {\beta_{1}}_0({\psi}, {\varphi}) = &{A}^{-1}\Big[\frac{1}{{s}^{{p}}}\Big(\sum\limits_{k = 0}^{m-1}\frac{{\beta_{1}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}\Big)\Big]\\& = {A}^{-1}\Big[\frac{{\beta_{1}}({{\psi}},0)}{{s}^2}\Big]\\& = {l_1} e^{-\psi^2},\\ \end{split} \end{equation*}$

$\begin{equation*} \begin{split} {\beta_{2}}_0({\psi}, {\varphi}) = &{A}^{-1}\Big[\frac{1}{{s}^{{p}}}\Big(\sum\limits_{k = 0}^{m-1}\frac{{\beta_{2}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}\Big)\Big]\\& = {A}^{-1}\Big[\frac{{\beta_{2}}({{\psi}},0)}{{s}^2}\Big]\\& = {l_2} e^{-\psi^2}. \end{split} \end{equation*}$

We replaced the RL integral in Eq (4.11) to get the equivalent variant.

$\begin{equation} \begin{split} &{\beta_{1}}({\psi}, {\varphi}) = {l_1} e^{-\psi^2}-{A}\Big[a\frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2}\Big],\\ &{\beta_{2}}({\psi}, {\varphi}) = {l_2} e^{-\psi^2}-{A}\Big[b\frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}+c{\beta_{1}}({\psi}, {\varphi})-d{\beta_{2}}({\psi}, {\varphi})\Big]. \end{split} \end{equation}$

(4.15)

Putting $a = 0.5, b = 3, c = 1$ , and $d = 0.8$ and taking the values of $l_1 = 160$ and $l_2 = 120$ , the following terms are then acquired by using the ATIM procedure:

$\begin{equation} \begin{split} &{\beta_{1}}_{0}({\psi},{\varphi}) = 160 e^{-\psi^2},\\ &{\beta_{2}}_{0}({\psi},{\varphi}) = 120 e^{-\psi^2},\\ &{\beta_{1}}_{1}({\psi},{\varphi}) = \frac{e^{-\psi^2} \left(320 \psi^2-160\right) \varphi^p}{\Gamma (p+1)},\\ &{\beta_{2}}_{1}({\psi},{\varphi}) = \frac{e^{-\psi^2} \left(1440 \psi^2-656\right) \varphi^p}{\Gamma (p+1)},\\ &{\beta_{1}}_{2}({\psi},{\varphi}) = \frac{e^{-\psi^2} \left(640 \psi^4-1920 \psi^2+478\right) \varphi^{2 p}}{\Gamma (2 p+1)},\\ &{\beta_{2}}_{2}({\psi},{\varphi}) = \frac{e^{-\psi^2} zeta^{2 p} \left(\left(364.8\, -832 \psi^2\right) \Gamma (p+1)+p \left(17280 \psi^4-51072 \psi^2+12576\right) \Gamma (p)\right)}{\Gamma (p+1) \Gamma (2 p+1)}. \end{split} \end{equation}$

(4.16)

The following is the final ATIM solution:

$\begin{equation} \begin{split} &{\beta_{1}}({\psi},{\varphi}) = {\beta_{1}}_{0}({\psi},{\varphi})+{\beta_{1}}_{1}({\psi},{\varphi})+{\beta_{1}}_{2}({\psi},{\varphi})+\cdots,\\ &{\beta_{2}}({\psi},{\varphi}) = {\beta_{2}}_{0}({\psi},{\varphi})+{\beta_{2}}_{1}({\psi},{\varphi})+{\beta_{2}}_{2}({\psi},{\varphi})+\cdots. \end{split} \end{equation}$

(4.17)

$\begin{equation} \begin{split} {\beta_{1}}({\psi},{\varphi}) = &160 e^{-\psi^2}+\frac{e^{-\psi^2} \left(320 \psi^2-160\right) \varphi^p}{\Gamma (p+1)}+\frac{e^{-\psi^2} \left(640 \psi^4-1920 \psi^2+478\right) \varphi^{2 p}}{\Gamma (2 p+1)}+\cdots,\\ {\beta_{2}}({\psi},{\varphi}) = &120 e^{-\psi^2}+\frac{e^{-\psi^2} \left(1440 \psi^2-656\right) \varphi^p}{\Gamma (p+1)}\\&+\frac{e^{-\psi^2} zeta^{2 p} \left(\left(364.8\, -832 \psi^2\right) \Gamma (p+1)+p \left(17280 \psi^4-51072 \psi^2+12576\right) \Gamma (p)\right)}{\Gamma (p+1) \Gamma (2 p+1)}+\cdots. \end{split} \end{equation}$

(4.18)

5. Problem 2

5.1. Problem 2 with ARPSM

Examine the KS model of fractional order as stated in [] with sensitivity term $\varpi({\beta_{2}}) = {\beta_{2}}(\psi, \varphi)$ . Then, the function $\frac{\partial}{\partial \psi}\Big({\beta_{1}}(\psi, \varphi) \frac{\partial\varpi({\beta_{2}})}{\partial \psi}\Big) = {\beta_{1}}(\psi, \varphi) \frac{\partial ^2{\beta_{2}}(\psi, \varphi)}{\partial \psi^2}+\frac{\partial {\beta_{1}}(\psi, \varphi)}{\partial \psi} \frac{\partial {\beta_{2}}(\psi, \varphi)}{\partial \psi}$ ,

$\begin{equation} \begin{split} &D_{{\varphi}}^{p} {\beta_{1}}({\psi}, {\varphi})-a \frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2}+{\beta_{1}}(\psi,\varphi) \frac{\partial ^2{\beta_{2}}(\psi,\varphi)}{\partial \psi^2}+\frac{\partial {\beta_{1}}(\psi,\varphi)}{\partial \psi} \frac{\partial {\beta_{2}}(\psi,\varphi)}{\partial \psi} = 0,\\ &D_{{\varphi}}^{p} {\beta_{2}}({\psi}, {\varphi})-b \frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}-c {\beta_{1}}({\psi}, {\varphi})+d {\beta_{2}}({\psi}, {\varphi}) = 0, \ \ \text { where }\ \ 0 < p \leq 1, \end{split} \end{equation}$

(5.1)

having IC's:

$\begin{equation} \begin{split} {\beta_{1}}({\psi}, 0) = {l_1} e^{-\psi^2},\\ {\beta_{2}}({\psi}, 0) = {l_2} e^{-\psi^2}. \end{split} \end{equation}$

(5.2)

AT is applied to Eq (5.1), the following results are obtained using Eq (5.2):

$\begin{equation} \begin{split} &{\beta_{1}}({\psi}, s)-\frac{{l_1} e^{-\psi^2}}{s^2}-\frac{a}{s^p}\Big[ \frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2}\Big]+\frac{1}{s^p}A_{{\varphi}}\Big[A_{{\varphi}}^{-1}{\beta_{1}}(\psi,\varphi)\times A_{{\varphi}}^{-1}\frac{\partial ^2{\beta_{2}}(\psi,\varphi)}{\partial \psi^2}\Big]\\&+\frac{1}{s^p}A_{{\varphi}}\Big[\frac{\partial A_{{\varphi}}^{-1}{\beta_{1}}(\psi,\varphi)}{\partial \psi}\times \frac{\partial A_{{\varphi}}^{-1}{\beta_{2}}(\psi,\varphi)}{\partial \psi}\Big] = 0, \\ &{\beta_{2}}({\psi}, s)-\frac{{l_2} e^{-\psi^2}}{s^2}-\frac{b}{s^p}\Big[\frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}\Big]-\frac{c}{s^p}\Big[{\beta_{1}}({\psi}, {\varphi})\Big]+\frac{d}{s^p}\Big[{\beta_{2}}({\psi}, {\varphi})\Big] = 0. \end{split} \end{equation}$

(5.3)

The $k^{th}$ truncated term series are

(5.4)

Aboodh residual functions (ARFs) are

$\begin{equation} \begin{split} {A}_{{{\varphi}}}Res({{\psi}},s) = &{\beta_{1}}({\psi}, s)-\frac{{l_1} e^{-\psi^2}}{s^2}-\frac{a}{s^p}\Big[ \frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2}\Big]+\frac{1}{s^p}A_{{\varphi}}\Big[A_{{\varphi}}^{-1}{\beta_{1}}(\psi,\varphi)\times A_{{\varphi}}^{-1}\frac{\partial ^2{\beta_{2}}(\psi,\varphi)}{\partial \psi^2}\Big]\\&+\frac{1}{s^p}A_{{\varphi}}\Big[\frac{\partial A_{{\varphi}}^{-1}{\beta_{1}}(\psi,\varphi)}{\partial \psi}\times \frac{\partial A_{{\varphi}}^{-1}{\beta_{2}}(\psi,\varphi)}{\partial \psi}\Big] = 0,\\ {A}_{{{\varphi}}}Res({{\psi}},s) = &{\beta_{2}}({\psi}, s)-\frac{{l_2} e^{-\psi^2}}{s^2}-\frac{b}{s^p}\Big[\frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}\Big]-\frac{c}{s^p}\Big[{\beta_{1}}({\psi}, {\varphi})\Big]+\frac{d}{s^p}\Big[{\beta_{2}}({\psi}, {\varphi})\Big] = 0, \end{split} \end{equation}$

(5.5)

and the ${k}^{th}$ -LRFs as:

$\begin{equation} \begin{split} {A}_{{{\varphi}}}Res_{k}({{\psi}},s) = &{\beta_{1}}_{k}({\psi}, s)+\frac{{l_1} e^{-\psi^2}}{s^2}-\frac{a}{s^p}\Big[ \frac{\partial ^2{\beta_{1}}_{k}(\psi,\varphi)}{\partial \psi^2}\Big]+\frac{1}{s^p}A_{{\varphi}}\Big[A_{{\varphi}}^{-1}{\beta_{1}}_{k}(\psi,\varphi)\times A_{{\varphi}}^{-1}\frac{\partial ^2{\beta_{2}}_{k}(\psi,\varphi)}{\partial \psi^2}\Big]\\&+\frac{1}{s^p}A_{{\varphi}}\Big[\frac{\partial A_{{\varphi}}^{-1}{\beta_{1}}_{k}(\psi,\varphi)}{\partial \psi}\times \frac{\partial A_{{\varphi}}^{-1}{\beta_{2}}_{k}(\psi,\varphi)}{\partial \psi}\Big] = 0,\\ {A}_{{{\varphi}}}Res_{k}({{\psi}},s) = &{\beta_{2}}_{k}({\psi}, s)-\frac{{l_2} e^{-\psi^2}}{s^2}-\frac{b}{s^p}\Big[\frac{\partial ^2{\beta_{2}}_{k}({\psi}, {\varphi})}{\partial \psi^2}\Big]-\frac{c}{s^p}\Big[{\beta_{1}}_{k}({\psi}, {\varphi})\Big]+\frac{d}{s^p}\Big[{\beta_{2}}_{k}({\psi}, {\varphi})\Big] = 0. \end{split} \end{equation}$

(5.6)

To determine $f_{r}({{\psi}}, s)$ and $j_{r}({{\psi}}, s)$ , for $r = 1, 2, 3, ...$ . Then, we iteratively solve $\lim_{s\rightarrow \infty}(s^{r{p}+1})$ by multiplying the resulting equation by $s^{r{p}+1}$ , substituting the $r^{th}$ -Aboodh residual function Eq (5.6) for the $r^{th}$ -truncated series Eq (5.4). ${A}_{{{\varphi}}}Res_{{\beta_{1}}, r}({{\psi}}, s)) = 0$ and ${A}_{{{\varphi}}}Res_{{\beta_{2}}, r}({{\psi}}, s)) = 0$ , and $r = 1, 2, 3, \cdots$ . Putting $a = 0.5, b = 3, c = 1,$ and $d = 0.8$ and taking the values of $l_1 = 160$ and $l_2 = 120$ , we find the first few terms as:

$\begin{equation} \begin{split} & f_{1}({{\psi}},s) = e^{-\psi^2} \left(320 \psi^2-160\right)-38400 e^{-2 \psi^2},\\ & j_{1}({{\psi}},s) = e^{-\psi^2} \left(1440 \psi^2-656\right), \end{split} \end{equation}$

(5.7)

$\begin{equation} \begin{split} f_{2}({{\psi}},s) = &e^{-3 \psi^2} \left(18432000 \psi^2+9216000\right)+e^{-2 \psi^2} \left(785920-1612800 \psi^2\right)\\&+e^{-\psi^2} \left(640 \psi^4-1920 \psi^2+480\right),\\ j_{2}({{\psi}},s) = &e^{-\psi^2} \left(17280 \psi^4-51904 \psi^2+12941\right)-38400 e^{-2 \psi^2}, \end{split} \end{equation}$

(5.8)

and so on.

Equation (5.4) is used to obtain $f_{r}({{\psi}}, s)$ for $r = 1, 2, 3, \cdots,$

$\begin{equation} \begin{split} &{\beta_{1}}({\psi},s) = \frac{160 e^{-\psi^2}}{s^2}+\frac{e^{-\psi^2} \left(320 \psi^2-160\right)-38400 e^{-2 \psi^2}}{s^{p+1}}\\&+ \frac{e^{-3 \psi^2} \left(18432000 \psi^2+9216000\right)+e^{-2 \psi^2} \left(785920-1612800 \psi^2\right)+e^{-\psi^2} \left(640 \psi^4-1920 \psi^2+480\right)}{s^{2 p+1}}+\cdots,\\ &{\beta_{2}}({\psi},s) = \frac{120 e^{-\psi^2}}{s^2}+\frac{e^{-\psi^2} \left(1440 \psi^2-656\right)}{s^{p+1}}+\frac{e^{-\psi^2} \left(17280 \psi^4-51904 \psi^2+12941\right)-38400 e^{-2 \psi^2}}{s^{2 p+1}}+\cdots. \end{split} \end{equation}$

(5.9)

AIT is applied to get

$\begin{equation} \begin{split} &{\beta_{1}}({\psi},\varphi) = 160 e^{-\psi^2}+\frac{e^{-\psi^2}\varphi^p \left(320 \psi^2-160\right)-38400 e^{-2 \psi^2}}{\Gamma(p+1)}\\&+ \frac{e^{-3 \psi^2}\varphi^{2p} \left(18432000 \psi^2+9216000\right)+e^{-2 \psi^2} \left(785920-1612800 \psi^2\right)+e^{-\psi^2} \left(640 \psi^4-1920 \psi^2+480\right)}{\Gamma(2p+1)}+\cdots,\\ &{\beta_{2}}({\psi},\varphi) = \frac{120 e^{-\psi^2}}{s^2}+\frac{e^{-\psi^2}\varphi^p \left(1440 \psi^2-656\right)}{\Gamma(p+1)}+\frac{e^{-\psi^2}\varphi^{2p} \left(17280 \psi^4-51904 \psi^2+12941\right)-38400 e^{-2 \psi^2}}{\Gamma(2p+1)}+\cdots. \end{split} \end{equation}$

(5.10)

5.2. Problem 2 with ATIM

$\begin{equation} \begin{split} &D_{{\varphi}}^{p} {\beta_{1}}({\psi}, {\varphi}) = a \frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2}-{\beta_{1}}(\psi,\varphi) \frac{\partial ^2{\beta_{2}}(\psi,\varphi)}{\partial \psi}-\frac{\partial {\beta_{1}}(\psi,\varphi)}{\partial \psi} \frac{\partial {\beta_{2}}(\psi,\varphi)}{\partial \psi},\\ &D_{{\varphi}}^{p} {\beta_{2}}({\psi}, {\varphi}) = b \frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}+c {\beta_{1}}({\psi}, {\varphi})-d {\beta_{2}}({\psi}, {\varphi}), \ \ \text { where }\ \ 0 < p \leq 1, \end{split} \end{equation}$

(5.11)

having IC's:

$\begin{equation} \begin{split} {\beta_{1}}({\psi}, 0) = {l_1} e^{-\psi^2},\\ {\beta_{2}}({\psi}, 0) = {l_2} e^{-\psi^2}, \end{split} \end{equation}$

(5.12)

when applying the AT to both sides of Eq (5.11), we get the following result:

$\begin{equation} \begin{split} &{A} [D_{{\varphi}}^{p} {\beta_{1}}({\psi},{\varphi})] = \frac{1}{{s}^{{p}}}\Big(\sum\limits_{k = 0}^{m-1}\frac{{\beta_{1}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}+{A} \Big[a \frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2}-{\beta_{1}}(\psi,\varphi) \frac{\partial ^2{\beta_{2}}(\psi,\varphi)}{\partial \psi}-\frac{\partial {\beta_{1}}(\psi,\varphi)}{\partial \psi} \frac{\partial {\beta_{2}}(\psi,\varphi)}{\partial \psi}\Big]\Big),\\ &{A} [D_{{\varphi}}^{p} {\beta_{2}}({\psi}, {\varphi})] = \frac{1}{{s}^{{p}}}\Big(\sum\limits_{k = 0}^{m-1}\frac{{\beta_{2}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}+{A} \Big[b\frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}+c{\beta_{1}}({\psi}, {\varphi})-d{\beta_{2}}({\psi}, {\varphi})\Big]\Big), \end{split} \end{equation}$

(5.13)

applying the AIT to both sides of Eq (5.13) yields the following result:

$\begin{equation} \begin{split} &{\beta_{1}}({\psi}, {\varphi}) = {A}^{-1}\Big[\frac{1}{{s}^{{p}}}\Big(\sum\limits_{k = 0}^{m-1}\frac{{\beta_{1}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}+{A} \Big[a \frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2}-{\beta_{1}}(\psi,\varphi) \frac{\partial ^2{\beta_{2}}(\psi,\varphi)}{\partial \psi}-\frac{\partial {\beta_{1}}(\psi,\varphi)}{\partial \psi} \frac{\partial {\beta_{2}}(\psi,\varphi)}{\partial \psi}\Big]\Big)\Big],\\ &{\beta_{2}}({\psi}, {\varphi}) = {A}^{-1}\Big[\frac{1}{{s}^{{p}}}\Big(\sum\limits_{k = 0}^{m-1}\frac{{\beta_{2}}^{(k)}({{\psi}},0)}{{s}^{2-{p}+k}}+{A} \Big[b\frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}+c{\beta_{1}}({\psi}, {\varphi})-d{\beta_{2}}({\psi}, {\varphi})\Big]\Big)\Big]. \end{split} \end{equation}$

(5.14)

This equation is obtained by using the AT's iterative procedure:

The RL integral is applied to Eq (5.1) to yield the equivalent form.

$\begin{equation} \begin{split} &{\beta_{1}}({\psi},{\varphi}) = {l_1} e^{-\psi^2}-{A}\Big[a \frac{\partial ^2{\beta_{1}}(\psi,\varphi)}{\partial \psi^2}-{\beta_{1}}(\psi,\varphi) \frac{\partial ^2{\beta_{2}}(\psi,\varphi)}{\partial \psi}-\frac{\partial {\beta_{1}}(\psi,\varphi)}{\partial \psi} \frac{\partial {\beta_{2}}(\psi,\varphi)}{\partial \psi}\Big],\\ &{\beta_{2}}({\psi}, {\varphi}) = {l_2} e^{-\psi^2}-{A}\Big[b\frac{\partial ^2{\beta_{2}}({\psi}, {\varphi})}{\partial \psi^2}+c{\beta_{1}}({\psi}, {\varphi})-d{\beta_{2}}({\psi}, {\varphi})\Big]. \end{split} \end{equation}$

(5.15)

Putting $a = 0.5, b = 3, c = 1,$ and $d = 0.8$ and taking the values of $l_1 = 160$ and $l_2 = 120$ , the following terms are then acquired by using the ATIM procedure: These terms are obtained using the ATIM process,

$\begin{align} {\beta_{1}}_{0}({\psi},{\varphi}) = &160 e^{-\psi^2},\quad\quad {\beta_{2}}_{0}({\psi},{\varphi}) = 120 e^{-\psi^2},\\ {\beta_{1}}_{1}({\psi},{\varphi}) = &\frac{e^{-2 {\beta_{2}}^2} \left(e^{{\beta_{2}}^2} \left(320 {\beta_{2}}^2-160\right)-38400\right) \varphi^p}{\Gamma (p+1)},\quad\quad {\beta_{2}}_{1}({\psi},{\varphi}) = \frac{e^{-x^2} \left(1440 x^2-656\right) t^p}{\Gamma (p+1)},\\ {\beta_{1}}_{2}({\psi},{\varphi}) = &-\Big(225\ 4^{6-p} e^{-2 \psi ^2} \varphi ^p \Big(p \Gamma (p) \Big(4^p \Gamma \Big(p+\frac{1}{2}\Big) \Big(e^{-\psi ^2} \Big(-240 \psi ^4-10.66 \psi ^2\\&+e^{\psi ^2} \Big(\psi ^4-0.044 \psi ^2+0.727\Big)+174.66\Big) \varphi ^{2 p} \Gamma (2 p+1)+0.0416 \Gamma (p+1) \Gamma (3 p+1)\Big)\\&+\varphi ^p \Gamma (p+1) \Gamma (3 p+1) \Big(3.10179 \psi ^2+\Big(-0.0012 \psi ^4\\&+35.4528 \psi ^2+17.7236\Big) \sinh \Big(\psi ^2\Big)+\Big(-0.00123 \psi ^4-35.4453 \psi ^2-17.7254\Big) \cosh \Big(\psi ^2\Big)\\&-1.5115\Big)\Big)-0.0738 \Gamma (p+1) \Gamma (2 p+1) \Gamma (3 p+1)\Big)\Big)/\Big(p \Gamma (p) \Gamma \Big(p+\frac{1}{2}\Big) \Gamma (p+1)^2 \Gamma (3 p+1)\Big),\\ {\beta_{2}}_{2}({\psi},{\varphi}) = &\Big(e^{-2 \psi ^2} \varphi ^{2 p} \Big(e^{\psi ^2} \Big(\Big(524.8\, -1152 \psi ^2\Big) \Gamma (p+1)\\&+p \Big(17280 \psi ^4-50752 \psi ^2+12416\Big) \Gamma (p)\Big)-38400 p \Gamma (p)\Big)\Big)/\Big(\Gamma (p+1) \Gamma (2 p+1)\Big). \end{align}$

(5.16)

The following is the ATIM procedure's ultimate solution:

(5.17)

$\begin{equation} \begin{split} {\beta_{1}}({\psi},{\varphi}) = &160 e^{-\psi^2}+\Big(e^{-2 {\beta_{2}}^2} \Big(e^{{\beta_{2}}^2} \Big(320 {\beta_{2}}^2-160\Big)-38400\Big) \varphi^p\Big)/\Big(\Gamma (p+1)\Big)\\&-\Big(225\ 4^{6-p} e^{-2 \psi ^2} \varphi ^p\times \Big(p \Gamma (p) \Big(4^p \Gamma \Big(p+\frac{1}{2}\Big) \Big(e^{-\psi ^2} \Big(-240 \psi ^4-10.66 \psi ^2\\&+e^{\psi ^2} \Big(\psi ^4-0.044 \psi ^2+0.727\Big)+174.66\Big) \varphi ^{2 p} \Gamma (2 p+1)\\&+0.0416 \Gamma (p+1) \Gamma (3 p+1)\Big)+\varphi ^p \Gamma (p+1) \Gamma (3 p+1) \Big(3.10179 \psi ^2\\&+\Big(-0.0012 \psi ^4+35.4528 \psi ^2+17.7236\Big) \sinh \Big(\psi ^2\Big)\\&+\Big(-0.00123 \psi ^4-35.4453 \psi ^2-17.7254\Big) \cosh \Big(\psi ^2\Big)-1.5115\Big)\Big)\\&-0.0738 \Gamma (p+1) \Gamma (2 p+1) \Gamma (3 p+1)\Big)\Big)/\Big(p \Gamma (p) \Gamma \Big(p+\frac{1}{2}\Big) \Gamma (p+1)^2 \Gamma (3 p+1)\Big)+\cdots, \end{split} \end{equation}$

(5.18)

$\begin{equation} \begin{split} {\beta_{2}}({\psi},{\varphi}) = &120 e^{-\psi^2}+\Big(e^{-x^2} \left(1440 x^2-656\right) t^p\Big)/\Big(\Gamma (p+1)\Big)+\Big(e^{-2 \psi ^2} \varphi ^{2 p} \Big(e^{\psi ^2} \Big(\Big(524.8\, -1152 \psi ^2\Big) \Gamma (p+1)\\&+p \Big(17280 \psi ^4-50752 \psi ^2+12416\Big) \Gamma (p)\Big)-38400 p \Gamma (p)\Big)\Big)/\Big(\Gamma (p+1) \Gamma (2 p+1)\Big)+\cdots. \end{split} \end{equation}$

(5.19)

6. Graphics and tables discussion

In Problem 1, we embark on a comprehensive exploration of the solutions ${\beta_{1}}(\psi, \varphi)$ and ${\beta_{2}}(\psi, \varphi)$ through both graphical and numerical analyses employing two distinct methodologies: the Aboodh residual power series method (ARPSM) and the Aboodh transform iteration method (ATIM). Beginning with ${\beta_{1}}(\psi, \varphi)$ , offers an insightful depiction of the approximate solution obtained via ARPSM for a specific value of $p = 1$ . Building upon this foundation, extends the analysis, providing both 3D and 2D representations to elucidate the influence of varying $p$ on the solution when $\varphi = 0.1$ . Similarly, and delve into the corresponding analyses for ${\beta_{2}}(\psi, \varphi)$ . These visualizations offer a nuanced understanding of how changes in the parameter $p$ affect the behavior of the solutions across different dimensions. In conjunction with the graphical exploration, and complement our investigation by presenting detailed fractional order analyses for ARPSM applied to ${\beta_{1}}(\psi, \varphi)$ and ${\beta_{2}}(\psi, \varphi)$ , respectively. These tables provide valuable insights into the fractional characteristics of the solutions and contribute to a comprehensive understanding of their properties.

Figure 1. Approximate solution of

${\beta_{1}}({\psi}, {\varphi})$ via ARPSM for

$p = 1.$ .

DownLoad: Full-Size Img PowerPoint

Figure 2.

${\beta_{1}}({\psi}, {\varphi})$ , (a) shows three-dimensional analysis of different values of

$p$ ; (b) shows two-dimensional analysis of different values of

$p$ at

$\varphi = 0.1$ via ARPSM.

DownLoad: Full-Size Img PowerPoint

Figure 3. Approximate solution of

${\beta_{2}}({\psi}, {\varphi})$ via ARPSM for

$p = 1$ .

DownLoad: Full-Size Img PowerPoint

Figure 4.

${\beta_{2}}({\psi}, {\varphi})$ : (a) shows three-dimensional analysis of different values of

$p$ ; (b) shows two-dimensional analysis of different values of

$p$ at

$\varphi = 0.1$ via ARPSM.

DownLoad: Full-Size Img PowerPoint

Table 1. Analysis of various values of fractional order of ARPSM of Problem 1 of

$\beta_{1}(\psi, \varphi)$ .

$\varphi$	${p=0.40}$	${p=0.60}$	${p=0.80}$	$p=1.00$
0.2	0.0048028	0.0022338	0.00098417	0.00043033
0.4	0.0082096	0.0049227	0.00273349	0.00144405
0.6	0.0112581	0.0078699	0.00506050	0.00305916
0.8	0.0140966	0.0110044	0.00787829	0.00527565
1.0	0.0167894	0.0142882	0.0111338	0.00809353

| Show Table

DownLoad: CSV

Table 2. Analysis of various values of fractional order of ARPSM of Problem 1 of

$\beta_{2}(\psi, \varphi)$ .

$\varphi$	${p=0.40}$	${p=0.60}$	${p=0.80}$	$p=1.00$
0.2	0.121739	0.054477	0.022375	0.0086336
0.4	0.211321	0.124300	0.066837	0.0334867
0.6	0.291879	0.201620	0.127169	0.0745724
0.8	0.367092	0.284277	0.200911	0.1318910
1.0	0.438569	0.371156	0.286582	0.2054420

| Show Table

DownLoad: CSV

Shifting focus to the ATIM method, and display the approximate solutions of ${\beta_{1}}(\psi, \varphi)$ and ${\beta_{2}}(\psi, \varphi)$ for $p = 1$ , respectively. and further extend the analysis, offering insights into the impact of varying $p$ at $\varphi = 0.1$ . The fractional order sensitivity is examined through and for ${\beta_{1}}(\psi, \varphi)$ and ${\beta_{2}}(\psi, \varphi)$ under ATIM.

Figure 5. Approximate solution of

$\beta_{1}(\psi, \varphi)$ via ATIM for

$p = 1.$ .

DownLoad: Full-Size Img PowerPoint

Figure 6.

$\beta_{1}(\psi, \varphi)$ : (a) shows three-dimensional analysis of different values of

$p$ ; (b) shows two dimensional analysis of different values of

$p$ at

$\varphi = 0.1$ via ATIM.

DownLoad: Full-Size Img PowerPoint

Figure 7. Approximate solution of

$\beta_{2}(\psi, \varphi)$ via ATIM for

$p = 1.$ .

DownLoad: Full-Size Img PowerPoint

Figure 8.

$\beta_{2}(\psi, \varphi)$ : (a) shows three-dimensional analysis of different values of

$p$ ; (b) shows two-dimensional analysis of different values of

$p$ at

$\varphi = 0.1$ via ATIM.

DownLoad: Full-Size Img PowerPoint

Table 3. Analysis of various values of fractional order of ATIM of Problem 1 of

$\beta_{1}(\psi, \varphi)$ .

$\varphi$	${p=0.40}$	${p=0.60}$	${p=0.80}$	$p=1.00$
0.2	0.0048028	0.0022338	0.00098416	0.00043033
0.4	0.0082095	0.0049226	0.00273345	0.00144403
0.6	0.0112579	0.0078698	0.00506043	0.00305912
0.8	0.0140964	0.0110042	0.00787818	0.00527558
1.0	0.0167892	0.0142880	0.01113370	0.00809342

| Show Table

DownLoad: CSV

Table 4. Analysis of various values of fractional order of ATIM of Problem 1 of

$\beta_{2}(\psi, \varphi)$ .

$\varphi$	${p=0.40}$	${p=0.60}$	${p=0.80}$	$p=1.00$
0.2	0.121748	0.054481	0.022376	0.0086342
0.4	0.211336	0.124308	0.066842	0.0334890
0.6	0.291899	0.201634	0.127178	0.0745776
0.8	0.367118	0.284297	0.200924	0.1319000
1.0	0.438600	0.371182	0.286602	0.2054560

| Show Table

DownLoad: CSV

In order to facilitate a comprehensive comparison, and juxtapose the results obtained from both ARPSM and ATIM for ${\beta_{1}}(\psi, \varphi)$ and ${\beta_{2}}(\psi, \varphi)$ in Problem 1. These tables provide a nuanced understanding of the performance of each method, shedding light on their strengths and limitations in solving the given fractional-order equations. The two most basic approaches to solving fractional differential equations are the ATIM and the ARPSM, as stated in ^[52,53] and ^[49,50,51], respectively. These techniques provide numerical solutions to PDEs that do not need discretization or linearization, making the symbolic terms in analytical solutions instantly visible. The primary objective of this study is to compare and contrast the performance of ARPSM and ATIM in solving the Keller-Segel (KS) model. It is worth mentioning that several linear and nonlinear fractional differential problems have been solved using these two approaches.

Table 5. Problem 1: comparison of both methods for

$\beta_{1}(\psi, \varphi)$ .

	$p=1.00$		$p=0.80$		$p=0.60$		$p=0.40$
$\varphi$	ARPSM	ATIM	ARPSM	ATIM	ARPSM	ATIM	ARPSM	ATIM
0.2	0.00043033	0.00043033	0.00098417	0.00098416	0.0022338	0.0022338	0.0048028	0.0048028
0.4	0.00144405	0.00144403	0.00273349	0.00273345	0.0049227	0.0049226	0.0082096	0.0082095
0.6	0.00305916	0.00305912	0.00506050	0.00506043	0.0078699	0.0078698	0.0112581	0.0112579
0.8	0.00527565	0.00527558	0.00787829	0.00787818	0.0110044	0.0110042	0.0140966	0.0140964
1.0	0.00809353	0.00809342	0.01113380	0.01113370	0.0142882	0.0142880	0.0167894	0.0167892

| Show Table

DownLoad: CSV

Table 6. Problem 1: comparison of both methods for

$\beta_{2}(\psi, \varphi)$ .

	$p=1.00$		$p=0.80$		$p=0.60$		$p=0.40$
$\varphi$	ARPSM	ATIM	ARPSM	ATIM	ARPSM	ATIM	ARPSM	ATIM
0.2	0.0086336	0.0086342	0.022375	0.022376	0.054477	0.054481	0.121739	0.121748
0.4	0.0334867	0.0334890	0.066837	0.066842	0.124300	0.124308	0.211321	0.211336
0.6	0.0745724	0.0745776	0.127169	0.127178	0.201620	0.201634	0.291879	0.291899
0.8	0.1318910	0.1319000	0.200911	0.200924	0.284277	0.284297	0.367092	0.367118
1.0	0.2054420	0.2054560	0.286582	0.286602	0.371156	0.371182	0.438569	0.438600

| Show Table

DownLoad: CSV

In Problem 2, the analysis of solutions ${\beta_{1}}(\psi, \varphi)$ and ${\beta_{2}}(\psi, \varphi)$ is carried out using the ARPSM and the ATIM. For ${\beta_{1}}(\psi, \varphi)$ , illustrates the approximate solution via ARPSM for $p = 1$ . Subsequently, presents 3D and 2D analyses, demonstrating the influence of varying $p$ on the solution at $\varphi = 0.1$ . Analogously, and provide the corresponding analyses for ${\beta_{2}}(\psi, \varphi)$ . Complementing the graphical exploration, details the fractional order analysis for ARPSM of ${\beta_{1}}(\psi, \varphi)$ , and does the same for ${\beta_{2}}(\psi, \varphi)$ .

Figure 9. Approximate solution of

$\beta_{1}(\psi, \varphi)$ via ARPSM for

$p = 1.$ .

DownLoad: Full-Size Img PowerPoint

Figure 10. 2D analysis of different values of

$p$ at

$\varphi = 0.1$ .

DownLoad: Full-Size Img PowerPoint

Figure 11. Approximate solution of

$\beta_{2}(\psi, \varphi)$ via ARPSM for

$p = 1.$ .

DownLoad: Full-Size Img PowerPoint

Figure 12. 2D analysis of different values of

$p$ at

$\varphi = 0.1$ .

DownLoad: Full-Size Img PowerPoint

Table 7. Analysis of various values of fractional order of ARPSM of Problem 2 of

$\beta_{1}(\psi, \varphi)$ .

$\varphi$	${p=0.40}$	${p=0.60}$	${p=0.80}$	$p=1.00$
0.2	0.0048027	0.0022338	0.00098415	0.00043032
0.4	0.0082094	0.0049226	0.00273344	0.00144403
0.6	0.0112578	0.0078698	0.00506040	0.00305910
0.8	0.0140963	0.0110042	0.00787813	0.00527555
1.0	0.0167891	0.0142880	0.01113360	0.00809337

| Show Table

DownLoad: CSV

Table 8. Analysis of various values of fractional order of ARPSM of Problem 2 of

$\beta_{2}(\psi, \varphi)$ .

$\varphi$	${p=0.40}$	${p=0.60}$	${p=0.80}$	$p=1.00$
0.2	0.121739	0.054477	0.022375	0.0086336
0.4	0.211321	0.124300	0.066837	0.0334867
0.6	0.291879	0.201620	0.127169	0.0745724
0.8	0.367092	0.284277	0.200911	0.1318910
1.0	0.438569	0.371156	0.286582	0.2054420

| Show Table

DownLoad: CSV

Figure 13. Approximate solution of

$\beta_{1}(\psi, \varphi)$ via ATIM for

$p = 1.$ .

DownLoad: Full-Size Img PowerPoint

Figure 14. 2D analysis of different values of

$p$ at

$\varphi = 0.1$ .

DownLoad: Full-Size Img PowerPoint

Figure 15. Approximate solution of

$\beta_{2}(\psi, \varphi)$ via ATIM for p = 1.

DownLoad: Full-Size Img PowerPoint

Figure 16. 2D analysis of different values of

$p$ at

$\varphi = 0.1$ .

DownLoad: Full-Size Img PowerPoint

Table 9. Analysis of various values of fractional order of ATIM of Problem 2 of

$\beta_{1}(\psi, \varphi)$ .

$\varphi$	${p=0.40}$	${p=0.60}$	${p=0.80}$	$p=1.00$
0.2	0.0047363	0.00220436	0.00097225	0.00042586
0.4	0.0080934	0.00485484	0.00269728	0.00142614
0.6	0.0110971	0.00775928	0.00499109	0.00301877
0.8	0.0138937	0.01084780	0.00776810	0.00520373
1.0	0.0165466	0.01408330	0.01097610	0.00798096

| Show Table

DownLoad: CSV

Table 10. Analysis of various values of fractional order of ATIM of Problem 2 of

$\beta_{2}(\psi, \varphi)$ .

$\varphi$	${p=0.40}$	${p=0.60}$	${p=0.80}$	$p=1.00$
0.2	0.121748	0.054481	0.0223765	0.0086342
0.4	0.211336	0.124308	0.0668425	0.0334890
0.6	0.291899	0.201634	0.1271780	0.0745776
0.8	0.367118	0.284297	0.2009240	0.1319000
1.0	0.438600	0.371182	0.2866020	0.2054560

| Show Table

DownLoad: CSV

To facilitate a comprehensive comparison, and juxtapose the results obtained from both ARPSM and ATIM for ${\beta_{1}}(\psi, \varphi)$ and ${\beta_{2}}(\psi, \varphi)$ in Problem 2. These tables provide a nuanced understanding of the performance of each method, shedding light on their strengths and limitations in solving the given fractional-order equations in the context of Problem 2.

Table 11. Problem 2 comparison of both methods for

$\beta_{1}(\psi, \varphi)$ .

	$p=1.00$		$p=0.80$		$p=0.60$		$p=0.40$
$\varphi$	ARPSM	ATIM	ARPSM	ATIM	ARPSM	ATIM	ARPSM	ATIM
0.2	0.0004303	0.0004258	0.0009841	0.0009722	0.0022338	0.0022043	0.004802	0.0047363
0.4	0.0014440	0.0014261	0.0027334	0.0026972	0.0049226	0.0048548	0.008209	0.0082096
0.6	0.0030591	0.0030187	0.0050604	0.0049910	0.0078698	0.0077592	0.011257	0.0110971
0.8	0.0052755	0.0052037	0.0078781	0.0077681	0.0110042	0.0108478	0.014096	0.0138937
1.0	0.0080933	0.0079809	0.0111336	0.0109761	0.0142880	0.0140833	0.016789	0.0165466

| Show Table

DownLoad: CSV

Table 12. Problem 2 comparison of both methods for

$\beta_{2}(\psi, \varphi)$ .

	$p=1.00$		$p=0.80$		$p=0.60$		$p=0.40$
$\varphi$	ARPSM	ATIM	ARPSM	ATIM	ARPSM	ATIM	ARPSM	ATIM
0.2	0.0086336	0.0086342	0.022375	0.022376	0.054477	0.054481	0.121739	0.121748
0.4	0.0334867	0.0334890	0.066837	0.066842	0.124300	0.124308	0.211321	0.211336
0.6	0.0745724	0.0745776	0.127169	0.127178	0.201620	0.201634	0.291879	0.291899
0.8	0.1318910	0.1319000	0.200911	0.200924	0.284277	0.284297	0.367092	0.367118
1.0	0.2054420	0.2054560	0.286582	0.286602	0.371156	0.371182	0.438569	0.438600

| Show Table

DownLoad: CSV

7. Conclusions

In summary, this research has focused on improving numerical methods designed for solving the fractional Keller-Segel (KS) model, which is a crucial framework for studying chemotaxis phenomena. By utilizing the Caputo operator framework, we have employed two distinct methodologies: the Aboodh residual power series method (ARPSM) and the Aboodh transform iteration method (ATIM). These approaches have enabled us to obtain accurate solutions to the fractional KS equation, contributing to a better understanding of chemotactic behavior in biological systems. Through a comparative analysis of ARPSM and ATIM, we have revealed their individual strengths and applications in addressing complex fractional models. This work not only advances numerical techniques tailored for fractional differential equations but also improves our understanding of chemotaxis dynamics through precise modeling.

Use of AI tools declaration

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

Acknowledgments

The researcher would like to thank the Deanship of Scientific Research, Shaqra University, for funding the publication of this project.

Conflict of interest

The authors declare that they have no competing interests.

References

[1]	O. P. Agrawal, Fractional variational calculus in terms of Riesz fractional derivatives, J. Phys. A, 40 (2007), 6287-6303. doi: 10.1088/1751-8113/40/24/003
[2]	B. C. Carlson, Lauricella's hypergeometric function F_D, J. Math. Anal. Appl., 7 (1963), 452-470. doi: 10.1016/0022-247X(63)90067-2
[3]	B. C. Carlson, A connection between elementary and higher transcendental functions, SIAM J. Appl. Math., 17 (1969), 116-148. doi: 10.1137/0117013
[4]	B. C. Carlson, Special Functions of Applied Mathematics, Academic Press, New York, 1977.
[5]	B. C. Carlson, B-splines, hypergeometric functions and Dirichlet average, J. Approx. Theory, 67 (1991), 311-325. doi: 10.1016/0021-9045(91)90006-V
[6]	W. Zu Castell, Dirichlet splines as fractional integrals of B-splines, Rocky Mountain J. Math., 32 (2002), 545-559. doi: 10.1216/rmjm/1030539686
[7]	J. Daiya, Representing double Dirichlet average in term of K-Mittag-Leffler function associated with fractional derivative, Journal of Chemical, Biological and Physical Sciences, Section C, 6 (2016), 1034-1045.
[8]	J. Daiya, J. Ram, Dirichlet averages of generalized Hurwitz-Lerch zeta function, Asian J. Math. Comput. Res., 7 (2015), 54-67.
[9]	R. Diaz, E. Pariguan, On hypergeometric functions and k-Pochhammer symbol, Divulg. Mat., 15 (2007), 179-192.
[10]	A. Erdélyi, W. Magnus, F. Oberhettinger, et al. Tables of Integral Transforms, Vol. II, McGrawHill, New York-Toronto-London, 1954.
[11]	S. C. Gupta and B. M. Agrawal, Double Dirichlet average and fractional derivative, Ganita Sandesh, 5 (1991), 47-52.
[12]	R. K. Gupta, B. S. Shaktawat, D. Kumar, Certain relation of generalized fractional calculus associated with the generalized Mittag-Leffler function, J. Raj. Acad. Phy. Sci., 15 (2016), 117-126.
[13]	R. K. Gupta, B. S. Shaktawat, D. Kumar, A study of Saigo-Maeda fractional calculus operators associated with the multiparameter K-Mittag-Leffler function, Asian J. Math. Comput. Res., 12 (2016), 243-251.
[14]	R. K. Gupta, B. S. Shaktawat, D. Kumar, Generalized fractional differintegral operators of the K-series, Honam Math. J., 39 (2017), 61-71. doi: 10.5831/HMJ.2017.39.1.61
[15]	V. Kiryakova, A brief story about the operators of the generalized fractional calculus, Fract. Calc. Appl. Anal., 11 (2008), 203-220.
[16]	D. Kumar, On certain fractional calculus operators involving generalized Mittag-Leffler function, Sahand Communication in Mathematical Analysis, 3 (2016), 33-45.
[17]	D. Kumar, R. K. Gupta, D. S. Rawat, et al. Hypergeometric fractional integrals of multiparameter K-Mittag-Leffler function, Nonlinear Science Letters A: Mathematics, Physics and Mechanics, 9 (2018), 17-26.
[18]	D. Kumar, S. D. Purohit, Fractional differintegral operators of the generalized Mittag-Leffler type function, Malaya J. Mat., 2 (2014), 419-425.
[19]	D. Kumar, S. Kumar, Fractional integrals and derivatives of the generalized Mittag-Leffler type function, Internat. Sch. Res. Not., 2014 (2014), 1-5. doi: 10.1093/imrn/rns215
[20]	Shy-Der Lin, H. M. Srivastava, Some miscellaneous properties and applications of certain operators of fractional calculus, Taiwanese J. Math., 14 (2010), 2469-2495. doi: 10.11650/twjm/1500406085
[21]	P. Massopust, B. Forster, Multivariate complex B-splines and Dirichlet averages, J. Approx. Theory, 162 (2010), 252-269. doi: 10.1016/j.jat.2009.05.002
[22]	E. Neuman, Stolarsky means of several variables, J. Inequal. Pure Appl. Math., 6 (2005), 1-10.
[23]	E. Neuman, P. J. Van Fleet, Moments of Dirichlet splines and their applications to hypergeometric functions, J. Comput. Appl. Math., 53 (1994), 225-241. doi: 10.1016/0377-0427(94)90047-7
[24]	S. G. Samko, A. A. Kilbas, O. I. Marichev, Fractional Integrals and Derivatives: Theory and Applications, Translated from the Russian: Integrals and Derivatives of Fractional Order and Some of Their Applications ("Nauka i Tekhnika", Minsk, 1987), Gordon and Breach Science Publishers, UK, 1993.
[25]	R. K. Saxena, J. Daiya, Integral transforms of the S-function, Le Mathematiche, 70 (2015), 147-159.
[26]	R. K. Saxena, J. Daiya, A. Singh, Integral transforms of the k-Mittag-Leffler function $E_{k,\alpha ,\beta }^\gamma \left( z \right)$ , Le Mathematiche, 69 (2014), 7-16.
[27]	R. K. Saxena, T. K. Pogány, J. Ram, et al. Dirichlet averages of generalized multi-index MittagLeffler functions, Armen. J. Math., 3 (2010), 174-187.
[28]	T. Zhang, L. Xiong, Periodic motion for impulsive fractional functional differential equations with piecewise Caputo derivative, Appl. Math. Lett., 101 (2020), 106072.

Reader Comments

Your name:*

Email:*
© 2020 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.1

Metrics

Article views(4154) PDF downloads(342) Cited by(3)

Preview PDF

Download XML

Export Citation

Article outline

Show full outline

AIMS Mathematics

S-function associated with fractional derivative and double Dirichlet average

Related Papers:

Abstract

1. Introduction

2. Basic concepts

3. Proposed methodologies

3.1. Aboodh residual power series method (ARPSM)

3.2. Aboodh transform iteration method (ATIM)

4. Problem 1

4.1. Problem 1 with ARPSM

4.2. Problem 1 with ATIM

5. Problem 2

5.1. Problem 2 with ARPSM

5.2. Problem 2 with ATIM

6. Graphics and tables discussion

7. Conclusions

Use of AI tools declaration

Acknowledgments

Conflict of interest

References

Reader Comments

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

Metrics

Other Articles By Authors

Catalog

AIMS Mathematics

S-function associated with fractional derivative and double Dirichlet average

Related Papers:

Abstract

1. Introduction

2. Basic concepts

3. Proposed methodologies

3.1. Aboodh residual power series method (ARPSM)

3.2. Aboodh transform iteration method (ATIM)

4. Problem 1

4.1. Problem 1 with ARPSM

4.2. Problem 1 with ATIM

5. Problem 2

5.1. Problem 2 with ARPSM

5.2. Problem 2 with ATIM

6. Graphics and tables discussion

7. Conclusions

Use of AI tools declaration

Acknowledgments

Conflict of interest

References

Reader Comments

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

Metrics

Other Articles By Authors

Related pages

Tools

Export File

Citation

Format

Content

Catalog