In the present research work, we designed a hybrid stochastic numerical solver to investigate nonlinear singular two-point boundary value problems with Neumann and Robin boundary conditions arising in various physical models. In this method, we hybridized Harris Hawks Optimizer with Interior Point Algorithm named HHO-IPA. We construct artificial neural networks (ANNs) model for the problem, and this model is tuned with the proposed scheme. This scheme overcomes the singular behavior of problems. The accuracy and applicability of the method are illustrated by finding absolute errors in the solution. The outcomes are compared with the results present in the literature to demonstrate the effectiveness and robustness of the scheme by considering four different nonlinear singular boundary value problems. Further, the convergence of the scheme is proved by performing computational complexity analysis. Moreover, the graphical overview of statistical analysis is added to our investigation to elaborate further on the scheme's stability, accuracy, and consistency.