Evaluation of Scipy.ode Integrators in Solving the Lane-Emden Equation for Polytropes as a Boundary Value Problem with a Fitting Method
Keywords:Lane-Emden equation, Two point BVP with fitting method, Scipy ode solvers
The use of Scipy integrators like dopri5 and others in accurately solving the Lane-Emden equation of a polytrope as a two-point BVP with fitting is investigated by comparing the Emden radius with the extended precision reference value obtained by Boyd's Chebyshev spectral method. It is found that both dopri5 and dop853 integrators provide acceptable accuracy upto 14 decimal digits.
 J.P. Boyd, "Chebyshev Spectral Methods and the Lane-Emden Problem", Numer. Math. Theor. Meth. Appl., Vol. 4, No. 2, pp. 142-157, 2011.