Numerical study of singular fractional Lane–Emden type equations arising in astrophysics

The well-known Lane–Emden equation plays an important role in describing some phenomena in mathematical physics and astrophysics. Recently, a new type of this equation with fractional order derivative in the Caputo sense has been introduced. In this paper, two computational schemes based on collocation method with operational matrices of orthonormal Bernstein polynomials are presented to obtain numerical approximate solutions of singular Lane–Emden equations of fractional order. Four illustrative examples are implemented in order to verify the efficiency and demonstrate solution accuracy.