The effect of temperature dependent variable viscosity on magnetohydrodynamic (MHD) natural convection flow of viscous incompressible fluid along a uniformly heated vertical wavy surface has been investigated. The governing boundary layer equations are first transformed into a nondimensional form using suitable set of dimensionless variables. The resulting nonlinear system of partial differential equations are mapped into the domain of a vertical flat plate and then solved numerically employing the implicit finite difference method, known as Keller-box scheme. The numerical results of the surface shear stress in terms of skin friction coefficient and the rate of heat transfer in terms of local Nusselt number, the stream lines and the isotherms are shown graphically for a selection of parameters set consisting of viscosity parameter (ε), magnetic parameter (M), and Prandtl number (Pr). Numerical results of the local skin friction coefficient and the rate of heat transfer for different values are also presented in tabular form.