In the present study, a novel fractional Meyer neuro-evolution-based intelligent computing solver (FMNEICS) is presented for numerical treatment of doubly singular multi-fractional Lane-Emden system (DSMF-LES) using combined heuristics of Meyer wavelet neural networks (MWNN) optimized with global search efficacy of genetic algorithms (GAs) and sequential quadratic programming (SQP), i.e., MWNN-GASQP. The design of novel FMNE-ICS for DSMF-LES is presented after derivation from standard Lane-Emden equation, and the singular points and shape factors along with fractional-order terms are analyzed. The MWNN modeling strength is used to represent the system model DSMF-LES in the meansquared error-based merit function and optimization of the networks is carried out with integrated optimization ability of GASQP. The verification, validation, and perfection of the FMNEICS for three different cases of DSMF-LES are established through comparative studies from reference solutions on convergence, robustness, accuracy, and stability measures. Moreover, the observations through the statistical analysis further authenticate the worth of proposed fractional MWNN-GASQP-based stochastic solver. Keywords Multi-fractional Lane-Emden model • Multi-singular systems • Artificial neural networks • Meyer wavelet neural networks • Sequential quadratic programming • Genetic algorithms Mathematics Subject Classification 34-XX • 34A08 • 34A34 • 82B31