On the Application of Adomian Decomposition Method to Special Equations in Physical Sciences

Aishah Alsulami, Mariam Al-Mazmumy, Huda Bakodah, Nawal Alzaid^*

Department of Mathematics, College of Science, University of Jeddah, Jeddah, Saudi Arabia.

DOI: 10.4236/ajcm.2023.133022                

Abstract

The current manuscript makes use of the prominent iterative procedure, called the Adomian Decomposition Method (ADM), to tackle some important special differential equations. The equations of curiosity in this study are the singular equations that arise in many physical science applications. Thus, through the application of the ADM, a generalized recursive scheme was successfully derived and further utilized to obtain closed-form solutions for the models under consideration. The method is, indeed, fascinating as respective exact analytical solutions are accurately acquired with only a small number of iterations.

Keywords

Iterative Scheme, Adomian Decomposition Method, Initial-Value Problems, Singular Ordinary Differential Equations

1. Introduction

Special differential equations of second-order are held with utmost esteem in many fields of physical sciences due to their fascinating properties [1] [2] [3] . Among others, they are characterized by the possession of the well-known special orthogonal functions, which are widely used in science and engineering domains. Additionally, these functions are greatly utilized in numerical methods, and approximation theories, just to mention a few. Besides, solutions of many differential equations are recast to the form of these special functions in mathematical physics [4] - [9] .

However, the method of interest in the present study is the celebrated Adomian Decomposition Method (ADM) [10] [11] . This method was devised in the 1980s by George Adomian to treat differential and integral equations. Various reformations and extensions of this method have been reported in the literature; in addition to different reliable modifications of the method among others [12] [13] [14] . Indeed, there exists a huge number of related literature with regard to the development of the ADM that is associated with various forms of Initial-Value Problems (IVPs) of ordinary and partial differential equations [15] [16] [17] [18] . Both the IVPs and boundary-value problems have been effectively treated via the application of the ADM to obtain their resulting closed-form (series) solutions or even the exact analytical solutions in many cases. This is, however, related to the fact that the procedure devised by Adomian tends to rapidly converge to the exact analytical solution, whenever obtainable, see [19] and the references therein.

So far, there have been few works on the solution of special functions. In [9] , the series solutions of some special functions are presented. Also, the authors in [20] studied the solution of Legendre differential equation and Chebyshev’s differential equation by multiplying them by the singular coefficient and expressed it as based on ADM. In addition, the exact solutions for the Hypergeometric equation and Legendre equation, by writing the equation in the general operator form and finding its inverse, are given in [21] . Next, classical results on the Gauss Hypergeometric equation, Confluent Hypergeometric equation and Bessel equation are discussed in [22] . Additionally, differential transform method for Gauss Hypergeometric equation and Laguerre equations is applied in [23] . As well, the Frobenius method around all regular singular points for κ-Hypergeometric equation is employed in [24] , Recently, [25] [26] thoroughly discussed the technique for obtaining an exact solution based on the modified Adomian decomposition method for the Laguerre equations and Chebyshev’s differential equation.

What’s more, the present study aims to tackle some important special second-order differential equations that arise in many physical models of physical sciences. More specifically, the singular models to be examined in the present study are mentioned as follows: Legendre’s differential equation

$\left(1-x^2\right){v^{\prime }}^{\prime }-2x{v}^{\prime }+n\left(n+1\right)v=0,-1<x<1.$ (1)

Chebyshev’s differential equation

$\left(1-x^2\right){v^{\prime }}^{\prime }-x{v}^{\prime }+n^{2}v=0,-1<x<1.$ (2)

Hermite’s differential equation

${v^{\prime }}^{\prime }-2x{v}^{\prime }+2nv=0,$ (3)

where in the above equations, n is a non-negative integer. In doing so, we will utilize the aforementioned classical ADM procedure to effectively acquire their respective recursive relations, leading to their closed forms or even exact analytical solutions when the value of n is specified. Additionally, we arrange the paper as follows: Section 2 gives the analysis of the approach of interest—the ADM. Section 3 gives the application of the ADM procedure to the models of concern; while Section 4 gives some concluding remarks.

2. Analysis of the Method

The present section gives a generalized derivation procedure for tackling nonlinear Initial-Value Problems (IVPs) based on the ADM. To do so, let us consider the following differential equation [16]

$G\left(v\left(x\right)\right)=g\left(x\right)\mathrm{,}$ (4)

with G representing a generalized ordinary (or partial) differential operator. This operator being general, it can equally be expressed to involve both linear and nonlinear operators. Thus, we decompose the operator further, and rewrite the above equation as follows

$Lv+Rv+Nv=g,$ (5)

where L is the highest linear operator that is invertible, with $R<L$ ; while N is specifically the nonlinear operator. More so, we rewrite the latter equation as follows

$Lv=g-Rv-Nv,$

such that applying the inverse operator $L^{-1}$ to both sides of the above equation yields

$v=\varphi \left(x\right)+L^{-1}g-L^{-1}Rv-L^{-1}Nv\mathrm{.}$ (6)

where $\varphi \left(x\right)$ is the function emanating from the prescribed initial data.

Further, the iterative procedure by the name ADM decomposes the solution $v\left(x\right)$ using an infinite series of the following form

$v\left(x\right)={\displaystyle \underset{m=0}{\overset{\infty }{\sum }}}{v}_m\left(x\right),$ (7)

while the nonlinear component Nv is equally decomposed using the following infinite series

$N\left(v\right)={\displaystyle \underset{m=0}{\overset{\infty }{\sum }}}{B}_m,$ (8)

where $B_m$ ’s are polynomials devised by Adomian, and recursively determined using the following scheme [18]

$B_m=\frac{1}{m!}\frac{{\text{d}}^m}{\text{d}{\lambda }^m}{\left[N\left({\displaystyle \underset{j=0}{\overset{m}{\sum }}}{\lambda }^j{v}_j\right)\right]}_{\lambda =0},m=0,1,2,\cdots$ (9)

Therefore, upon substituting Equations (7) and (8) into Equation (6), one gets

${\displaystyle \underset{m=0}{\overset{\infty }{\sum }}}{v}_m\left(x\right)=\varphi \left(x\right)+L^{-1}g\left(x\right)-L^{-1}R{\displaystyle \underset{m=0}{\overset{\infty }{\sum }}}{v}_m(x)-L^{-1}{\displaystyle \underset{m=0}{\overset{\infty }{\sum }}}{B}_m.$

Furthermore, the ADM procedure swiftly reveals the generalized recursive solution for the problem from the above equation as follows

$\{\begin{array}{l}{v}_0=\varphi \left(x\right)+L^{-1}g\left(x\right),\hfill \\ v_1=-L^{-1}R{v}_0-L^{-1}{B}_0,\hfill \\ v_2=-L^{-1}R{v}_1-L^{-1}{B}_1,\hfill \\ v_3=-L^{-1}R{v}_2-L^{-1}{B}_2,\hfill \\ ⋮\hfill \\ v_{m+1}=-L^{-1}R{v}_m-L^{-1}{B}_m,m\ge 0,\hfill \end{array}$ (10)

where $B_m$ ’s are the Adomian polynomials computed from Equation (9). Expressing few of these terms, we get

$\begin{array}{l}{B}_0=N\left(v_0\right),\\ B_1=\frac{\text{d}N\left(v_0\right)}{d{v}_0}{v}_1,\\ B_2=\frac{\text{d}N\left(v_0\right)}{\text{d}{v}_0}{v}_2+\frac{1}{2}\frac{{\text{d}}^{2}N\left(v_0\right)}{\text{d}{v}_0^2}{v}_1^2,\\ B_3=\frac{\text{d}N\left(v_0\right)}{\text{d}{v}_0}{v}_3+\frac{{\text{d}}^{2}N\left(v_0\right)}{\text{d}{v}_0^2}{v}_1{v}_2+\frac{1}{3!}\frac{{\text{d}}^{3}N\left(v_0\right)}{\text{d}{v}_0^3}{v}_1^3,\\ ⋮\end{array}$ (11)

Remarkable, it is obvious that the Adomian polynomials $B_m$ ’s depend on the solution components $v_m$ . For instance, $B_0$ relies merely on $v_0$ ; $B_1$ relies merely on $v_0$ and $v_1$ ; $B_2$ relies merely on $v_0$ , $v_1$ and $v_2$ , and so on.

Finally, a realistic solution is obtained by considering the following m-term approximations as

${\Psi }_m={\displaystyle \underset{j=0}{\overset{m-1}{\sum }}}{v}_j,$ (12)

where

$v\left(x\right)=\underset{m\to \infty }{\lim }{\Psi }_m\left(x\right)={\displaystyle \underset{j=0}{\overset{\infty }{\sum }}}{v}_j\left(x\right).$ (13)

3. Application to Standard Equations

This section presents the application of the ADM procedure presented earlier in the above section to some selected singular special differential equations’ IVPs. More specifically, we make consideration to Legendre’s differential equation, Chebyshev’s differential equation and Hermite’s differential equation.

3.1. Legendre’s Differential Equation

Consider the Legendre’s differential equation

$\left(1-x^2\right){v^{\prime }}^{\prime }-2x{v}^{\prime }+n\left(n+1\right)v=0,-1<x<1,$ (14)

or alternatively rewritten as follows

${v^{\prime }}^{\prime }=x^2{v^{\prime }}^{\prime }+2x{v}^{\prime }-n\left(n+1\right)v,-1<x<1.$ (15)

So, we consider the right-hand side of the equation as a normal nonhomogeneous term, where the differential operator L is defined by $\frac{{\text{d}}^2}{\text{d}{x}^2}$ .

More so, we consider the inverse operator $L^{-1}$ as a two-fold integral operator defined by

$L^{-1}\left(\mathrm{.}\right)={\displaystyle {\int }_0^x}{\displaystyle {\int }_0^x}\left(\right)\text{d}x\text{d}x\mathrm{.}$

Thus, applying the inverse operator $L^{-1}$ to both side of Equation (15), one gets

$v=\varphi \left(x\right)+L^{-1}\left[x^2{v^{\prime }}^{\prime }+2x{v}^{\prime }-n\left(n+1\right)v\right]\mathrm{,}$

such that

$L\varphi \left(x\right)=0.$

Thus, based on the ADM procedure, the solution $v\left(x\right)$ is introduced through an infinite summation of components $v_m\left(x\right)$ earlier discussed in the methodology. Hence, the recursive solution of Legendre’s differential equation is as follows

$\begin{array}{l}{v}_0\left(x\right)=\varphi \left(x\right),\\ v_{k+1}\left(x\right)=L^{-1}\left[x^2{{v^{\prime }}^{\prime }}_k+2x{v^{\prime }}_k-n\left(n+1\right)v_k\right],k\ge 0\end{array}$

where the overall solution $v\left(x\right)$ follows immediately by summing the above components as follows

$v\left(x\right)=\underset{m\to \infty }{\lim }{\Psi }_m\left(x\right)={\displaystyle \underset{j=0}{\overset{\infty }{\sum }}}{v}_j\left(x\right).$ (16)

Case 1: $n=1$

Let us consider an IVP featuring Legendre’s differential equation with $n=1$ as follows

$\begin{array}{l}\left(1-x^2\right){v^{\prime }}^{\prime }-2x{v}^{\prime }+2v=0,\\ v\left(0\right)=0,v^{\prime }\left(0\right)=1.\end{array}$ (17)

Therefore, we re-express the governing differential equation using operator notation as follows

$Lv=x^2{v^{\prime }}^{\prime }+2x{v}^{\prime }-2v.$ (18)

What’s more, we apply the inverse operator $L^{-1}$ to the both sides of the above equation to obtain

$v=x+L^{-1}\left[x^2{v^{\prime }}^{\prime }+2x{v}^{\prime }-2v\right]\mathrm{.}$

Lastly, through the application of the ADM, the overall recursive relation takes the following form

$\begin{array}{l}{v}_0\left(x\right)=x,\\ v_1\left(x\right)=L^{-1}\left[x^2{{v^{\prime }}^{\prime }}_0+2x{v^{\prime }}_0-2v_0\right]=0,\\ v_{k+1}\left(x\right)=L^{-1}\left[x^2{{v^{\prime }}^{\prime }}_k+2x{v^{\prime }}_k-2v_k\right]=0,k\ge 1,\end{array}$

and it implies that

$v\left(x\right)=x.$ (19)

In fact, this is a well-known exact analytical solution for Legendre’s differential equation when $n=1$ .

Case 2: $n=2$

Let us consider an IVP featuring Legendre’s differential equation with $n=2$ as follows

$\begin{array}{l}\left(1-x^2\right){v^{\prime }}^{\prime }-2x{v}^{\prime }+6v=0,\\ v\left(0\right)=-\frac{1}{2},v^{\prime }\left(0\right)=0.\end{array}$ (20)

Thus, upon expressing the governing differential equation using operator notation as follows

$Lv=x^2{v^{\prime }}^{\prime }+2x{v}^{\prime }-6v,$ (21)

we further apply the inverse operator $L^{-1}$ to both sides of the latter equation to obtain

$v=-\frac{1}{2}+L^{-1}\left[x^2{v^{\prime }}^{\prime }+2x{v}^{\prime }-6v\right].$

Then, the application of the ADM yields the following overall recursive relation

$\begin{array}{l}{v}_0\left(x\right)=-\frac{1}{2},\\ v_1\left(x\right)=L^{-1}\left[x^2{{v^{\prime }}^{\prime }}_0+2x{v^{\prime }}_0-6v_0\right]=\frac{3}{2}{x}^2,\\ v_2\left(x\right)=L^{-1}\left[x^2{{v^{\prime }}^{\prime }}_1+2x{v^{\prime }}_1-6v_1\right]=0,\\ v_{k+1}\left(x\right)=L^{-1}\left[x^2{{v^{\prime }}^{\prime }}_k+2x{v^{\prime }}_k-6v_k\right]=0,k\ge 2,\end{array}$

and it implies,

$v\left(x\right)=-\frac{1}{2}+\frac{3}{2}{x}^2,$ (22)

which is a known exact analytical solution for Legendre’s differential equation when $n=2$ .

3.2. Chebyshev’s Differential Equation

Consider the Chebyshev’s differential equation

$\left(1-x^2\right){v^{\prime }}^{\prime }-x{v}^{\prime }+n^{2}v=0,-1<x<1.$ (23)

Equally, using the ADM procedure described above, the corresponding recursive scheme is obtained by

$\begin{array}{l}{v}_0\left(x\right)=\varphi \left(x\right),\\ v_{k+1}\left(x\right)=L^{-1}\left[x^2{{v^{\prime }}^{\prime }}_k+x{v^{\prime }}_k-n^2{v}_k\right],k\ge 0.\end{array}$

Case 1: $n=1$

Let us consider an IVP featuring Chebyshev’s differential equation with $n=1$ as follows

$\begin{array}{l}\left(1-x^2\right){v^{\prime }}^{\prime }-x{v}^{\prime }+v=0,\\ v\left(0\right)=0,v^{\prime }\left(0\right)=1.\end{array}$ (24)

Now, expressing the governing differential equation using operator notation becomes

$Lv=x^2{v^{\prime }}^{\prime }+x{v}^{\prime }-v,$ (25)

such that upon applying $L^{-1}$ to both sides of Equation (25) reveals

$v=x+L^{-1}\left[x^2{v^{\prime }}^{\prime }+x{v}^{\prime }-v\right]\mathrm{.}$

What’s more, without loss of generality, the resulting recursive scheme for the IVP is found as follows

$\begin{array}{l}{v}_0\left(x\right)=x,\\ v_{k+1}\left(x\right)=L^{-1}\left[x^2{{v^{\prime }}^{\prime }}_k+x{v^{\prime }}_k-v_k\right]=0,k\ge 0,\end{array}$

implying

$v\left(x\right)=x,$ (26)

which is indeed the exact analytical solution for Chebyshev’s differential equation when $n=1$ .

Case 2: $n=2$

Let us consider an IVP featuring Chebyshev’s differential equation with $n=2$ as follows

$\begin{array}{l}\left(1-x^2\right){v^{\prime }}^{\prime }-x{v}^{\prime }+4v=0,\\ v\left(0\right)=-1,v^{\prime }\left(0\right)=0.\end{array}$ (27)

In the same manner, we express the governing differential equation above using operator notation as follows

$Lv=\left[x^2{v^{\prime }}^{\prime }+x{v}^{\prime }-4v\right],$ (28)

such that after applying the inverse linear differential operator $L^{-1}$ to both sides of the latter equation reveals

$v=-1+L^{-1}\left[x^2{v^{\prime }}^{\prime }+x{v}^{\prime }-4v\right]\mathrm{.}$

As in the preceding problem, the resulting recursive relation is thus obtained as follows

$\begin{array}{l}{v}_0\left(x\right)=-1,\\ v_1\left(x\right)=L^{-1}\left[x^2{{v^{\prime }}^{\prime }}_0+x{v^{\prime }}_0-4v_0\right]=2x^2,\\ v_2\left(x\right)=L^{-1}\left[x^2{{v^{\prime }}^{\prime }}_1+x{v^{\prime }}_1-4v_1\right]=0,\\ v_{k+1}\left(x\right)=L^{-1}\left[x^2{{v^{\prime }}^{\prime }}_k+x{v^{\prime }}_k-4v_k\right]=0,k\ge 2,\end{array}$

which upon taking the sum gives

$v\left(x\right)=-1+2x^2,$ (29)

that is also the exact analytical solution for Chebyshev’s differential equation when $n=2$ .

3.3. Hermite’s Differential Equation

Consider the Hermite’s differential equation

${v^{\prime }}^{\prime }-2x{v}^{\prime }+2nv=0.$ (30)

Equally, using the ADM procedure described above, the corresponding recursive scheme is obtained by

$\begin{array}{l}{v}_0\left(x\right)=\varphi \left(x\right),\\ v_{k+1}\left(x\right)=L^{-1}\left[2x{v^{\prime }}_k-2n{v}_k\right],k\ge 0.\end{array}$

Case 1: $n=1$

Let us consider an IVP of Hermite’s differential equation with $n=1$ as follows

$\begin{array}{l}{v^{\prime }}^{\prime }-2x{v}^{\prime }+2v=0,\\ v\left(0\right)=0,v^{\prime }\left(0\right)=2.\end{array}$ (31)

In a similar way, we express the governing differential equation above via an operator notation as follows

$Lv=2x{v}^{\prime }-2v,$ (32)

such that after applying the inverse linear differential operator $L^{-1}$ to both sides of the latter equation reveals

$v=2x+L^{-1}\left[2x{v}^{\prime }-2v\right]\mathrm{.}$

Then, as in the preceding problem, the resulting recursive relation is thus obtained as follows

$\begin{array}{l}{v}_0\left(x\right)=2x,\\ v_{k+1}\left(x\right)=L^{-1}\left[2x{v^{\prime }}_k-2v_k\right]=0,k\ge 0,\end{array}$

which upon summing the components $v_n$ gives

$v\left(x\right)=2x,$ (33)

which is a well-known exact analytical solution for Hermite’s differential equation when $n=1$ .

Case 2: $n=2$

Let us consider an IVP of Hermite’s differential equation with $n=2$ as follows

$\begin{array}{l}{v^{\prime }}^{\prime }-2x{v}^{\prime }+4v=0,\\ v\left(0\right)=-2,v^{\prime }\left(0\right)=0.\end{array}$ (34)

Therefore, expressing the governing differential equation above via an operator notation, one gets

$Lv=\left[2x{v}^{\prime }-4v\right]\mathrm{,}$ (35)

such that after applying the inverse linear differential operator $L^{-1}$ to both sides of the latter equation reveals

$v=-2+L^{-1}\left[2x{v}^{\prime }-4v\right].$

Lastly, the resulting recursive relation is thus obtained via the application of the ADM procedure as follows

$\begin{array}{l}{v}_0=-2,\\ v_1\left(x\right)=L^{-1}\left[2x{v^{\prime }}_0-4v_0\right]=4x^2,\\ v_2\left(x\right)=L^{-1}\left[2x{v^{\prime }}_1-4v_1\right]=0,\\ v_{k+1}\left(x\right)=L^{-1}\left[2x{v^{\prime }}_k-4v_k\right]=0,k\ge 2,\end{array}$

that sums to,

$v\left(x\right)=-2+4x^2,$ (36)

which is the exact analytical solution for Hermite’s differential equation when $n=2$ .

Notably, it is worth mentioning here that, as the present study considers only the Legendre’s, Chebyshev’s, and Hermite’s differential equations for $n=1$ and $n=2$ , the same ADM procedure can equally be extended to other classes of special differential equations in mathematical physics for any value of n, in general. Equations like Laguerre, Gegenbauer, Jacobi differential equations, and others (both singular and nonsingular) could, in the same way, be treated via the ADM for any given n. More so, other special equations that have no orthogonal functions like Bessel, modified Bessel, Riccati, and Euler to state a few could also be examined by the method; in addition to the higher-order differential equations that are also evenly covered by the technique.

4. Conclusion

In conclusion, the present study makes use of the prominent Adomian iterative procedure to recurrently tackle some important special differential equations in physical sciences. As a particular interest, the known singular equations including Legendre’s equation, Chebyshev’s equation, and Hermite’s equation have been successfully treated as test problems. Respective exact closed-form solutions for the models have been fruitfully acquired when $n=1$ and $n=2$ . In fact, the acquired solutions happen to be precisely the same as the available results in the literature. Lastly, the used methodology is highly recommended for the study of both singular and nonsingular functional equations.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Alkarawi, A.H. and Al-Saiq, I.R. (2020) Applications of Modified Adomian Decomposition Method for Solving the Second-Order Ordinary Differential Equations. Journal of Physics: Conference Series, 1530, Article ID: 012155. https://doi.org/10.1088/1742-6596/1530/1/012155
[2]	Hasan, Y.Q. and Zhu, L.M. (2007) Solving Singular Initial Value Problems in the Second-Order Ordinary Differential Equations. Journal of Applied Sciences, 7, 2505-2508. https://doi.org/10.3923/jas.2007.2505.2508
[3]	Wazwaz, A. (2002) A New Method for Solving Singular Initial Value Problems in the Second-Order Ordinary Differential Equations. Applied Mathematics and Computation, 128, 45-57. https://doi.org/10.1016/S0096-3003(01)00021-2
[4]	Bildik, N. and Deniz, S. (2015) Modified Adomian Decomposition Method for Solving Riccati Differential Equations. Review of the Air Force Academy, 3, 21-26. https://doi.org/10.19062/1842-9238.2015.13.3.3
[5]	Liu, Y.C. (2012) Application of Legendre Polynomials in Adomian Decomposition Method. Proceedings of 2012 International Conference on Computer, Electrical, and Systems Sciences, Amsterdam, 13-14 May 2012, 567-571.
[6]	Hosseini, M.M. (2006) Adomian Decomposition Method with Chebyshev Polynomials. Applied Mathematics and Computation, 175, 1685-1693. https://doi.org/10.1016/j.amc.2005.09.014
[7]	Çenesiz, Y. and Kurnaz, A. (2011) Adomian Decomposition Method by Gegenbauer and Jacobi Polynomials. International Journal of Computer Mathematics, 88, 3666-3676. https://doi.org/10.1080/00207160.2011.611503
[8]	Mahmoudi, Y., Abdollahi, M., Karimian, N. and Khalili, H. (2012) Adomian Decomposition Method with Laguerre Polynomials for Solving Ordinary Differential Equation. Journal of Basic and Applied Scientific Research, 2, 2236-12241.
[9]	Bell, W.W. (2004) Special Functions for Scientists and Engineers. Courier Corporation, Chelmsford.
[10]	Adomian, G. (1990) A Review of the Decomposition Method and Some Recent Results for Nonlinear Equations. Mathematical and Computer Modelling, 13, 17-43. https://doi.org/10.1016/0895-7177(90)90125-7
[11]	Adomian, G. (1994) Solving Frontier Problems of Physics: The Decomposition Method. Fundamental Theories of Physics. Springer, Berlin.
[12]	Wazwaz, A. (2002) A New Algorithm for Calculating Adomian Polynomials for Nonlinear Operators. Applied Mathematics and Computation, 111, 33-51. https://doi.org/10.1016/S0096-3003(99)00063-6
[13]	Wazwaz, A. and El-Sayed, S.M. (2001) A New Modification of the Adomian Decomposition Method for Linear and Nonlinear Operators. Applied Mathematics and Computation, 122, 393-405. https://doi.org/10.1016/S0096-3003(00)00060-6
[14]	Wazwaz, A. (1999) A Reliable Modification of Adomian Decomposition Method. Applied Mathematics and Computation, 102, 77-86. https://doi.org/10.1016/S0096-3003(98)10024-3
[15]	Liu, Y. (2009) Adomian Decomposition Method with Second Kind Chebyshev Polynomials. Proceedings of the Jangjeon Mathematical Society, 12, 57-67.
[16]	Tomaizeh, H. (2017) Modified Adomian Decomposition Method for Differential Equations. Hebron University, Hebron.
[17]	Hosseini, M.M. and Nasabzadeh, H. (2007) Modified Adomian Decomposition Method for Specific Second Order Ordinary Differential Equations. Applied Mathematics and Computation, 186, 17-123. https://doi.org/10.1016/j.amc.2006.07.094
[18]	Nuruddeen, R.I., Muhammad, L., Nass, A.M. and Sulaiman, T.A. (2018) A Review of the Integral Transforms-Based Decomposition Methods and Their Applications in Solving Nonlinear PDEs. Palestine Journal of Mathematics, 7, 262-280.
[19]	Bakodah, H.O., Banaja, M.A., Alrigi, B.A., Ebaid, A. and Rach, R. (2019) An Efficient Modification of the Decomposition Method with a Convergence Parameter for Solving Korteweg de Vries Equations. Journal of King Saud University-Science, 31, 1424-1430. https://doi.org/10.1016/j.jksus.2018.11.010
[20]	Rach, R., Baghdasarian, A. and Adomian, G. (1992) Differential Equations with Singular Coefficients. Applied Mathematics and Computation, 47, 179-184. https://doi.org/10.1016/0096-3003(92)90045-3
[21]	Bougoffa, L. (2009) On the Exact Solutions for Initial Value Problems of Second-Order Differential Equations. Applied Mathematics Letters, 22, 1248-1251. https://doi.org/10.1016/j.aml.2009.01.038
[22]	Dita, P. and Grama, N. (1997) On Adomian’s Decomposition Method for Solving Differential Equations. ArXiv: Solv-int/9705008.
[23]	Suresh, P.L., Krishna, G.V., Maheswari, K.U. and Ramanaiah, J.V. (2017) Solutions of Gauss’s Hypergeometric Equation, Laguerre’s Equation by Differential Transform Method. International Journal of Current Engineering and Scientific Research, 4, 32-38.
[24]	Mubeen, S., Naz, M., Rehman, A. and Rahman, G. (2014) Solutions of κ-Hypergeometric Differential Equations. Journal of Applied Mathematics, 2014, Article ID: 128787. https://doi.org/10.1155/2014/128787
[25]	Al-Mazmumy, M. and Alsulami, A.A. (2023) Solution of Laguerre’s Differential Equations via Modified Adomian Decomposition Method. Journal of Applied Mathematics and Physics, 11, 85-100. https://doi.org/10.4236/jamp.2023.111007
[26]	Al Mazmumy, M., Alsulami, A., Bakodah, H. and Alzaid, N. (2023) Adomian Modification Methods for the Solution of Chebyshev’s Differential Equations. Applied Mathematics, 14, 512-530.