An efficient numerical approach for solving a general class of nonlinear singular boundary value problems

This paper is concerned with the development of a collocation method based on the Bessel polynomials for numerical solution of a general class of nonlinear singular boundary value problems (SBVPs). Due to the existence of singularity at the point x=0,\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x=0,$$\end{document} we first modify the problem at the singular point. The proposed method is then developed for solving the resulting regular boundary value problem. To demonstrate the effectiveness and accuracy of the method, we apply it on several numerical examples. The numerical results obtained confirm that the present method has an advantage in terms of numerical accuracy over the uniform mesh cubic B-spline collocation (UCS) method (Roul and Goura in Appl Math Comput 341:428–450, 2019), non-standard finite difference (NSFD) method (Verma and Kayenat in J Math Chem 56:1667–1706, 2018), three-point finite difference methods (FDMs) (Pandey and Singh in Int J Comput Math 80:1323–1331, 2003; Pandey and Singh in J Comput Appl Math 205:469–478, 2007) and the cubic B-spline collocation (CBSC) method (Caglar et al. in Chaos Solitons Fractals 39:1232–1237, 2009)