A new high-order numerical method for solving singular two-point boundary value problems

Recently, Goh et al. (2012) [28] proposed a numerical technique based on quartic B-spline collocation for solving a class of singular boundary value problems (SBVP) with Neumann and Dirichlet boundary conditions (BC). This method is only fourth-order accurate. In this paper, we propose an optimal numerical technique for solving a more general class of nonlinear SBVP subject to Neumann and Robin BC. The method is based on high order perturbation of the problem under consideration. The convergence of the proposed method is analyzed. To demonstrate the applicability and efficiency of the method, we consider four numerical examples, three of which arise in various physical models in applied science and engineering. A comparison with other available numerical solutions has been carried out to justify the advantage of the proposed technique. Numerical result reveals that the proposed method is sixth order convergent, which in turn is two orders of magnitude larger than in Goh et al. (2012) [28].