“…The overview, importance and applications of classical numerical approaches for TFE can be seen (Kirzhnits 1957 ; Bush and Caldwell 1931 ). The research community has shown great interest in the accurate and reliable calculation of the solution for the TFE including sinc-collocation method (Parand et al 2013a ), Laguerre pseudospectral approximation (Liu and Zhu 2015 ), Chebyshev pseudospectral method (Kılıçman et al 2014 ), rational Chebyshev pseudospectral approach (Parand and Shahini 2009 ), hermite collocation method (Bayatbabolghani and Parand 2014 ), rational approximation (Fernández 2011 ), Homotopy Analysis Method (HAM) (Yao 2008 ) improved HAM (Zhao et al 2012 ), rational Bessel functions collocation method (Parand et al 2016a ), rational Euler functions based methods (Parand et al 2016b ), methods based on Jacobi rational functions with Gauss quadrature formula (Bhrawy and El-Soubhy 2015 ) and Padé–Hankel method (Amore et al 2014 ), Beside these there are many other studies for solving Thomas–Fermi models, see (Fernández 2008 ; Liao 2003a ; Filobello-Nino et al 2015 ; Dahmani and Anber 2015 ; Feng et al 2015 ) and the references therein.…”