3. Numerical scheme

Numerical solutions of ordinary differential Eqs. (8) and (9) subject to boundary conditions (10) are obtained using a shooting method. First we have converted the boundary value problem (BVP) into initial value problem (IVP) and assumed a suitable finite value for the far field boundary condition, i.e.η! ∞, say η∞. To solve the IVP, the values for f00(0) and θ<sup>0</sup> (0) are needed but no such values are given prior to the computation. The initial guess values of f00(0) and θ<sup>0</sup> (0) are chosen and fourth order Runge-Kutta method is applied to obtain a solution. We compared the calculated values of f 0 (η) and θ(η) at the far field boundary condition η∞(=20) with the given boundary conditions (10) and the values of f00(0) and θ<sup>0</sup> (0), are adjusted using Secant method for better approximation. The step-size is taken as Δη = 0.01 and accuracy to the fifth decimal place as the criterion of convergence. It is important to note that the dual solutions are obtained by setting two different initial guesses for the values of f00(0).
