Interaction between nonuniform heat source/sink, magnetic field and thermal radiation with heat flux through the flow of non-Newtonian power-law fluid due to a linearly stretching sheet was studied numerically using an implicit finite difference method (FDM). The heat flux is assumed to depend on both the thermal conductivity and the thermal radiation. Besides, the effects of all governing parameters, such as the magnetic parameter, thermal conductivity parameter, the power-law index parameter, Prandtl number, the space-dependent heat source/sink parameter, the temperature-dependent heat source/sink parameter, and the radiation parameter, on the profiles of velocities and temperature are studied and discussed. In particular, thermal radiation was found to play a key role in the heat transfer characteristics and in the formation of thermal boundary layer. Generally, our numerical results reveal that both the velocity and temperature distributions are marginally influenced by both the magnetic parameter and power-law index. A good agreement is observed between our results via finite difference method and the previously published numerical results for some special case.
