Abstract
This paper introduces an efficient approach for solving Lane–Emden–Fowler problems. Our method utilizes two Nyström schemes to perform the integration. To overcome the singularity at the left end of the interval, we combine an optimized scheme of Nyström type with a set of Nyström formulas that are used at the fist subinterval. The optimized technique is obtained after imposing the vanishing of some of the local truncation errors, which results in a set of symmetric hybrid points. By solving an algebraic system of equations, our proposed approach generates simultaneous approximations at all grid points, resulting in a highly effective technique that outperforms several existing numerical methods in the literature. To assess the efficiency and accuracy of our approach, we perform some numerical tests on diverse real-world problems, including singular boundary value problems (SBVPs) from chemical kinetics.