This paper seeks to develop an efficient method for solving natural convection of a non-Newtonian nanofluid flow between two vertical flat plates. Sodium alginate is considered as the non-Newtonian fluid, and then two distinct types of nanoparticles, namely silver and copper, are added to it. To do so, the governing boundary layer and temperature equations are reduced to a set of ordinary differential equations. This approach is based on a global collocation method using Sinc basis functions, and the resulting set of ordinary differential equations are replaced by a system of algebraic equations. It is well known that the Sinc procedure converges to the solution at an exponential rate. Numerical results are included to demonstrate the validity and applicability of the method, and a comparison is made with the existing results. Also, the effect of various parameters such as Prandtl number (Pr), dimensionless non-Newtonian viscosity number () and nanoparticle volume fraction () on non-dimensional velocity and temperature profiles are discussed. It was concluded from this study that velocity and temperature increased with increasing Pr. Moreover, our results indicate that both the velocity and temperature decrease as increases. Finally, the results demonstrated that, when increases, the velocity increases but the temperature values decrease. Keywords Natural convection • Sinc-collocation • Non-Newtonian fluid • Nanoparticles List of symbols 2b Distance between the plates (m) C p Specific heat at constant pressure (J kg K −1) Ec Eckert number k Thermal conductivity (W m −1 K −1) Pr Prandtl number T Temperature (K) Greek symbols Dynamic viscosity (N s m −2) Density (kg m −3) Dimensionless temperature Solid volume fraction Dimensionless non-Newtonian viscosity Subscripts f Pure fluid s Nanoparticle nf Nanofluid