The Adomian Decomposition Method for Solving Nonlinear Partial Differential Equation Using Maple ()

Dalal Adnan Maturi^{}, Honaida Mohammed Malaikah^{}

Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia.

**DOI: **10.4236/apm.2021.116038
PDF
HTML XML
455
Downloads
2,894
Views
Citations

Department of Mathematics, Faculty of Science, King Abdulaziz University, Jeddah, Saudi Arabia.

The nonlinear partial differential equation is solved using the Adomian decomposition method (ADM) in this article. A number of examples have been provided to illustrate the numerical results, which is the comparison of the exact and numerical solutions, and it has been discovered through the tables that the amount of error between the exact and numerical solutions is very small and almost non-existent, and the graph also shows how the exact solution of absolutely applies to the numerical solution. This demonstrates the precision of the Adomian decomposition method (ADM) for solving the nonlinear partial differential equation with Maple18. And that in terms of obtaining numerical results, this approach is characterized by ease, speed, and high accuracy.

Share and Cite:

Maturi, D. and Malaikah, H. (2021) The Adomian Decomposition Method for Solving Nonlinear Partial Differential Equation Using Maple. *Advances in Pure Mathematics*, **11**, 595-603. doi: 10.4236/apm.2021.116038.

1. Introduction

The aim of this study is to use Maple18 to solve the Volterra-Fredholm integral equation with the Adomian decomposition process. Integral equations are fundamental sciences in our everyday lives; they describe physical, chemical, engineering, and medical phenomena, and they also help us find analytical and numerical solutions to these problems.

A reliable Modification of Adomain Decomposition Method [1]. A new Modification of the Adomain Decomposition Method for Linear and Nonlinear Operators [2]. Lorenz equations are solved using a decomposition method [3]. Using Adomian’s decomposition method to solve the Riccati differential equation [4].

Numeric-analytic integration of strongly nonlinear and chaotic oscillators using Adomian decomposition [5]. For fourth-order boundary value problems, the extended Adomian decomposition method [6]. The Adomian decomposition approach has been used to solve multipoint boundary value problems [7]. A new algorithm for evaluating Adomian polynomials has been developed [8]. An efficient algorithm for the multivariable Adomian polynomials [9]. Convenient analytic recurrence algorithms for the Adomian polynomials [10]. A review of the Adomian decomposition method and its applications to fractional differential equations are discussed in this paper [11]. A bibliography of the Adomian decomposition method's theory and applications [12]. Nonlinear integral equation solutions are more difficult to solve than linear integral equation solutions [13], and there are several analytical and computational methods for solving both linear and nonlinear integral equations [14] - [19].

MATLAB and Maple were used to implement the Adomian decomposition method for the Fredholm integral equation of the second kind. To solve the Fredholm integral equation of the second kind [20], the Adomian decomposition method was employed. In addition, using Maple, a Modified research approach for solving the Volterra integral equation of the second kind [21]. The Adomian Decomposition Method for Solving Volterra-Fredholm Integral Equation Using Maple [22].

In this article, we used the Maple algorithm to apply the Adomian decomposition method to various cases, such as finding the approximate solution, comparing it to the exact solution, and determining the amount of error between the approximate solution and the exact solution.

2. Adomian Decomposition Method

Consider the nonlinear partial differential equation given in an operator form

${L}_{x}u\left(x,y\right)+{L}_{y}u\left(x,y\right)+R\left(u\left(x,y\right)\right)+F\left(u\left(x,y\right)\right)=g\left(x,y\right)$ (1)

where *L _{x}* is the highest order differential in

Assuming that the operator *L _{x}* meets the two bases of selection, therefore we set

${L}_{x}u\left(x,y\right)=g\left(x,y\right)-{L}_{y}u\left(x,y\right)-R\left(u\left(x,y\right)\right)-F\left(u\left(x,y\right)\right)$ (2)

Applying ${L}_{x}^{-1}$ to both sides of (1) gives

$u\left(x,y\right)={\varnothing}_{0}-{L}_{x}^{-1}g\left(x,y\right)-{L}_{x}^{-1}{L}_{y}u\left(x,y\right)-{L}_{x}^{-1}R\left(u\left(x,y\right)\right)-{L}_{x}^{-1}R\left(u\left(x,y\right)\right),$ (3)

${\varnothing}_{0}=\{\begin{array}{l}u\left(0,y\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}L=\frac{\partial}{\partial x}\\ u\left(0,y\right)+x{u}_{x}\left(0,y\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}L=\frac{{\partial}^{2}}{\partial {x}^{2}}\\ u\left(0,y\right)+x{u}_{x}\left(0,y\right)+\frac{1}{2!}{x}^{2}{u}_{xx}\left(0,y\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}L=\frac{{\partial}^{3}}{\partial {x}^{3}}\\ u\left(0,y\right)+x{u}_{x}\left(0,y\right)+\frac{1}{2!}{x}^{2}{u}_{xx}\left(0,y\right)+\frac{1}{3!}{x}^{3}{u}_{xxx}\left(0,y\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{for}\text{\hspace{0.17em}}L=\frac{{\partial}^{4}}{\partial {x}^{4}}\end{array}$

We proceed in exactly the same manner by calculating the solution $u\left(x,y\right)$.

In a series form

$u\left(x,y\right)={\displaystyle {\sum}_{n=0}^{\infty}{u}_{n}}\left(x,y\right),$ (4)

And the nonlinear term

$F\left(u\left(x,y\right)\right)={\displaystyle {\sum}_{n=0}^{\infty}{A}_{n}},$

where *A _{n}* are Adomian polynomials that can be generated for all forms of nonlinearity. Based on these assumptions, Equation (2) becomes

$\begin{array}{c}{\displaystyle {\sum}_{n=0}^{\infty}{u}_{n}\left(x,y\right)}={\varnothing}_{0}-{L}_{x}^{-1}g\left(x,y\right)-{L}_{x}^{-1}{L}_{y}\left({\displaystyle {\sum}_{n=0}^{\infty}{u}_{n}\left(x,y\right)}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.05em}}-{L}_{x}^{-1}R\left({\displaystyle {\sum}_{n=0}^{\infty}{u}_{n}\left(x,y\right)}\right)-{L}_{x}^{-1}\left({\displaystyle {\sum}_{n=0}^{\infty}{A}_{n}}\right),\end{array}$ (5)

The components, ${u}_{n}\left(x,y\right),n\ge 0$ of the solution $u\left(x,y\right)$ can be recursively determined by using the relation

${u}_{0}\left(x,y\right)={\varnothing}_{0}-{L}_{x}^{-1}g\left(x,y\right),$

${u}_{k+1}\left(x,y\right)=-{L}_{x}^{-1}{L}_{y}{u}_{k}\left(x,y\right)-{L}_{x}^{-1}R\left({u}_{k}\left(x,y\right)\right)-{L}_{x}^{-1}{A}_{k},\text{\hspace{0.17em}}\text{\hspace{0.17em}}k\ge 0$

Using the algorithms described before for calculating *A _{n}* for the nonlinear term
$F\left(u\right)$.

The first few components can be identified by

${u}_{0}\left(x,y\right)={\varnothing}_{0}-{L}_{x}^{-1}g\left(x,y\right),$

${u}_{1}\left(x,y\right)=-{L}_{x}^{-1}{L}_{y}{u}_{0}\left(x,y\right)-{L}_{x}^{-1}R\left({u}_{0}\left(x,y\right)\right)-{L}_{x}^{-1}{A}_{0},$

${u}_{2}\left(x,y\right)=-{L}_{x}^{-1}{L}_{y}{u}_{1}\left(x,y\right)-{L}_{x}^{-1}R\left({u}_{1}\left(x,y\right)\right)-{L}_{x}^{-1}{A}_{1},$

${u}_{3}\left(x,y\right)=-{L}_{x}^{-1}{L}_{y}{u}_{2}\left(x,y\right)-{L}_{x}^{-1}R\left({u}_{2}\left(x,y\right)\right)-{L}_{x}^{-1}{A}_{2},$

${u}_{4}\left(x,y\right)=-{L}_{x}^{-1}{L}_{y}{u}_{3}\left(x,y\right)-{L}_{x}^{-1}R\left({u}_{3}\left(x,y\right)\right)-{L}_{x}^{-1}{A}_{3},$

where each components can be determined by using the preceding component. Having calculated the components ${u}_{n}\left(x,y\right)$.

3. Numerical Examples

In this section, we solve some examples, and we can compare the numerical results with the exact solution.

Example 1. Consider the nonlinear partial differential equation

${u}_{t}+u{u}_{x}=0,\text{\hspace{0.17em}}\text{\hspace{0.17em}}u\left(x,0\right)=x,\text{\hspace{0.17em}}\text{\hspace{0.17em}}t>0,$ (6)

the exact Solution $u\left(x\right)=\frac{x}{1+t},\left|t\right|<1$.

Applying the Adomian decomposition method using Maple18.

Example 2. Consider the nonlinear partial differential equation.

${u}_{t}+u{u}_{x}=1+t\cdot \mathrm{cos}x+\mathrm{sin}x\mathrm{cos}x,\text{\hspace{0.17em}}\text{\hspace{0.17em}}u\left(x,0\right)=\mathrm{sin}x$ (7)

the exact Solution $u\left(x\right)=t+\mathrm{sin}x$.

Applying the Adomian Decomposition Method using Maple18.

Example 3. Consider the nonlinear partial differential equation

${u}_{t}={x}^{2}+\frac{1}{4}{u}_{x}^{2},\text{\hspace{0.17em}}\text{\hspace{0.17em}}u\left(x,0\right)=0$ (8)

the exact Solution $u\left(x\right)={x}^{2}\mathrm{tan}t$.

Applying the Adomian Decomposition Method using Maple18.

Example 4. Consider the nonlinear partial differential equation

${u}_{t}+\frac{1}{36}x{u}_{xx}^{2}={x}^{3},\text{\hspace{0.17em}}\text{\hspace{0.17em}}u\left(x,0\right)=0,$ (9)

the exact Solution $u\left(x,t\right)={x}^{3}\mathrm{tanh}t$.

Applying the Adomian Decomposition Method using Maple.

Figure 1. Plot of the solutions of nonlinear partial differential equation for example 1.

Figure 2. Plot of the solutions of nonlinear partial differential equation for example 2.

Figure 3. Plot of the solutions of nonlinear partial differential equation for example 3.

Table 1. Approximation solution and exact solution of nonlinear partial differential equation for example 1.

Table 2. Approximation solution and exact solution of nonlinear partial differential equation for example 2.

Table 3. Approximation solution and exact solution of Volterra Fredholm integral equations for example 3.

Figure 4. Plot of the solutions of of nonlinear partial differential equation for example 4.

Table 4. Approximation solution and exact solution of nonlinear partial differential equation for example 4.

4. Conclusion

In this paper, the Adomian decomposition method is applied to solve the nonlinear partial differential equation using Maple18 software. Results were obtained by tables and drawing with figures. The exact solution and numerical solution are shown in Tables 1-4 and Figures 1-4. By comparing the numerical results, we find that the numerical solution is widely applied to the precise solution, which proves the efficiency of the method used and the ability to obtain the numerical solution corresponding to the precise solution easily and conveniently using Maple 18 software. Moreover, high accuracy of the results obtained.

Acknowledgements

This project was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University. The authors, therefore, acknowledge with thanks DSR technical and financial support.

Conflicts of Interest

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

[1] |
Wazwa, A.-M. (1999) A Reliable Modification of Adomain Decomposition Method. Applied Mathematics and Computation, 102, 77-86. https://doi.org/10.1016/S0096-3003(98)10024-3 |

[2] |
Wazwa, A.-M. and El-Sayed, S.A. (2001) Anew Modification of the Adomain Decomposition Method for Linear and Nonlinear Operators. Applied Mathematics and Computation, 122, 393-405. https://doi.org/10.1016/S0096-3003(00)00060-6 |

[3] |
Vadasz, P. and Olek, S. (2000) Convergence and Accuracy of Adomian’s Decomposition Method for the Solution of Lorenz Equations. International Journal of Heat and Mass Transfer, 43, 1715-1734. https://doi.org/10.1016/S0017-9310(99)00260-4 |

[4] |
El-Tawil, M.A., Bahnasawi, A.A. and Abdel-Naby, A. (2004) Solving Riccati Differential Equation Using Adomian’s Decomposition Method. Applied Mathematics and Computation, 157, 503-514. https://doi.org/10.1016/j.amc.2003.08.049 |

[5] |
Ghosh, S., Roy, A. and Roy, D. (2007) An Adaptation of Adomian Decomposition for Numeric-Analytic Integration of Strongly Nonlinear and Chaotic Oscillators. Computer Methods in Applied Mechanics and Engineering, 196, 1133-1153. https://doi.org/10.1016/j.cma.2006.08.010 |

[6] | Ebadi, G. and Rashedi, S. (2010) The Extended Adomian Decomposition Method for Fourth Order Boundary Value Problems. Acta Universitatis Apulensis, 22, 65-78. |

[7] |
Tatari, M. and Dehghan, M. (2006) The Use of the Adomian Decomposition Method for Solving Multipoint Boundary Value Problems. Physica Scripta, 73, 672-676. https://doi.org/10.1088/0031-8949/73/6/023 |

[8] | Azreg-Ainou, M. (2009) A Developed New Algorithm for Evaluating Adomian Polynomials. Computer Modeling in Engineering & Sciences, 42, 1-18. |

[9] |
Duan, J.S. (2010) An Efficient Algorithm for the Multivariable Adomian Polynomials. Applied Mathematics and Computation, 217, 2456-2467. https://doi.org/10.1016/j.amc.2010.07.046 |

[10] |
Duan, J.S. (2011) Convenient Analytic Recurrence Algorithms for the Adomian Polynomials. Applied Mathematics and Computation, 217, 6337-6348. https://doi.org/10.1016/j.amc.2011.01.007 |

[11] | Duan, J.S., Rach, R., Baleanu, D. and Wazwaz, A.M. (2012) A Review of the Adomian Decomposition Method and Its Applications to Fractional Differential Equations. Commun. Frac. Calc., 3, 73-99. |

[12] |
Rach, R. (2012) A Bibliography of the Theory and Applications of the Adomian Decomposition Method, 1961-2011. Kybernetes, 41, 1087-1148. https://doi.org/10.1108/k.2012.06741gaa.008 |

[13] | Abdul-Majid, W. (2011) Linear and Nonlinear Integral Equations. Higher Education Press, Beijing. |

[14] | Linz, P. (1985) Analytical and Numerical Methods for Volterra Equations (Studies in Applied Mathematics, Series Number 7). SIAM, Philadelphia. |

[15] |
Doddrell, D.M., Forbes, L.K. and Crozier, S. (1997) Calculating Current Densities and Fields Produced by Shielded Magnetic Resonance Imaging Probes. SIAM Journal on Applied Mathematics, 57, 401-425. https://doi.org/10.1137/S0036139995283110 |

[16] | Dalal, A.M., Amani, Z.B. and Badreeh, M.G. (2014) Numerical Solution of Volterra Integral Equation of Second Kind Using Implicit Trapezoidal. Journal of Advances in Mathematics, 8, 1540-1553. |

[17] | Dalal, A.M. (2014) Adomian Decomposition Method of Fredholm Integral Equation of the Second Kind Using Maple. Journal of Advances in Mathematics, 9, 1868-1875. |

[18] | Dalal, A.M. (2014) Application of Adomian Decomposition Method for Solving of Fredholm Integral Equation of the Second Kind. European Journal of Science and Engineering, 9, 1-9. |

[19] | Dalal, A.M. (2016) Adomian Decomposition Method for Solving of Fredholm Integral Equation of the Second Kind Using Matlab. International Journal of GEOMATE, 11, 2830-2833. |

[20] |
Maturi, D.A. and Malaikah, H. (2018) Numerical Solution of System of Three Nonlinear Volterra Integral Equation Using Implicit Trapezoidal. Journal of Mathematics Research, 10, 44-58. https://doi.org/10.5539/jmr.v10n1p44 |

[21] |
Maturi, D.A. (2019) The Modified Decomposition Method for Solving Volterra Integral Equation of the Second Kind Using Maple. International Journal of GEOMATE, 17, 23-28. https://doi.org/10.21660/2019.62.4590 |

[22] |
Malaikah, H.H. (2020) The Adomian Decomposition Method for Solving Volterra-Fredholm Integral Equation Using Maple. Applied Mathematics, 11, 779-787. https://doi.org/10.4236/am.2020.118052 |

Journals Menu

Contact us

customer@scirp.org | |

+86 18163351462(WhatsApp) | |

1655362766 | |

Paper Publishing WeChat |

Copyright © 2023 by authors and Scientific Research Publishing Inc.

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.