The main purpose of this work is to investigate the influence of multi slip implications on steady-state MHD fluid flow in the occurrence of Soret and Dufour through a suction/ injection over a non-Isothermal stretching surface. The governing equations are converted into non-linear ODEs using similarity variables. The Keller box technique is used to achieve the numerical solution to ODEs. The influence of the numerous pertinent flow factors on the fluid velocity, temperature, and species concentration fields, For engineering interest we also calculated the physical quantities the skin friction coefficient, Nusselt and Sherwood number are derived. And discussed in detail with slips and without slips (Velocity, Thermal, and Concentration) with the aid of graphs and tables. It is observed that The velocity profile increases as an increase in Gr, Gc and reverse pattern for an increase in magnetic parameter (M), permeability parameter (K). It is inferred from the graphs that the temperature of the plate increases with increasing the values of M, Ec and Df and reverse trend for increasing suction and Soret number. The rate of mass transfer upsurge with an increase in Dufour number, and Suction parameter and coincide with previously published results with fast convergence. To ensure the validity of our findings, we compared them to previously published work and found significant agreement.