A mathematical model of 2D-double diffusive layer flow model of boundary in MHD Maxwell fluid created by a sloping slope surface is constructed in this paper. The numerical findings of non-Newtonian fluid are important to the chemical processing industry, mining industry, plastics processing industry, as well as lubrication and biomedical flows. The diversity of regulatory parameters like buoyancy rate, magnetic field, mixed convection, absorption, Brownian motion, thermophoretic diffusion, Deborah number, Lewis number, Prandtl number, Soret number, as well as Dufour number contributes significant impact on the current model. The steps of research methodology are as followed: a) conversion from a separate matrix (PDE) to standard divisive calculations (ODEs), b) Final ODEs are solved in bvp4c program, which developed in MATLAB software, c) The stability analysis part also being developed in bvp4c program, to select the most effective solution in the real liquid state. Lastly, the numerical findings are built on a system of tables and diagrams. As a result, the profiles of velocity, temperature, and concentration are depicted due to the regulatory parameters, as mentioned above. In addition, the characteristics of the local Nusselt, coefficient of skin-friction as well as Sherwood numbers on the Maxwell fluid are described in detail.