“…The coupled six harmonized high order ordinary differential equations (13)- (18), subjected to the boundary conditions, equation (19), are explained numerically using the function bvp4c from matlab software for different values of physical parameters, namely slip parameters (A, B, D, H), Prandtl number Pr, Dufour number Df, Soret number Sr, Brownian motion parameter Nb, thermophoresis parameter Nt, Schmidt number Sc, stretching ratio parameter λ, material parameter β and micro rotation parameter n. Numerical results are found using Matlab BVP solver bvp4c from matlab which is a finite difference code that applies the three-stage Lobatto IIIa formula. To apply bvp4c from matlab: equations (13)-(18) are converted into a system of first-order equations.…”