Numerical solution for the system of Lane-Emden type equations using cubic B-spline method arising in engineering

In this study, we develop a collocation method based on cubic B-spline functions for effectively solving the system of Lane-Emden type equations arising in physics, star structure, and astrophysics. To overcome the singularity behavior of the considered system at tau=0, we apply the L'Hopital rule. Furthermore, we have carried out a convergence analysis of the proposed method and have demonstrated that it has a second-order convergence. To demonstrate the effectiveness, accuracy, simplicity, and practicality of the method, five test problems are solved numerically and the maximum absolute errors of the proposed method are compared with those of some existing methods.