- About this Journal
- Abstracting and Indexing
- Aims and Scope
- Annual Issues
- Article Processing Charges
- Articles in Press
- Author Guidelines
- Bibliographic Information
- Citations to this Journal
- Contact Information
- Editorial Board
- Editorial Workflow
- Free eTOC Alerts
- Publication Ethics
- Reviewers Acknowledgment
- Submit a Manuscript
- Subscription Information
- Table of Contents
Journal of Applied Mathematics
Volume 2012 (2012), Article ID 316534, 9 pages
A Numerical Method for Lane-Emden Equations Using Hybrid Functions and the Collocation Method
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.
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.
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  and Wazwaz [5, 6] using the Adomian decomposition method which provides a convergent series solution. Ramos  presented a series approach to the Lane-Emden equation and made comparisons with He’s homotopy perturbation method. Dehghan and Shakeri  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  applied an integral operator, converted Lane-Emden equations to integral equations, and solved these integral equations using Legendre wavelets. Recently, Parand et al.  proposed an approximation algorithm for the solution of the nonlinear Lane-Emden type equation using Hermite functions collocation method. Moreover, Adibi and Rismani  introduced a modified Legendre-spectral method for solving (1.1). Wazwaz  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.
2.3. Operational Matrix of the Integration
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 ).
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  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 , .
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  with is shown in Figure 1. It is easily found that the present approximations are more efficient.
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.
C. Yang is grateful to the National Natural Science Foundation of China which provided support through Grant No. 40806011.
- 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.