This analysis discuss the influence of heat and mass transfer on natural convection boundary layer flow of thermally radiating wavy surface. To make the surface radiating; Stephan-Boltzmann law is used in the boundary conditions. Therefore, convection and thermal radiation interact simultaneously with the aid of Stephan-Boltzmann law and ultimately producing highly nonlinear boundary conditions. The governing coupled boundary layer equations are switched into suitable form by using primitive variable formulations over which two-point iterative finite difference scheme is applied to obtain the unknown quantities. Physical quantities like wall shear stress, rate of heat transfer and rate of mass transfer are expressed graphically by varying significant emerging parameters: surface radiation (R), radiative length parameter (0, buoyancy ratio parameter (N) and amplitude of the wavy surface (alpha). Comparison of numerical results is also done in tabular form with the earlier study of Siddiqa et al. (2013) in order to validate the results. (C) 2017 Elsevier Ltd. All rights reserved.