Numerical resolution of Emden's equation using Adomian polynomials

Purpose: In this paper the authors aim to show the advantages of using the decomposition method introduced by Adomian to solve Emden's equation, a classical non‐linear equation that appears in the study of the thermal behaviour of a spherical cloud and of the gravitational potential of a polytropic fluid at hydrostatic equilibrium. Design/methodology/approach: In their work, the authors first review Emden's equation and its possible solutions using the Frobenius and power series methods; then, Adomian polynomials are introduced. Afterwards, Emden's equation is solved using Adomian's decomposition method and, finally, they conclude with a comparison of the solution given by Adomian's method with the solution obtained by the other methods, for certain cases where the exact solution is known. Findings: Solving Emden's equation for n in the interval [0, 5] is very interesting for several scientific applications, such as astronomy. However, the exact solution is known only for n=0, n=1 and n=5. The experiments show that Adomian's method achieves an approximate solution which overlaps with the exact solution when n=0, and that coincides with the Taylor expansion of the exact solutions for n=1 and n=5. As a result, the authors obtained quite satisfactory results from their proposal. Originality/value: The main classical methods for obtaining approximate solutions of Emden's equation have serious computational drawbacks. The authors make a new, efficient numerical implementation for solving this equation, constructing iteratively the Adomian polynomials, which leads to a solution of Emden's equation that extends the range of variation of parameter n compared to the solutions given by both the Frobenius and the power series methods.

Assessment of rural nonpoint source pollution : a model based on the Universal soil loss equation /

Mode of access: Internet.

The development of a numerical solution to transport equation /

"Prepared for U.S. Army Engineer District, Mobile."

Crash padding research. Volume II: constitutive equation models. Final report.

National Highway Traffic Safety Administration, Office of Research and Development, Washington, D.C.

An algorithm for the solution of a quadratic equation using continued fractions /

Thesis (M. S.)--University of Illinois at Urbana-Champaign, 1972.

Iterative and direct methods for solving Poisson's equation and their adaptability to ILLIAC IV,

Includes bibliography.

A multirate ordinary differential equation integrator /

Originally presented as the author's thesis, University of Illinois at Urbana-Champaign.

Generation and comparison of equivalent equation sets in a general purpose simulation and modeling package /

"C00-2383-0018."

Automatic integration of the heat equation /

Thesis (M. Sc.)--University of Illinois at Urbana-Champaign.

The finite range Wiener-Hopf integral equation and a boundary value problem in a waveguide /

"Contract No. AF33(616)-6079 Project No. 9-(13-6278) Task 40572. Sponsored by: Wright Air Development Center"