In this article, with the help of fractal quintic spline, numerical solutions of singularly perturbed boundary value problem of ordinary differential equations are obtained. To check the accuracy of the proposed method, convergence analysis is derived. The developed method has fourth order convergence. To check the correctness of the proposed method, numerical examples are provided.