In this study, a high-order compact finite difference method is used to solve Lane–Emden equations with various boundary conditions. The norm use first-order scheme approximate Neumann and Robin conditions, but that compromises the accuracy of entire scheme. As result, new higher-order schemes for approximating conditions are developed in work. We test applicability performance using different ...