The current work explores the intelligent computational strength of neural networks based on the Levenberg–Marquardt backpropagation (LMBP-NNs) neural networks technique for simulation of Maxwell nanofluid flow past a linear stretchable surface model. The fluid flow is incorporated Rosseland’s thermal radiation, and Darcy’s Forchheimer law. The Maxwell nanofluid model gives more relaxing time to momentum boundary layer. For the nanofluid phenomena that concentrate on thermophoresis and Brownian motion, Buongiorno’s model is used. The procedure transforms partial differential equations arising in nanofluidics systems with an appropriate degree of similarity into nonlinear differential equation systems. For the nonlinear nanofluid problem with accuracy having order 4–5, the (FDM) finite difference method (Lobatto IIIA) is implemented via various selections of collocation points. The strong aspect of Lobatto IIIA is its ability to handle very nonlinear couple differential equations in an easy manner. The precise results of (FDM) are used to build the reference datasets for LMBP-NNs technique for the various factors of fluid problem. The design scheme for various factors of fluid problem carries out a series of operations based on training, testing, and authentication on reference dataset. The accuracy of LMBP-NNs is checked through statistical based neural network tools such that mean square error, regression plot, curve fitting graphs, and error histogram. Furthermore, the investigation of flow model parameters for momentum, energy, and concentration profiles is described via visual representation.