In this thesis, the problem of the magnetohydrodynamics boundary layer flow of an incompressible viscous fluid over a non-linear stretching sheet is studied. Using a suitable similarity transformation, the governing partial differential equations are transformed into a system of nonlinear higher order ordinary differential equations. The resulting equations are solved numerically using implicit finite difference scheme known as Keller box method by implementing in matlab. The effects of flow parameters like magnetic field parameter M, thermal radiation R, the suction/injection parameter s, Prandtl number Pr, Eckert number Ec and power index parameter n are demonstrated graphically for velocity and temperature profiles, skin friction coefficient and surface heat transfer rate are presented numerically as well. Numerical results obtained on the values of surface heat transfer rate are compared with previously reported cases available from literature and they are found to be in a very good agreement.