Journal of Applied MathematicsVolume 2012 (2012), Article ID 316534, 9 pagesdoi:10.1155/2012/316534

Research Article

A Numerical Method for Lane-Emden Equations Using Hybrid Functions and the Collocation Method

Changqing Yang and Jianhua Hou

Department of Science, Huaihai Institute of Technology, Jiangsu Lianyungang, 222005, China

Received 26 September 2011; Accepted 18 November 2011

Academic Editor: Song Cen

Copyright © 2012 Changqing Yang and Jianhua Hou. 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.

A numerical method to solve Lane-Emden equations as singular initial value problems is presented in this work. This method is based on the replacement of unknown functions through a truncated series of hybrid of block-pulse functions and Chebyshev polynomials. The collocation method transforms the differential equation into a system of algebraic equations. It also has application in a wide area of differential equations. Corresponding numerical examples are presented to demonstrate the accuracy of the proposed method.

1. Introduction

In recent years, studies on singular initial value problems in second-order ordinary differential equations have attracted much attention among mathematicians and physicists. One of the equations describing this type of problem is the Lane-Emden equation which can be written in the following form: where and are constants, and is a real-valued continuous function. This equation has been used in the fields of mathematical physics and astrophysics to model several problems such as the stellar structure theory, the thermal behavior of a spherical gas cloud, isothermal gas spheres, and the thermionic current theory [1–3].

On the other hand, another type of singular initial value problems of the Lane-Emden equation can also be given in the following form: where and are constants, is a continuous real-valued function, and .

Many methods have been used to solve Lane-Emden equations. However, most algorithms in current use to deal with Lane-Emden problems are based on either series solutions or perturbation techniques. For instance, approximate solutions to (1.1) were presented by Shawagfeh [4] and Wazwaz [5, 6] using the Adomian decomposition method which provides a convergent series solution. Ramos [7] presented a series approach to the Lane-Emden equation and made comparisons with He’s homotopy perturbation method. Dehghan and Shakeri [8] first applied exponential transformation to the Lane-Emden equation in order to address the difficulty of a singular point at and solve the resulting nonsingular problem using the variational iteration method. Yousefi [9] applied an integral operator, converted Lane-Emden equations to integral equations, and solved these integral equations using Legendre wavelets. Recently, Parand et al. [10] proposed an approximation algorithm for the solution of the nonlinear Lane-Emden type equation using Hermite functions collocation method. Moreover, Adibi and Rismani [11] introduced a modified Legendre-spectral method for solving (1.1). Wazwaz [12] applied the variational iteration method to solve the singular initial problems or boundary value problems of Emden-Fowler type of equations.

In this work, we present a reliable algorithm based on hybrid functions and the collocation method in order to obtain approximate solutions for (1.3).

The main advantage of our proposed method is its direct application to all types of differential equations whether linear or nonlinear, homogeneous or inhomogeneous, and with constant coefficients or with variable coefficients. Furthermore, the method can greatly reduce the size of computational work while still maintaining the numerical solution’s high accuracy because it does not employ numerical integration [9, 11, 13].

2. Properties of Hybrid Functions

2.1. Hybrid Functions of Block-Pulse and Chebyshev Polynomials

Hybrid functions , , , are defined on the interval as Here, are the well-known Chebyshev polynomials of order which satisfy the following recursive formula:

2.2. Function Approximation

A function defined over the interval 0 to may be expanded as where in which denotes the inner product.

If the infinite series in (2.3) is truncated, then (2.3) can be written as where In (2.5) and (2.6), , , , are the coefficients expansions of the function and , , , are defined in (2.1).

2.3. Operational Matrix of the Integration

The integration of the vector defined in (2.7) can be obtained as where is matrix, that is called the operational matrix for integration and is given in [14].

3. Numerical Solution Using the Hybrid Functions and the Collocation Method

We then consider the Lane-Emden equations given in (1.3). First, we assume that the unknown function is given by Then using Condition (1.4), we obtain Substituting (3.1) and (3.2) into (1.3), we obtain Rewriting (3.3) in the form We now collocate (3.4) at points as Suitable collocation points are zeros of Chebyshev polynomials System (3.5) consists of nonlinear equations which can be solved using conventional iterative methods such as Newton’s method or the simplex method.

4. Accuracy of Solution

We can easily verify the accuracy of the method. Given that the truncated hybrid functions series in (2.5) are an approximate solution of (1.3), it must have approximately satisfied these equations. Thus, for each If max ( is any positive integer) is prescribed, then the truncation limit is increased until the difference at each of the points becomes smaller than the prescribed .

Proposition 4.1. Let (Sobolev space) and be the best approximation polynomial of in -norm. Thus, the truncation error is where is a positive constant, which is dependent on the selected norm and independent of and (proof [15]).

Proposition 4.2. Let , then By using of Proposition 4.1, it is obvious [16].

5. Numerical Examples

In the section, we demonstrate the accuracy of the proposed method in solving the Lane-Emden equation. In all the examples, we use the relations shown in Section 3 to convert a differential equation to a system of nonlinear equations. The numerical results are presented in tabular form to illustrate the efficiency of the method. All computations were carried out using Matlab7.0.

Example 5.1. We consider the linear singular initial value problem [6] We applied the method presented in this paper and solved (5.1) with , . For this equation, we find Using (2.6), we obtain which is the exact solution.

Example 5.2. We consider the linear singular initial value problem [6, 13] The exact solution for this problem is . We assume that Using the boundary condition in (5.7)), we derive Now, by inserting (5.7)–(5.9) into (5.5), we obtain The collocation points are Using these points, we obtain By solving (5.12), we obtain the coefficients ’s. Using (5.9), we then derive the approximate solution. Table 1 shows the absolute error of the exact and approximate solution for , and , .

Table 1: Comparison of the absolute errors for Example 5.2.

Method in [10]
Present method
Present method

(, )
(, )

0
0
0
0

0.1
1.78e−08
5.73e−006
2.52e−009

0.2
2.09e−08
7.22e−006
3.05e−010

0.5
2.62e−08
1.17e−005
2.61e−009

0.7
3.27e−08
1.22e−006
1.28e−009

0.8
3.79e−08
1.88e−006
6.19e−009

0.9
5.48e−08
1.31e−006
2.02e−010

1.0
2.51e−09
8.49e−006
1.30e−010

Example 5.3. We consider the nonlinear singular initial value problem [6, 10, 13] with the exact solution .

We solved (5.13) using the method presented in this paper. Table 2 shows the absolute errors of exact and approximate solutions in some points of the interval [] obtained by presented method for , . The comparison of absolute errors using present method and the method in [10] with is shown in Figure 1. It is easily found that the present approximations are more efficient.

Table 2: Absolute error of Example 5.3′s exact and approximated solution.

Approximated solution
Exact solution
Absolute error

0
0
0
0

0.1
−1.990065975846963e−002
−1.990066170633618e−002
1.947866552082234e−009

0.2
−7.844142227675365e−002
−7.844142630656266e−002
4.029809014260977e−009

0.3
−1.723553877014332e−001
−1.723553924821048e−001
4.780671636339662e−009

0.4
2.968400304512400e−001
−2.968400102365468e−001
2.021469319712921e−008

0.5
−4.462871382464894e−001
−4.462871026284195e−001
3.561806982688154e−008

0.6
−6.149693761315079e−001
−6.149693994959211e−001
2.336441318195881e−008

0.7
−7.975522016254086e−001
−7.975522399147356e−001
3.828932693128451e−008

0.8
−9.893925027250091e−001
−9.893924836722143e−001
1.905279478631172e−008

0.9
−1.186653691831187e+000
−1.186653690555469e+000
1.275718419790906e−009

1.0
−1.386294358347456e+000
−1.386294361119891e+000
2.772434948283831e−009

Figure 1: Comparison of absolute errors using the present method and method in [10] for Example 5.3.

6. Conclusion

We have solved the Lane-Emden equation by using hybrid functions and the collocation method. The properties of the hybrid functions are used to reduce the equation to the solution of nonlinear algebraic equations. Illustrative examples are given to demonstrate the validity and applicability of the proposed method. We believe that the proposed method is applicable to higher-order differential equations and to other sets of orthogonal polynomials.

Acknowledgment

C. Yang is grateful to the National Natural Science Foundation of China which provided support through Grant No. 40806011.

References

S. Chandrasekhr, Introduction to Study of Stellar Structure, Dover, New York, NY, USA, 1967.

H. T. Davis, Introduction to Nonlinear Differential and Integral Equations, Dover, New York, NY, USA, 1962.

O. U. Richardson, The Emission of Electricity from Hot Bodies, Longman, Green and Co., London, UK, 1921.

N. T. Shawagfeh, “Nonperturbative approximate solution for Lane-Emden equation,” Journal of Mathematical Physics, vol. 34, no. 9, pp. 4364–4369, 1993.

A.-M. Wazwaz, “A new algorithm for solving differential equations of Lane-Emden type,” Applied Mathematics and Computation, vol. 118, no. 2-3, pp. 287–310, 2001.

A.-M. Wazwaz, “A new method for solving singular initial value problems in the second-order ordinary differential equations,” Applied Mathematics and Computation, vol. 128, no. 1, pp. 45–57, 2002.

J. I. Ramos, “Series approach to the Lane-Emden equation and comparison with the homotopy perturbation method,” Chaos, Solitons and Fractals, vol. 38, no. 2, pp. 400–408, 2008.

M. Dehghan and F. Shakeri, “Approximate solution of differential equation arising in astrophysics using the variation method,” New Astronomy, vol. 13, pp. 53–59, 2008.

S. A. Yousefi, “Legendre wavelets method for solving differential equations of Lane-Emden type,” Applied Mathematics and Computation, vol. 181, no. 2, pp. 1417–1422, 2006.

K. Parand, M. Dehghan, A. R. Rezaei, and S. M. Ghaderi, “An approximation algorithm for the solution of the nonlinear Lane-Emden type equations arising in astrophysics using Hermite functions collocation method,” Computer Physics Communications, vol. 181, no. 6, pp. 1096–1108, 2010.

H. Adibi and A. M. Rismani, “On using a modified Legendre-spectral method for solving singular IVPs of Lane-Emden type,” Computers & Mathematics with Applications, vol. 60, no. 7, pp. 2126–2130, 2010.

A.-M. Wazwaz, “A reliable treatment of singular Emden-Fowler initial value problems and boundary value problems,” Applied Mathematics and Computation, vol. 217, no. 24, pp. 10387–10395, 2011.

S. Karimi Vanani and A. Aminataei, “On the numerical solution of differential equations of Lane-Emden type,” Computers & Mathematics with Applications, vol. 59, no. 8, pp. 2815–2820, 2010.

M. Razzaghi and H.-R. Marzban, “Direct method for variational problems via hybrid of block-pulse and Chebyshev functions,” Mathematical Problems in Engineering, vol. 6, no. 1, pp. 85–97, 2000.

C. Canuto, M. Y. Hussaini, A. Quarteroni, and T. A. Zhang, Spectral Methods on Fluid Dynamics, Springer, New York, NY, USA, 1988.

K. Maleknejad and M. Tavassoli Kajani, “Solving linear integro-differential equation system by Galerkin methods with hydrid functions,” Applied Mathematics and Computation, vol. 159, no. 3, pp. 603–612, 2004.