Abstract
This work aims to find the numerical solution to a class of nonlinear singular boundary value problems (SBVPs). The considered problem has a singularity at x = 0. We introduce a computational technique comprising an optimal homotopy analysis (OHA) approach and exponential compact finite difference method (ECFDM) to solve this SBVPs. In this technique, the domain of the problem is divided into two subintervals as
(the point
is chosen sufficiently close to the singularity). In interval
, we employ the OHA scheme to overcome the singularity. In interval
, an ECFDM is designed to solve the resultant boundary value problem (BVP). Convergence analysis of the ECFDM is discussed. Furthermore, numerical experiments are performed to confirm the theoretical claims. The proposed ECFDM is shown to be fourth-order convergence.
Acknowledgments
The authors are very grateful to the anonymous referee for his insightful comments and suggestions.
Disclosure statement
No potential conflict of interest was reported by the author(s).