Volume 2012 (2012), Article ID 541083, 22 pages
http://dx.doi.org/10.1155/2012/541083
Research Article

## New Approach for Solving a Class of Doubly Singular Two-Point Boundary Value Problems Using Adomian Decomposition Method

Department of Mathematics, Indian Institute of Technology Kharagpur, Kharagpur 721302, India

Received 13 October 2012; Revised 3 December 2012; Accepted 5 December 2012

Copyright © 2012 Randhir Singh et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We propose two new modified recursive schemes for solving a class of doubly singular two-point boundary value problems. These schemes are based on Adomian decomposition method (ADM) and new proposed integral operators. We use all the boundary conditions to derive an integral equation before establishing the recursive schemes for the solution components. Thus we develop recursive schemes without any undetermined coefficients while computing successive solution components, whereas several previous recursive schemes have done so. This modification also avoids solving a sequence of nonlinear algebraic or transcendental equations for the undetermined coefficients with multiple roots, which is required to complete calculation of the solution by several earlier modified recursion schemes using the ADM. The approximate solution is computed in the form of series with easily calculable components. The effectiveness of the proposed approach is tested by considering four examples and results are compared with previous known results.

#### 1. Introduction

Consider the following class of doubly singular two-point boundary value problems: with boundary conditions where , , , , and are any finite constants. The condition characterizes that problem (1.1) is singular and in addition to this if is allowed to be discontinuous at , then problem (1.1) is called doubly singular Bobisud [1]. Consider problem (1.1) with the following conditions on , , and .

Type 1. Dirichlet boundary conditions: , , with in and ; in and ; (i) ; (ii)   exists and continuous, for all and all real .

Type 2. Mixed type boundary conditions: ,  ,   with in ; in (0,1], and ;  the same as ().

In recent years, the study of such singular boundary value problems (1.1) has attracted the attention of several researchers [29]. In particular, if ,   and problem (1.1) with Type 1 boundary conditions is known as Thomas-Fermi equation. Thomas [10] and Fermi [11] independently derived a boundary value problem for determining the electrical potential in an atom. The analysis leads to the nonlinear SBVP with boundary conditions given by ,  .

Chan and Hon [12] considered the generalized Thomas-Fermi equation: , , with , where , , , , which is doubly singular problem. Problem (1.1) with where arises in the study of the distribution of heat sources in the human head [13] with , , .

There is a huge literature available on numerical methods for problem (1.1) with , but very few numerical methods are available to tackle doubly singular boundary value problems. Reddien [14] studied the linear form of problem (1.1) and derived numerical methods for which is stronger assumption than ().

Chawla and Shivakumar [15] established the existence as well as uniqueness of solution for problem (1.1) where . The existence and uniqueness for problem (1.1) have also been discussed by Dunninger and Kurtz [16] and Bobisud [1]. Later Pandey and Verma [17] extended the results on the existence uniqueness for problem (1.1) with Types 1 and 2 boundary conditions.

In this subsection, we briefly describe standard ADM for nonlinear second-order equation.

Recently, many researchers [4, 1826] have shown interest to the study of ADM for different scientific models. Adomian [18] asserted that the ADM provides an efficient and computationally suitable method for generating approximate series solution for a large class of differential equations. Let us consider nonlinear second-order ordinary differential equation where is the second-order linear derivative operator, is the linear remainder operator, represents the nonlinear term, and is a source term. The above equation can be rewritten as The inverse operator of is defined as .

Operating the inverse linear operator on both the sides of (1.4) yields Next, we decompose the solution and the nonlinear function by an infinite series as where are Adomian polynomials that can be constructed for various classes of nonlinear functions with the formula given by Adomian and Rach [19] Substituting the series (1.6) into (1.5), we obtain From (1.8), the various components of the solution can be determined by using the recursive relation For numerical purpose, the -term truncated series may be used to give the approximate solution.

The ADM has been used to solve nonlinear boundary value problems (BVPs) for ordinary differential equations by several researchers [4, 2026]. Solving nonlinear BVP by standard ADM is always a computationally involved task as it requires the computation of undetermined coefficients in (1.9) in a sequence of nonlinear equations which increases computational complexity.

It is important to note that the standard ADM (1.9) can not be applied directly to solve two-point boundary value problem (1.1) as the component in scheme (1.9) is not independent from undetermined coefficient. Many researchers [2325] have proposed modified ADM to overcome the difficulty by setting undermined coefficient , and then it will be determined by using second boundary condition satisfying . In this case, it requires additional computational work to solve nonlinear equation for and may not be uniquely determined.

Benabidallah and Cherruault [24] have considered the following boundary value problem: The inverse operator was proposed, then unidentified constant is set, as , and each solution component is obtained by the following scheme: In order to determine the unidentified constants , it is also required -term approximate solution satisfying all the boundary conditions. It is clearly that . For boundary condition , the following scheme was proposed However, since each constant is determined by solving (1.12), this scheme also require additional computational work.

In order to avoid solving such nonlinear algebraic equations for a two-point boundary value problems, Jang [26] introduced extended ADM for nonsingular problems with Dirichlet boundary condition. Khuri and Sayfy [27] applied a novel approach based on the mixed decomposition spline for solving singular problems arising in physiology. The method proposed by Ebaid [4] is based on the modification of Lesnic’s work [28]. Later, Duan and Rach [29] introduced special modified inverse linear integral operators for higher order boundary value problem.

To the best of our knowledge, no one has applied ADM to solve a class of doubly singular two-point boundary value problem (1.1) with Types 1 and 2 boundary conditions.

In this work, a new approach based on ADM and new integral operator are proposed for solving doubly singular two-point boundary value problems (1.1) with Types 1 and 2. To set up the modified scheme, we first consider DSBVP (1.1) and use all the boundary conditions and Adomian decomposition method to establish the recursion scheme. Thus a modified recursion scheme is developed which does not require the computation of undermined coefficients, whereas most of previous recursive schemes do require the computation of undermined coefficients (see [2325]). In fact, the proposed recursion scheme is useful for solving problem (1.1) whether they are linear or nonlinear singular (nonsingular). The main advantage of the method is that it provides a direct scheme for solving the singular boundary value problem, that is, without linearization and discretization. Numerical results are presented to demonstrate the effectiveness of proposed recursive scheme. The symbolic and numerical computations have been performed using “MATHEMATICA” software.

The organization of the paper is as follows. Section 2 presents the inverse integral operator with modified Adomian decomposition method for doubly singular boundary value problems of (1.1) with Types 1 and 2 boundary conditions. In Section 3, the convergence of the method is discussed. In Section 4, we illustrate our method with numerical results along with graphical representation.

#### 2. Modified Adomian Decomposition Method

In this section, we establish two modified recursive schemes for solving doubly singular boundary value problems (1.1) with Types 1 and 2 boundary conditions.

We again write (1.1) with Type 1 boundary conditions as: which can be rewritten as where is the linear differential operator to be inverted and is an analytic nonlinear operator. It is also assumed that the solution of the problem (2.1) exits and unique. A twofold integral operator , regarded as the inverse operator of , is proposed as To establish the recursive scheme, we operate on the left hand side of (2.1) and use initial condition ; we have where .

We again operate the inverse operator on both sides of (2.1) and use (2.5), yielding For simplicity, we set and

Then (2.6) can be written as Using the boundary condition in (2.7), we get By using into (2.7), we obtain It is important to note that the right hand side of (2.9) does not contain any unknown constants.

The equation (2.9) can be rewritten in the following operator equation form as where Next, the solution and the nonlinear function are decomposed by infinite series of the form respectively, where are Adomian polynomials [19].

Substituting the series (2.12) into (2.9) gives Using the above equation (2.13), the solution components are determined by the following recursive scheme: We further modify the above algorithm (2.14) to get a more efficient and economic algorithm. To obtain modified algorithm, the zeroth component is divided into the sum of two parts, namely, , where ,  . The first part, , is kept in and the rest part, , is added to . Thus, the zeroth component will be in modified algorithm as Note that new modified recursion scheme (2.15) does not require the computation of undermined coefficients, since this algorithm does not contain any undermined coefficient. This slight change plays a major role for minimizing the size of calculations.

Now, we again consider the same problem (1.1) with Type 2 boundary conditions: which can be rewritten as where is the linear differential operator to be inverted and is an analytic nonlinear operator. It is also assumed that the solution of problem (2.16) exits and is unique. In this case, we propose the following inverse integral operator: To setup an algorithm, we operate the inverse operator on the left hand side of (2.16) and using the boundary condition , we have where .

By operating the inverse operator on both sides of (2.16) and using (2.20) we have Thus (2.21) can be written as Using boundary condition into (2.22), we obtain The equation (2.23) may be rewritten in operator equation form where Substituting the series (2.12) into (2.23) gives Upon matching both sides of (2.26), the solution components can be determined by the following modified recursive scheme: Thus, we have established two modified recursive schemes which give the complete determination of solution components , and hence the approximate series solution can be obtained for the singular boundary value problems (1.1) with Types 1 and 2 boundary conditions. For numerical purpose, the truncated -term approximate series solution is given by Note that the truncated -term approximate series solution is obtained by adding solutions components , that is, .

#### 3. Convergence of Method

In this section, we discuss the convergence analysis of modified ADM for doubly singular boundary value problem (1.1).

Many authors [3033] established the convergence of ADM for differential and integral equations. The first proof of convergence of ADM for the general functional equation was given by Cherruault [31]. Cherruault and Adomian [32] proposed a new convergence proof of Adomian decomposition method for the general nonlinear functional equation based on the properties of convergent series. Recently, Hosseini and Nasabzadeh [33] introduced a simple technique to determine the rate of convergence of ADM for initial value problem.

Now, we discuss the convergence of doubly singular two-point boundary value problems (1.1) with Types 1 and 2 boundary conditions. To do so, note that (2.9) and (2.23) may be written in the operator equation form where is an operator from a Banach space to , is a given function in , and we are looking for satisfying (3.1).

For Type 1 boundary conditions, we have For Type 2 boundary conditions, and are given as Let the solution and the nonlinear function be considered as sum of series where are Adomian polynomials [19].

Substituting the series (3.6) into the operator equation (3.1) yields where for Types 1 and 2 boundary conditions, respectively.

Upon matching both sides of (3.7), we obtain the following scheme:

Let be approximate solution, then the modified ADM for (3.1) is equivalent to the following problem: Finding the solution of (3.1) is equivalent to finding the sequence such that satisfies (3.10).

Theorem 3.1. Let be the nonlinear operator defined by (3.3) or (3.5) which satisfies the Lipschitz condition , for all with Lipschitz constant , . If , then there holds , for all and the sequence defined by (3.10) converges to .

Proof. Since we have , .
We now show that the sequence is convergent. To prove this, it is sufficient to show that is cauchy sequence in Banach space . Now using the nonlinearity of , we have Thus we obtain Now for every , and , we have As we know is geometric series with common ration , . Therefore, we have which converges to zero, that is, , as . This implies that there exits such that . But, we have , that is, which is solution of (3.1). Hence, we have convergence in norm.

#### 4. Numerical Examples and Discussion

In this section, we consider four examples to demonstrate the effectiveness of proposed recursive schemes (2.15) and (2.27). The numerical results are compared with known results and maximum absolute error is also calculated. We also plot approximate and exact solutions to show how approximate solutions converge to exact solution.

Example 4.1. Consider linear two-point boundary value problem [6] with exact solution . For every and , (4.1) is called singular.

In particular, if ,  , that is, and then problem (4.1) is called doubly singular boundary value problem. Applying modified recursive scheme (2.15) to (4.1), the solution components are obtained as For numerical purpose, the approximate series solution is Now, we define error function as and the maximum absolute errors as For different values of and , the maximum absolute error , for , and 8, is shown in Tables 1 and 2, and we observe that when increases the error decreases. On the other hand, Chawla and Katti [6] applied method based on three evaluations of to solve the same problem in which a huge amount of computation work is needed to obtain the numerical solution. The results are compared with our approximate series solution in Table 3. Furthermore, the approximate series solutions ,   and exact solution are plotted in Figures 1 and 2. It can be seen from the figures that series solution is very close to exact solution.

Table 1: Maximum error of Example 4.1, when .
Table 2: Maximum error of Example 4.1, when .
Table 3: Maximum error of Example 4.1.
Figure 1: Comparison of exact and approximate ,   solutions of Example 4.1, ,  .
Figure 2: Comparison of exact and approximate ,   solutions of Example 4.1, ,  .

Example 4.2. Consider nonlinear singular two-point boundary value problem with exact solution and for any real and , this equation is called singular.

If ,  , that is, , , then problem (4.5) is called doubly singular boundary value problem.

Now we apply recursive scheme (2.15) to obtain the solution components as To compare the results numerically, the 6-term truncated series solution is given by Similarly, for ,  , that is, and , the solution components are obtained using the scheme (2.15) as The maximum absolute error , for , and , is listed in Tables 4 and 5. As the value of increases, the maximum absolute error decreases. In addition, to verify how close is our approximate series solutions ,   to exact solution , we have plotted approximate solution and exact solution in Figures 3 and 4. It is shown from figures that the approximate series solution is very near to the exact solution .

Table 4: Maximum error for Example 4.2, when .
Table 5: Maximum error for Example 4.2, when .
Figure 3: Comparison of exact and approximate ,   solutions of Example 4.2, when ,  .
Figure 4: Comparison of exact and approximate ,   solutions of Example 4.2, when ,  .

Example 4.3. The heat conduction model in human head problem Duggan and Goodman [13] is We apply the recursive scheme (2.27) to (4.9), where , , , , and then scheme (2.27) for (4.9) becomes as

We can easily determine the solution components from above scheme as

The truncated series solution is obtained by adding , given by The comparison of numerical results of solutions with known solution is given in Table 6. It can be noted that the approximate solution is very close to known solution [13]. However, Duggan and Goodman [13] utilized the maximum principle based on lower and an upper solutions technique in which a lot of calculation is required as we need to solve sequence of linear differential equation analytically.

Table 6: Comparison of numerical results with known results of Example 4.3.

Example 4.4. The nonlinear singular boundary problem considered by Chawla et al. [34] is with exact solution , where , .

In the above problem, we have and ,  ,  ,  , by applying recursive scheme (2.27), we have following recursive scheme: We have successive solution components as and then the truncated series is obtained as In order to check whether the approximate solution converges to exact solution , we plot approximate solution , for and exact solution in Figure 5. It is clear from the figure that the approximate solution converges rapidly to exact solution. The maximum absolute error , for , and , is presented in Table 7. It can be noted that Chawla et al. [34] employed fourth-order finite difference method to obtain the numerical solution in which a lot of computation work is required. However, the maximum absolute error obtained was using .

Table 7: Numerical results and maximum error of Example 4.4.
Figure 5: Comparison of exact and approximate ,   solutions of Example 4.4.

#### 5. Conclusion

In this paper, we we have illustrated how the modified recursive schemes (2.15) and (2.27) can be used to solve a class of doubly singular two-point boundary value problems (1.1) with Types 1 and 2 boundary conditions. The accuracy of the numerical results indicates that the method is well suited for the solution of this type problem. The main advantage of this approach is that it provides a direct scheme to obtain approximate solutions, and we have also shown graphically that these approximate solutions are almost identical to the analytic solution. Another advantage of a modified recursion scheme is that this does not require the computation of undermined coefficients, whereas most of previous recursive schemes do require the computation of undermined coefficients (see [2325]). The method provides a reliable technique which requires less work compared to the traditional techniques such as finite difference method, Cubic spline method, and standard ADM. The numerical results of the examples are presented and only a few terms are required to obtain accurate solutions. By comparing the results with other existing methods, it has been proved that proposed modified ADM is a more powerful method for solving the singular problems (1.1).

#### References

1. L. E. Bobisud, “Existence of solutions for nonlinear singular boundary value problems,” Applicable Analysis, vol. 35, no. 1–4, pp. 43–57, 1990.
2. T. Aziz and M. Kumar, “A fourth-order finite-difference method based on non-uniform mesh for a class of singular two-point boundary value problems,” Journal of Computational and Applied Mathematics, vol. 136, no. 1-2, pp. 337–342, 2001.
3. M. M. Chawla and C. P. Katti, “A uniform mesh finite difference method for a class of singular two-point boundary value problems,” SIAM Journal on Numerical Analysis, vol. 22, no. 3, pp. 561–565, 1985.
4. A. Ebaid, “A new analytical and numerical treatment for singular two-point boundary value problems via the Adomian decomposition method,” Journal of Computational and Applied Mathematics, vol. 235, no. 8, pp. 1914–1924, 2011.
5. M. Kumar, “A new finite difference method for a class of singular two-point boundary value problems,” Applied Mathematics and Computation, vol. 143, no. 2-3, pp. 551–557, 2003.
6. M. M. Chawla and C. P. Katti, “Finite difference methods and their convergence for a class of singular two-point boundary value problems,” Numerische Mathematik, vol. 39, no. 3, pp. 341–350, 1982.
7. M. M. Chawla and C. P. Katti, “A finite-difference method for a class of singular two-point boundary value problems,” IMA Journal of Numerical Analysis, vol. 4, no. 4, pp. 457–466, 1984.
8. M. M. Chawla, S. McKee, and G. Shaw, “Order 2 method for a singular two-point boundary value problem,” BIT Numerical Mathematics, vol. 26, no. 3, pp. 318–326, 1986.
9. S. R. K. Iyengar and P. Jain, “Spline finite difference methods for singular two point boundary value problems,” Numerische Mathematik, vol. 50, no. 3, pp. 363–376, 1987.
10. L. Thomas, “The calculation of atomic fields,” in Mathematical Proceedings of the Cambridge Philosophical Society, vol. 23, pp. 542–548, Cambridge University Press, 1927.
11. E. Fermi, “Un metodo statistico per la determinazione di alcune priorieta dell’atome,” Rendicondi Accademia Nazionale de Lincei, vol. 32, no. 6, pp. 602–607, 1927.
12. C. Y. Chan and Y. C. Hon, “A constructive solution for a generalized Thomas-Fermi theory of ionized atoms,” Quarterly of Applied Mathematics, vol. 45, no. 3, pp. 591–599, 1987.
13. R. Duggan and A. Goodman, “Pointwise bounds for a nonlinear heat conduction model of the human head,” Bulletin of Mathematical Biology, vol. 48, no. 2, pp. 229–236, 1986.
14. G. W. Reddien, “Projection methods and singular two point boundary value problems,” Numerische Mathematik, vol. 21, pp. 193–205, 1974.
15. M. M. Chawla and P. N. Shivakumar, “On the existence of solutions of a class of singular nonlinear two-point boundary value problems,” Journal of Computational and Applied Mathematics, vol. 19, no. 3, pp. 379–388, 1987.
16. D. R. Dunninger and J. C. Kurtz, “Existence of solutions for some nonlinear singular boundary value problems,” Journal of Mathematical Analysis and Applications, vol. 115, no. 2, pp. 396–405, 1986.
17. R. K. Pandey and A. K. Verma, “A note on existence-uniqueness results for a class of doubly singular boundary value problems,” Nonlinear Analysis: Theory, Methods & Applications, vol. 71, no. 7-8, pp. 3477–3487, 2009.
18. G. Adomian, Solving Frontier Problems of Physics: The Decomposition Method, vol. 60 of Fundamental Theories of Physics, Kluwer Academic, Dordrecht, The Netherlands, 1994.
19. G. Adomian and R. Rach, “Inversion of nonlinear stochastic operators,” Journal of Mathematical Analysis and Applications, vol. 91, no. 1, pp. 39–46, 1983.
20. G. Adomian and R. Rach, “A new algorithm for matching boundary conditions in decomposition solutions,” Applied Mathematics and Computation, vol. 58, no. 1, pp. 61–68, 1993.
21. G. Adomian and R. Rach, “Modified decomposition solution of linear and nonlinear boundary-value problems,” Nonlinear Analysis: Theory, Methods & Applications, vol. 23, no. 5, pp. 615–619, 1994.
22. A.-M. Wazwaz, “Approximate solutions to boundary value problems of higher order by the modified decomposition method,” Computers & Mathematics with Applications, vol. 40, no. 6-7, pp. 679–691, 2000.
23. A. M. Wazwaz, “A reliable algorithm for obtaining positive solutions for nonlinear boundary value problems,” Computers & Mathematics with Applications, vol. 41, no. 10-11, pp. 1237–1244, 2001.
24. M. Benabidallah and Y. Cherruault, “Application of the Adomian method for solving a class of boundary problems,” Kybernetes, vol. 33, no. 1, pp. 118–132, 2004.
25. M. Inc and D. J. Evans, “The decomposition method for solving of a class of singular two-point boundary value problems,” International Journal of Computer Mathematics, vol. 80, no. 7, pp. 869–882, 2003.
26. B. Jang, “Two-point boundary value problems by the extended Adomian decomposition method,” Journal of Computational and Applied Mathematics, vol. 219, no. 1, pp. 253–262, 2008.
27. S. A. Khuri and A. Sayfy, “A novel approach for the solution of a class of singular boundary value problems arising in physiology,” Mathematical and Computer Modelling, vol. 52, no. 3-4, pp. 626–636, 2010.
28. D. Lesnic, “A computational algebraic investigation of the decomposition method for time-dependent problems,” Applied Mathematics and Computation, vol. 119, no. 2-3, pp. 197–206, 2001.
29. J.-S. Duan and R. Rach, “A new modification of the Adomian decomposition method for solving boundary value problems for higher order nonlinear differential equations,” Applied Mathematics and Computation, vol. 218, no. 8, pp. 4090–4118, 2011.
30. K. Abbaoui and Y. Cherruault, “Convergence of Adomian's method applied to differential equations,” Computers & Mathematics with Applications, vol. 28, no. 5, pp. 103–109, 1994.
31. Y. Cherruault, “Convergence of Adomian's method,” Kybernetes, vol. 18, no. 2, pp. 31–38, 1989.
32. Y. Cherruault and G. Adomian, “Decomposition methods: a new proof of convergence,” Mathematical and Computer Modelling, vol. 18, no. 12, pp. 103–106, 1993.
33. M. M. Hosseini and H. Nasabzadeh, “On the convergence of Adomian decomposition method,” Applied Mathematics and Computation, vol. 182, no. 1, pp. 536–543, 2006.
34. M. M. Chawla, R. Subramanian, and H. L. Sathi, “A fourth order method for a singular two-point boundary value problem,” BIT Numerical Mathematics, vol. 28, no. 1, pp. 88–97, 1988.