“…In recent years, Many powerful methods have been presented for solving of Lane-Emden type equations. For instance, the homotopy perturbation method [4,6], the Legendre wavelets [7], the variational iteration method [8,9], the B-spline method [10], the Adomian decomposition method [11], the Bessel collocation method [12], the Pade series method [13], the rational Legendre pseudospectral method [14], the nonperturbative approximate method [15], the Hermite functions collocation method [16], and the variational approach method [17]. Continuous or piecewise polynomials are incredibly useful mathematical tools as they are precisely defined, calculated rapidly on a modern computer system and can represent a great variety of functions.…”