This paper examines the solutions of nonlinear higher-order singular Emden-Fowler type equations arising in various physical models. Generally, it becomes difficult to obtain the solution near the point of singularity. To overcome this problem, an iterative technique is introduced that depends on the variational iteration method (VIM) and the homotopy perturbation method (HPM). Such a technique generates the solution in terms of a series, which is highly practical from computing perspective. An equivalent recursive integral representation (involving Lagrange's multiplier) for the higher order nonlinear singular Emden-Fowler type (SEFT) equations with initial conditions (ICs) is established with the support of the variational iteration method (VIM). Making use of the concept of homotopy, a system of integral equations is established, which helps to deal with nonlinearity. Some numerical examples are studied through the proposed iterative technique to show the applicability and efficiency of the technique.
[1] V. Ananthaswamy and S. Punitha, Mathematical study on infinite boundary value problem for MHD flow of a micropolar nanofluid, Comput. Methods Differ. Equ., 14(2) (2026), 701–720.
[2] S. Aydinlik and A. Kiris, A high-order numerical method for solving nonlinear Lane-Emden type equations arising in astrophysics, Astrophys Space Sci., 363 (2018), 1-12.
[3] S. Chandrasekhar, An introduction to the study of Stellar structure, Dover Publications, New York, 1967.
[4] P. L. Chambre, On the solution of the Poisson-Boltzmann equation with application to the theory of thermal explosions, J. Chem. Phys., 20 (1952), 1795-1797.
[5] H. T. Davis, Introduction to nonlinear differential and integral equations, US Atomic Energy Commission, 1960.
[6] A. Dezhbord, T. Lotfi, and K. Mahdiani, A numerical approach for solving the high-order nonlinear singular Emden-Fowler type equations, Adv. Difference Equ., 2018(1) (2018), 1-17.
[7] R. C. Duggan and A. M. Goodman, Pointwise bounds for a nonlinear heat conduction model of the human head, Bull. Math. Biol., 48(2) (1986), 229-236.
[8] R. Emden, Gaskugeln Anwendungen der Mechan, Warmtheorie. Teubner, Leipzig/Berlin. (1907).
[9] R. H. Fowler, Further studies of Emden’s and similar differential equations, Q. J. Math., 2(1) (1931), 259-288.
[10] J. H. He, Homotopy perturbation technique, Comput. Methods Appl. Mech. Engrg., 178 (1999), 257-262.
[11] J. H. He, Variational iteration method-a kind of nonlinear analytical technique, Some examples, Int. J. Non-Linear Mech., 34(4) (1999), 699-708.
[12] Y. J. Huang and H. K. Liu, A new modification of the variational iteration method for van der Pol equations, Appl. Math. Model., 37(16-17) (2013), 8118-8130.
[13] M. K. Iqbal, M. Abbas, and I. Wasim, New cubic B-spline approximation for solving third order Emden-Fowler type equations, Appl. Math. Comput., 331 (2018), 319-333.
[14] J. B. Keller, Electrohydrodynamics I. The equilibrium of a charged gas in a container, J. Ration. Mech. Anal., 5(4) (1956), 715-724.
[15] J. H. Lane, On the theoretical temperature of the sun; under the hypothesis of a gaseous mass maintaining its volume by its internal heat, and depending on the laws of gases as known to terrestial experiment, Am. J. Sci., 2(148) (1870), 57-74.
[16] S. Liao, An optimal homotopy-analysis approach for strongly nonlinear differential equations, Commun. Nonlinear Sci. Numer. Simulat., 15(8) (2010), 2003-2016.
[17] B. Lin, A new numerical scheme for third-order singularly Emden-Fowler equations using quintic B-spline function, Int. J. Comput. Math., 98(12) (2021), 2406-2422.
[18] V. Marinca and N. Herisanu, Application of optimal homotopy asymptotic method for solving nonlinear equations arising in heat transfer, Int. Commun. Heat Mass Transf., 35(6) (2008), 710-715.
[19] M. Merafina, G. S. Bisnovatyi-Kogan, and S. O. Tarasov, A brief analysis of self-graviating polytropic models with a non-zero cosmological constant, Astron. Astrophys., 541 (2012), A84.
[20] B. Muatjetjeja and C. M. Khalique, First integrals for a generalized coupled Lane-Emden system, Nonlinear Anal. Real World Appl., 12(2) (2011), 1202-1212.
[21] R. K. Pandey and A. K. Verma, Existence-uniqueness results for a class of singular boundary value problems arising in physiology, Nonlinear Anal. Real World Appl., 9 (2008), 40-52.
[22] K. Parand, M. Shahini, and M. Dehghan, Rational Legendre pseudospectral approach for solving nonlinear differential equations of Lane-Emden type, J. Comput. Phys., 228 (2009), 8830-8840.
[23] Z. Perveen, Z. Fatima, A. H. Majeed, and A. Refaie ALi, Exact and iterative solutions for DEs, including FokkerPlanck and Newell-Whitehead-Segel equations using Shehu Transform and HPM, Comput. Methods Differ. Equ., (1-22) (2024).
[24] O. W. Richardson, The Emission of Electricity from Hot Bodies, Longmans Green and Company, 4 (1921).
[25] M. A. Rufai and H. Ramos, Numerical integration of third-order singular boundary value problems of EmdenFowler type using hybrid block techniques, Commun. Nonlinear Sci. Numer. Simul., 105 (2022), 106069.
[26] J. Shahni and R. Singh, An efficient numerical technique for Lane-Emden-Fowler boundary value problems, Bernstein collocation method, Eur. Phys. J. Plus, 135(6) (2020), 1-21.
[27] M. Singh and A. K. Verma, An effective computational technique for a class of Lane-Emden equations, J. Math. Chem., 54(1) (2016), 231-251.
[28] M. Singh, A. K. Verma, and R. P. Agarwal, On iterative method for class of 2 point & 3 point nonlinear SBVPs, J. Appl. Anal. Comput., 9 (2019), 1242-1260.
[29] R. Singh, H. Garg, and V. Guleria, Haar wavelet collocation method for Lane-Emden equations with Dirichlet, Neumann and Neumann-Robin boundary conditions, J. Comput. Appl. Math., 346 (2018), 150-161.
[30] Swati, M. Singh, and K. Singh, An efficient technique based on higher order Haar wavelet method for Lane-Emden equations, Math Comput. Simul., 206 (2023), 21-39.
[31] Swati, K. Singh, A. K. Verma, and M. Singh, Higher order Emden-Fowler type equations via uniform Haar Wavelet resolution technique, J. Comput. Appl. Math., 376 (2020), 112836.
[32] A. K. Verma and S. Kayenat, Applications of modified Mickens-type NSFD schemes to Lane-Emden equations, Comp. Appl. Math, 39 (2020), 1-25.
[33] A. K. Verma, M. Singh, and R. P. Agarwal, Regions of existence for a class of nonlinear diffusion type problems, Appl. Anal. Discret. Math., 14(1) (2020), 106-121.
[34] A. K. Verma, B. Pandit, and R. P. Agarwal, On approximate stationary radial solutions for a class of boundary value problems arising in epitaxial growth theory, J. Appl. Comput. Mech., 6(4) (2020), 713-734.
[35] A. K. Verma, B. Pandit, and C. Escudero, Numerical solutions for a class of singular boundary value problems arising in the theory of epitaxial growth, Eng. Comput., 37(7) (2020), 2539-2560.
[36] A. Verma and M. Kumar, Numerical solution of third-order Emden-Fowler type equations using artificial neural network technique, Eur. Phys. J. Plus, 135(9) (2020), 1-14.
[37] A. M. Wazwaz, R. Rach, L. Bougoffa, and J. S. Duan, Solving the Lane-Emden-Fowler type equations of higher orders by the Adomian decomposition method, Comput. Model. Eng. Sci. (CMES), 100(6) (2014), 507-529.
[38] A. M. Wazwaz, The variational iteration method for solving new fourth-order Emden-Fowler type equations, Chem. Eng. Commun., 202(11) (2015), 1425-1437.
[39] A. M. Wazwaz, Solving two Emden-Fowler type equations of third order by the variational iteration method, Appl. Math. Inf. Sci., 9(5) (2015), 2429-2436.
[40] A. M. Wazwaz, Adomian decomposition method for a reliable treatment of the Emden-Fowler equation, Appl. Math. Comput., 161 (2005), 543-560.
[41] J. S. W. Wong, On the generalized Emden-Fowler Equation, SIAM Rev., 17 (1975), 339-360.
[42] H. Zou, A priori estimates for a semilinear elliptic system without variational structure and their applications, Math. Ann., 323(4) (2002), 713-735.
-, J. and Singh, M. (2026). A computational iterative technique for a kind of nonlinear higher-order singular Emden-Fowler type equations. Computational Methods for Differential Equations, 14(2), 852-868. doi: 10.22034/cmde.2025.65131.2979
MLA
-, J. , and Singh, M. . "A computational iterative technique for a kind of nonlinear higher-order singular Emden-Fowler type equations", Computational Methods for Differential Equations, 14, 2, 2026, 852-868. doi: 10.22034/cmde.2025.65131.2979
HARVARD
-, J., Singh, M. (2026). 'A computational iterative technique for a kind of nonlinear higher-order singular Emden-Fowler type equations', Computational Methods for Differential Equations, 14(2), pp. 852-868. doi: 10.22034/cmde.2025.65131.2979
CHICAGO
J. - and M. Singh, "A computational iterative technique for a kind of nonlinear higher-order singular Emden-Fowler type equations," Computational Methods for Differential Equations, 14 2 (2026): 852-868, doi: 10.22034/cmde.2025.65131.2979
VANCOUVER
-, J., Singh, M. A computational iterative technique for a kind of nonlinear higher-order singular Emden-Fowler type equations. Computational Methods for Differential Equations, 2026; 14(2): 852-868. doi: 10.22034/cmde.2025.65131.2979