The present study is related to design a novel multi-fractional multi-singular Lane–Emden model (MFMS-LEM) by keeping the ideas of the literature LEM and by extension of the work of doubly singular multi-fractional LEM. This mathematical novel MFMS-LEM is numerically treated by applying the fractional Meyer neuro-evolution intelligent solver (FMNEICS). The optimization is performed using the mutual heuristics of fractional Mayer wavelet neural networks (FMW-NN), the global search aptitude of genetic algorithms (GAs) and interior-point algorithm (IPA), i.e., FMW-NN-GAIPA. The derivation steps, details of the singular points, fractional terms, shape factors and singular points are also provided. The modeling strength of MW-NN is implemented to characterize the novel model in the sagacity of mean squared error of objective function and network optimization is performed with the integrated capability of GAIPA. The authentication, perfection and verification of FMNEICS is checked for three diverse cases of the novel model which are conventional via relative studies through the reference solutions based on accuracy, stability, robustness and convergence procedures. Furthermore, the explanations via the statistical measures validate the value of the designed stochastic solver FMW-NN-GAIPA.