Next Article in Journal
Asymptotic Diffusion Analysis of Retrial Queueing System M/M/1 with Impatient Customers, Collisions and Unreliable Servers
Previous Article in Journal
Chromatic Number of Fuzzy Graphs: Operations, Fuzzy Graph Coloring, and Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modified Adomian Method through Efficient Inverse Integral Operators to Solve Nonlinear Initial-Value Problems for Ordinary Differential Equations

Department of Mathematics, Faculty of Science, University of Jeddah, Jeddah P. O. Box 80327, Saudi Arabia
*
Author to whom correspondence should be addressed.
Axioms 2022, 11(12), 698; https://doi.org/10.3390/axioms11120698
Submission received: 20 July 2022 / Revised: 9 November 2022 / Accepted: 28 November 2022 / Published: 6 December 2022

Abstract

:
The present manuscript examines different forms of Initial-Value Problems (IVPs) featuring various types of Ordinary Differential Equations (ODEs) by proposing a proficient modification to the famous standard Adomian decomposition method (ADM). The present paper collected different forms of inverse integral operators and further successfully demonstrated their applicability on dissimilar nonlinear singular and nonsingular ODEs. Furthermore, we surveyed most cases in this very new method, and it was found to have a fast convergence rate and, on the other hand, have high precision whenever exact analytical solutions are reachable.

1. Introduction

Many real-life problems are modeled mathematically using Initial-Value Problems (IVPs) featuring various types of Ordinary Differential Equations (ODEs). Solutions to these problems are very important as they improve human lives. Thus, this importance is what necessitates various researchers to devise different mathematical methods to tackle countless models arising in, for instance, physics, chemistry, engineering, biology, economics, and social sciences, to mention a few. In particular, the literature is full of diverse studies of the competent method called the Adomian decomposition method (ADM) and its various modifications and extensions [1,2,3]. In [4,5] some modifications of the Adomian decomposition method are presented for solving initial value problems in Ordinary Differential Equations. However, in [6], a dependable semi-analytical method via the application of a modified Adomian decomposition method (ADM) to tackle the coupled system of Emden–Fowler-type equations has been proposed, and an effective differential operator together with its corresponding inverse is successfully constructed. The present study [7] investigates certain singular Initial-Value Problems (IVPs) featuring the classical and generalized inhomogeneous LaneEmden-type equations. This study proposes different forms of inverse integral operators that are based on the Adomian method to accelerate the convergence rate of the standard Adomian decomposition method (ADM), which includes some cases from the survey that we present in this work. This method and its variants have been comprehensively utilized to treat different forms of linear and nonlinear ODEs, including integral equations and together with the combination of the two [8,9,10,11].
However, it is the aim of the present study to examine different forms of IVPs portraying different types of ODEs by proposing a proficient modification to the famous standards of ADM. The method devises different forms of inverse integral operators based on the available literature and further successfully demonstrates their applicability to a class of ODEs of physical relevance. Furthermore, we will assess this new method by establishing a comparative examination with the standard ADM and, on the other hand, with the exact analytical solutions whenever they are reachable.

2. Standard Adomian Decomposition Method

To present the standard ADM methodology, we take into consideration the following generalized ODE
L w ( t ) + R w ( t ) + F w ( t ) = g ( t ) ,
where L and R are linear operators with R < L , L is the highest linear operator, and R is an operator with a degree less than L, while F is a nonlinear operator from a Hilbert H. g ( t ) is a given function in H, and we are looking for w H satisfying (1). We assume that (1) has a unique solution for g H [12], and g ( t ) is an inhomogeneous or source term. Next, we rewrite the above equation as follows
L w ( t ) = g ( t ) R w ( t ) F w ( t ) ,
such that when the inverse operator L 1 of L is applied to both sides of the later equation it yields
w ( t ) = ψ ( t ) + L 1 g ( t ) L 1 R w ( t ) L 1 F w ( t ) ,
where function ψ ( t ) emanates from the prescribed initial data.
Therefore, the ADM decomposes the solution w ( x ) and the nonlinear term F ( w ) as series forms as follows
w ( t ) = n = 0 w n ( t ) , F w ( t ) = n = 0 A n ,
where A n is the Adomian polynomials that are recurrently computed using the following relation [13,14]
A n = 1 n ! d n d λ n F i = 0 n λ i w i λ = 0 , n = 0 , 1 , 2 , . . .
Additionally, from the above Adomian polynomials relation, we express some of these components as follows
A 0 = F w 0 , A 1 = F w 0 w 1 , A 2 = F w 0 w 2 + 1 2 F w 0 w 1 2 , A 3 = F w 0 w 3 + F w 0 w 1 w 2 + 1 3 ! F w 0 w 1 3 ,
Therefore, substituting Equation (4) into Equation (3) yields the following
n = 0 w n ( t ) = ψ ( t ) + L 1 g ( t ) L 1 R n = 0 w n ( t ) L 1 n = 0 A n .
Hence, from the above equation, the recurrent solution is thus obtained via the ADM process as follows
w 0 = ψ ( t ) + L 1 g ( t ) , w k + 1 = L 1 R w k L 1 A k , k 0 ,
such that the n-term approximation takes the form
Φ n = k = 0 n 1 w k ,
where the closed-form solution is finally revealed as
w ( t ) = lim n Φ n ( t ) = k = 0 w k ( t ) .
The theoretical treatment of the convergence of ADM has been considered in [12,15,16]. Cherruault [12] has given the first proof of convergence of the Adomian decomposition method, and he used fixed-point theorems for abstract functional equations. Abbaoui and Cherruault [15] have given new formulae that easily calculate the Adomian’s polynomials used in the decomposition methods. A simple proof of convergence of Adomian’s technique is presented in [16].

3. Modified Adomian Decomposition Method

Let us now present a modification methodology based on the standard ADM to solve certain classes of nonsingular and singular ODEs featuring IVPs, including, for instance, low- and high-order and systems of inhomogeneous ODEs.

3.1. First-Order IVPs

To present this procedure on ODEs of the first order, let us consider the following first-order IVP [17]
w + p ( t ) w + F ( t , w ) = g ( t ) , w ( 0 ) = A ,
where g ( t ) and p ( t ) are given functions, F ( t , w ) is the general nonlinear real term, and A is a supplied real constant.
Next, we rewrite Equation (11) using the operator denotation as follows
L w = g ( t ) F ( t , w ) ,
where L is the linear differential operator, together with its corresponding one-fold inverse integral operator L 1 defined as follows
L ( . ) = e p ( t ) d t d d t e p ( t ) d t ( . ) , L 1 ( . ) = e p ( t ) d t 0 t e p ( t ) d t ( . ) d t .
Now, we apply the inverse operator L 1 expressed in Equation (13) to the first two terms of Equation (11) as follows
L 1 ( w + p ( t ) w ) = e p ( t ) d t 0 t e p ( t ) d t ( w + p ( t ) w ) d t ,
= w w ( 0 ) ψ ( 0 ) e p ( t ) d t ,
where ψ ( t ) = e p ( t ) d t . Furthermore, applying the inverse operator L 1 to Equation (12) yields the following
w ( t ) = w ( 0 ) ψ ( 0 ) e p ( t ) d t + L 1 g ( t ) L 1 F ( t , w ) .
Therefore, we make use of the modification of ADM by first decomposing w ( t ) and F ( t , w ) as suggested in Equation (4), and then we obtain
n = 0 w n ( t ) = w ( 0 ) ψ ( 0 ) e p ( t ) d t + L 1 g ( t ) L 1 n = 0 A n ,
such that the overall recursive relation is acquired as follows
w 0 = w ( 0 ) ψ ( 0 ) e p ( t ) d t + L 1 g ( t ) , w k + 1 = L 1 A k , k 0 .
Lastly, on using Equation (18) via Equation (4), the closed-form solution is finally revealed for computational purposes as follows
Φ n = k = 0 n 1 w k ,
where the closed-form solution is finally revealed as
w ( t ) = lim n Φ n ( t ) = k = 0 w k ( t ) .

3.2. Second-Order IVPs

To present an efficient method based on the standard ADM to solve IVPs featuring singular ODEs of the second-order we refer to the well-known modification of ADM as suggested in [10]. In doing so, we take into consideration the following generalized second-order IVP
w + p ( t ) w + F ( t , w ) = g ( t ) , w ( 0 ) = A 1 , w ( 0 ) = A 2 ,
where g ( t ) and p ( t ) are given functions, F ( t , w ) is the general nonlinear real term, and A 1 and A 2 are supplied real constants.
What is more, we rewrite the ODE given in Equation (19) using the operator denotation as given in Equation (12) and further employ the following differential linear operator L together with its corresponding two-fold integral inverse L 1 as [10]
L ( . ) = e p ( t ) d t d d t e p ( t ) d t d ( . ) d t , L 1 ( . ) = 0 t e p ( t ) d t 0 t e p ( t ) d t ( . ) d t d t .
Therefore, applying the inverse operator L 1 given above to the resulting operator equation gives
w ( t ) = ψ ( t ) + L 1 g ( t ) L 1 F ( t , w ) ,
such that
L ψ ( t ) = 0 .
Hence, on decomposing the solution w ( t ) and the nonlinear term F ( t , w ) via infinite series earlier defined in Equation (4 ), the recurrent solution is thus given as follows
w 0 = ψ ( t ) + L 1 g ( t ) , w n + 1 = L 1 A n , n 0 ,
such that the n-term approximation takes the form
Φ n = k = 0 n 1 w k ,
where the closed-form solution is finally revealed as
w ( t ) = lim n Φ n ( t ) = k = 0 w k ( t ) .

3.3. Second-Order Singular IVPs

More importantly, we mention here that the method presented in the above subsection for the second-order IVPs was generalized by Hosseini and Jafari [18] for singular IVPs. This generalization is very powerful as it tackles different forms of second-order IVPs, including, for instance, linear, nonlinear, singular, and nonsingular ODEs.
However, considering a nonlinear singular second-order IVP of the form given in Equation (19), we suppose that the function p ( t ) is of the following singular form
p ( t ) = 1 t a h ( t ) ,
where Taylor’s series expansion of h ( t ) exists at t = a . Now, having already considered the differential linear operator and its inverse in Equation (20) based on the suggestion in [10], it will be very difficult to obtain a closed-form solution in the presence of such a singularity in the above equation. Thus, it is pertinent to make use of polynomials to approximate e p ( t ) d t and e p ( t ) d t in order to swiftly obtain the components w i ’s. Therefore, we further obtain Taylor’s series expansion of h ( t ) at t = a (for m N ) and re-express Equation (23) as follows
p ( t ) = 1 t a k = 0 m ( t a ) k k ! h ( k ) ( a ) .
Thus, we have
p ( t ) d t = ln ( t a ) h ( a ) + ( t a ) h ( a ) + . . . + ( t a ) m m × m ! h ( m ) ( a ) ,
and
e p ( t ) d t = ( t a ) h ( a ) S ( t ) ,
where
S ( t ) = e ( t a ) h ( a ) + + ( t a ) m m × m ! h ( m ) ( a ) .
Additionally, for any v N , we substitute Taylor’s series expansion of S ( t ) into (26) to yield
e p ( t ) d t = ( t a ) h ( a ) S ( a ) + ( t a ) S ( a ) + + ( t a ) v v ! S ( v ) ( a ) ,
such that in the same manner, we obtain
e p ( t ) d t = ( t a ) h ( a ) S ¯ ( a ) + ( t a ) S ¯ ( a ) + + ( t a ) v v ! S ¯ ( v ) ( a ) ,
where
S ¯ ( t ) = e ( t a ) h ( a ) ( t a ) m m × m ! h ( m ) ( a ) .
Finally, the difficulty associated with the singular function p ( t ) with regards to the operators given in Equation (20) is thus solved in line with the present development presented above. Thus, the recurrent solution follows by easily computing the components w i ’s.

3.4. Higher-Order IVPs

As higher-order IVPs arise in many real-life applications, we present here a promising technique based on the standard ADM to solve higher-order IVPs, as asserted in [19]. Thus, we take into consideration the following generalized n-order IVP [19]
w ( n ) + p ( t ) w ( n 1 ) + F ( w ) = g ( t ) , w ( 0 ) = β 0 , w ( 0 ) = β 1 , , w ( n 1 ) ( 0 ) = β n 1 ,
where g ( t ) and p ( t ) are given functions, F is a nonlinear differential operator of the order less than ( n 1 ) , and β 0 , β 1 , , β n 1 are prescribed real constants.
Furthermore, we equally express the ODE in the above system using differential operator denotation as follows
L w = g ( t ) F ( w ) ,
such that the differential operator L and its corresponding n-fold inverse integral operator L 1 are defined by
L ( . ) = e p ( t ) d t d d t e p ( t ) d t d n 1 ( . ) d t n 1 , L 1 ( . ) = 0 t 0 t 0 t e p ( t ) d t 0 t e p ( t ) d t ( . ) d t d t .
Thus, the application L 1 on Equation (31) transforms the equation to the following
w ( t ) = ψ ( t ) + L 1 g ( t ) L 1 F ( w ) ,
such that
L ψ ( t ) = 0 .
As we proceed through the use of the ADM procedure, we receive the following equation
n = 0 w n = ψ ( t ) + L 1 g ( t ) L 1 n = 0 A n ,
which yields the following recurrent solution
w 0 = ψ ( t ) + L 1 g ( t ) , w k + 1 = L 1 A k , k 0 ,
and a closed-form solution of
w ( t ) = lim n Φ n ( t ) = k = 0 w k ( t ) , Φ n = k = 0 n 1 w k .

3.5. Nonlinear System of IVPs

Let us take into consideration the following generalized system of nonlinear IVPs of ODEs,
w 1 ( n ) + p ( t ) w 1 ( n 1 ) + F 1 ( t , w 1 , , w 1 ( n 2 ) , w 2 , , w 2 ( n 2 ) , w n , , w n ( n 2 ) ) = g 1 ( t ) , w 2 ( n ) + p ( t ) w 2 ( n 1 ) + F 2 ( t , w 1 , , w 1 ( n 2 ) , w 2 , , w 2 ( n 2 ) , w n , , w n ( n 2 ) ) = g 2 ( t ) , w n ( n ) + p ( t ) w n ( n 1 ) + F n ( t , w 1 , , w 1 ( n 2 ) , w 2 , , w 2 ( n 2 ) , w n , , w n ( n 2 ) ) = g n ( t ) , w 1 ( 0 ) = β 1 , w 2 ( 0 ) = β 2 , , w n ( 0 ) = β n ,
where g ( t ) and p ( t ) are prescribed nice functions, F 1 , F 2 , , F n are nonlinear real functions, and β 1 , β 2 , , β n are supplied real constants.
Thus, without a loss of generalization, the modification of ADM personated in the above subsections can be equally extended to successfully tackle the system of nonlinear IVPs given above. This is, of course, can be performed by suitably constructing a generalized differential operator, L, together with its corresponding n-fold-generalized integral operator.

4. Numerical Illustrations

The present section demonstrates the application of the proposed methods on a number of test problems featuring different forms of ODEs.
Example 1.
Let us consider the following inhomogeneous first-order nonlinear IVP [17]
w + 2 t w = 1 + t 2 + w 2 , w ( 0 ) = 1 .
Standard Adomian decomposition method
First, we define a differential operator L together with its corresponding one-fold inverse integral operator L 1 as follows
L = d d t , L 1 ( . ) = 0 t ( . ) d t .
Next, we express Equation (37) in operator form as follows
L w = 2 t w + 1 + t 2 + w 2 ,
such that after applying L 1 to both sides of Equation (39) yields
w = w ( 0 ) 2 L 1 ( t w ) + L 1 ( 1 + t 2 ) + L 1 ( w 2 ) .
Therefore, without a loss in generality, we obtain the following recurrent relation
w 0 = w ( 0 ) + L 1 ( 1 + t 2 ) , w n + 1 = 2 L 1 ( t w n ) + L 1 ( A n ) , n 0 ,
where A n is the Adomian polynomial corresponding to the nonlinear term w 2 with a few components as follows
A 0 = w 0 2 , A 1 = 2 w 0 w 1 , A 2 = w 1 2 + 2 w 0 w 2 , A 3 = 2 w 1 w 2 + 2 w 0 w 3 ,
Therefore, substituting the above polynomial components into the recurrent relation determined in Equation (41) gives
w 0 = w ( 0 ) + L 1 ( 1 + t 2 ) = 1 + t + 1 3 t 3 , w 1 = 2 L 1 ( t w 0 ) + L 1 ( A 0 ) = t 1 3 t 3 + 1 6 t 4 + 1 63 t 7 , w 2 = 2 L 1 ( t w 1 ) + L 1 ( A 1 ) = t 2 1 6 t 4 + 1 5 t 5 2 63 t 7 + ,
Finally, from the above iterates, we obtain the following series solution
w ( t ) = 1 + 2 t + t 2 + t 3 + t 4 + t 5 + ,
whose closed-form solution is
w ( t ) = t + 1 1 t .
Modified Adomian decomposition method
Let us define a differential operator L together with its corresponding one-fold inverse integral operator L 1 as follows
L ( . ) = e t 2 d d t ( e t 2 ( . ) ) , L 1 ( . ) = e t 2 0 t e t 2 ( . ) d t .
Then, Equation (37) in operator form becomes
L w = 1 + t 2 + w 2 ,
such that after operating L 1 in the later equation reveals Equation (47)
w ( t ) = e t 2 + L 1 ( 1 + t 2 ) + L 1 ( w 2 ) ,
with
w 0 ( t ) = e t 2 + e t 2 0 t e t 2 ( 1 + t 2 ) d t .
Therefore, on making use of Taylor’s series expansion on e t 2 and e t 2 of order 6 , the following solution iterates are obtained
w 0 = 1 + t t 2 t 3 3 + t 4 2 + , w 1 = t + t 2 t 3 7 t 4 6 + , w 2 = t 2 + 4 t 3 3 t 4 + ,
Hence, we obtain the following series solution of the form given by
w ( t ) = 1 + 2 t + t 2 + t 3 + t 4 + t 5 + ,
whose closed-form solution is
w ( t ) = t + 1 1 t .
Example 2.
Let us consider the following inhomogeneous second-order linear singular IVP [10]
w + cos t sin t w = 2 cos t , w ( 0 ) = 1 , w ( 0 ) = 0 .
Standard Adomian decomposition method
We define a differential operator L together with its corresponding two-fold inverse integral operator L 1 as follows
L ( . ) = d 2 ( . ) d t 2 , L 1 ( . ) = 0 t 0 t ( . ) d t d t .
Expressing Equation (53) in an operator form becomes
L w = cos t sin t w 2 cos t ,
such that after taking L 1 of the later equation yields
w = w ( 0 ) + t w ( 0 ) L 1 cos t sin t w L 1 ( 2 cos t ) .
Accordingly, we obtain the following recurrent relation
{ w 0 = w ( 0 ) + t w ( 0 ) L 1 ( 2 cos t ) w n + 1 = L 1 ( cos t sin t w n ) , n 0 ,
where some of its iterates are expressed as follows
w 0 = 2 cos t 1 , w 1 = 2 cos t + 2 , w 2 = 2 cos t 2 , w 3 = 2 cos t + 2 ,
We, therefore, conclude from the above components that the standard ADM fails as the obtained series solution is divergent.
Modified Adomian decomposition method
Let us define a differential operator L together with its corresponding two-fold inverse integral operator L 1 as follows
L ( . ) = 1 sin t d ( . ) d t sin t d d t ( . ) , L 1 ( . ) = 0 t 1 sin t 0 t sin t ( . ) d t d t .
Expressing Equation (53) in an operator form becomes
L w = 2 cos t ,
while making use of L 1 on the above equations gives
L 1 L w = 2 0 t 1 sin t 0 t sin t ( cos t ) d t d t .
Without a loss in generality, the proposed modified ADM reveals the following exact solution
w ( t ) = w ( 0 ) + t w ( 0 ) + cos t 1 = cos t .
In fact, this shows the power of the proposed method over the standard ADM.
Example 3.
Let us consider the following inhomogeneous second-order nonlinear IVP [10]
w + t w + t 2 w 3 = ( 2 + 6 t 2 ) e t 2 + t 2 e 3 t 2 , w ( 0 ) = 1 , w ( 0 ) = 0 ,
with the exact solution
w ( t ) = e t 2 .
Standard Adomian decomposition method
Let us define a differential operator L together with its corresponding two-fold inverse integral operator L 1 as follows
L ( . ) = d 2 d t 2 , L 1 ( . ) = 0 t 0 t ( . ) d t d t .
Accordingly, we obtain the following recurrent relation
w 0 = w ( 0 ) + t w ( 0 ) + L 1 ( g ( t ) ) , w n + 1 = L 1 ( t w n ) L 1 ( A n ) , n 0 ,
where
g ( t ) = ( 2 + 6 t 2 ) e t 2 + t 2 e 3 t 2 ,
and the Adomian polynomials, A n , of nonlinear term t 2 w 3 are given as follows
A 0 = t 2 w 0 3 , A 1 = t 2 ( 3 w 0 2 w 1 ) , A 2 = t 2 ( 3 w 0 2 w 2 + 3 w 0 w 1 2 ) , A 3 = t 2 ( 3 w 0 2 w 3 + 6 w 0 w 1 w 2 + w 1 3 ) ,
We mention here that Taylor’s series of order 10 was utilized on g ( t ) for the computation of w 0 . Thus, the solution becomes
w 0 = 1 + t 2 + 3 4 t 4 + 1 3 t 6 + , w 0 + w 1 = 1 + t 2 + 1 2 t 4 + 2 15 t 6 + 1 96 t 8 + , w 0 + w 1 + w 2 = 1 + t 2 + 1 2 t 4 + 1 6 t 6 + 19 420 t 8 + , w 0 + w 1 + w 2 + w 3 = 1 + t 2 + 1 2 t 4 + 1 6 t 6 + 1 24 t 8 + 101 12600 t 10 + ,
Clearly, this solution converges to the exact solution w ( t ) = e t 2 , as Taylor’s series expansion of order 10 of e t 2 is expressed as
e t 2 = 1 + t 2 + 1 2 t 4 + 1 6 t 6 + 1 24 t 8 +
Modified Adomian decomposition method
Let us define a differential operator L together with its corresponding two-fold inverse integral operator L 1 as follows
L ( . ) = e t 2 / 2 d d t e t 2 / 2 d ( . ) d t , L 1 ( . ) = 0 t e t 2 / 2 0 t e t 2 / 2 ( . ) d t d t .
Without a loss in generality, we obtain the following recurrent relation
{ w 0 = w ( 0 ) + t w ( 0 ) + L 1 ( g ( t ) ) , w n + 1 = L 1 ( A n ) , n 0 .
Further, making use of Taylor’s series expansion of order 10 on g ( t ) , e t 2 / 2 and e t 2 / 2 , and coupling with obtaining the Adomian polynomials of the given nonlinearity terms in the original equation expressed in Equation (63), we receive the following solution
w 0 = 1 + t 2 + 7 12 t 4 + 23 90 t 6 + , w 0 + w 1 = 1 + t 2 + 1 2 t 4 + 1 6 t 6 + 25 672 t 8 + , w 0 + w 1 + w 2 = 1 + t 2 + 1 2 t 4 + 1 6 t 6 + 1 24 t 8 + 1 120 t 10 + ,
It is also obvious that the proposed modified ADM converges faster than the standard ADM; this can clearly be observed by comparing the two solutions.
Example 4.
Consider the following inhomogeneous second-order nonlinear singular IVP [18]
w + 1 + 5 t 2 t ( t + 1 ) w + w 3 = 5 t + 11 t 2 + 8 27 t 9 + 2 3 t 10 + 1 2 t 11 + 1 8 t 12 , w ( 0 ) = 0 = w ( 0 ) ,
which admits the following exact solution
w ( t ) = 2 3 t 3 + 1 2 t 4 .
Standard Adomian decomposition method
Let us define a differential operator L together with its corresponding two-fold inverse integral operator L 1 as follows
L ( . ) = d 2 ( . ) d t 2 , L 1 ( . ) = 0 t 0 t ( . ) d t d t .
As preceded, the following recurrent relation is obtained
w 0 = w ( 0 ) + t w ( 0 ) + L 1 ( g ( t ) ) , w n + 1 = L 1 ( f ( t ) w n ) L 1 ( A n ) , n 0 ,
where
g ( t ) = 5 t + 11 t 2 + 8 27 t 9 + 2 3 t 10 + 1 2 t 11 + 1 8 t 12 ,
f ( t ) is determined through the application of Taylor’s series by expanding 1 + 5 t 2 t ( t + 1 ) at t = 0 with order 9 to yield
f ( t ) = 1 t 1 2 + 2 t 2 t 2 + 2 t 3 + 2 t 8 ,
while the Adomian polynomials, A n , of the nonlinear term w 3 are expressed for some terms as follows
A 0 = w 0 3 , A 1 = 3 w 0 2 w 1 , A 2 = 3 w 0 2 w 2 + 3 w 0 w 1 2 , A 3 = 3 w 0 2 w 3 + 6 w 0 w 1 w 2 + w 1 3 ,
In this case, we obtain
w 0 = 5 6 t 3 + 11 12 t 4 + 4 1485 t 11 + w 0 + w 1 = 5 8 t 3 + 25 72 t 4 7 60 t 5 + w 0 + w 1 + w 2 + + w 5 = 1365 2048 t 3 + 372575 746496 t 4 41797 19906560 t 5 +
It is easy to see that the standard Adomian decomposition method converges to the exact solution (72) very slowly.
Modified Adomian decomposition method
We consider the following differential operator L together with its corresponding two-fold inverse integral operator L 1 as follows
L ( . ) = e p ( t ) d t d d t e p ( t ) d t d ( . ) d t , L 1 ( . ) = 0 t e p ( t ) d t 0 t e p ( t ) d t ( . ) d t d t ,
such that
p ( t ) d t = 1 + 5 t 2 t ( t + 1 ) d t = ln t ( 1 + t ) 2 ,
and
e p ( t ) d t = t ( 1 + t ) 2 = t 1 2 + 2 t 3 2 + t 5 2 , e p ( t ) d t = 1 t ( 1 + t ) 2 .
Therefore, with the application of Taylor’s series expansion on 1 ( 1 + t ) 2 of order 9 with regards to Equation (80), we acquire
e p ( t ) d t = t 1 2 1 2 t + 3 t 2 + + 9 t 8 ,
of which the recurrent relation is finally obtained as follows
w 0 = L 1 ( g ( t ) ) , w n + 1 = L 1 A n , n 0 .
Thus, through substituting Equations (79) and (80) into Equation (78), we have
w 0 = 2 3 t 3 + 1 2 t 4 + 16 6237 t 11 + , w 0 + w 1 = 2 3 t 3 + 1 2 t 4 + O ( t 12 ) ,
where w w 0 + w 1 is pretty close to the exact analytical solution earlier stated.
Example 5.
Let us consider the following inhomogeneous third-order nonlinear IVP [19]
w + e t w + 4 t 2 w + t 2 w 3 = g ( t ) , w ( 0 ) = w ( 0 ) = w ( 0 ) = 0 ,
where g ( t ) is compatible with the following exact solution
w ( t ) = t 3 e t .
What is more, expressing function g ( t ) using Taylor’s series expansion of order 9 yields
g ( t ) = 6 + 30 t + 48 t 2 + 45 t 3 + 171 4 t 4 + 164 5 t 5 + 529 30 t 6 + 243 35 t 7 + 2881 1344 t 8 .
Standard Adomian decomposition method
We consider the following differential operator L together with its corresponding three-fold inverse integral operator L 1 as follows
L ( . ) = d 3 ( . ) d t 3 , L 1 ( . ) = 0 t 0 t 0 t ( . ) d t d t d t .
As preceded, the following recurrent relation is obtained
w 0 = L 1 ( g ( t ) ) , w n + 1 = L 1 ( f ( t ) w n ) 4 L 1 ( t 2 w n ) L 1 ( t 2 A n ) , n 0 ,
where f ( t ) is determined through the application of Taylor’s series by expanding e t at t = 0 as follows
f ( t ) 1 + t + t 2 2 + + t 8 8 ! ,
while the Adomian polynomials, A n , corresponding to the nonlinear term w 3 are given for some components as follows
A 0 = w 0 3 , A 1 = 3 w 1 w 0 2 , A 2 = 3 w 2 w 0 2 + 3 w 0 w 1 2 ,
Hence, substituting Equations (88) and (89) into Equation (87) gives
w 0 = t 3 + 5 4 t 4 + , w 0 + w 1 = t 3 + t 4 + , w 0 + w 1 + w 2 + + w 7 = t 3 + t 4 + 1 2 t 5 + 1 6 t 6 + ,
Additionally, we further affirm the obtained series solution in the above equation by applying Taylor’s series expansion of order 9 to the exact solution as follows
w ( t ) = t 3 + t 4 + 1 2 t 5 + 1 6 t 6 +
Certainly, the obtained series solution gradually progresses to the exact closed-form solution, but slowly. Thus, the convergence rate of standard ADM is slow; this can clearly be seen in the proposed scheme.
Modified Adomian decomposition method
Accordingly, we consider the following differential operator L together with its corresponding three-fold inverse integral operator L 1 as follows
L ( . ) = e p ( t ) d t d d t e p ( t ) d t d 2 ( . ) d t 2 , L 1 ( . ) = 0 t 0 t e p ( t ) d t 0 t e p ( t ) d t ( . ) d t d t d t
such that
p ( t ) d t = e t d t = e t
and
e p ( t ) d t = e e t , e p ( t ) d t = e e t .
Therefore, with the application of Taylor’s series expansion of order 9 on e e t and e e t with regards to Equation (92), we obtain
a = e e t = e 1 1 t + 1 6 t 3 + + 5 4032 t 8 , b = e e t = e 1 + t + t 2 + 5 6 t 3 + + 23 224 t 8 .
In addition, we rewrite the above inverse operator in terms of a and b as follows
L 1 ( . ) = 0 t 0 t ( a ) 0 t ( b ) ( . ) d t d t d t .
Thus, the resulting recurrent relation is obtained based on Equation (35) as follows
w 0 = L 1 g ( t ) , w n + 1 = L 1 ( t 2 A n ) , n 0 ,
such that some of its component sums are as follows
w 0 = t 3 + t 4 + 1 2 t 5 + 1 6 t 6 + 83 840 t 7 + , w 0 + w 1 = t 3 + t 4 + 1 2 t 5 + 1 6 t 6 +
Therefore, it is obvious that the proposed modified ADM converges faster than the standard ADM, as earlier demonstrated; this can clearly be seen by comparing the two series solutions in Equations (90) and (98), respectively. The modified ADM solution of Φ 2 equals the standard ADM solution of Φ 8 ; in fact, the convergence rate of this method is higher by far.
Example 6.
Consider the following nonlinear system of inhomogeneous second-order IVPs [19]
u + tan ( t ) u + v 2 = g ( t ) , u ( 0 ) = 0 , u ( 0 ) = 0 , v + 100 v + u 2 = h ( t ) , v ( 0 ) = 0 , v ( 0 ) = 0 ,
where g ( t ) and h ( t ) are compatible with the following exact solution set
u ( t ) = t sin t , v ( t ) = t tan t .
More so, expressing functions g ( t ) and h ( t ) using Taylor’s series expansion order 9 yields
g ( t ) = 2 + 5 4 t 4 + 3 4 t 6 + , h ( t ) = 2 + 200 t + 4 t 2 + 400 3 t 3 + 5 t 4 80 t 5 + 121 45 t 6 + .
Standard Adomian decomposition method
Without a loss in generality, the system admits the following recurrent relation
u 0 = L 1 ( g ( t ) ) = t 2 + 1 24 t 6 + 3 224 t 8 , v 0 = L 1 ( h ( t ) ) = t 2 + 100 3 t 3 + 1 3 t 4 + 20 3 t 5 + 1 6 t 6 40 21 t 7 + 121 2520 t 8 + ,
and
u n + 1 = L 1 ( f ( t ) u n ) L 1 A n , n 0 , v n + 1 = L 1 ( 100 v n ) L 1 B n , n 0 ,
where A n and B n are the Adomian polynomials corresponding to the nonlinear terms v 2 and u 2 , respectively. Additionally, function f ( t ) represents Taylor’s series expansion of tan t of order 9. Thus, by considering Equation (101) and (102), we obtain
u 0 = t 2 + 1 24 t 6 + 3 224 t 8 u 0 + u 1 = t 2 1 6 t 4 1 72 t 6 + , u 0 + u 1 + u 2 + + u 6 = t 2 1 6 t 4 + 1 120 t 6 + ,
and
v 0 = t 2 + 100 3 t 3 + 1 3 t 4 + , v 0 + v 1 = t 2 883 t 4 4994 45 t 6 + , v 0 + v 1 + v 2 + + v 6 = t 2 + 1 3 t 4 + 2 15 t 6 +
Therefore, the obtained standard ADM solution converges to the exact solution; this could be seen clearly by expanding the exact solution using Taylor’s expansion and thereafter comparing the two solutions.
Modified Adomian decomposition method
Accordingly, we consider the following differential operator L together with its corresponding two-fold inverse integral operator L 1 as follows
L ( . ) = e p ( t ) d t d d t e p ( t ) d t d ( . ) d t , L 1 ( . ) = 0 t e p ( t ) d t 0 t e p ( t ) d t ( . ) d t d t
such that
e p ( t ) d t = cos ( t ) , e p ( t ) d t = 1 cos ( t ) .
Therefore, with the application of Taylor’s series expansion of order 9 on cos ( t ) and 1 cos ( t ) with regards to Equation (105), we obtain
a = cos ( t ) = 1 1 2 t 2 + 1 24 t 4 + + 1 40320 t 8 , b = 1 cos ( t ) = 1 + 1 2 t 2 + 5 24 t 4 + + 277 8064 t 8 .
In addition, we rewrite the above inverse operator in terms of a and b as follows
L 1 ( . ) = 0 t ( a ) 0 t ( b ) ( . ) d t d t .
Thus, the resulting recurrent relation is obtained as follows
u 0 = L 1 ( g ( t ) ) = t 2 1 6 t 4 + 1 24 t 6 + , u n + 1 = L 1 A n , n 0 ,
Now, we consider the second differential operator L together with its corresponding two-fold inverse integral operator L 1 as follows
L ( . ) = e p ( t ) d t d d t e p ( t ) d t d ( . ) d t , L 1 ( . ) = 0 t e p ( t ) d t 0 t e p ( t ) d t ( . ) d t d t
such that
e p ( t ) d t = e 100 d t = e 100 t , e p ( t ) d t = e 100 d t = e 100 t .
Therefore, with the application of Taylor’s series expansion of order 9 on e 100 t and e 100 t with regards to Equation (105), we obtain
c = e 100 t = 1 100 t + 5000 t 2 + + 15625000000000 63 t 8 , d = e 100 t = 1 + 100 t + 5000 t 2 + + 15625000000000 63 t 8 .
In addition, we rewrite the above inverse operator in terms of c and d as follows
L 1 ( . ) = 0 t ( c ) 0 t ( d ) ( . ) d t d t .
Thus, the resulting recurrent relation is obtained as follows
v 0 = L 1 ( h ( t ) ) = t 2 + 1 3 t 4 + 1 6 t 6 + v n + 1 = L 1 B n , n 0 ,
such that some components are expressed as follows
u 0 = t 2 1 6 t 4 + 1 24 t 6 + , u 0 + u 1 = t 2 1 6 t 4 + 1 120 t 6 + ,
and
v 0 = t 2 + 1 3 t 4 + 1 6 t 6 + , v 0 + v 1 = t 2 + 1 3 t 4 + 2 15 t 6 + ,
where the present method also outperforms the standard ADM. This could undoubtedly be noted by comparing the two solutions where the modified ADM solution of Φ 2 matches the standard ADM solution of Φ 7 . This further affirms the higher convergence rate of the proposed method.

5. Conclusions

In conclusion, the present manuscript examined various forms of IVPs by proposing a proficient modification to the famous standard ADM. The proposed method collected different forms of inverse integral operators and successfully applied them to dissimilar inhomogeneous nonlinear singular and nonsingular ODEs. The efficiency of the method was further assessed, taking into account its faster convergence rate and, on the other hand, its higher precision with the available exact analytical solutions. Thus, we finally recommend that the proposed method should be utilized to solve physical models arising in different nonlinear sciences.

Author Contributions

Data curation, A.A.A.; Investigation, A.A.A.; Methodology, M.A.-M.; Project administration, M.A.-M.; Software, N.A.; Validation, N.A.;Writing—review editing, H.O.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Adomian, G. A review of the decomposition method and some recent results for nonlinear equations. Math. Comput. Model. 1990, 13, 17–43. [Google Scholar] [CrossRef]
  2. Wazwaz, A.M. A reliable modification of Adomian decomposition method. Appl. Math. Comput. 1999, 102, 77–86. [Google Scholar] [CrossRef]
  3. Hosseini, M.M. Adomian decomposition method with Chebyshev polynomials. Appl. Math. Comput. 2006, 175, 1685–1693. [Google Scholar] [CrossRef]
  4. Almazmumy, M.; Hendi, F.; Bakodah, H.; Alzumi, H. Recent modifications of Adomian decomposition method for initial value problem in ordinary differential equations. Am. J. Comput. Math. 2012, 12, 23204. [Google Scholar] [CrossRef] [Green Version]
  5. Bakodah, H.; Banaja, M.; Alrigi, B.; Ebaid, A.; Rach, R. An efficient modification of the decomposition method with a convergence parameter for solving Korteweg de Vries equations. J. King Saud Univ.-Sci. 2019, 31, 1424–1430. [Google Scholar] [CrossRef]
  6. AL-Mazmumy, M.; Alsulami, A.; Bakodah, H.; Alzaid, N. Modified Adomian Method for the Generalized Inhomogeneous Lane-Emden-Type Equations. Nonlinear Anal. Differ. Equ 2022, 10, 15–35. [Google Scholar] [CrossRef]
  7. Alsulami, A.A.; AL-Mazmumy, M.; Bakodah, H.O.; Alzaid, N. A Method for the Solution of Coupled System of Emden–Fowler–Type Equations. Symmetry 2022, 14, 843. [Google Scholar] [CrossRef]
  8. Wazwaz, A.M. A new method for solving singular initial value problems in the second-order ordinary differential equations. Appl. Math. Comput. 2002, 128, 45–57. [Google Scholar] [CrossRef]
  9. Wazwaz, A.M. Adomian decomposition method for a reliable treatment of the Emden–Fowler equation. Appl. Math. Comput. 2005, 161, 543–560. [Google Scholar] [CrossRef]
  10. Hosseini, M.M.; Nasabzadeh, H. Modified Adomian decomposition method for specific second order ordinary differential equations. Appl. Math. Comput. 2007, 186, 117–123. [Google Scholar] [CrossRef]
  11. Hasan, Y.Q.; Zhu, L.M. Modified Adomian decomposition method for singular initial value problems in the second-order ordinary differential equations. Surv. Math. Its Appl. 2008, 3, 183–193. [Google Scholar]
  12. Cherruault, Y. Convergence of Adomian’s method. Math. Comput. Model. 1990, 14, 83–86. [Google Scholar] [CrossRef]
  13. Wazwaz, A. Partial Differential Equations and Solitary Waves Theory; Higher Education Press: Beijing, China; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  14. Wazwaz, A.M. A new algorithm for calculating Adomian polynomials for nonlinear operators. Appl. Math. Comput. 2000, 111, 33–51. [Google Scholar] [CrossRef]
  15. Abbaoui, K.; Cherruault, Y. New ideas for proving convergence of decomposition methods. Comput. Math. Appl. 1995, 29, 103–108. [Google Scholar] [CrossRef] [Green Version]
  16. Hosseini, M.M.; Nasabzadeh, H. On the convergence of Adomian decomposition method. Appl. Math. Comput. 2006, 182, 536–543. [Google Scholar] [CrossRef]
  17. Hasan, Y.Q. Solving first-order ordinary differential equations by Modified Adomian decomposition method. Adv. Intell. Transp. Syst. 2012, 1, 86–89. [Google Scholar]
  18. Hosseini, M.; Jafari, M. An efficient method for solving nonlinear singular initial value problems. Int. J. Comput. Math. 2009, 86, 1657–1666. [Google Scholar] [CrossRef]
  19. Hosseini, M.; Jafari, M. A note on the use of Adomian decomposition method for high-order and system of nonlinear differential equations. Commun. Nonlinear Sci. Numer. Simul. 2009, 14, 1952–1957. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

AL-Mazmumy, M.; Alsulami, A.A.; Bakodah, H.O.; Alzaid, N. Modified Adomian Method through Efficient Inverse Integral Operators to Solve Nonlinear Initial-Value Problems for Ordinary Differential Equations. Axioms 2022, 11, 698. https://doi.org/10.3390/axioms11120698

AMA Style

AL-Mazmumy M, Alsulami AA, Bakodah HO, Alzaid N. Modified Adomian Method through Efficient Inverse Integral Operators to Solve Nonlinear Initial-Value Problems for Ordinary Differential Equations. Axioms. 2022; 11(12):698. https://doi.org/10.3390/axioms11120698

Chicago/Turabian Style

AL-Mazmumy, Mariam, Aishah A. Alsulami, Huda O. Bakodah, and Nawal Alzaid. 2022. "Modified Adomian Method through Efficient Inverse Integral Operators to Solve Nonlinear Initial-Value Problems for Ordinary Differential Equations" Axioms 11, no. 12: 698. https://doi.org/10.3390/axioms11120698

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop