In the present study the magnetohydrodynamic (MHD) liquid flow and heat transfer over non-linear permeable stretching surface has been presented in the presence of chemical reactions and partial slip. By means of proper similarity variables, the fundamental equations of the boundary layer are transformed to ordinary differential equations which for the fixed values of the x-coordinate along the plate local similarity solution would be valid appropriately. The ordinary differential equations are solved numerically using an explicit Runge-Kutta (4, 5) formula, the Dormand-Prince pair and shooting method. As a result, the velocity profiles, the concentration profiles, temperature profiles, the wall shear stress, the local Sherwood number and the local Nusselt number for the various values of the involved parameters of the problem are presented and discussed in details. (C) 2011 Elsevier Ltd. All rights reserved.