In this paper, a neuro-evolution based numerical computing approach is presented for the solution of nonlinear third order multi-singular Emden–Fowler system of differential equations (MS-EF-SDEs) by manipulating the proficiency of continuous mapping through exploitation of feed-forward artificial neural networks (ANN). The weights or decision variables of these networks are optimized with genetic algorithms (GAs) and sequential quadratic programming (SQP), i.e., ANN-GA-SQP. An error based figure of merit is introduced using the differential model of MS EF-SDE along with corresponding boundary conditions. The objective/cost function is optimized by integrating capability of global and local search with GA and SQP, respectively. The competency of the designed ANN-GA-SQP approach in terms of significance, efficiency and consistency is perceived by solving MS-EF-SDEs. Moreover, statistical based investigations are implemented to validate the correctness of ANN-GA-SQP.