The present study aims to design stochastic intelligent computational heuristics for the numerical treatment of a nonlinear SITR system representing the dynamics of novel coronavirus disease 2019 (COVID-19). The mathematical SITR system using fractal parameters for COVID-19 dynamics is divided into four classes; that is, susceptible (S), infected (I), treatment (T), and recovered (R). The comprehensive details of each class along with the explanation of every parameter are provided, and the dynamics of novel COVID-19 are represented by calculating the solution of the mathematical SITR system using feed-forward artificial neural networks (FF-ANNs) trained with global search genetic algorithms (GAs) and speedy fine tuning by sequential quadratic programming (SQP)—that is, an FF-ANN-GASQP scheme. In the proposed FF-ANN-GASQP method, the objective function is formulated in the mean squared error sense using the approximate differential mapping of FF-ANNs for the SITR model, and learning of the networks is proficiently conducted with the integrated capabilities of GA and SQP. The correctness, stability, and potential of the proposed FF-ANN-GASQP scheme for the four different cases are established through comparative assessment study from the results of numerical computing with Adams solver for single as well as multiple autonomous trials. The results of statistical evaluations further authenticate the convergence and prospective accuracy of the FF-ANN-GASQP method.

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 FMNEICS 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 mean-squared 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.