“…As the Equations (36) to (41) form coupled, the solutions can be found by solving the equations separately pairwise and successively. In this study, implicit Runge-Kutta-Butcher 28 processes together with Nachtsheim-Swigert 29 iterative technique is applied to find the solutions of (36) to (41) up to Ο ξ ( ) 10 . Also, the Pade's approximant 30 is employed to get the better approximation to the amplitude and phase angles of the undulating part of the shear stress, wall temperature, and wall species concentration.…”