An analysis is carried out to thoroughly understand the characteristics of heat and mass transfer for the natural convection boundary layer flow along a triangular horizontal wavy surface. Combine buoyancy driven boundary layer equations for the flow are switched into convenient form via co-ordinate transformations. Full non-linear equations are integrated numerically for Pr = 0.051. Interesting results for the uneven surface are found which are expressed in the form of wall shear stress, rate of heat transfer and rate of mass transfer. Solutions are also visualized via streamlines, isotherms, and isolines for concentration. Computational results certify that, shear stress, temperature gradient and concentration gradient enhances as soon as the amplitude of the wavy surface, a, increases, but complex geometry do not allow to carry simulations for a > 1.5. This factor probably ensures that sinusoidal waveform is better than triangular waveform. <br><br><font color="red"><b> This article has been corrected. Link to the correction <u><a href="http://dx.doi.org/10.2298/TSCI170525126E">10.2298/TSCI170525126E</a><u></b></font>