Natural convection of non-Newtonian fluids along a vertical wavy surface with uniform surface temperature has been investigated using a modified power-law viscosity model. An important parameter of the problem is the ratio of the length scale introduced by the power-law and the wavelength of the wavy surface. In this model there are no physically unrealistic limits in the boundary-layer formulation for power-law, non-Newtonian fluids. The governing equations are transformed into parabolic coordinates and the singularity of the leading edge removed; hence, the boundary-layer equations can be solved straightforwardly by marching downstream from the leading edge. Numerical results are presented for the case of shear-thinning as well as shear-thickening fluid in terms of the viscosity, velocity, and temperature distribution, and for important physical properties, namely, the wall shear stress and heat transfer rates in terms of the local skin-friction coefficient and the local Nusselt number, respectively. Also results are presented for the variation in surface amplitude and the ratio of length scale to surface wavelength. The numerical results demonstrate that a Newtonian-like solution for natural convection exists near the leading edge where the shear-rate is not large enough to trigger non-Newtonian effects. After the shear-rate increases beyond a threshold value, non-Newtonian effects start to develop.