In this article, we propose a new numerical method based on fractal cubic spline for the numerical approximate solutions of linear and nonlinear singular two-point boundary-value problems. Error analysis of the proposed method is examined and it is shown that proposed method has second-order convergence. Numerical examples are provided to validate the theoretical results. © 2020, Springer Nature India Private Limited.