The main aim of this article is to introduce the approximate solution for MHD flow of an electrically conducting Newtonian fluid over an impermeable stretching sheet with a power law surface velocity and variable thickness in the presence of thermal-radiation and internal heat generation/absorption. The flow is caused by the non-linear stretching of a sheet. Thermal conductivity of the fluid is assumed to vary linearly with temperature. The obtaining PDEs are transformed into non-linear system of ODEs using suitable boundary conditions for various physical parameters. We use the Chebyshev spectral method to solve numerically the resulting system of ODEs. We present the effects of more parameters in the proposed model, such as the magnetic parameter, the wall thickness parameter, the radiation parameter, the thermal conductivity parameter and the Prandtl number on the flow and temperature profiles are presented, moreover, the local skin-friction and Nusselt numbers. A comparison of obtained numerical results is made with previously published results in some special cases, and excellent agreement is noted. The obtained numerical results confirm that the introduced technique is powerful mathematical tool and it can be implemented to a wide class of non-linear systems appearing in more branches in science and engineering.