Symbolic Generation of Adomian Polynomials for Different Nonlinearities by Python

The Adomian decomposition method (ADM) is a powerful mathematical technique to find closed-form solutions to nonlinear functional equations including ODEs, PDEs, differential-difference, integral, integro-differential, algebraic, and transcendental equations or systems of such equations. It features a particular infinite series for the representation of nonlinear terms of the equation under study, referred to as the Adomian polynomials. Nevertheless, the computation of such polynomials manually, devoid of any assistance from computational resources, can often be a laborious and protracted endeavor. In this paper, an innovative Python code is proposed, which exploits the SymPy library to perform the involved symbolic calculus operations to generate the Adomian polynomials of any given nonlinear expressions. The use of the code would substantially facilitate the implementation of the ADM to the equations arising in various branches of science and engineering. A number of nonlinear expressions are decomposed to their relevant Adomian polynomials for the sake of demonstration.