A mathematical model will be analyzed in order to study the effects of viscous dissipation and Ohmic heating (Joule heating) on magnetohydrodynamic (MHD) natural convection flow of a temperature dependent viscosity from heated vertical wavy surface. The present physical problem is studied numerically by using the appropriate variables, which reduce the complex wavy surface into a flat one. An implicit marching Chebyshev collocation scheme is employed for the analysis. Numerical solutions are obtained for velocity, temperature, local skin friction, and Nusselt number for a selection of parameter sets consisting of Eckert number, Prandtl number, MHD variation, and amplitude-wavelength ratio parameter. Numerical results show that these parameters have significant influences on the velocity and the temperature profiles as well as for the local skin friction and Nusselt number